Changeset 1720 in lmdz_wrf

Timestamp:

Dec 14, 2017, 2:12:51 AM (7 years ago)

Author:

lfita

Message:

Adding:

• `var_td_Uhus': Function to compute dew-point air temperature from temperature and pressure values using hus

Location:

trunk/tools

Files:

• diag_tools.py (modified) (4 diffs)
• get_stations.py (modified) (1 diff)

trunk/tools/diag_tools.py

@@ -40 +40 @@
 # var_mslp: Fcuntion to compute mean sea-level pressure
 # var_td: Function to compute dew-point air temperature from temperature and pressure values
+# var_td_Uhus: Function to compute dew-point air temperature from temperature and pressure values using hus
 # var_virtualTemp: This function returns virtual temperature in K, 
 # var_WRFtime: Function to copmute CFtimes from WRFtime variable
@@ -1358 +1359 @@
     return td
 
+def var_td_Uhus(t, p, hus):
+    """ Function to compute dew-point air temperature from temperature and pressure values using hus
+      t= temperature [K]
+      hus= specific humidity [1]
+      formula:
+          temp = theta*(p/p0)**(R/Cp)
+
+    """
+    fname = 'compute_td'
+
+    tk = (t)*(p/fdef.module_definitions.p0ref)**(fdef.module_definitions.rcp)
+    data1 = 10.*0.6112*np.exp(17.67*(tk-273.16)/(tk-29.65))
+    data2 = 0.622*data1/(0.01*p-(1.-0.622)*data1)
+
+    qv = hus/(1.+hus)
+
+    rh = qv/data2
+                
+    pa = rh * data1
+    td = 257.44*np.log(pa/6.1121)/(18.678-np.log(pa/6.1121))
+
+    return td
+
 def var_wd(u, v):
     """ Function to compute the wind direction
@@ -1474 +1498 @@
             """
             vn = 'td'
-            CF3Dvars = ['ta', 'plev', 'hus']
+            CF3Dvars = ['ta', 'plev', 'q']
             for v3D in CF3Dvars:
                 if not vars3D.has_key(v3D):
@@ -1494 +1518 @@
             ta = ncobj.variables[vars3D['ta']][:]
             p = ncobj.variables[vars3D['plev']][:]
-            hus = ncobj.variables[vars3D['hus']][:]
+            q = ncobj.variables[vars3D['q']][:]
 
             if len(ta.shape) != len(p.shape):
                 p = gen.fill_Narray(p, ta*0., filldim=[0,2,3])
 
-            self.values = var_td(ta, p, hus)
+            self.values = var_td(ta, p, q)
+            self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
+              dictdims['lon']]
+            self.units = 'K'
+
+        elif Cdiag == 'td_Uhus':
+            """ Computing dew-point temperature 
+            """
+            vn = 'td'
+            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_td_Uhus(ta, p, hus)
             self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
               dictdims['lon']]

trunk/tools/get_stations.py

@@ -21 +21 @@
 
 # Gneric variables prepared to be computed
-Cvars = ['C_hur', 'C_hur_Uhus', 'C_td', 'C_wd', 'C_ws']
+Cvars = ['C_hur', 'C_hur_Uhus', 'C_td', 'C_td_hus', 'C_wd', 'C_ws']
 
 # WRF specific diagnostics