Changeset 619 in lmdz_wrf for trunk/tools/TS_ASCII_netCDF.py

Timestamp:

Aug 21, 2015, 1:06:47 PM (9 years ago)

Author:

lfita

Message:

Adding right names for the variables more CF like!

File:

• trunk/tools/TS_ASCII_netCDF.py (modified) (2 diffs)

trunk/tools/TS_ASCII_netCDF.py

@@ -149 +149 @@
 ncvar.set_attribute(newvar, 'calendar', 'standard')
 
+dt = time_stepv[1] - time_stepv[0]
+
 # time-series variables
 for iv in range(Ntsvariables):
-    newvar = objofile.createVariable(tsvn[iv], 'f4', ('time'))
+    if tsvn[iv] == 't' or tsvn[iv] == 'u' or tsvn[iv] == 'v' or tsvn[iv] == 'q': 
+      varname, stdname, minvar, maxvar, longname, unitsvar, cbarvar =                \
+      ncvar.variables_values('TS' + tsvn[iv])
+      tsu = unitsvar
+    else:
+      varname, stdname, minvar, maxvar, longname, unitsvar, cbarvar =                \
+      ncvar.variables_values(tsvn[iv])
+      tsu = tsvu[iv]
+
+    newvar = objofile.createVariable(varname, 'f4', ('time'))
     newvar[:] = tsvaluesv[:,iv]
-    newattr = ncvar.basicvardef(newvar, tsvn[iv], tsvln[iv], tsvu[iv] )
+
+    newattr = ncvar.basicvardef(newvar, stdname, longname.replace('|',' '), tsu)
+    newattr = ncvar.set_attribute(newvar, 'wrfTSname', tsvn[iv])
+    newattr = ncvar.set_attribute(newvar, 'wrfTSdesc', tsvln[iv])
 
 # Extra vars
@@ -160 +174 @@
 varvals = np.zeros((dimt), dtype=np.float)
 varvals[1:dimt] = pracc[1:dimt] - pracc[0:dimt-1]
+varname, stdname, minvar, maxvar, longname, unitsvar, cbarvar =                      \
+      ncvar.variables_values('RAINTOT')
 
-newvar = objofile.createVariable('pr', 'f4', ('time'))
-newvar[:] = varvals
-newattr = ncvar.basicvardef(newvar, 'pr', 'precipitation', 'mm' )
+newvar = objofile.createVariable(varname, 'f4', ('time'))
+newvar[:] = varvals / dt
+newattr = ncvar.basicvardef(newvar, stdname, longname.replace('|',' '), unitsvar )
 
 objofile.sync()