使用r将多个.nc文件合并成一个多文件netCDF？

install.packages("ncdf4")
library(ncdf4)

# 创建一个新的netCDF文件
nc <- nc_create("merged.nc")

# 定义维度
nc_def_dim(nc, "time", unlimited = TRUE)
nc_def_dim(nc, "lat", length = nlat)
nc_def_dim(nc, "lon", length = nlon)

# 定义变量
nc_def_var(nc, "temperature", "double", c("time", "lat", "lon"))
nc_def_var(nc, "precipitation", "double", c("time", "lat", "lon"))

# 定义变量属性
ncatt_put(nc, "temperature", "units", "Celsius")
ncatt_put(nc, "precipitation", "units", "mm")

# 打开第一个.nc文件
file1 <- nc_open("file1.nc")

# 读取数据
data1 <- ncvar_get(file1, "temperature")

# 关闭文件
nc_close(file1)

# 打开第二个.nc文件
file2 <- nc_open("file2.nc")

# 读取数据
data2 <- ncvar_get(file2, "temperature")

# 关闭文件
nc_close(file2)

# 写入数据
ncvar_put(nc, "temperature", data1, start = c(1, 1, 1), count = c(length(data1), nlat, nlon))
ncvar_put(nc, "temperature", data2, start = c(length(data1)+1, 1, 1), count = c(length(data2), nlat, nlon))

# 关闭文件
nc_close(nc)