我在下面有这段代码,使用xarray将一堆不同的变量和不同的lat.lon合并为netcdfs。它永远停滞不休,永无止境,但没有错误消息。
import xarray as xr
import glob
import numpy as np
datasets_WFHS= (xr.open_mfdataset(fname) for fname in
glob.glob(r'Z:\travelers\shp\test\*WFHS.nc'))
datasets_WFHS=xr.align(*datasets_WFHS, join='outer')
dataset_WFHS =np.maximum.reduce([d['AREASCORE']
for d in datasets_WFHS])
datasets_AF= (xr.open_mfdataset(fname)
for fname in glob.glob(r'Z:\travelers\shp\test\*AF.nc'))
datasets_AF=xr.align(*datasets_AF, join='outer')
dataset_AF=np.maximum.reduce([d['AREAFUEL']
for d in datasets_AF])
datasets = xr.merge([dataset_WFHS, dataset_AF], join='outer')
datasets.to_netcdf(r"Z:\travelers\shp\test\results.nc")