[gdal-dev] Problem reading WAFS GRIB2 file with gdal

[Even Rouault]

Pierre,

> A bug seems to be already opened for this error :
> https://trac.osgeo.org/gdal/ticket/5768

I had a look and just added support for that kind of products in GDAL trunk. 
The issue was that Template 4.15 ( 
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_temp4-15.shtml ) used for 
bands 4-37 wasn't supported before.

I've also refreshed the GRIB tables used to get parameter names from latest 
degrib v2.0.3

Now you get:

$ gdalinfo WAFS_blended_2016090700f06.grib2
Driver: GRIB/GRIdded Binary (.grb)
Files: WAFS_blended_2016090700f06.grib2
Size is 288, 145
Coordinate System is:
GEOGCS["Coordinate System imported from GRIB file",
    DATUM["unknown",
        SPHEROID["Sphere",6371229,0]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]
Origin = (-0.625000000000000,90.625000000000000)
Pixel Size = (1.250000000000000,-1.250000000000000)
Corner Coordinates:
Upper Left  (  -0.6250000,  90.6250000) (  0d37'30.00"W, 90d37'30.00"N)
Lower Left  (  -0.6250000, -90.6250000) (  0d37'30.00"W, 90d37'30.00"S)
Upper Right (     359.375,      90.625) (359d22'30.00"E, 90d37'30.00"N)
Lower Right (     359.375,     -90.625) (359d22'30.00"E, 90d37'30.00"S)
Center      ( 179.3750000,   0.0000000) (179d22'30.00"E,  0d 0' 0.01"N)
Band 1 Block=288x1 Type=Float64, ColorInterp=Undefined
  Description = 0[-] RESERVED(11) (Reserved)
  Metadata:
    GRIB_COMMENT=ICAO Standard Atmosphere Reference Height [m]
    GRIB_ELEMENT=ICAHT
    GRIB_FORECAST_SECONDS=21600 sec
    GRIB_PDS_PDTN=0
    GRIB_PDS_TEMPLATE_NUMBERS=3 3 2 0 96 0 0 0 1 0 0 0 6 11 0 0 0 0 0 255 0 0 
0 0 0
    GRIB_REF_TIME=  1473206400 sec UTC
    GRIB_SHORT_NAME=0-RESERVED(11)
    GRIB_UNIT=[m]
    GRIB_VALID_TIME=  1473228000 sec UTC
[...snip...]
Band 4 Block=288x1 Type=Float64, ColorInterp=Undefined
  Description = 70000[Pa] ISBL="Isobaric surface"
  Metadata:
    GRIB_COMMENT=In-Cloud Turbulence [%]
    GRIB_ELEMENT=CTP
    GRIB_FORECAST_SECONDS=21600 sec
    GRIB_PDS_PDTN=15
    GRIB_PDS_TEMPLATE_NUMBERS=19 21 2 0 96 0 0 0 1 0 0 0 6 100 0 0 1 17 112 
255 0 0 0 0 0 0 3 1
    GRIB_REF_TIME=  1473206400 sec UTC
    GRIB_SHORT_NAME=70000-ISBL
    GRIB_UNIT=[%]
    GRIB_VALID_TIME=  1473228000 sec UTC
[... snip ... ]
Band 37 Block=288x1 Type=Float64, ColorInterp=Undefined
  Description = 30000[Pa] ISBL="Isobaric surface"
  Metadata:
    GRIB_COMMENT=Icing [%]
    GRIB_ELEMENT=ICIP
    GRIB_FORECAST_SECONDS=21600 sec
    GRIB_PDS_PDTN=15
    GRIB_PDS_TEMPLATE_NUMBERS=19 20 2 0 96 0 0 0 1 0 0 0 6 100 0 0 0 117 48 
255 0 0 0 0 0 2 3 1
    GRIB_REF_TIME=  1473206400 sec UTC
    GRIB_SHORT_NAME=30000-ISBL
    GRIB_UNIT=[%]
    GRIB_VALID_TIME=  1473228000 sec UTC


Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com