thinned gaussian grib1

[Jennifer Adams]

On Feb 13, 2008, at 5:06 PM, Wesley Ebisuzaki wrote:
> ...  Not everything works
> with 2.0.a1.  For example I couldn't get ecwmf thinned gaussian
> grids to
> work.

Hi, Wesley --
I made sure that I could read that gaussian_reduced.grib1 sample file
you gave me. I used grib2ctl  v0.9.12.5p39a to write a descriptor
file and a pdef file, then I used GrADS v2 gribmap:

# .gribmap -e -v -i gaussian_reduced.grib1.ctl
grib1map:  opening GRIB file: /Users/jma/Desktop/gaussian_reduced.grib1
!!!!! MATCH:     1    1688688   0  255     1    0  129 100
500       1703   0 btim: 2006081712:00 tau:     12 dtim: 2006081800:00
grib1map:  reached end of files
grib1map:  writing the map...


And then I opened and displayed it with GrADS:


Config: v2.0.a1 little-endian readline printim grib2 netcdf hdf4-sds
Issue 'q config' command for more information.
GX Package Initialization: Size = 11 8.5
Scanning description file:  /home/Desktop/gaussian_reduced.grib1.ctl
Data file /home/Desktop/gaussian_reduced.grib1 is open as file 1
LON set to 0 360
LAT set to -89.828 89.828
LEV set to 1 1
Time values set: 2006:8:18:0 2006:8:18:0
E set to 1 1
Notice: Implied interpolation for file /home/Desktop/
gaussian_reduced.grib1.ctl
  Interpolation will be performed on any data displayed from this file
ga-> q file
File 1 : gaussian_reduced.grib1
   Descriptor: /home/Desktop/gaussian_reduced.grib1.ctl
   Binary: /home/Desktop/gaussian_reduced.grib1
   Type = Gridded
   Xsize = 1600  Ysize = 800  Zsize = 1  Tsize = 1  Esize = 1
   Number of Variables = 1
     z500mb 0 129 ** 500 mb Geopotential [m**2 s**-2]
ga-> d z500mb
Notice:  Automatic Grid Interpolation Taking Place
Contouring: 47000 to 58000 interval 1000
ga->


What part didn't work on your end?
--Jennifer
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080214/2c61ea2d/attachment.html