Gribmap for GrADS 2.0.a2 with template option issue? (grib2 files)

Uwe Gramann weather at UNISERVE.COM
Fri Jun 6 16:48:05 EDT 2008


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 
  To: 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
      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










  --
  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/3eac5c30/attachment.html 


More information about the gradsusr mailing list