[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-0003.html
More information about the gradsusr
mailing list