[gradsusr] Opening a netCDF file using xdfopen/sdfopen
Charles Seman
Charles.Seman at noaa.gov
Tue Apr 26 13:29:04 EDT 2011
Jeff,
Sorry to hear that adding the "t" didn't solve the problem... Are the
data "shifted on a diagonal"? (suggesting a dataset size specification
issue: is the dataset global, and if so 248 degrees doesn't cover the
globe in the xdef record...)
Hope this helps,
Chuck
Jeffrey Duda wrote:
> 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
> <mailto: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.6514998928
> <tel:1.6514998928>4377,
> > -1.17649989284377, -0.7014998928 <tel:0.7014998928>43798,
> -0.2264998928 <tel:0.2264998928>43742,
> > 0.2485001071 <tel:0.2485001071>56228, 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>
> <mailto: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>
> <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>
> > <http://www.meteor.iastate.edu/%7Ejdduda>>
> >
> >
> >
> ------------------------------------------------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> <mailto: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>
> > <mailto: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> <mailto:cjs at gfdl.noaa.gov
> <mailto:cjs at gfdl.noaa.gov>>.
> >
> >
> ********************************************************************
> > Charles Seman
> > Charles.Seman at noaa.gov <mailto:Charles.Seman at noaa.gov>
> <mailto: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/> <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>
> <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>
> > <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>
> > <http://atmos.1981-2000.01_subvars.nc>
> > atmos.1981-2000.02_subvars.nc
> <http://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>
> > <http://atmos.1981-2000.03_subvars.nc>
> > atmos.1981-2000.04_subvars.nc
> <http://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>
> > <http://atmos.1981-2000.05_subvars.nc>
> > atmos.1981-2000.06_subvars.nc
> <http://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>
> > <http://atmos.1981-2000.07_subvars.nc>
> > atmos.1981-2000.08_subvars.nc
> <http://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>
> > <http://atmos.1981-2000.09_subvars.nc>
> > atmos.1981-2000.10_subvars.nc
> <http://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>
> > <http://atmos.1981-2000.11_subvars.nc>
> > atmos.1981-2000.12_subvars.nc
> <http://atmos.1981-2000.12_subvars.nc>
> > <http://atmos.1981-2000.12_subvars.nc> climatology.nc
> <http://climatology.nc>
> > <http://climatology.nc>\n",
> > "/home/z1l/bin/tools/fregrid_parallel
> > --input_mosaic atmos_mosaic.nc <http://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>
> <mailto: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 <mailto:gradsusr at gradsusr.org>
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> _______________________________________________
> 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/~jdduda>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
--
Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.
********************************************************************
Charles Seman Charles.Seman at noaa.gov
U.S. Department of Commerce / NOAA / OAR
Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
201 Forrestal Road fax: (609) 987-5063
Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
********************************************************************
"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."
More information about the gradsusr
mailing list