Changeset 1711 in lmdz_wrf

Timestamp:

Dec 12, 2017, 11:48:07 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `var_rh': Function to compute relative humidity following 'Tetens' equation (T,P) ...'

File:

• trunk/tools/diag_tools.py (modified) (4 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_mslp: Fcuntion to compute mean sea-level pressure
 # var_td: Function to compute dew-point air temperature from temperature and pressure values
@@ -1264 +1265 @@
     return wd
 
+####### Variables (as they arrive without dimensions staff)
+
+def var_rh(p, t, q):
+    """ Function to compute relative humidity following 'Tetens' equation (T,P) ...'
+      [t]= temperature (assuming [[t],z,y,x] in [K])
+      [p] = pressure field (assuming in [hPa])
+      [q] = mixing ratio in [kgkg-1]
+    """
+    fname = 'var_rh'
+
+    data1 = 10.*0.6112*np.exp(17.67*(t-273.16)/(t-29.65))
+    data2 = 0.622*data1/(0.01*p-(1.-0.622)*data1)
+
+    rh = q/data2
+
+    return rh
+
 def var_td(t, p, qv):
     """ Function to compute dew-point air temperature from temperature and pressure values
@@ -1329 +1347 @@
         self.units = None
 
-        if Cdiag == 'td':
+        if Cdiag == 'hur':
+            """ Computing dew-point temperature 
+            """
+            vn = 'hur'
+            CF3Dvars = ['ta', 'plev', 'hus']
+            for v3D in CF3Dvars:
+                if not vars3D.has_key(v3D):
+                    print gen.errormsg
+                    print '  ' + fname + ": missing variable '" + v3D +              \
+                      "' attribution to compute '" + vn + "' !!"
+                    print '  Equivalence of 3D variables provided _______'
+                    gen.printing_dictionary(vars3D)
+                    quit(-1)
+                if not self.incvars.has_key(vars3D[v3D]):
+                    print gen.errormsg
+                    print '  ' + fname + ": missing variable '" + vars3D[v3D] +      \
+                      "' in input file to compute '" + vn + "' !!"
+                    print '  available variables:', self.incvars.keys()
+                    print '  looking for variables _______'  
+                    gen.printing_dictionary(vars3D)
+                    quit(-1)
+
+            ta = ncobj.variables[vars3D['ta']][:]
+            p = ncobj.variables[vars3D['plev']][:]
+            hus = ncobj.variables[vars3D['hus']][:]
+
+            if len(ta.shape) != len(p.shape):
+                p = gen.fill_Narray(p, ta*0., filldim=[0,2,3])
+
+            self.values = var_rh(ta, p, hus)
+            self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
+              dictdims['lon']]
+            self.units = '1'
+
+        elif Cdiag == 'td':
             """ Computing dew-point temperature 
             """
@@ -1353 +1405 @@
             ta = ncobj.variables[vars3D['ta']][:]
             p = ncobj.variables[vars3D['plev']][:]
-            hur = ncobj.variables[vars3D['hus']][:]
+            hus = ncobj.variables[vars3D['hus']][:]
 
             if len(ta.shape) != len(p.shape):
                 p = gen.fill_Narray(p, ta*0., filldim=[0,2,3])
 
-            self.values = var_td(ta, p, hur)
+            self.values = var_td(ta, p, hus)
             self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
               dictdims['lon']]