NOAA-EMC / EMC_verif-global

Global Forecast System (GFS) verification package using MET and METplus
8 stars 15 forks source link

Can not convert grib2 to grib1 data above 1mb #116

Open XiaqiongZhou-NOAA opened 1 year ago

XiaqiongZhou-NOAA commented 1 year ago

The "get_data_files.py" script utilizes "cvngrib -g21" to convert grib2 data for analysis and forecasting into grib1 format. However, if the data level is below 1 mb, "cnvgrib" is unable to perform the conversion. An example of a grib1 file record is:

"11:1043396:d=19120300:HGT:kpds5=7:kpds6=100:kpds7=0:TR=10:P1=0:P2=102:TimeU=1:0 mb:102hr fcst:NAve=0"

Due to this limitation, the model level must be an integer, and decimals are not allowed. Consequently, levels above 1mb cannot be verified or plotted.

malloryprow commented 1 year ago

Is this a problem or EMC_verif-global to resolve? It sounds like something isn't working right with cnvgrib.

XiaqiongZhou-NOAA commented 1 year ago

This issue impacts the capability of the EMC_verif-global. It is better to find a way to solve this. change the model level in grib2 to "Pa" from mb? any other utilities to convert? Can this be solved with cnvgrib?

malloryprow commented 1 year ago

Can you elaborate more on "change the model level in grib2 to "Pa" from mb"?

To me, this sounds like an issue with cnvgrib as if you did the same command outside EMC_verif-global it would also fail.