Changeset 1962 in lmdz_wrf

Timestamp:

Jul 24, 2018, 4:24:09 PM (7 years ago)

Author:

lfita

Message:

Adding `WRFtds': 2m dew-point temperature

Location:

trunk/tools

Files:

• diagnostics.inf (modified) (2 diffs)
• diagnostics.py (modified) (2 diffs)

trunk/tools/diagnostics.inf

@@ -5 +5 @@
 bils, WRFbils, HFX@LH
 cape, WRFcape_afwa, WRFt@WRFrh@WRFp@WRFgeop@HGT
-clt, clt, CLDFRA
+clt, clt, CLDFRA # WRF
 cll, cllmh, CLDFRA@WRFp
 clm, cllmh, CLDFRA@WRFp
@@ -48 +48 @@
 wa, OMEGAw, vitw@pres@temp
 wds, TSwds, u@v
-wds, wds, U10@V10
+wds, wds, U10@V10 # WRF
 ws, ws, U@V
 #wss, SPDUV10, U10@V10
 #wss, TSwss, u@v
-wss, wss, U10@V10
+wss, wss, U10@V10 # WRF
 wss, wss, u10m@v10m
 wssturb, turbulence, wss

trunk/tools/diagnostics.py

@@ -90 +90 @@
   'WRFmrso', 'WRFpotevap_orPM', 'WRFp', 'WRFpsl_ecmwf',                              \
   'WRFpsl_ptarget', 'WRFrvors', 'WRFslw', 'ws', 'wds', 'wss', 'WRFheight',           \
-  'WRFheightrel', 'WRFtda', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log',            \
+  'WRFheightrel', 'WRFtda', 'WRFtds', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log',  \
   'WRFzwindMO']
 
@@ -1228 +1228 @@
         ncvar.insert_variable(ncobj, 'tda', td, dnames, diagoutvd, newnc)
 
+# WRFtds (PSFC, T2, Q2)
+    elif diagn == 'WRFtds':
+        ARM2 = fdef.module_definitions.arm2
+        ARM3 = fdef.module_definitions.arm3
+
+        var0 = ncobj.variables[depvars[0]][:]
+        var1 = ncobj.variables[depvars[1]][:]
+        var2 = ncobj.variables[depvars[2]][:]
+
+        dnamesvar = list(ncobj.variables[depvars[1]].dimensions)
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        rhs, diagoutd, diagoutvd = diag.compute_rh(var0,var1,var2,dnamesvar,dvnamesvar)
+
+        gammatarhs = np.log(rhs) + ARM2*(var1-273.15)/((var1-273.15)+ARM3)
+        tds = ARM3*gammatarhs/(ARM2-gammatarhs)
+
+        # 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, 'tds', tds, dnames, diagoutvd, newnc)
+
 # WRFua (U, V, SINALPHA, COSALPHA) to be rotated !!
     elif diagn == 'WRFua':