Gribmap for GrADS 2.0.a2 with template option issue? (grib2 files)
Uwe Gramann
weather at UNISERVE.COM
Mon Jun 9 17:36:49 EDT 2008
Thanks for that, Wesley,
I have written a template ctl file for now, that is updated by a simple
script ... not necessarily the smoothest solution, but it will be okay until
there is more time to go into depth!
Thanks both of you again ... my scripts are happily running!
Cheers,
Uwe
----- Original Message -----
From: "Wesley Ebisuzaki" <Wesley.Ebisuzaki at NOAA.GOV>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Monday, June 09, 2008 7:55 AM
Subject: Re: Gribmap for GrADS 2.0.a2 with template option issue? (grib2
files)
> Uwe,
>
> I haven't thought much about adding EDEF support to g2ctl. Writing
> a perl script
> would be the best. If the perl script is general enough, we could try
> integrating
> it into g2ctl later.
>
> Wesley
>
>
> Jennifer Adams wrote:
>> Wesley would have to provide the last word, but I do not believe he
>> has added support for EDEF statements in g2ctl yet.
>> I generally compose them manually using output from grib2scan, and
>> gribmap with the -v option.
>> Jennifer
>>
>> On Jun 6, 2008, at 4:48 PM, Uwe Gramann wrote:
>>
>>> Do I understand you correctly then:
>>> I cannot get g2ctl.pl to create a control file with EDEF statements,
>>> since none of the following can find my grib files:
>>>
>>> g2ctl.pl -verf gep%e.t00z.pgrb2af%f3.grib > uwe.ctl
>>> g2ctl.pl -verf gep%e.t00z.pgrb2af000.grib > uwe.ctl
>>>
>>> So I should write a little perl script that will do the job for me?
>>>
>>> Cheers,
>>> Uwe
>>>
>>>
>>>
>>>
>>> ----- Original Message -----
>>> *From:* Jennifer Adams <mailto:jma at COLA.IGES.ORG>
>>> *To:* GRADSUSR at LIST.CINECA.IT <mailto:GRADSUSR at LIST.CINECA.IT>
>>> *Sent:* Friday, June 06, 2008 1:28 PM
>>> *Subject:* Re: Gribmap for GrADS 2.0.a2 with template option
>>> issue? (grib2 files)
>>>
>>>
>>> On Jun 6, 2008, at 3:19 PM, Uwe Gramann wrote:
>>>
>>>> Hmm that's interesting, I was not aware that I HAD to use the
>>>> ensemble descriptor.
>>> You need to include the ensemble grib codes when the Product
>>> Definition Template Number is 1, 2, 11, or 12. There are examples
>>> (such as ecmwf tigge data) in which the PDT number is the only
>>> way to distinguish two fields in the same file. The GrADS
>>> interface doesn't have a place to put the PDT number, it is
>>> assumed to be 0 or 8 unless specified otherwise in the extended
>>> EDEF entry.
>>>
>>> The grib2scan utility will tell you what the PDT number is, and
>>> so will the gribmap -v output, and so will wgrib2 -Sec4.
>>>
>>>
>>>> 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??
>>>
>>> The %e template substitution marker doesn't require a number the
>>> way the time template markers do. The ensemble name (as given in
>>> your EDEF statement) is substituted completely where the %e
>>> appears in your DSET entry.
>>>
>>> Jennifer
>>>
>>>
>>>>
>>>> Thanks so much and Cheers,
>>>> Uwe
>>>>
>>>>
>>>> ----- Original Message -----
>>>> *From:* Jennifer Adams <mailto:jma at COLA.IGES.ORG>
>>>> *To:* GRADSUSR at LIST.CINECA.IT <mailto: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 <mailto:jma at COLA.IGES.ORG>
>>>>> *To:* GRADSUSR at LIST.CINECA.IT
>>>>> <mailto: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
>>>>>> <mailto:weather at uniserve.com>
>>>>>> web: www.mountainweatherservices.com
>>>>>> <http://www.mountainweatherservices.com>
>>>>>> ---------------------------------------
>>>>>>
>>>>>>
>>>>>> ----- Original Message -----
>>>>>> *From:* Jennifer Adams <mailto:jma at COLA.IGES.ORG>
>>>>>> *To:* GRADSUSR at LIST.CINECA.IT
>>>>>> <mailto: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 <mailto:jma at cola.iges.org>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>> --
>>>>> Jennifer M. Adams
>>>>> IGES/COLA
>>>>> 4041 Powder Mill Road, Suite 302
>>>>> Calverton, MD 20705
>>>>> jma at cola.iges.org <mailto: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 <mailto:jma at cola.iges.org>
>>>>
>>>>
>>>>
>>>>
>>>
>>> --
>>> Jennifer M. Adams
>>> IGES/COLA
>>> 4041 Powder Mill Road, Suite 302
>>> Calverton, MD 20705
>>> jma at cola.iges.org <mailto:jma at cola.iges.org>
>>>
>>>
>>>
>>>
>>
>> --
>> Jennifer M. Adams
>> IGES/COLA
>> 4041 Powder Mill Road, Suite 302
>> Calverton, MD 20705
>> jma at cola.iges.org <mailto:jma at cola.iges.org>
>>
>>
>>
>
More information about the gradsusr
mailing list