<p>TRMM文件没有时间变量或维度。我无法得到军士的工具来处理这种情况,但我不是专家。我的解决方案是使用python创建文件的副本,同时根据输入文件名添加时间维度和时间值。然后按照@CharlieZender的建议使用ncrcat。对不起,我的剧本太长了。在</p>
<pre><code>#! /usr/bin/env python
import os
import glob
import netCDF4
import datetime
# get all the TRMM files with file names like, add a time dimension and time variable value from the file name.
# E.G. 3B43.19980101.7.HDF.nc, 3B43.19980201.7.HDF.nc, 3B43.20160701.7.HDF.nc
# Creates a list of monthly files which can be concatenated by NCO Tools ncrcat command.
# assumes all the 3B43.*.nc files are in the subdir files/
# creates out_3B43.*.nc in the current directory
infiles = glob.iglob( "files/*.nc" )
# add time dimension to these gridded variables. we only care about precipiation.
#wanted_vars = ['precipitation', 'relativeError', 'gaugeRelativeWeighting']
wanted_vars = ['precipitation']
for infile in sorted(infiles):
try:
nci = netCDF4.Dataset(infile, 'r')
except RuntimeError:
print "ERROR: Could not open " + infile
exit()
ofile = 'out_' + os.path.basename(infile)
try:
nco = netCDF4.Dataset(ofile,'w',clobber=True, format=nci.data_model)
except RuntimeError:
print "ERROR: Could not open " + ofile
exit()
# copy global attributes
in_g_attdict = nci.__dict__
# copy dimensions
nco.setncatts(in_g_attdict)
for dname, d in nci.dimensions.iteritems():
if d.isunlimited():
d_size = None
else:
d_size = len(d)
nco.createDimension(dname, d_size)
# Add a time dimension, as unlimited, None means unlimited.
nco.createDimension('time', None)
# Get YYYYMMDD from file name you could use dy = 15 to reflect that these are monthly files.
file_date = infile.split('.')[1]
yr = file_date[0:4]
mo = file_date[4:6]
dy = file_date[6:]
odt = datetime.datetime(year=int(yr), month=int(mo), day=int(dy))
# create the time variable. Note: the time variable must go first to be the outer or record dimension for ncrcat
time_v = nco.createVariable('time', 'int32', ('time'))
time_v.setncattr('standard_name', 'time')
time_v.setncattr('long_name', 'reference time of data field')
time_v.setncattr('axis', 'T')
# for time_urnits, I've used seconds since the epoch. But I'm guessing that 'days since XXX' should work as well. But harded to
time_v.setncattr('units', 'seconds since 1970-01-01 00:00:00 UTC')
# calculate the number.
ntime = netCDF4.date2num(odt, time_v.units)
time_v[:] = ntime
# Copy variables, skip the wanted_vars, for now
for iname, ivar in nci.variables.iteritems():
if iname in wanted_vars:
continue
ovar = nco.createVariable(iname, ivar.dtype, ivar.dimensions)
iattdict = ivar.__dict__
ovar.setncatts(iattdict)
ovar[:] = ivar[:]
# now copy the wanted 2-D gridded variables, adding a time dimension, as axis 0
for data_var_name in wanted_vars:
ivar = nci.variables[data_var_name]
# original precipation variable is dimensioned by nlon,nlat
ovar = nco.createVariable(data_var_name, ivar.dtype, ('time', 'nlon', 'nlat') )
iattdict = ivar.__dict__
ovar.setncatts(iattdict)
ovar[0,:,:] = ivar[:]
nci.close()
nco.close()
exit()
</code></pre>
<p>这会让你离开_xxxxx.nc公司可以使用当前目录和ncrcat中的文件。在</p>
^{pr2}$