[RegCNET] Post-processing Omega

[Graziano Giuliani]

On 01/16/2013 03:04 PM, Rana Samuels wrote:
> Hello All,
> We are having difficulties postprocessing the omega variable. For some
> reason the code used for the rest of parameters (and works) is providing
> only a field of zeros for the Omega.
[...]

Rana,

which post-processing software are you using? One possible reason is that the 
values for omega are small. This is a vertical profile in one point in a small 
domain:

ncks -H -v omega -d iy,10 -d jx,10 -d time,3 output/verysmall_ATM.1990010100.nc

iy[10]=-300000

jx[10]=-240000

time[3]=351408 kz[0] iy[10]=-300000 jx[10]=-240000 omega[48846]=-3.37515e-05
time[3]=351408 kz[1] iy[10]=-300000 jx[10]=-240000 omega[49745]=0.000180258
time[3]=351408 kz[2] iy[10]=-300000 jx[10]=-240000 omega[50644]=0.000122793
time[3]=351408 kz[3] iy[10]=-300000 jx[10]=-240000 omega[51543]=0.000136113
time[3]=351408 kz[4] iy[10]=-300000 jx[10]=-240000 omega[52442]=4.12693e-05
time[3]=351408 kz[5] iy[10]=-300000 jx[10]=-240000 omega[53341]=0.000179324
time[3]=351408 kz[6] iy[10]=-300000 jx[10]=-240000 omega[54240]=0.00049728
time[3]=351408 kz[7] iy[10]=-300000 jx[10]=-240000 omega[55139]=0.000432289
time[3]=351408 kz[8] iy[10]=-300000 jx[10]=-240000 omega[56038]=0.000627681
time[3]=351408 kz[9] iy[10]=-300000 jx[10]=-240000 omega[56937]=0.000570556
time[3]=351408 kz[10] iy[10]=-300000 jx[10]=-240000 omega[57836]=0.000242735
time[3]=351408 kz[11] iy[10]=-300000 jx[10]=-240000 omega[58735]=-0.000960659
time[3]=351408 kz[12] iy[10]=-300000 jx[10]=-240000 omega[59634]=-0.00143808
time[3]=351408 kz[13] iy[10]=-300000 jx[10]=-240000 omega[60533]=-0.000733171
time[3]=351408 kz[14] iy[10]=-300000 jx[10]=-240000 omega[61432]=-0.000476119
time[3]=351408 kz[15] iy[10]=-300000 jx[10]=-240000 omega[62331]=-0.000337746
time[3]=351408 kz[16] iy[10]=-300000 jx[10]=-240000 omega[63230]=-0.000250138
time[3]=351408 kz[17] iy[10]=-300000 jx[10]=-240000 omega[64129]=-0.000187176

time[3]=351408

If the post-processing has a fixed format in the printout in a fortran program, 
check if the number of significant digits accomodate for the values.

Graziano.

-- 
Graziano Giuliani - Earth System Physics Section
The Abdus Salam International Centre for Theoretical Physics
Strada Costiera, 11 - I - 34151 Trieste Italy