使用 cfgrib 加载具有 "unknown" typeOfLevel 的 GRIB 变量

Question

我正在尝试使用xarray和cfgrib从DWD的ICON模型加载GRIB2文件。大多数变量工作良好，但对于某些变量(例如CLCL)，xarray.open_dataset会引发以下错误:

Traceback (most recent call last):
File "/shared/conda/envs/devtools/lib/python3.8/site-packages/xarray/core/dataset.py", line 1398, in _construct_dataarray
variable = self._variables[name]
KeyError: 'CLCL'
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "DAT-2634/grib_file_loaders.py", line 276, in <module>
df_cfgrib_sf = read_cfgrib_sf(grib_fns, nodes, params)
File "DAT-2634/grib_file_loaders.py", line 115, in read_cfgrib_sf
x = xr.open_dataset(
File "/shared/conda/envs/devtools/lib/python3.8/site-packages/xarray/core/dataset.py", line 1502, in __getitem__
return self._construct_dataarray(key)
File "/shared/conda/envs/devtools/lib/python3.8/site-packages/xarray/core/dataset.py", line 1400, in _construct_dataarray
_, name, variable = _get_virtual_variable(
File "/shared/conda/envs/devtools/lib/python3.8/site-packages/xarray/core/dataset.py", line 173, in _get_virtual_variable
ref_var = variables[ref_name]
KeyError: 'CLCL'

当我通过grib_ls检查文件时，我得到

$ grib_ls ICON_europe_reg_0.125x0.125_2021101900_f000.grib2 | grep CLCL
2            edzw         20211019     fc           regular_ll   0            unknown      800          CLCL         grid_simple 

我在cfgrib问题#195和#213中读到，这可能是由于使用了本地代码，我可能能够使用适当的代码表来解决这个问题。DWD确实为ecCodes提供了GRIB表，但我不确定如何与cfgrib结合使用它们。

我如何使用自定义代码表与cfgrib?或者是否有其他方法来加载这些变量?

Answer 1

在没有得到任何进一步的我已经发布了相同的问题作为cfgribGitHub项目的问题，并得到了一个回应:解决方案是将自定义代码表的路径添加到ECCODES_DEFINITION_PATH环境变量:

import os
from pathlib import Path
import xarray
ECCODES_DEFINITION_PATHS = [
# The custom codes I want to use
Path('path/to/the/custom/codes'),
# The codes from my system's ecCodes installation
Path('/usr/share/eccodes/definitions'),
]
os.environ['ECCODES_DEFINITION_PATH'] = ':'.join(
str(p.resolve())
for p in ECCODES_DEFINITION_PATHS
)
x = xarray.open_dataset(...)

ECMWF wiki中有详细的文档。