Changeset 1716 in lmdz_wrf

Timestamp:

Dec 13, 2017, 10:39:05 PM (7 years ago)

Author:

lfita

Message:

Fixing variable issues from 'var_hur'

File:

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

trunk/tools/diag_tools.py

@@ -36 +36 @@
 # var_cllmh: Fcuntion to compute cllmh on a 1D column
 # var_clt: Function to compute the total cloud fraction following 'newmicro.F90' from LMDZ using 1D vertical column values
-# var_rh: Function to compute relative humidity following 'Tetens' equation (T,P) ...'
+# var_hur: Function to compute relative humidity following 'August - Roche - Magnus' formula
 # var_mslp: Fcuntion to compute mean sea-level pressure
 # var_td: Function to compute dew-point air temperature from temperature and pressure values
@@ -1267 +1267 @@
 ####### Variables (as they arrive without dimensions staff)
 
-def var_rh(p, t, q):
-    """ Function to compute relative humidity following 'Tetens' equation (T,P) ...'
+def var_hur(p, t, q):
+    """ Function to compute relative humidity following 'August - Roche - Magnus' formula
       [t]= temperature (assuming [[t],z,y,x] in [K])
       [p] = pressure field (assuming in [Pa])
@@ -1276 +1276 @@
     0.250250256174
     """
-    fname = 'var_rh'
+    fname = 'var_hur'
 
     # Enthalpy of vaporization [Jkg-1]
@@ -1293 +1293 @@
     ws = 0.622*es/(0.01*p-es)
 
-    rh = q/(ws*1000.)
-
-    return rh
+    hur = q/(ws*1000.)
+
+    return hur
 
 def var_td(t, p, qv):
@@ -1391 +1391 @@
                 p = gen.fill_Narray(p, ta*0., filldim=[0,2,3])
 
-            self.values = var_rh(ta, p, hus)
+            self.values = var_hur(ta, p, hus)
             self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
               dictdims['lon']]