-
Notifications
You must be signed in to change notification settings - Fork 179
Support GRIB2 PDS Template 32 #846
Comments
We are seeking a sample grib2 file that can be distributed to public. |
The HWRF data from NCEP could be a good test: from OSGeo/gdal#249
BTW, the template 4.32 is quite similar to 4.31 . @lesserwhirls I think I can add the 4.32 very easily do you want me to try? |
Sure, that would be awesome! |
@lesserwhirls I have hit a wall. The data contains info about satellite id, sensor id and band been use, but it's been defined as additional octets in the PDS template. The problem it's how to map those information... to a parameter id? to a new level type? to a new band coordinate? all of them? In the above example all data corresponds to the same grib2 parameter you can see that 18 GRIB2 messages have the same parameter and forecast step:
but they differ on platform/source, sensor/instrument, band and polarization. I have seen some examples where the band it's been used as a new coordinate ... any ideas how to map this? there is any examples for templates 4.30 or 4.31? |
HI @cofinoa! Thanks for taking a look at this. Supporting this kind of GRIB2 collection will be a quite a bit of work in the code. It's not just supporting the PDS, but we will need to extend the GRID code to include a wavelength dimension, much like we do with the ensemble dimension. Then, we will need to upgrade the GRIB variable naming code to take into account, say, platform name. Between those two enhancements, I think we will be able to uniquely identify the grids in this kind of file. This might also require a change to the GRIB index files, which means we would need to tackle this in v5.0 of the code, rather than 4.6, as we've already bumped the version of the index files there. |
@lesserwhirls Is a wavelength dimension required? Couldn't the different wavelengths be kept as separate variables? Not as elegant as a new dimension but much less work. |
True - the wavelengths could be included in the variable name to separate the grids. |
@lesserwhirls here is a sample file like discussed in Unidata/netcdf-java#285 (comment) It has 3 grib messages. Two have the same wavelength but different satellite identifiers (message 1 and 2). Two have the same satellite identifiers but different wavelengths (message 1 and 3). The values for message 2 and 3 are identical. It is a 22 by 34 grid (approximately 3.7 by 3.7 degrees). Forecast hour 43 for the 12Z 12 May 2020 model basetime. |
Excellent - thank you @ColemanTom! I'll take a look and see where things stand. |
We are missing Product Definition Section (PDS) Template 32:
ucar.nc2.grib.grib2.Grib2Pds
will need to be extended. The sample file is not licensed for redistribution (thus not attached), so will need to find a public dataset to use for tests.The text was updated successfully, but these errors were encountered: