xarray使用及問題

1、讀取grib2

ds = xr.open_dataset('xxx.grib2', engine='cfgrib', backend_kwargs={'filter_by_keys':{'typeOfLevel': 'isobaricInhPa'}})

isobaricInhPa 氣壓座標系、垂直座標系，一般是多層高空數據
surface是地面數據
讀取會生成.idx文件，讀取報錯時記得刪除.idx文件

ds = xr.open_dataset('hrrr.t18z.wrfnatf00.grib2', engine='cfgrib', backend_kwargs={'filter_by_keys':{'stepType': 'instant', 'typeOfLevel': 'surface'}}) #地面實時數據
ds = cfgrib.open_datasets(r'hrrr.t18z.wrfnatf00.grib2') #所有數據
dat = xr.open_dataset(r'hrrr.t18z.wrfnatf00.grib2', engine='cfgrib',
                          backend_kwargs={'filter_by_keys': {'paramId': 168}}) #paramId唯一纔可

instant是實時的
accum是累計

[{'typeOfLevel': 'hybrid'}, {'typeOfLevel': 'depthBelowLandLayer'}, {'typeOfLevel': 'atmosphere'}, {'typeOfLevel': 'cloudTop'}, {'typeOfLevel': 'surface'}, {'typeOfLevel': 'heightAboveGround'}, {'typeOfLevel': 'pressureFromGroundLayer'}, {'typeOfLevel': 'sigmaLayer'}, {'typeOfLevel': 'meanSea'}, {'typeOfLevel': 'isobaricInhPa'}, {'typeOfLevel': 'heightAboveGroundLayer'}, {'typeOfLevel': 'unknown'}, {'typeOfLevel': 'isobaricLayer'}, {'typeOfLevel': 'cloudBase'}, {'typeOfLevel': 'nominalTop'}, {'typeOfLevel': 'adiabaticCondensation'}]