[gradsusr] Opening a netCDF file using xdfopen/sdfopen

Jeffrey Duda jdduda at iastate.edu
Mon Apr 25 21:32:16 EDT 2011


Jim,
Thanks for pointing that out.  Unfortunately, adding the "t" in the variable
list did not solve the problem.  Perhaps I am wrong about the vertical level
the field is on?  I didn't see anywhere in the ncdump output that would tell
me where vertically the variable is located, so I'm not sure.  I'm using
1000 mb as a default right now, but I could be wrong.

Jeff

On Mon, Apr 25, 2011 at 8:20 PM, James T. Potemra <jimp at hawaii.edu> wrote:

> Hi Jeff:
>
> I think since PRECLIM is a function of time, lat, lon, you need to
> include "t" in the variable definition, e.g.,
> PRECLIM 0 t,y,x PRECIPITATION
>
> Jim
>
> Jeffrey Duda wrote:
> > Chuck,
> > Thanks for the help.  I have gotten Grads to recognize the data and
> > open the .ctl file, but now the data doesn't look right.  What have I
> > done wrong?
> >
> > SGS2.ctl text:
> > DSET ^SGS.nc
> > UNDEF -1e.+34 _FillValue
> > OPTIONS 365_day_calendar
> > XDEF 248 LINEAR 1 1.0
> > YDEF 144 LEVELS
> > -88, -84, -80, -76, -72, -68, -65.4684377765723, -63.4528127641699,
> >     -61.4371877517676, -59.4215627393652, -57.4059377269629,
> >     -55.3903127145605, -53.3746877021582, -51.3590626897559,
> >     -49.3434376773535, -47.3278126649512, -45.3121876525488,
> >     -43.2965626401465, -41.2809376277441, -39.2653126153418,
> >     -37.2496876029395, -35.2340625905371, -33.2184375781348,
> >     -31.2028125657324, -29.1871875533301, -27.1715625409277,
> >     -25.1559375285254, -23.140312516123, -21.1246875037207,
> >     -19.1090624913184, -17.093437478916, -15.0778124665137,
> >     -13.0621874541113, -11.046562441709, -9.03093742930666,
> >     -7.01531241690431, -4.99968740450198, -3.07649989284378,
> >     -2.60149989284378, -2.12649989284376, -1.65149989284377,
> >     -1.17649989284377, -0.701499892843798, -0.226499892843742,
> >     0.248500107156228, 0.723500107156225, 1.19850010715621,
> > 1.6735001071562,
> >     2.14850010715623, 2.62350010715626, 3.09850010715625,
> > 3.57350010715625,
> >     4.04850010715625, 4.52350010715621, 4.99850010715626,
> > 5.47350010715624,
> >     5.94850010715622, 6.42350010715623, 6.89850010715621,
> > 7.37350010715624,
> >     7.84850010715624, 8.32350010715629, 8.79850010715625,
> > 9.27350010715626,
> >     9.74850010715623, 10.2235001071563, 10.6985001071563,
> > 11.1735001071562,
> >     11.6485001071562, 12.1235001071562, 12.5985001071563,
> > 13.0735001071562,
> >     13.5485001071562, 14.0235001071562, 14.4985001071562,
> > 14.9735001071563,
> >     15.4485001071563, 15.9235001071563, 16.3985001071563,
> > 16.8735001071562,
> >     17.3485001071563, 17.8235001071562, 18.2985001071562,
> > 18.7735001071562,
> >     19.2485001071563, 19.7235001071563, 20.1985001071563,
> > 20.6735001071563,
> >     21.1485001071562, 21.6235001071562, 22.0985001071562,
> > 22.5735001071562,
> >     23.0485001071562, 23.5235001071562, 23.9985001071562,
> > 24.4735001071562,
> >     24.9485001071563, 25.4235001071563, 25.8985001071563,
> > 26.3735001071563,
> >     26.8485001071563, 27.3235001071563, 27.7985001071563,
> > 28.2735001071563,
> >     28.7485001071563, 29.2235001071562, 29.6985001071562,
> > 30.1735001071563,
> >     30.6485001071562, 31.1235001071562, 31.5985001071562,
> > 32.0735001071562,
> >     32.5485001071562, 33.0235001071562, 33.4985001071562,
> > 33.9735001071562,
> >     34.4485001071563, 34.9235001071563, 35.3985001071563,
> > 35.8735001071563,
> >     36.3485001071563, 36.8235001071562, 37.2985001071563,
> > 37.7735001071562,
> >     39.2653126153417, 41.2809376277441, 43.2965626401464,
> > 45.3121876525488,
> >     47.3278126649512, 49.3434376773535, 51.3590626897559,
> > 53.3746877021582,
> >     55.3903127145606, 57.4059377269629, 59.4215627393652,
> > 61.4371877517676,
> >     63.4528127641699, 65.4684377765723, 68, 71.9999999999999, 76,
> >     79.9999999999999, 83.9999999999999, 87.9999999999999
> > ZDEF 1 LEVELS 1000
> > TDEF 12 LINEAR 00Z01SEP1996 1mo
> > VARS 1
> > PRECLIM 0 y,x PRECIPITATION
> > ENDVARS
> >
> > Jeff Duda
> >
> > On Mon, Apr 25, 2011 at 11:24 AM, Charles Seman
> > <Charles.Seman at noaa.gov <mailto:Charles.Seman at noaa.gov>> wrote:
> >
> >     Jeff,
> >
> >     Try this (http://iges.org/grads/gadoc/gradcomdxdfopen.html):
> >
> >     TDEF TIME 12 LINEAR 00Z01SEP1996 1mo
> >
> >     Note the use of "TIME 12" and dt = "1mo"...
> >     For reference, please find attached a GrADS ctl file which opens a
> >     NOLEAP climatology.nc <http://climatology.nc> file in GrADS v2 and
> >     an "ncdump -c" file showing the time coordinates, etc...
> >
> >     Hope this helps,
> >     Chuck
> >
> >     Jeffrey Duda wrote:
> >
> >         Hello,
> >         I need help opening a netcdf file to be viewed in Grads.  I've
> >         attached it to this email.  I tried sdfopen , but I got this
> >         error:
> >
> >         SDF Error: 365 day calendars are no longer supported by sdfopen.
> >          To open this file with GrADS, use a descriptor file with
> >          a complete TDEF entry and OPTIONS 365_day_calendar.
> >          Documentation is at
> >         http://iges.org/grads/gadoc/SDFdescriptorfile.html
> >
> >         So I tried writing a very simple descriptor file:
> >
> >         DSET ^SGS.nc
> >         OPTIONS 365_day_calendar
> >         TDEF 365 LINEAR 00Z01SEP1996 1dy
> >
> >         When I tried to open this using xdfopen, I got this error:
> >
> >         gadxdf error: Missing or invalid dimension size.
> >          --> The invalid description file record is:
> >          --> tdef 365 linear 00z01sep1996 1dy
> >          The data file was not opened.
> >
> >         So I'm not really sure what I'm doing wrong...or what I'm
> >         doing at all, for that matter.  Please help.  The output from
> >         ncdump on the file is:
> >
> >         netcdf SGS {
> >         dimensions:
> >                LON = 248 ;
> >                LAT = 144 ;
> >                bnds = 2 ;
> >                TIME = UNLIMITED ; // (12 currently)
> >         variables:
> >                double LON(LON) ;
> >                        LON:units = "degrees_east" ;
> >                        LON:modulo = 360. ;
> >                        LON:point_spacing = "even" ;
> >                        LON:axis = "X" ;
> >                double LAT(LAT) ;
> >                        LAT:units = "degrees_north" ;
> >                        LAT:point_spacing = "uneven" ;
> >                        LAT:axis = "Y" ;
> >                        LAT:bounds = "LAT_bnds" ;
> >                double LAT_bnds(LAT, bnds) ;
> >                double TIME(TIME) ;
> >                        TIME:units = "days since 1996-09-01 00:00:00" ;
> >                        TIME:time_origin = "01-SEP-1996 00:00:00" ;
> >                        TIME:calendar = "NOLEAP" ;
> >                        TIME:axis = "T" ;
> >                        TIME:bounds = "TIME_bnds" ;
> >                double TIME_bnds(TIME, bnds) ;
> >                float PRECLIM(TIME, LAT, LON) ;
> >                        PRECLIM:missing_value = -1.e+34f ;
> >                        PRECLIM:_FillValue = -1.e+34f ;
> >                        PRECLIM:long_name = "(PRECIP[D=1,L=1])" ;
> >
> >         // global attributes:
> >                        :history = "FERRET V6.65   30-Mar-11" ;
> >                        :Conventions = "CF-1.0" ;
> >         }
> >
> >         Thanks.
> >
> >         Jeff Duda
> >
> >         --
> >         Jeff Duda
> >         Iowa State University
> >         Meteorology Graduate Student
> >         3134 Agronomy Hall
> >         www.meteor.iastate.edu/~jdduda<http://www.meteor.iastate.edu/%7Ejdduda>
> >         <http://www.meteor.iastate.edu/%7Ejdduda>
> >         <http://www.meteor.iastate.edu/~jdduda
> >         <http://www.meteor.iastate.edu/%7Ejdduda>>
> >
> >
> >
> ------------------------------------------------------------------------
> >
> >         _______________________________________________
> >         gradsusr mailing list
> >         gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> >         http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >     --
> >
> >     Please note that Charles.Seman at noaa.gov
> >     <mailto:Charles.Seman at noaa.gov> should be considered my NOAA
> >     email address, not cjs at gfdl.noaa.gov <mailto:cjs at gfdl.noaa.gov>.
> >
> >     ********************************************************************
> >      Charles Seman
> >      Charles.Seman at noaa.gov <mailto:Charles.Seman at noaa.gov>
> >      U.S. Department of Commerce / NOAA / OAR
> >      Geophysical Fluid Dynamics Laboratory         voice: (609)
> >     452-6547 <tel:%28609%29%20452-6547>
> >      201 Forrestal Road                              fax: (609)
> >     987-5063 <tel:%28609%29%20987-5063>
> >      Princeton, NJ  08540-6649
> >      http://www.gfdl.noaa.gov/~cjs/ <http://www.gfdl.noaa.gov/%7Ecjs/>
> >     ********************************************************************
> >
> >     "The contents of this message are mine personally and do not
> >     reflect any
> >     official or unofficial position of the United States Federal
> >     Government,
> >     the United States Department of Commerce, or NOAA."
> >
> >     dset ^climatology.nc <http://climatology.nc>
> >     dtype netcdf
> >     options 365_day_calendar
> >     undef -999.
> >     xdef 144 linear 1.25 2.5
> >     ydef 90 linear -89 2
> >     zdef 23 levels 1000 925 850 775 700 600 500 400
> >      300 250 200 150 100 70 50 30 20 10
> >      7 5 3 2 1
> >     tdef 12 linear 12Z11JAN2000 1mo
> >     vars 5
> >     olr=>olr 0 t,y,x outgoing longwave radiation
> >     swdn_toa=>swdn_toa 0 t,y,x SW flux down at TOA
> >     swup_toa=>swup_toa 0 t,y,x SW flux up at TOA
> >     olr_clr=>olr_clr 0 t,y,x clear skyoutgoing longwave radiation
> >     swup_toa_clr=>swup_toa_clr 0 t,y,x clear skySW flux up at TOA
> >     endvars
> >
> >     netcdf climatology { // format variant: 64bit
> >     dimensions:
> >            time = UNLIMITED ; // (12 currently)
> >            level = 23 ;
> >            lat = 90 ;
> >            lon = 144 ;
> >     variables:
> >            double time(time) ;
> >                    time:long_name = "time" ;
> >                    time:units = "days since 1980-01-01 00:00:00" ;
> >                    time:cartesian_axis = "T" ;
> >                    time:calendar_type = "NOLEAP" ;
> >                    time:calendar = "NOLEAP" ;
> >                    time:climatology = "climatology_bounds" ;
> >            float level(level) ;
> >                    level:units = "hPa" ;
> >                    level:long_name = "pressure" ;
> >                    level:cartesian_axis = "Z" ;
> >                    level:axis = "Z" ;
> >                    level:positive = "down" ;
> >            float lat(lat) ;
> >                    lat:long_name = "latitude" ;
> >                    lat:units = "degrees_N" ;
> >                    lat:cartesian_axis = "Y" ;
> >            float lon(lon) ;
> >                    lon:long_name = "longitude" ;
> >                    lon:units = "degrees_E" ;
> >                    lon:cartesian_axis = "X" ;
> >            float olr(time, lat, lon) ;
> >                    olr:long_name = "outgoing longwave radiation" ;
> >                    olr:units = "watts/m2" ;
> >                    olr:missing_value = -999.f ;
> >                    olr:interp_method = "conserve_order2" ;
> >                    olr:cell_methods = "time: mean within months time:
> >     mean over years" ;
> >                    olr:time_avg_info =
> >     "average_T1,average_T2,average_DT" ;
> >            float swdn_toa(time, lat, lon) ;
> >                    swdn_toa:long_name = "SW flux down at TOA" ;
> >                    swdn_toa:units = "watts/m2" ;
> >                    swdn_toa:missing_value = -999.f ;
> >                    swdn_toa:interp_method = "conserve_order2" ;
> >                    swdn_toa:cell_methods = "time: mean within months
> >     time: mean over years" ;
> >                    swdn_toa:time_avg_info =
> >     "average_T1,average_T2,average_DT" ;
> >            float swup_toa(time, lat, lon) ;
> >                    swup_toa:long_name = "SW flux up at TOA" ;
> >                    swup_toa:units = "watts/m2" ;
> >                    swup_toa:missing_value = -999.f ;
> >                    swup_toa:interp_method = "conserve_order2" ;
> >                    swup_toa:cell_methods = "time: mean within months
> >     time: mean over years" ;
> >                    swup_toa:time_avg_info =
> >     "average_T1,average_T2,average_DT" ;
> >            float olr_clr(time, lat, lon) ;
> >                    olr_clr:long_name = "clear skyoutgoing longwave
> >     radiation" ;
> >                    olr_clr:units = "watts/m2" ;
> >                    olr_clr:missing_value = -999.f ;
> >                    olr_clr:interp_method = "conserve_order2" ;
> >                    olr_clr:cell_methods = "time: mean within months
> >     time: mean over years" ;
> >                    olr_clr:time_avg_info =
> >     "average_T1,average_T2,average_DT" ;
> >            float swup_toa_clr(time, lat, lon) ;
> >                    swup_toa_clr:long_name = "clear skySW flux up at TOA"
> ;
> >                    swup_toa_clr:units = "watts/m2" ;
> >                    swup_toa_clr:missing_value = -999.f ;
> >                    swup_toa_clr:interp_method = "conserve_order2" ;
> >                    swup_toa_clr:cell_methods = "time: mean within
> >     months time: mean over years" ;
> >                    swup_toa_clr:time_avg_info =
> >     "average_T1,average_T2,average_DT" ;
> >
> >     // global attributes:
> >                    :filename = "atmos.1981-2000.01.nc
> >     <http://atmos.1981-2000.01.nc>" ;
> >                    :title = "c48L48_am3p8" ;
> >                    :grid_type = "mosaic" ;
> >                    :grid_tile = "1" ;
> >                    :history = "Thu Mar  4 15:44:48 2010: ncrcat -O
> >     atmos.1981-2000.01_subvars.nc
> >     <http://atmos.1981-2000.01_subvars.nc>
> >     atmos.1981-2000.02_subvars.nc
> >     <http://atmos.1981-2000.02_subvars.nc>
> >     atmos.1981-2000.03_subvars.nc
> >     <http://atmos.1981-2000.03_subvars.nc>
> >     atmos.1981-2000.04_subvars.nc
> >     <http://atmos.1981-2000.04_subvars.nc>
> >     atmos.1981-2000.05_subvars.nc
> >     <http://atmos.1981-2000.05_subvars.nc>
> >     atmos.1981-2000.06_subvars.nc
> >     <http://atmos.1981-2000.06_subvars.nc>
> >     atmos.1981-2000.07_subvars.nc
> >     <http://atmos.1981-2000.07_subvars.nc>
> >     atmos.1981-2000.08_subvars.nc
> >     <http://atmos.1981-2000.08_subvars.nc>
> >     atmos.1981-2000.09_subvars.nc
> >     <http://atmos.1981-2000.09_subvars.nc>
> >     atmos.1981-2000.10_subvars.nc
> >     <http://atmos.1981-2000.10_subvars.nc>
> >     atmos.1981-2000.11_subvars.nc
> >     <http://atmos.1981-2000.11_subvars.nc>
> >     atmos.1981-2000.12_subvars.nc
> >     <http://atmos.1981-2000.12_subvars.nc> climatology.nc
> >     <http://climatology.nc>\n",
> >                            "/home/z1l/bin/tools/fregrid_parallel
> >     --input_mosaic atmos_mosaic.nc <http://atmos_mosaic.nc>
> >     --input_file 19810101.atmos_month --interp_method conserve_order2
> >     --nlon 144 --nlat 90 --scalar_field (**please see the field list
> >     in this file**)" ;
> >                    :code_version = "$Name: preQ_z1l $" ;
> >                    :comment = "FMS time averaging, version 3.0,
> >     precision=double" ;
> >                    :nco_openmp_thread_number = 1 ;
> >     data:
> >
> >      time = 7315.5, 7345, 7374.5, 7405, 7435.5, 7466, 7496.5, 7527.5,
> >     7558,
> >        7588.5, 7619, 7649.5 ;
> >
> >      level = 1000, 925, 850, 775, 700, 600, 500, 400, 300, 250, 200,
> >     150, 100,
> >        70, 50, 30, 20, 10, 7, 5, 3, 2, 1 ;
> >
> >      lat = -89, -87, -85, -83, -81, -79, -77, -75, -73, -71, -69, -67,
> >     -65, -63,
> >        -61, -59, -57, -55, -53, -51, -49, -47, -45, -43, -41, -39,
> >     -37, -35,
> >        -33, -31, -29, -27, -25, -23, -21, -19, -17, -15, -13, -11, -9,
> >     -7, -5,
> >        -3, -1, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29,
> >     31, 33,
> >        35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65,
> >     67, 69,
> >        71, 73, 75, 77, 79, 81, 83, 85, 87, 89 ;
> >
> >      lon = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75, 21.25,
> >     23.75,
> >        26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25,
> >     48.75,
> >        51.25, 53.75, 56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25,
> >     73.75,
> >        76.25, 78.75, 81.25, 83.75, 86.25, 88.75, 91.25, 93.75, 96.25,
> >     98.75,
> >        101.25, 103.75, 106.25, 108.75, 111.25, 113.75, 116.25, 118.75,
> >     121.25,
> >        123.75, 126.25, 128.75, 131.25, 133.75, 136.25, 138.75, 141.25,
> >     143.75,
> >        146.25, 148.75, 151.25, 153.75, 156.25, 158.75, 161.25, 163.75,
> >     166.25,
> >        168.75, 171.25, 173.75, 176.25, 178.75, 181.25, 183.75, 186.25,
> >     188.75,
> >        191.25, 193.75, 196.25, 198.75, 201.25, 203.75, 206.25, 208.75,
> >     211.25,
> >        213.75, 216.25, 218.75, 221.25, 223.75, 226.25, 228.75, 231.25,
> >     233.75,
> >        236.25, 238.75, 241.25, 243.75, 246.25, 248.75, 251.25, 253.75,
> >     256.25,
> >        258.75, 261.25, 263.75, 266.25, 268.75, 271.25, 273.75, 276.25,
> >     278.75,
> >        281.25, 283.75, 286.25, 288.75, 291.25, 293.75, 296.25, 298.75,
> >     301.25,
> >        303.75, 306.25, 308.75, 311.25, 313.75, 316.25, 318.75, 321.25,
> >     323.75,
> >        326.25, 328.75, 331.25, 333.75, 336.25, 338.75, 341.25, 343.75,
> >     346.25,
> >        348.75, 351.25, 353.75, 356.25, 358.75 ;
> >     }
> >
> >     _______________________________________________
> >     gradsusr mailing list
> >     gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> >     http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> >
> > --
> > Jeff Duda
> > Iowa State University
> > Meteorology Graduate Student
> > 3134 Agronomy Hall
> > www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda><
> http://www.meteor.iastate.edu/%7Ejdduda>
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110425/8303fd84/attachment-0001.html 


More information about the gradsusr mailing list