从多个Netcdf文件中提取给定变量并连接到单个文件

import netCDF4 as nc
import numpy as np
import xarray as xr

# 打开和读取文件
file_paths = ['file1.nc', 'file2.nc', 'file3.nc']
datasets = [nc.Dataset(file_path) for file_path in file_paths]

# 提取变量
variable_name = 'temperature'
variables = [dataset.variables[variable_name] for dataset in datasets]

# 连接数据
combined_data = xr.concat(variables, dim='time')

# 创建新文件
output_file = nc.Dataset('combined_file.nc', 'w', format='NETCDF4')

# 配置文件维度
for dim_name, dim in datasets[0].dimensions.items():
    output_file.createDimension(dim_name, len(dim))

# 配置文件变量
output_variable = output_file.createVariable(variable_name, np.float32, datasets[0].variables[variable_name].dimensions)

# 将数据写入新文件
output_variable[:] = combined_data.values

# 关闭文件
output_file.close()
for dataset in datasets:
    dataset.close()