def ncf_write_by_varinfo(ncflnm,
dimname,
dimtype,
dimvar,
varinfos,
dimattr=None):
dimTypes = dimtype
dimVars = dimvar
# print size(dimVars[0])
# print size(dimVars[1])
dimNames = dimname
nf.netCDF_def_dims(ncflnm, dimNames, dimTypes, dimVars, dimattr=dimattr)
tc = list()
if (type(tc) == type(varinfos)):
nvar = len(varinfos)
for varinfo in varinfos:
varName = varinfo.name
varType = varinfo.type
varData = varinfo.var
vardimnames = varinfo.dimname
varattr = varinfo.varattr
# print 'winfo', varName
# print shape(varData)
# print vardimnames
nf.netCDF_var_write(ncflnm, vardimnames, varName, varType, varData,
varattr)
def ncf_write(ncflnm, dimname, dimtype, dimvar, varnames, vartypes, vars):
dimTypes = dimtype
dimVars = dimvar
dimNames = dimname
nf.netCDF_def_dims(ncflnm, dimNames, dimTypes, dimVars)
tc = list()
if (type(tc) == type(varnames)):
nvar = len(varnames)
for ivar in range(nvar):
varName = varnames[ivar]
varType = vartypes[ivar]
varData = vars[ivar]
print varName, size(varData), size(dimvar)
nf.netCDF_var_write(ncflnm, dimNames, varName, varType, varData)
def save_to_netcdf(self, gpfile, gpname):
""" dump the data to the netcdf file
"""
import netCDF_gen as nf
dimTypes = list()
dimVars = list()
dimNames = list()
for ax in self.axis_set:
dimTypes.append('f')
axvar = array(ax[:])
axname = ax.getprop('name')
dimNames.append(axname)
dimVars.append(axvar)
nf.netCDF_def_dims(gpfile, dimNames, dimTypes, dimVars)
varName = gpname
varType = 'f'
varData = self.gpdict[gpname]
nf.netCDF_var_write(gpfile, dimNames, varName, varType, varData)
rnd_err=rnd.normal(scale=err_val)
rnd_obs_err[iobs]=rnd_err
cflag=obs.obs_cflag
otime=obs.obs_time
olat=obs.obs_lat
olon=obs.obs_lon
od=obs.obs_od
lwi=obs.obs_lwi
dimTypes=['f']
dimVars=[otime]
dimNames=['time']
nf.netCDF_def_dims(ncflnm,dimNames,dimTypes, dimVars)
# lon
varName='lon'
varType='f'
varData=olon
nf.netCDF_var_write(ncflnm,dimNames,varName, varType, varData)
# lat
varName='lat'
varType='f'