Python netCDF4在shap发生变化的情况下向现有文件添加新变量

Question

我有两个netCDF文件：file1.nc和{} 唯一的区别是file1.nc包含一个变量'rho'，我想把它附加到file2.nc上，但是要修改这些变量。原始的file2.nc中不包含“rho”。 我使用的是Python模块netCDF4。在

import netCDF4 as ncd  
file1data=ncd.dataset('file1.nc')

file1data.variables['rho']

<class 'netCDF4._netCDF4.Variable'> float64 rho(ocean_time, s_rho, eta_rho, xi_rho)
    long_name: density anomaly
    units: kilogram meter-3
    time: ocean_time
    grid: grid
    location: face
    coordinates: lon_rho lat_rho s_rho ocean_time
    field: density, scalar, series
    _FillValue: 1e+37 
unlimited dimensions: ocean_time 
current shape = (2, 15, 1100, 1000) 
filling on

所以rho的形状是[2,1511001000]，但在添加到file2.nc时，我只想添加rho[1,1511001000]，即只添加第二个时间步的数据。这将导致file2.nc中的“rho”的形状为[1511001000]。但我一直做不到。在

我一直在尝试这样的代码：

^{pr2}$

我错过了什么？在

Answer 1

您尚未在第二个netCDF文件中指定变量rho的大小。在

你要做的是：

file2data.createVariable('rho','float64')

当它应该是

^{pr2}$