Plotting RTOFS grib2 data in Grads - Problems
Pablo Santos
Pablo.Santos at NOAA.GOV
Wed Sep 30 11:04:47 EDT 2009
Hi:
I am attempting to plot data from the following file:
http://www.ftp.ncep.noaa.gov/data/nccf/com/ofs/prod/ofs.20090930/ofs_atl.t00z.n000.grb.grib2
When I run it through the g2ctl.pl utility I get the output below.
Notice the pdef line. Further below I take this output from a .ctl file
and run it through the gribmap -i command and get the error message
shown complaining about the pdef line. Does anyone know what I need to
do to get passed this?
Thanks, Pablo
[ifps at MFL-LW-WAVES ncep_hourly]$ g2ctl.pl ofs_atl.t00z.n000.grb.grib2
dset ^ofs_atl.t00z.n000.grb.grib2
index ^ofs_atl.t00z.n000.grb.grib2.idx
undef 9.999E+20
title ofs_atl.t00z.n000.grb.grib2
* produced by g2ctl v0.0.4n
* griddef=1:0:(1200 x 1684):grid_template=204:Curvilinear Orthongonal
grid: see lat lon fields in this grib file:winds(grid)
dtype grib2
pdef 1200 1684 bilin sequential binary-big pdef_ncep_12.atl
xdef 1300 linear -105.04 0.1325
ydef 1500 linear -25.2 0.0678
tdef 1 linear 00Z29sep2009 1mo
zdef 1 linear 1 1
vars 15
LAPPsfc 0,1,0 0,2,202 ** surface Latitude of Presure Point [deg]
LAUVsfc 0,1,0 0,2,198 ** surface Latitude of U Wind Component of
Velocity [deg]
LAVVsfc 0,1,0 0,2,200 ** surface Latitude of V Wind Component of
Velocity [deg]
LOPPsfc 0,1,0 0,2,203 ** surface Longitude of Presure Point [deg]
LOUVsfc 0,1,0 0,2,199 ** surface Longitude of U Wind Component of
Velocity [deg]
LOVVsfc 0,1,0 0,2,201 ** surface Longitude of V Wind Component of
Velocity [deg]
MIXHTsfc 0,1,0 0,19,3 ** surface Mixed Layer Depth [m]
MNTSFsfc 0,1,0 0,2,6 ** surface Montgomery Stream Function [m^2/s]
SALINhl1p5 0,105,1,2 10,4,193 ** 1 hybrid level - 2 hybrid level 3-D
Salinity [-]
SSHGsfc 0,1,0 10,3,195 ** surface Sea Surface Height Relative to
Geoid [m]
UBAROocn 0,201,0 10,1,194 ** entire ocean (considered as a single
layer) Barotropic U velocity [m/s]
UOGRDhl1p5 0,105,1,2 10,1,2 ** 1 hybrid level - 2 hybrid level
U-Component of Current [m/s]
VBAROocn 0,201,0 10,1,195 ** entire ocean (considered as a single
layer) Barotropic V velocity [m/s]
VOGRDhl1p5 0,105,1,2 10,1,3 ** 1 hybrid level - 2 hybrid level
V-Component of Current [m/s]
WTMPChl1p5 0,105,1,2 10,4,192 ** 1 hybrid level - 2 hybrid level 3-D
Temperature [degc]
ENDVARS
[ifps at MFL-LW-WAVES ncep_hourly]$
[ifps at MFL-LW-WAVES ncep_hourly]$ gribmap -i test2.ctl
Error opening pdef file: pdef_ncep_12.atl
Open Error: Invalid pdef record.
--> The invalid description file record is:
--> pdef 1200 1684 bilin sequential binary-big pdef_ncep_12.atl
The data file was not opened.
[ifps at MFL-LW-WAVES ncep_hourly]$
More information about the gradsusr
mailing list