[gradsusr] netcdf data: Getting wrong

[M.Mahakur]

Hi, 

I have a small netCDF file (attached with this mail). It has only 42 data values (1-D), prepared a ctl to display it (pasted below). But, the values displayed (plot attached) and the values from 'gxout print' is wrong as well. 

If I read the file and also by ncdump ; the output date values seems to be okay. 

What wrong I am doing with GrADS? I am not able to guess, kindly suggest. 


ga->open test_1D.ctl 
ga-> set gxout print 
ga-> d qi 
Printing Grid -- 42 Values -- Undef = -9.99e+08 
12562.1 7.41672e-33 7.37557e-40 0 0 -9.99e+08 -9.99e+08 3.19692e-41 
0 6.72623e-44 0 9.31128e-08 2.10757e+11 5.62216e-09 5.59809e-09 5.6272e-09 
5.59809e-09 5.55654e-36 1.21997e-41 3.76158e-37 2.8026e-45 0 -9.99e+08 -9.99e+08 
-9.99e+08 -9.99e+08 -9.99e+08 -9.99e+08 2.35113e-38 0 1.50592e-36 3.76158e-37 
8.65165e-36 0 1.8996e-35 0 2.25695e-35 0 4.64906e-16 0 
1.43633e-42 0 
ga-> 

test_1D.ctl 

dset test_1D.nc 
title 1-D field 
#dtype netcdf 
undef 1.0E+15 
xdef 1 linear 90.0 0.6 
ydef 1 linear 33.5 0.5 
zdef 42 levels 1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 40 30 20 10 7 5 4 3 2 1 0.699999988079071 0.5 0.400000005960464 0.300000011920929 0.100000001490116 
tdef 1 linear 1dec2014 1yr 
vars 1 
M2IMNPASM_5_12_4_QI=>qi 42 t,z,y,x data (kg/kg) 
endvars 


$ ncdump test_1D.nc 
netcdf test_1D { 
dimensions: 
lev = 42 ; 
variables: 
float M2IMNPASM_5_12_4_QI(lev) ; 
M2IMNPASM_5_12_4_QI:long_name = "Mass fraction of cloud ice water" ; 
M2IMNPASM_5_12_4_QI:units = "kg/kg" ; 
M2IMNPASM_5_12_4_QI:_FillValue = 1.e+15f ; 
M2IMNPASM_5_12_4_QI:missing_value = 1.e+15f ; 
M2IMNPASM_5_12_4_QI:fmissing_value = 1.e+15f ; 
M2IMNPASM_5_12_4_QI:vmax = 1.e+15f ; 
M2IMNPASM_5_12_4_QI:vmin = -1.e+15f ; 
M2IMNPASM_5_12_4_QI:origname = "QI" ; 
M2IMNPASM_5_12_4_QI:fullnamepath = "/QI" ; 
M2IMNPASM_5_12_4_QI:quantity_type = "Atmospheric Moisture" ; 
M2IMNPASM_5_12_4_QI:product_short_name = "M2IMNPASM" ; 
M2IMNPASM_5_12_4_QI:coordinates = "lev" ; 
double lev(lev) ; 
lev:long_name = "vertical level" ; 
lev:units = "hPa" ; 
lev:positive = "down" ; 
lev:vmax = 1.e+15f ; 
lev:vmin = -1.e+15f ; 
lev:origname = "lev" ; 
lev:fullnamepath = "/lev" ; 
lev:standard_name = "lev" ; 

// global attributes: 
:NCO = "\"4.5.3\"" ; 
:nco_openmp_thread_number = 1 ; 
:Conventions = "CF-1.4" ; 
... 
... 
:plot_hint_x_axis_label = "kg/kg" ; 
data: 

M2IMNPASM_5_12_4_QI = 1.085931e-34, 2.425139e-15, 2.551117e-15, 
2.137763e-10, 3.473309e-10, 5.428896e-10, 1.224933e-09, 2.271552e-09, 
3.657936e-09, 4.386702e-09, 3.968771e-09, 3.949917e-09, 5.202131e-09, 
1.241103e-08, 3.866537e-08, 1.015603e-07, 3.370216e-07, 1.236549e-06, 
3.431156e-06, 6.231193e-06, 7.820233e-06, 7.548194e-06, 4.589851e-06, 
1.440135e-06, 2.245211e-08, 1.667856e-11, 7.0466e-27, 5.366316e-36, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; 

lev = 1000, 975, 950, 925, 900, 875, 850, 825, 800, 775, 750, 725, 700, 650, 
600, 550, 500, 450, 400, 350, 300, 250, 200, 150, 100, 70, 50, 40, 30, 
20, 10, 7, 5, 4, 3, 2, 1, 0.699999988079071, 0.5, 0.400000005960464, 
0.300000011920929, 0.100000001490116 ; 
} 


Thank you. 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170831/928d68e2/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_1D.nc
Type: application/x-netcdf
Size: 22814 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20170831/928d68e2/attachment-0001.nc 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.png
Type: image/png
Size: 47320 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20170831/928d68e2/attachment-0001.png