Changeset 1183 in lmdz_wrf for trunk/tools
- Timestamp:
- Oct 12, 2016, 3:11:31 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/tools/nc_var_tools.py
r1182 r1183 11015 11015 11016 11016 if varn == 'all': 11017 varstats = list(ncobj.variables )11017 varstats = list(ncobj.variables.keys()) 11018 11018 else: 11019 11019 varstats = [varn] … … 11038 11038 # Dimension variables 11039 11039 add_vars(ncobj,onewnc,vardims) 11040 for vardn in vardims: 11041 varstats.remove(vardn) 11040 11042 11041 11043 # Statistics var … … 11045 11047 onewnc.sync() 11046 11048 11047 for vn in varstats: 11049 for Vn in varstats: 11050 print ' ' + fname + ": statistics of '" + Vn + "' ..." 11048 11051 statsvariable = {} 11049 11052 Ntotrun = 1 … … 11052 11055 print ' ' + fname + '______ ______ _____ ____ ___ __ _' 11053 11056 for vn in varstats: 11054 if not ncobj.variables.has_key( vn):11057 if not ncobj.variables.has_key(Vn): 11055 11058 print errormsg 11056 print ' ' + fname + ": file do not have variable '" + vn + "' !!"11059 print ' ' + fname + ": file do not have variable '" + Vn + "' !!" 11057 11060 quit(-1) 11058 11061 11059 objfield = ncobj.variables[ vn]11062 objfield = ncobj.variables[Vn] 11060 11063 11061 11064 dtype = objfield.dtype … … 11071 11074 if len(runvLdims) == 0: 11072 11075 print warnmsg 11073 print ' ' + fname + ': variable without running dimension!'11074 print ' computing statistics for the whole variable'11076 print ' ' + fname + ': variable without running dimension!' 11077 print ' computing statistics for the whole variable' 11075 11078 Ntotrun = 1 11076 11079 … … 11112 11115 11113 11116 if stdout: 11114 print ' ' + vn + '[' + Srundim + ']... .. .'11117 print ' ' + Vn + '[' + Srundim + ']... .. .' 11115 11118 if stats == 'full': 11116 11119 minv = np.min(field) … … 11156 11159 print ' ' 11157 11160 print Sm + 'Matrix of the statisitcs _______' 11158 print Sm + ' '.join('{:<15}'.format( vn) for vn in statns)11161 print Sm + ' '.join('{:<15}'.format(Vn) for vn in statns) 11159 11162 for key in statsvariable.keys(): 11160 11163 vals = statsvariable[key] 11161 11164 print Sm+'{:<15}'.format(key)+' '+' '.join('{:<15g}'.format(v) for v in vals) 11162 11165 11163 newvar = onewnc.createVariable( vn + '_stats', 'f4', tuple(['stats'] + \11166 newvar = onewnc.createVariable(Vn + '_stats', 'f4', tuple(['stats'] + \ 11164 11167 runvdimns)) 11165 11168 newvar[:] = statsresults[:]
Note: See TracChangeset
for help on using the changeset viewer.