У меня есть этот код ниже, используйте 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")