[gradsusr] netcdf data: Getting wrong
Jennifer M Adams
jadams21 at gmu.edu
Thu Aug 31 10:53:52 EDT 2017
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
More information about the gradsusr
mailing list