Changeset 353 in lmdz_wrf

Timestamp:

Mar 10, 2015, 1:08:16 PM (10 years ago)

Author:

lfita

Message:

Adding de-accumulated total precipitation 'pr'

File:

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

trunk/tools/TS_ASCII_netCDF.py

@@ -118 +118 @@
 tsvaluesv = np.zeros( (dimt,Ntsvariables), dtype= np.float)
 
+pracc = np.zeros((dimt), dtype=np.float)
+
 itz = 0
 for it in range(dimt):
@@ -125 +127 @@
         tsvaluesv[it,iv] = np.float(tsvals[itz][iv+5])
 
+    pracc[it] = np.float(tsvals[it][16]) + np.float(tsvals[it][17])
     itz = itz + 1
-# Surface variables
+
+# time
 newvar = objofile.createVariable('time','f8',('time'))
 newvar[:] = time_stepv
@@ -139 +143 @@
     newattr = ncvar.basicvardef(newvar, tsvn[iv], tsvln[iv], tsvu[iv] )
 
+# Extra vars
+
+# pr
+varvals = np.zeros((dimt), dtype=np.float)
+varvals[1:dimt] = pracc[1:dimt] - pracc[0:dimt-1]
+
+newvar = objofile.createVariable('pr', 'f4', ('time'))
+newvar[:] = varvals
+newattr = ncvar.basicvardef(newvar, 'pr', 'precipitation', 'mm' )
+
 objofile.sync()
 objofile.close()