Changeset 1942 in lmdz_wrf for trunk/tools

Timestamp:

Jul 17, 2018, 9:50:29 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `WRFtda:' dew-point temperature

Location:

trunk/tools

Files:

• diagnostics.inf (modified) (1 diff)
• diagnostics.py (modified) (4 diffs)

trunk/tools/diagnostics.inf

@@ -36 +36 @@
 slw, WRFslw,  SH2O@DZS
 ta, WRFt, T@P@PB
-td, WRFtd, WRFp@WRFt2@QVAPOR
+td, WRFtd, WRFp@WRFt@QVAPOR
 tds, TStd, psfc@t@q
 tds, WRFtds, PSFC@T2@Q2

trunk/tools/diagnostics.py

@@ -57 +57 @@
 import datetime as dtime
 import module_ForDiag as fdin
+import module_ForDef as fdef
 import diag_tools as diag
 
@@ -89 +90 @@
   'WRFmrso', 'WRFpotevap_orPM', 'WRFp', 'WRFpsl_ecmwf',                              \
   'WRFpsl_ptarget', 'WRFrvors', 'WRFslw', 'ws', 'wds', 'wss', 'WRFheight',           \
-  'WRFheightrel', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log', 'WRFzwindMO']
+  'WRFheightrel', 'WRFtda', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log',            \
+  'WRFzwindMO']
 
 methods = ['accum', 'deaccum']
@@ -1204 +1206 @@
         ncvar.insert_variable(ncobj, 'ta', WRFt, dnames, diagoutvd, newnc)
 
+# WRFtda (WRFrh, WRFt)
+    elif diagn == 'WRFtda':
+        ARM2 = fdef.module_definitions.arm2
+        ARM3 = fdef.module_definitions.arm3
+
+        gammatarh = np.log(WRFrh) + ARM2*(WRFt-273.15)/((WRFt-273.15)+ARM3)
+        td = ARM3*gammatarh/(ARM2-gammatarh) + 273.15
+
+        dnamesvar = list(ncobj.variables['T'].dimensions)
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        # Removing the nonChecking variable-dimensions from the initial list
+        varsadd = []
+        diagoutvd = list(dvnames)
+        for nonvd in NONchkvardims:
+            if gen.searchInlist(dvnames,nonvd): diagoutvd.remove(nonvd)
+            varsadd.append(nonvd)
+
+        ncvar.insert_variable(ncobj, 'tda', td, dnames, diagoutvd, newnc)
+
 # WRFua (U, V, SINALPHA, COSALPHA) to be rotated !!
     elif diagn == 'WRFua':
@@ -1464 +1486 @@
     #   where not in the original file
     for vadd in varsadd:
-        if not gen.searchInlist(newnc.variables.keys(),vadd):
+        if not gen.searchInlist(newnc.variables.keys(),vadd) and                     \
+          dictcompvars.has_key(vadd):
             attrs = dictcompvars[vadd]
             vvn = attrs['name']