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

Diane Stokes diane.stokes at noaa.gov
Sat May 26 23:36:03 EDT 2012


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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120526/3fe0b9d4/attachment-0003.html 


More information about the gradsusr mailing list