Changeset 2834 in lmdz_wrf

Timestamp:

Apr 29, 2020, 7:23:56 PM (5 years ago)

Author:

lfita

Message:

Fixing `var_hur' computation

File:

• trunk/tools/diag_tools.py (modified) (2 diffs)

trunk/tools/diag_tools.py

@@ -2246 +2246 @@
       [q] = mixing ratio in [kgkg-1]
       ref.: M. G. Lawrence, BAMS, 2005, 225
-    >>> print var_rh(101300., 290., 3.)
+    >>> var_hur(101300., 290., 0.003)
     0.250250256174
+    >>> var_hur(100000., 284.8892, 0.004031838)
+    0.464525998553
     """
     fname = 'var_hur'
@@ -2260 +2262 @@
     #es = 6.11*np.exp(L/Rv*(1.-273./t)/273.)
 
+    tC = t - 273.15
     # August - Roche - Magnus formula
-    es = 6.1094*np.exp(17.625*(t-273.15)/((t-273.15)+243.04))
+    es = 6.1094*np.exp(17.625*tC/(tC+243.04))
 
     # Saturated mixing ratio [g/kg]
-    ws = 0.622*es/(0.01*p-es)
-
-    hur = q/(ws*1000.)
+    ws = 0.622*es/(p/100.-es)
+
+    hur = q/ws
 
     #print 'q:', q[5,5,5,5], 't:', t[5,5,5,5], 'p:', p[5,5,5,5]