[gradsusr] reading in a netCDF file from NARR
Don Hooper
hoop at colorado.edu
Fri Mar 12 16:09:51 EST 2010
Joseph Pollina,
As it says at the beginning of the narropen.gs script:
**args are:
* 1 dset = full path of netCDF data file
* 2 varname = variable name in netcdf file
* 3 ntimes = number of timesteps in netcdf file
* 4 start = first time/date in netcdf file : e.g. 00z01jan1979
* 5 delta = increment between timesteps : e.g 3hr,1dy, 1mo,or 1yr
* 6 out = full path of output data descriptor file (default is ./.ddf)
* 7 vert = vertical coordinate type (sfc or press or subsfc: sfc is default)
* 8 res = horizontal resolution to regrid to (0.33 or 1.0: 0.33 is default)
So, for this file, one would say:
narropen /path/to/file.nc mslet 1 0z01jan1979 1mo mslet.ddf press 0.33
-Hoop
gradsusr-request at gradsusr.org wrote:
> Send gradsusr mailing list submissions to
> gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
> gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
> gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
> 1. reading in a netCDF file from NARR (joseph.pollina at noaa.gov)
> 2. Re: reading in a netCDF file...ncdump (joseph.pollina at noaa.gov)
> 3. help (luis blacutt)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 12 Mar 2010 14:53:18 -0500
> From: joseph.pollina at noaa.gov
> Subject: [gradsusr] reading in a netCDF file from NARR
> To: gradsusr at gradsusr.org
> Message-ID: <2e809a8f45e64f26.4b9a555e at noaa.gov>
> Content-Type: text/plain; charset="us-ascii"
>
> I have been trying to read in a netCDF file that was creasted using the NARR monthly composite site. I've tried sdfopen and I tried using xdfopen with a descriptor file, however, I'm not sure if I configured the descriptor file correctly. I've even tried the narropen.gs script on the script library page, but again, not sure if I am inputting the correct variables. I was wondering if anyone has any clue as to how to get this file open.
> Also, for some reason, sdfopen and xdfopen is not recognized as a command in grads.2.0.a7.1...at least at school, but xdfopen seems to work ok at my office at work, not sure why that would be as well (they are both the same versions).
> Any help would be GREATLY appreciated.
> Joe
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100312/f2ee6d18/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 12 Mar 2010 15:06:01 -0500
> From: joseph.pollina at noaa.gov
> Subject: Re: [gradsusr] reading in a netCDF file...ncdump
> To: joseph.pollina at noaa.gov
> Cc: gradsusr at gradsusr.org
> Message-ID: <170f6a1e2befaef4.4b9a5859 at noaa.gov>
> Content-Type: text/plain; charset="iso-8859-1"
>
> The following is the output when I ncdump the netCDF file:
> dimensions:
> y = 277 ;
> x = 349 ;
> level = 29 ;
> time = UNLIMITED ; // (1 currently)
> variables:
> float level(level) ;
> level:units = "millibar" ;
> level:actual_range = 1000.f, 100.f ;
> level:long_name = "Level" ;
> level:positive = "down" ;
> level:GRIB_id = 100s ;
> level:GRIB_name = "hPa" ;
> level:axis = "Z" ;
> level:coordinate_defines = "point" ;
> float lat(y, x) ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:standard_name = "latitude" ;
> lat:axis = "Y" ;
> lat:coordinate_defines = "point" ;
> float lon(y, x) ;
> lon:units = "degrees_east" ;
> lon:long_name = "longitude" ;
> lon:standard_name = "longitude" ;
> lon:axis = "X" ;
> lon:coordinate_defines = "point" ;
> double time(time) ;
> time:units = "hours since 1800-1-1 00:00:0.0" ;
> time:long_name = "time" ;
> time:actual_range = 1569072., 1569069. ;
> time:delta_t = "0000-00-00 03:00:00" ;
> time:standard_name = "time" ;
> time:axis = "T" ;
> time:coordinate_defines = "point" ;
> float xdist(x) ;
> xdist:long_name = "eastward distance from southwest corner of domain in projection coordinates" ;
> xdist:units = "m" ;
> float ydist(y) ;
> ydist:long_name = "northward distance from southwest corner of domain in projection coordinates" ;
> ydist:units = "m" ;
> float mslet(time, level, y, x) ;
> mslet:long_name = "Mean Sea Level (ETA model) Pressure " ;
> mslet:valid_range = -9.e+16f, 9.e+16f ;
> mslet:actual_range = -9.e+15f, 9.e+15f ;
> mslet:units = "Pa" ;
> mslet:missing_value = -9.96921e+36f ;
> mslet:_FillValue = -32767.f ;
> mslet:precision = 2s ;
> mslet:least_significant_digit = 1s ;
> mslet:var_desc = "Mean Sea Level (ETA model) Pressure " ;
> mslet:dataset = "NARR" ;
> mslet:level_desc = "?^?\r???)??????Lf????\b" ;
> mslet:statistic = "Composite" ;
> mslet:parent_stat = "Other" ;
> mslet:standard_name = "Default" ;
> int Lambert_Conformal ;
> Lambert_Conformal:grid_mapping_name = "lambert_conformal_conic" ;
> Lambert_Conformal:standard_parallel = 50., 50. ;
> Lambert_Conformal:longitude_of_central_meridian = -107. ;
> Lambert_Conformal:latitude_of_projection_origin = 50. ;
> Lambert_Conformal:false_easting = 5632642.22547 ;
> Lambert_Conformal:false_northing = 4612545.65137 ;
> // global attributes:
> :Conventions = "CF-1.0" ;
> :title = "Composite Value" ;
> :history = "dummy file: created in web ap" ;
> :comments = "Data is from North American Regional Reanalysis\n",
> "(8x/day). It consists of most variables interpolated to\n",
> "pressure surfaces from model (sigma) surfaces." ;
> :institution = "National Centers for Environmental Prediction" ;
> :platform = "Model" ;
> :references = "http://wwwt.emc.ncep.noaa.gov/mmb/rreanl/index.html\n",
> "http://www.esrl.noaa.gov/psd/data/gridded/data.narr.html" ;
> :description = "Computed from the web page http://www.esrl.noaa.gov/psd/narr/ at NOAA/ESRL PSD: monthly mean composite" ;
> data:
> level = 1000, 975, 950, 925, 900, 875, 850, 825, 800, 775, 750, 725, 700,
> 650, 600, 550, 500, 450, 400, 350, 300, 275, 250, 225, 200, 175, 150,
> 125, 100 ;
> time = 51080.1711899228 ;
> }
>
> Joe
>
> ----- Original Message -----
> From: <joseph.pollina at noaa.gov>
> Date: Friday, March 12, 2010 2:53 pm
> Subject: reading in a netCDF file from NARR
>> I have been trying to read in a netCDF file that was creasted
>> using the NARR monthly composite site. I've tried sdfopen and I
>> tried using xdfopen with a descriptor file, however, I'm not sure
>> if I configured the descriptor file correctly. I've even tried
>> the narropen.gs script on the script library page, but again, not
>> sure if I am inputting the correct variables. I was wondering if
>> anyone has any clue as to how to get this file open.
>> Also, for some reason, sdfopen and xdfopen is not recognized as a
>> command in grads.2.0.a7.1...at least at school, but xdfopen seems
>> to work ok at my office at work, not sure why that would be as
>> well (they are both the same versions).
>> Any help would be GREATLY appreciated.
>> Joe
>>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100312/82036152/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Fri, 12 Mar 2010 17:58:57 -0300
> From: luis blacutt <luis.blacutt at gmail.com>
> Subject: [gradsusr] help
> To: gradsusr at gradsusr.org
> Message-ID:
> <5825ce3f1003121258qfda4395u24c5d147deb8cc46 at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear grads users I have two problems,
>
> I generated a netcdf file using xconv but now I can't read it on grads it
> gives me the following error:
> ga-> sdfopen mday_PRCDAY_AJ.nc
> Scanning self-describing file: mday_PRCDAY_AJ.nc
> NetCDF: Index exceeds dimension bound
> SDF Error: nc_get_vara_double failed to read coordinate axis value
> gadsdf: Error reading first time value in SDF file.
>
> The second problem is related to the first one, I have the same file in a
> different format. I can access it through a ctl file. When I try to convert
> it to netcdf format I have the following error:
>
> ga-> run dr2nc.gs
> No hardcopy metafile open
> All files closed; all defined objects released;
> All GrADS attributes have been reinitialized
> Define Error: Unable to allocate data memory
> Size of request was 469045200 grid elements
> violation segment
>
> Can anyone help me?
>
> Luis
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100312/c5f84a98/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 1, Issue 10
> ***************************************
More information about the gradsusr
mailing list