python 多层次grib2数据读取

以下是Python多层次grib2数据读取的示例：

1. 解压grib2文件

import bz2

pp = r'/mnt/f/2022/data/0601/W_NAFP_C_ECMF_20180601054641_P_C1D06010000060100001.bz2'
with open('/mnt/f/2022/data/test.grib', 'wb') as new_file, bz2.BZ2File(pp, 'rb') as file:
    for data in iter(lambda: file.read(100 * 1024), b''):
        new_file.write(data)

2. 读取grib2文件

import xarray as xr

data = xr.open_dataset('/mnt/f/2022/data/test.grib', engine='cfgrib')