[gradsusr] Error plotting full 120 hour GFS forecast, stops working after t30
Theo Carter
thibidottwo at hotmail.com
Sun May 27 00:51:27 EDT 2012
Hi Diane,
First, to reply to my own previous mailing - My Bad! When I read the instructions properly, I found that I should be setting all those initial YYMMDDHH's with an "i" prefix.
Then gribmap worked fine.
Plus - it was obviously some problem with the 2.2 Gig concatenated big file - because now, with the Options Template set I am able to plot right through the full 120 hours!
Thank you kindly for the help and suggestion. That was brilliant,
Thanks again,
Theo
From: thibidottwo at hotmail.com
To: gradsusr at gradsusr.org
Date: Sun, 27 May 2012 04:33:48 +0000
Subject: Re: [gradsusr] Error plotting full 120 hour GFS forecast, stops working after t30
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
_______________________________________________
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/8d562692/attachment-0003.html
More information about the gradsusr
mailing list