[gradsusr] netcdf data: Getting wrong

M.Mahakur mmahakur at tropmet.res.in
Thu Aug 31 11:22:46 EDT 2017


Thank you for your reply and suggestion. If I use 'dtype netcdf', it does not helped me. Please see plot and values as well. All gets almost equal to zero. 

ga-> d qi 
Printing Grid -- 42 Values -- Undef = -9.99e+08 
1.08593e-34 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 
2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 
2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 
2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 
2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 2.34474e-310 
2.34474e-310 2.34474e-310 
ga-> 

Best Regards, 
Mahakur 


From: "Jennifer M Adams" <jadams21 at gmu.edu> 
To: "GrADS Users Forum" <gradsusr at gradsusr.org> 
Sent: Thursday, August 31, 2017 8:23:52 PM 
Subject: Re: [gradsusr] netcdf data: Getting wrong 

Try uncommenting the ‘dytpe netcdf’ from your descriptor. —Jennifer 

> On Aug 31, 2017, at 10:22 AM, M.Mahakur <mmahakur at TROPMET.RES.IN> wrote: 
> 
> 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. 
> 
> <test_1D.nc><test.png>_______________________________________________ 
> gradsusr mailing list 
> gradsusr at gradsusr.org 
> http://gradsusr.org/mailman/listinfo/gradsusr 

-- 
Jennifer Miletta Adams 
Center for Ocean-Land-Atmosphere Studies (COLA) 
George Mason University 




_______________________________________________ 
gradsusr mailing list 
gradsusr at gradsusr.org 
http://gradsusr.org/mailman/listinfo/gradsusr 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170831/6a34179e/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_new.png
Type: image/png
Size: 45355 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20170831/6a34179e/attachment-0001.png 


More information about the gradsusr mailing list