2017-02-07 30 views
2

私は作成した多変数ncdfを持っており、各変数にデータを追加したいと考えています。緯度と経度のディメンションは同じままですが、新しいデータを各変数に追加することで時間ディメンションを拡張したいと考えています。新しいデータセットは同じ数の変数dim1とdim2を持ちますが、dim3は最初のデータセットのdim3の時間が終了するところから始まります。ここで をスタートアップコードです: 既存NCDFは(42変数を持っている)既存のnetcdfファイルに新しいデータを追加する方法

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) 

答えて

1

あなたは試してみてくださいNCOのncrcat

ncrcat in1.nc in2.nc out.nc 
1

あなたが使用して時間内にマージすることができ

cdo mergetime in1.nc in2.nc out.nc 
関連する問題