Changeset 1710 in lmdz_wrf

Timestamp:

Dec 12, 2017, 11:40:43 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `compute_WRFhur': Function to compute WRF relative humidity following Teten's equation

File:

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

trunk/tools/diag_tools.py

@@ -19 +19 @@
 # compute_wds: Function to compute the wind direction
 # compute_wss: Function to compute the wind speed
+# compute_WRFhur: Function to compute WRF relative humidity following Teten's equation
 # compute_WRFta: Function to compute WRF air temperature
 # compute_WRFtd: Function to compute WRF dew-point air temperature
@@ -982 +983 @@
     return times
 
+def compute_WRFhur(t, p, qv, dimns, dimvns):
+    """ Function to compute WRF relative humidity following Teten's equation
+      t= orginal WRF temperature
+      p= original WRF pressure (P + PB)
+      formula:
+          temp = theta*(p/p0)**(R/Cp)
+
+    """
+    fname = 'compute_WRFtd'
+
+    tk = (t+300.)*(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)
+
+    rh = qv/data2
+                
+    dnamesvar = ['Time','bottom_top','south_north','west_east']
+    dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dimns,dimvns)
+
+    return rh, dnamesvar, dvnamesvar
+
 def compute_WRFua(u, v, sina, cosa, dimns, dimvns):
     """ Function to compute geographical rotated WRF 3D winds
@@ -1429 +1451 @@
         self.units = None
 
-        if Wdiag == 'p':
+        if Wdiag == 'hur':
+            """ Computing relative humidity
+            """
+            vn = 'hur'
+            CF3Dvars = ['ta', '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)
+
+            ta = ncobj.variables['T'][:]
+            p = ncobj.variables['P'][:] + ncobj.variables['PB'][:]
+            hur = ncobj.variables['QVAPOR'][:]
+    
+            vals, dims, vdims = compute_WRFhur(ta, p, hur, indims, invardims)
+            self.values = vals
+            self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \
+              dictdims['lon']]
+            self.units = '1'
+
+        elif Wdiag == 'p':
             """ Computing air pressure
             """
@@ -1478 +1524 @@
             ta = ncobj.variables['T'][:]
             p = ncobj.variables['P'][:] + ncobj.variables['PB'][:]
-            hur = ncobj.variables['QVAPOR'][:]
-    
-            vals, dims, vdims = compute_WRFtd(ta, p, hur, indims, invardims)
+            hus = ncobj.variables['QVAPOR'][:]
+    
+            vals, dims, vdims = compute_WRFtd(ta, p, hus, indims, invardims)
             self.values = vals
             self.dims = [dictdims['time'], dictdims['plev'], dictdims['lat'],        \