search

ARPA-SIMC / libsim

Command-line utilities and fortran modules for meteorological data processing
GNU General Public License v2.0
7 stars 1 forks source link

wet bulb temperature from grib2 data #106

Closed brancomat closed 11 months ago

brancomat commented 12 months ago

In order to plot the Thom discomfort index in arkimaps, vg6d_transform is used to compute wet bulb temperature from t2m, 2d and p100 (see: https://github.com/ARPA-SIMC/arkimaps/blob/master/recipes/thomindex.yaml)

This used to work in COSMO (grib1) but not in ICON (grib2):

$ vg6d_transform --rounding --output-variable-list=B12102 pro_wb.grib wb.grib
[stderr] ERROR    _default - Timerange with 0>p1>p2 cannot be exported in grib2

example source grib: pro_wb.tar.gz

dcesari commented 11 months ago

There was an issue with "rounding" of levels that prevented the export in grib2, it has been fixed in 5f600a5e220fb686f010d2a5be0f4be285cd05f7 . BTW as a pressure field you could use surface pressure instead of 1000hPa, it still requires rounding because T/Td are at 2m, but it is more accurate.