创建一个新的nc文件，合并存储所有现存以一小时划分的nc文件

Question

现在有按小时划分的若干nc文件（即，所有nc文件的时间维度为1）；
想要利用Python创建一个新的时间维度为小时数的nc文件去整合存储所有单独的文件；

思路如下：
S1：生成指定文件夹内nc文件列表filelist;
S2：创建新的nc文件ncfile；
S3：用循环读取按小时划分的nc文件，并按时间先后存储其variable到新的ncfile中。

目前写到S2发生报错“AttributeError: NetCDF: Attribute not found”。
想知道需要如何修改代码去处理此报错。代码如下，谢谢大家的建议！

import os
from netCDF4 import Dataset
import numpy as np

fileList = []

##生成指定文件夹内nc文件列表fileList
path = 'D:\\Sort out\\ChinaRainfall py\\hourly_FiT_2008_output'
dirs = os.listdir(path)
for i in dirs:
    if os.path.splitext(i)[1] == '.nc':
        fileList.append(i)
#print(fileList )
filenumber = len(fileList)

##create the merged netcdf file to store the result.

#New netCDF file
ncfile = Dataset('Chinarain2008','w',format = 'NETCDF4_CLASSIC')

#Add dimensions
xdim = ncfile.createDimension('x',617)
ydim = ncfile.createDimension('y',386)
tdim = ncfile.createDimension('time',filenumber)

#Add global attributes
ncfile.addgroupsttr('title','Chinarainfall2008')
ncfile.addgroupsttr('geospatial_lat_min','15 degrees')
ncfile.addgroupsttr('geospatial_lat_max','60 degrees')
ncfile.addgroupsttr('geospatial_lon_min','70 degrees')
ncfile.addgroupsttr('geospatial_lon_min','140 degrees')

#Add variables
variables = []
var = ncfile.createVariable('x',np.float32,[xdim])
var.addattr('long_name','xgrid')
var.addattr('units','num')
variables.append(var)
var = ncfile.createVariable('y',np.float32,[ydim])
var.addattr('long_name','ygrid')
var.addattr('units','num')
variables.append(var)
tvar = ncfile.createVariable('time',np.int,[tdim])
tvar.addattr('long_name','time')
tvar.addattr('units','hours since 2008-01-01 00:00:0.0')
variables.append(tvar)
var = ncfile.createVariable('rainfall',np.float,[tdim,ydim,xdim])
var.addattr('long_name','rainfall')
var.addattr('units','mm')
variables.append(var)

#Create netCDF file
ncfile.close()

参考代码来源：https://www.cnblogs.com/yaqiang/p/4854415.html

Answer 1

不知道你这个问题是否已经解决, 如果还没有解决的话:

• 文章：nc格式文件的读取 中也许有你想要的答案，请看下吧

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 以帮助更多的人 ^-^