[gradsusr] wgrib and hirlam grib
Wesley Ebisuzaki
Wesley.Ebisuzaki at noaa.gov
Wed Dec 1 09:44:52 EST 2010
Fil,
wgrib doesn't support complex-packed grids. The error message "two
values of nxny"
is bogus. Sorry but wgrib/grib2ctl are more or less frozen with NCEP's
conversion to grib2.
AFAIK GrADS doesn't support grib1 complex-packed grids either.
Yours,
Wesley Ebisuzaki
Fil wrote:
> Hi all,
>
> I'm trying to do some plotting with hirlam rotated gribs but I'm
> getting an error with wgrib (v1.8.1.0d):
>
> # wgrib -V hirlam.grb
> grib header at record 1: two values of nxny 246768 50158
> LEN 171286 DataStart 14 UnusedBits 7 #Bits 9 nxny 246768
> rec 1:0:date 2010113018 PRES kpds5=1 kpds6=103 kpds7=0 levels=(0,0)
> grid=255 0 m above MSL 54hr fcst:
> PRES=Pressure [Pa]
> timerange 0 P1 54 P2 0 TimeU 1 nx 50158 ny 1 GDS grid 10 num_in_ave
> 0 missing 0
> center 214 subcenter 0 process 1 Table 1 scan: WE:SN winds(grid)
> rotated LatLon grid lat -35.500000 to 32.180000 lon -46.500000 to
> 46.460000
> nxny 50158 (582 x 424) dx 160 dy 160 scan 64 mode 136
> transform: south pole lat -35.000000 lon -15.000000 rot angle
> 0.000000
> *** Cannot decode complex packed fields n=50158***
>
>
> With grib2ctl.pl <http://grib2ctl.pl> (v0.9.12.5p46):
>
> # grib2ctl.pl <http://grib2ctl.pl> hirlam.grb > hirlam.ctl
> grib header at record 1: two values of nxny 246768 50158
> LEN 171286 DataStart 14 UnusedBits 7 #Bits 9 nxny 246768
> grib header at record 2: two values of nxny 246768 50158
> LEN 154832 DataStart 14 UnusedBits 6 #Bits 9 nxny 246768
> grib header at record 3: two values of nxny 246768 50158
> LEN 231734 DataStart 14 UnusedBits 0 #Bits 10 nxny 246768
> grib header at record 4: two values of nxny 246768 50158
> LEN 232468 DataStart 14 UnusedBits 14 #Bits 10 nxny 246768
> grib header at record 5: two values of nxny 246768 50158
> LEN 137810 DataStart 14 UnusedBits 13 #Bits 7 nxny 246768
> grib header at record 6: two values of nxny 246768 50158
> LEN 256406 DataStart 14 UnusedBits 7 #Bits 14 nxny 246768
> grib header at record 7: two values of nxny 246768 50160
> LEN 219698 DataStart 14 UnusedBits 0 #Bits 10 nxny 246768
> grib header at record 1: two values of nxny 246768 50158
> LEN 171286 DataStart 14 UnusedBits 7 #Bits 9 nxny 246768
> *** Cannot decode complex packed fields n=50158***
>
>
> # gribmap -i hirlam.ctl
> grib1map: opening GRIB file: hirlam.grb
> grib1map error: GRIB file format error (rc = 53)
>
> I have attached the grib file.
>
> Any help appreciated.
>
> Regards.
> ------------------------------------------------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
More information about the gradsusr
mailing list