Changeset 1966 in lmdz_wrf for trunk

Timestamp:

Jul 24, 2018, 6:44:03 PM (7 years ago)

Author:

lfita

Message:

Adopting 'tdas' as surface dew-point temperature

Location:

trunk/tools

Files:

• diagnostics.inf (modified) (1 diff)
• diagnostics.py (modified) (5 diffs)
• variables_values.dat (modified) (1 diff)

trunk/tools/diagnostics.inf

@@ -36 +36 @@
 slw, WRFslw,  SH2O@DZS
 ta, WRFt, T@P@PB
-td, WRFtd, WRFp@WRFt@QVAPOR
+tdas, WRFtda, WRFp@WRFt@QVAPOR
 tda, WRFtda, WRFrh@WRFt
-tds, TStd, psfc@t@q
-tds, WRFtds, PSFC@T2@Q2
+tdas, TStdas, psfc@t@q
+tdas, WRFtdas, PSFC@T2@Q2
 ua, WRFua, U@V@SINALPHA@COSALPHA
 va, WRFva, U@V@SINALPHA@COSALPHA

trunk/tools/diagnostics.py

@@ -90 +90 @@
   'WRFmrso', 'WRFpotevap_orPM', 'WRFp', 'WRFpsl_ecmwf',                              \
   'WRFpsl_ptarget', 'WRFrvors', 'WRFslw', 'ws', 'wds', 'wss', 'WRFheight',           \
-  'WRFheightrel', 'WRFtda', 'WRFtds', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log',  \
+  'WRFheightrel', 'WRFtda', 'WRFtdas', 'WRFua', 'WRFva', 'WRFzwind', 'WRFzwind_log', \
   'WRFzwindMO']
 
@@ -815 +815 @@
         diagout, diagoutd, diagoutvd = diag.compute_td(var0,var1,var2,dnamesvar,dvnamesvar)
 
-        ncvar.insert_variable(ncobj, 'tds', diagout, diagoutd, diagoutvd, newnc)
+        ncvar.insert_variable(ncobj, 'tdas', diagout, diagoutd, diagoutvd, newnc)
 
 # td (psfc, t, q) from TimeSeries files
@@ -830 +830 @@
         diagout, diagoutd, diagoutvd = diag.compute_td(var0,var1,var2,dnamesvar,dvnamesvar)
 
-        ncvar.insert_variable(ncobj, 'tds', diagout, diagoutd, diagoutvd, newnc)
+        ncvar.insert_variable(ncobj, 'tdas', diagout, diagoutd, diagoutvd, newnc)
 
 # wds (u, v)
@@ -1228 +1228 @@
         ncvar.insert_variable(ncobj, 'tda', td, dnames, diagoutvd, newnc)
 
-# WRFtds (PSFC, T2, Q2)
-    elif diagn == 'WRFtds':
+# WRFtdas (PSFC, T2, Q2)
+    elif diagn == 'WRFtdas':
         ARM2 = fdef.module_definitions.arm2
         ARM3 = fdef.module_definitions.arm3
@@ -1243 +1243 @@
 
         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)
+        tdas = 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, 'tdas', tdas, dnames, diagoutvd, newnc)
 
 # WRFua (U, V, SINALPHA, COSALPHA) to be rotated !!

trunk/tools/variables_values.dat

@@ -480 +480 @@
 tda, tda, air_dew_point_temperature, 240., 310., air|dew|point|temperature, C, YlGnBu, $tda$, tda
 DWPT, tda, air_dew_point_temperature, 240., 310., air|dew|point|temperature, C, YlGnBu, $tda$, tda
-tds, tds, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
-TD, tds, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
-TDS, tds, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
-TStd, tds, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
-TStdC, tds, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
+tdas, tdas, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
+TD, tdas, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
+TDS, tdas, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
+TStd, tdas, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
+TStdC, tdas, air_dew_point_temperature, 240., 310., air|dew|point|temperature|at|2m, C, YlGnBu, $tdas$, tdas
 th105, th105, thickness_1000_500, 0., 2500., atmosphere|thickness|1000|hPa|to|500|hPa, m, Reds, $m$, m
 1000 hPa to 500 hPa thickness, th105, thickness_1000_500, 0., 2500., atmosphere|thickness|1000|hPa|to|500|hPa, m, Reds, $m$, m