Как добавить новые данные в существующий файл netcdf

У меня есть многовариантный ncdf, который я создал и хотел добавить дополнительные данные к каждой переменной. Значения широты и долготы останутся прежними, но я хочу расширить измерение времени, добавив новые данные к каждой переменной. Новый набор данных имеет такое же количество переменных, dim1 и dim2, но его dim3 начинается там, где заканчивается время dim3 первого набора данных. Существующий ncdf (имеет 42 переменных): Вот код запуска:

library(ncdf4)
dim1 = ncdim_def("lat")
dim2 = ncdim_def( "long")
dim3 = ncdim_def( "time", "days since 2004-01-01", as.integer(time))
Var<-c("a","b","c","d",.....) # variables of existing "merged.nc" file
unit<-c("aa","ab","ac","ad",...)
mat<-(n by m data matrix)
mync = nc_open('merged.nc', write=TRUE)
for (k in 2:length(var)){
ncvar_put(mync,var[k],mat[,k])
}
nc_close(mync)

2 ответа

Вы можете попробовать ncrcat от NCO

ncrcat in1.nc in2.nc out.nc

Вы можете объединить во времени, используя

cdo mergetime in1.nc in2.nc out.nc
Другие вопросы по тегам