Changeset 884 in lmdz_wrf

Timestamp:

Jun 18, 2016, 8:05:16 PM (9 years ago)

Author:

lfita

Message:

Adding `WRFt': ta from WRF

File:

• trunk/tools/diagnostics.py (modified) (1 diff)

trunk/tools/diagnostics.py

@@ -1473 +1473 @@
         ncvar.insert_variable(ncobj, 'rvors', diagout, dnamesvar, dvnamesvar, newnc)
 
+# WRFt (T, P, PB)
+    elif diag == 'WRFt':
+        var0 = ncobj.variables[depvars[0]][:]
+        var1 = ncobj.variables[depvars[1]][:]
+        var2 = ncobj.variables[depvars[2]][:]
+
+        p0=100000.
+        p=var1 + var2
+
+        WRFt = (var0 + 300.)*(p/p0)**(2./7.)
+
+        dnamesvar = list(ncobj.variables[depvars[0]].dimensions)
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        ncvar.insert_variable(ncobj, 'ta', WRFt, dnames, dvnames, newnc)
 
 # WRFtime (u10, v10, WRFpos)