NetCDFファイルのサイズを小さくしてデータを変更するにはどうすればよいですか？

Question

dimensions: 
    i1 = 3 ; 
    x = 11 ; 
    s1 = 1 ; 
    mid1 = 8 ; 
    mid2 = 8 ; 
variables: 
    double Height(i1,x) ; 
    double Temp(s1, x) ; 
    short Soil(s1, x) ; 
    double Liq(mid1, x) ; 

は私が寸法mid1の1のサイズを小さくし、値を置き換えしたいのnetCDFファイルを持っている：

icond <- ncdf4::nc_open('dat.nc) 

#New dimensions for new file 
idim <- icond$dim[['i1']] 
xdim <- icond$dim[['x']] 
s1dim <- icond$dim[['s1']] 
mid1dim <- ncdim_def("mid1", "", 1:3) #3 layers 
mid2dim <- icond$dim[['mid2']] 

mv <- -9999 

#Get variable data 
Liqxdat <- ncvar_get(icond, 'Liq')[,1:3] #3 Layers 

#Define new variable 
Liqx = ncvar_def("Liq", "units", list(mid1dim, i1), mv, prec="double") 

#Create netCDF file 
nc = nc_create("test.nc", list(Height, Temp, Soil, Liqx) 

#Write data to the NetCDF file 
ncvar_put(nc, Liqx, Liqxdat) 

をしかし、これは私に、出力ファイル内のデータを与えていません。

Answer 1

残念ながら、コードを実行できません。だから私はちょうどそれが動作していない理由を推測することができます。

1. nc = nc_create（..）行にクロージングがありません。
2. 高さ、温度、および土壌は定義されていません。
3. nc_createを使用すると、新しいファイルを作成できますが、既存のファイルにデータを追加することはありません。新しい変数を追加する場合は、使用する必要があります。ncvar_add()