python xr.open_dataset 读取NCEP累计降水的bug

毕业两年的小白 · 发表于 2023-10-26 12:51:13

本帖最后由毕业两年的小白于 2023-10-26 12:53 编辑

从NCEP下载的标准grib2预测数据，想要读取其中累计降水
利用Panoply检查数据，发现tp有两个时间序列：

python xarray代码
import xarray as xr
fp = 'gdas.t00z.pgrb2.0p25.f009.grib2'
flt = {'shortName':'tp'}
d1 = xr.open_dataset(fp, engine='cfgrib', backend_kwargs={'indexpath':'', 'filter_by_keys': flt})

这样得出的d1只有第二个时间序列，如何获取第一个时间的数据呢？

毕业两年的小白 · 发表于 2023-10-27 11:37:23

已解决，在filter_keys中添加stepRange:'0-9'

毕业两年的小白 · 发表于 2023-10-27 11:49:59

公众号气象备忘录搜索利用python绘制GFS预报降水，有解决我这个问题，一模一样的

黄45 · 发表于 2023-10-28 08:59:29

毕业两年的小白发表于 2023-10-27 11:49
公众号气象备忘录搜索利用python绘制GFS预报降水，有解决我这个问题，一模一样的

感谢版主经验分享

		自动登录	找回密码
密码			立即注册