[gradsusr] Error plotting full 120 hour GFS forecast, stops working after t30

Sun May 27 00:33:48 EDT 2012

Hi Diane - thanks for the suggestion - perhaps it is a big file problem.

Trying dset ^%y2%m2%d2%h2.gfs.t%h2z.pgrb2f%f2 and then running grib2map for my gfs files which look like this:
12052418.gfs.t18z.pgrb2f00
12052418.gfs.t18z.pgrb2f03
etc..

only the first hour gets gribmapped - on the second file it tries to scan 
grib2map: scanning GRIB2 file: 12052421.gfs.t21z.pgrb2f03 

Where it is incrementing the 18 hour parts to 21 - but that should stay the same as with year, month and day?

I have never tried the templates option, is that %h2 the correct replacement to use

Regards,
Theo

Date: Sat, 26 May 2012 23:36:03 -0400
From: diane.stokes at noaa.gov
To: gradsusr at gradsusr.org
Subject: Re: [gradsusr] Error plotting full 120 hour GFS forecast, stops working after t30

    Theo,

    Have you tried using the template option in your descriptor file to
    process the individual files rather than concatenating them into one
    large file?  See:

      http://grads.iges.org/grads/gadoc/templates.html

    for more info.

         Diane

    On 5/26/2012 10:14 PM, Theo Carter wrote:

        Hi Diane,

        This is a new attempt at downloading gfs and creating text file
        output - so I cannot say it worked for me before: however, I
        have used a similar way (cat the files, create .ctl and idx
        files) to create graphical (not text) output using GrADS from
        wrf modeling. The one difference is that the wrf I use only
        forecasts out to 48 hours, not 120.

        I use gribmaster to download full gfs004grb2 3-hourly 0.5 deg
        files from f0 to f120 from the NCEP server.

        I then concatenate the files together using "cat" into one large
        file of 2.2 Gig:

        cat *pgrb2*>>gfsbig

        Then create .ctl and .idx using:

        ./g2ctl -verf gfsbig >gfsgrads.ctl

        gribmap -v -i gfsgrads.ctl

        (Running gribmap manually shows a "Match" in every instance)

        My created .ctl file shows:

        dtype grib2

        ydef 361 linear -90.000000 0.5

        xdef 720 linear 0.000000 0.500000

        tdef 41 linear 18Z24may2012 3hr

        And from within grads:

        "q file 1"  gives:  ! Xsize = 720  Ysize = 361  Zsize = 26 
        Tsize = 41  Esize = 1

        So all the times appear to be available.

        My hardware environment is 5.8GiB Ram and 4 processors Xeon
        E5520 at 2.27GHz allocated.

        My software environment is Scientific Linux 6.2_64 in a virtual
        machine (Basically same as Enterprise Linux). 

        Grads Version 2.0.a9 from the EPEL repository and wgrib2 from
        the same. Built Fri Sep  3 17:16:06 UTC 2010 for
        x86_64-redhat-linux-gnu/

        q config reveals:

        This version of GrADS has been configured with the following
        options:

          o Built on a LITTLE ENDIAN machine

          o Athena Widget GUI DISABLED

          o Command line editing ENABLED 

          o printim command for image output ENABLED 

          o GRIB2 interface ENABLED 

          o NetCDF interface ENABLED 

                netcdf 4.1.1 of Dec 27 2010 21:12:42 $  

          o OPeNDAP gridded data interface ENABLED

          o OPeNDAP! station data interface DISABLED

          o HDF4 and HDF5 interfaces ENABLED 

              HDF 4.2r5 

              HDF5 1.8.5 

          o GeoTIFF and KML/TIFF output ENABLED

          o KML contour output ENABLED

          o Shapefile interface ENABLED

          Date: Fri, 25 May 2012 09:50:56 -0400

          From: diane.stokes at noaa.gov

          To: gradsusr at gradsusr.org

          Subject: Re: [gradsusr] Error plotting full 120 hour GFS
          forecast, stops working after t30

          Theo,

          Answers to the following questions may help the forum point
          you in the right direction.

          What is the source of your data?

          Are you using the template option to read a series of files,
          or  are you reading one large file?

          Did the same request work for you for GFS runs prior to
          20120522?

          What version of GrADS are you running?  (Please provide output
          from 'q config')

            Regards,

              Diane Stokes

          On 5/25/2012 12:23 AM, Theo Carter wrote:

             Hi All,

              I am trying to plot a variable using 120 hours of a GFS
              run, but for some reason I cannot plot past T30 (T41 being
              120 hours). Everything plots fine up to that point, but
              anything past T30 crashes out with the GRIB2 I/O error
              specified below.

              Has anyone else come across this?

              Cheers,

              Theo

              ga-> set t 1 41

              Time values set: 2012:5:24:18 2012:5:29:18 

              ga-> d tmpsfc

              GRIB2 I/O error: fseeko failed 

              Data Request Error:  Error for variable 'tmpsfc'

                Error ocurred at column 1

              DISPLAY error:  Invalid expression 

                Expression = tmpsfc

              ga-> c

              ga-> set t 1 30

              Time values set: 2012:5:24:18 2012:5:28:9 

              ga-> d tmax2m

              ga-> d tmp2m

              ga-> 

          _______________________________________________
          gradsusr mailing list
          gradsusr at gradsusr.org
          http://gradsusr.org/mailman/listinfo/gradsusr

      _______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120527/e5d5e384/attachment-0003.html 