Ensemble GRIBs decoding
Mandelli Alessandro
mandelli.alessandro at TISCALI.IT
Sun Jan 23 10:50:53 EST 2005
I would like to make use of the gribs found in
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/mrf/prod/
I downloaded enspost.t00z.z500 and the index file enspost.t00z.z500i
Then I created the control file
ga-> !perl grib2ctl.pl -verf enspost.t00z.z500 enspost.t00z.z500i
>enspost.ctl
Result is
dset ^enspost.t00z.z500
index ^enspost.t00z.z500i
undef 9.999E+20
title enspost.t00z.z500
* produced by grib2ctl v0.9.12.5p33j beta11
dtype grib 2
options yrev
ydef 73 linear -90.000000 2.5
xdef 144 linear 0.000000 2.500000
tdef 65 linear 00Z22jan2005 6hr
zdef 1 linear 1 1
vars 1
HGTprs 0 7,100,500 ** Geopotential height [gpm]
ENDVARS
I just modified "zdef 1 linear 500 500" for consistency.
Then I mapped the grib
ga->!gribmap -i enspost.ctl
So far so good. I'm proudly plotting all timesteps of the variable.
Fact is that I'm pretty sure that in the grib all 12 members for the 65
timsteps are included.
In fact
ga-> !wgrib -v enspost.t00z.z500
gives fcst+0, -0, +1, -1, etc. for all the 12 members.
Unfortunately I'm not able to extract the members from the grib.
Anyone has a control file or a suggestion how to do it?
Thanks
Regards
More information about the gradsusr
mailing list