当我想用python在现有的netcdf文件上添加新的组时,我很困难。MWE在我的机器上工作。
分子量
import os
import numpy as np
from scipy.io import netcdf
from netCDF4 import Dataset
try:
os.remove('test.nc')
except OSError:
pass
mydata = np.linspace(1,30,30)
g=Dataset('test.nc','w',format='NETCDF4')
g.description = 'The data in this file was created for a test'
g.close()
#Do stuff
g=Dataset('test.nc','r+',format='NETCDF4')
g1 = g.createGroup('grp1')
g1.createDimension('dim', len(mydata))
ncdata = g1.createVariable('data1', 'f8', ('dim',)); ncdata[:] = mydata
ncdata = g1.createVariable('data2', 'f8', ('dim',)); ncdata[:] = mydata
g.close()
#Do stuff
g=Dataset('test.nc','r+',format='NETCDF4')
g1 = g.createGroup('grp2')
g1.createDimension('dim', len(mydata))
ncdata = g1.createVariable('data1', 'f8', ('dim',)); ncdata[:] = mydata
ncdata = g1.createVariable('data2', 'f8', ('dim',)); ncdata[:] = mydata
g.close()
但是我无法打开我的文件
ncview
相反,我明白了
Ncview 1.93g大卫·W·皮尔斯2009年2月24日
http://meteora.ucsd.edu:80/~pierce/ncview_home_page.html
版权(C)
1993年至2009年,David W.Pierce Ncview以绝对不
担保;详细信息类型
ncview -w'. This is free software licensed
under the Gnu General Public License version 3; type
ncview-c'用于
重新分发详细信息。
找不到可显示的变量!
上面写着
no displayable variables found!
,但通过
ncdump -h test.nc
我能够看到以下几行:;
netcdf测试{
//全局属性:
:description=“此文件中的数据是为测试创建的”;
组:grp1{尺寸:
dim=30;变量:
双数据1(dim);
双数据2(dim);}//组grp1
组:grp2{尺寸:
dim=30;变量:
双数据1(dim);
双数据2(dim);}//组grp2}
我是不是在我的
Python
剧本(我在想
ncview公司
是可以的,因为当我不将它们保存在
netcdf组
)
