Plotting RTOFS grib2 data in Grads - Problems

[Wesley Ebisuzaki]

Pablo,

   The "GDT Template ... not defined message" means that you are using a
old version
of g2clib.  Perhaps a new version of GrADS includes the latest g2clib.
If not recompile
GrADS with the current version of g2clib.  (Works for me.)

       Wesley Ebisuzaki




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