Megan输出nc文件合并，这是我写的python代码，有点问题，请指导

cmaq学习小白

import xarray as xr

# 定义文件的起始和结束日期
start_date = 2017153
end_date = 2017159

# 定义文件路径和名称的模板
file_template = 'MEGANv31.d01.mycase.CB6.{date}.BDSNP.ncf'

# 构建要合并的数据集列表
datasets = []
for date in range(start_date, end_date + 1):
    filename = file_template.format(date=str(date))
    ds = xr.open_dataset(filename)
    datasets.append(ds)

# 使用手动指定拼接的维度进行合并
merged_dataset = xr.concat(datasets, dim='TSTEP')

# 输出到新的NetCDF文件
merged_dataset.to_netcdf('merged_file.nc')

# 关闭数据集
for ds in datasets:
    ds.close()

DustforCMAQWRF

能用吗

                               
登录/注册后可看大图

来自苹果APP客户端