Plotting RTOFS grib2 data in Grads - Problems

[Pablo Santos]

Well:

I was directed to download pdef_ncep_12.atl. Now I got passed that but
am getting the following error when I run the control file through the
gribmap -i command:

gribmap -i test2.ctl
grib2map: scanning GRIB2 file: ofs_atl.t00z.n000.grb.grib2
getgridtemplate: GDT Template 3.204 not defined.
grib2map error: g2_getfld failed: ierr=10

Any suggestions on how to get passed this. I am looking into converting
the file to grib1 and see if I am luckier now. But thought I ask anyhow
in case someone knows.

Thanks,

Pablo

Pablo Santos wrote:
> 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]$
>