Gribmap for GrADS 2.0.a2 with template option issue? (grib2 files)
Uwe Gramann
weather at UNISERVE.COM
Fri Jun 6 15:19:22 EDT 2008
Hmm that's interesting, I was not aware that I HAD to use the ensemble descriptor.
This leads me then to a more basic question, though:
g2ctl.pl then would be able to generate the control file from the original files, right? The orginals are:
...
-rw-r--r-- 1 ugramann ugramann 64919 Jun 5 03:56 gep20.t00z.pgrb2af240.grib
-rw-r--r-- 1 ugramann ugramann 64866 Jun 5 03:56 gep20.t00z.pgrb2af252.grib
-rw-r--r-- 1 ugramann ugramann 65039 Jun 5 03:56 gep20.t00z.pgrb2af264.grib
-rw-r--r-- 1 ugramann ugramann 65079 Jun 5 03:56 gep20.t00z.pgrb2af276.grib
-rw-r--r-- 1 ugramann ugramann 64692 Jun 5 03:56 gep20.t00z.pgrb2af288.grib
-rw-r--r-- 1 ugramann ugramann 65104 Jun 5 03:56 gep20.t00z.pgrb2af300.grib
-rw-r--r-- 1 ugramann ugramann 64787 Jun 5 03:56 gep20.t00z.pgrb2af312.grib
-rw-r--r-- 1 ugramann ugramann 65517 Jun 5 03:56 gep20.t00z.pgrb2af324.grib
-rw-r--r-- 1 ugramann ugramann 65090 Jun 5 03:56 gep20.t00z.pgrb2af336.grib
-rw-r--r-- 1 ugramann ugramann 65445 Jun 5 03:56 gep20.t00z.pgrb2af348.grib
-rw-r--r-- 1 ugramann ugramann 64933 Jun 5 03:56 gep20.t00z.pgrb2af360.grib
-rw-r--r-- 1 ugramann ugramann 65708 Jun 5 03:56 gep20.t00z.pgrb2af372.grib
-rw-r--r-- 1 ugramann ugramann 65313 Jun 5 03:56 gep20.t00z.pgrb2af384.grib
but the following line does not find any files:
g2ctl.pl -verf gep%e.t00z.pgrb2af%f3.grib > uwe.ctl
(if I replace the %e template with '20' it finds the files, but not with '%e2')
I am still doing something wrong with the ensemble template??
Thanks so much and Cheers,
Uwe
----- Original Message -----
From: Jennifer Adams
To: GRADSUSR at LIST.CINECA.IT
Sent: Friday, June 06, 2008 11:17 AM
Subject: Re: Gribmap for GrADS 2.0.a2 with template option issue? (grib2 files)
Uwe,
There's no need to have a separate descriptor file for each member -- that's the whole purpose of the E dimension, to aggregate members into a single data set! An example of a GFS ensemble descriptor file with full EDEF is here:
http://iges.org/grads/gadoc/gfsens.ctl
Specifically, what your descriptor needs, is this:
edef 1
20 33 00z05jun2008 3,20
endedef
The name "20" is a guess from your filename. If you want to template over E, change DSET to
dset ^gens_%e_2008060500z_P%f3.grib
Then get the ensemble grib codes from the grib2scan output (grep on "Ens") or from the gribmap -v output (look for "ens=")
The EDEF documentation is at http://iges.org/grads/gadoc/descriptorfile.html#EDEF
Jennifer
On Jun 6, 2008, at 12:22 PM, Uwe Gramann wrote:
Hi Jennifer,
thanks for looking into this!
I have attached the output of the gribmap/grib2scan and the ctl file (which does not contain a edef declaration). I am not sure if I need an edef declaration, though, since I am treating each member as a separate model and create control files for each independently?
In case I will not get anywhere ... is there a server that has grib1 data of the gfs ensembles still?
Cheer and thanks again,
Uwe
----- Original Message -----
From: Jennifer Adams
To: GRADSUSR at LIST.CINECA.IT
Sent: Thursday, June 05, 2008 12:00 PM
Subject: Re: Gribmap for GrADS 2.0.a2 with template option issue? (grib2 files)
On Jun 4, 2008, at 10:13 PM, Uwe Gramann wrote:
Hello everybody,
I have not much luck trying to open GFS ensemble grib2 files with the template option. After downloading all members and trying to apply the time-template option to a member like:
g2ctl.pl -verf gens_20_2008060400z_P%f3.grib > uwe.ctl
gribmap -i uwe.ctl
gribmap reports no errors, the control file looks okay, but all variables are undef in Grads 2.0.a2!
What does your descriptor file look like? I don't think g2ctl handles the ensemble declarations for EDEF (yet), you may have to add those manually.
What is the output from gribmap -v ?
What is the output from grib2scan -v ?
Jennifer
I have no issues opening the individual files in Grads or combining each forecast-member's files into one large file and then opening those.
Jennifer just mentioned that a gribmap time-template bug has been fixed, but I still have one ....
Does anybody have any ideas ... suggestions ...
Cheers and thanks to all
Uwe
---------------------------------------
Uwe Gramann
Mountain Weather Services
P.O.Box 4191, 4262 Railway Ave.
Smithers, BC, V0J 2N0
Cell: +250.877.3785
e-mail: weather at uniserve.com
web: www.mountainweatherservices.com
---------------------------------------
----- Original Message -----
From: Jennifer Adams
To: GRADSUSR at LIST.CINECA.IT
Sent: Friday, May 02, 2008 9:40 AM
Subject: NEW RELEASE of GrADS 2.0.a2
Dear All,
I have posted a new version of GrADS: 2.0.a2. This version has several new features and important bug fixes. The GrADS downloads web page (http://iges.org/grads/downloads.html) has also been updated with links to the source code and two binary releases (more binaries to follow as soon as possible).
Features:
* ready for use with GDS-2.0 to access 5-dimensional data sets via OPeNDAP
* support for thinned grib2 grids
* gribmap has -0 option for grib2
* allows non-float data types for hdf coordinate axes
Bug Fixes:
* memory leak when replacing an existing defined variable
* multiple fixes for netcdf/hdf handling (templating, zrev, %nodim%, et al.)
* contour interval and label handling for double precision numbers
* changed 'query dims' output for ensembles
* fixed gribmap's handling of data sets templated over T but not E
* fixed 'set annot' command and other minor bugs
This release is a companion for GDS 2.0, which can now serve 5-D data sets. COLA is running a 2nd GDS for testing on port 9191 that is serving everything on the 9090 server plus GFS ensemble forecasts. Check 'em out at
http://monsoondata.org:9191/dods/gfsens/
A public release of GDS 2.0 is forthcoming -- we are still testing. If you'd like to help, please try using the 9191 server.
Please post comments and questions and bug reports here (low-tech bug tracking is still in use at COLA). As always, I must be able to reproduce a problem on my own systems in order to fix it, so keep that in mind when reporting bugs.
Jennifer
--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org
--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org
<gens_20_2008060500z_P000.ctl><grib2scan.with.gens_20_2008060500z_P012.grib><gribmap.of.gens_20_2008060500z_P000.ctl>
--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080606/682535c9/attachment.html
More information about the gradsusr
mailing list