Changeset 2839 in lmdz_wrf for trunk/tools

Timestamp:

May 15, 2020, 11:07:24 PM (5 years ago)

Author:

lfita

Message:

Fixing humidity calculations

File:

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

trunk/tools/diag_tools.py

@@ -2274 +2274 @@
     #print 'es:', es[5,5,5,5], 'ws:', ws[5,5,5,5], 'hur:', hur[5,5,5,5]
 
+
+    # Complementary formula from where?
+    svp1=611.2
+    svp2=17.67
+    svp3=29.65
+    svpt0=273.15
+    eps = 0.622
+
+#    rh = 1.E2 * (p*q/(q*(1.-eps) + eps))/(svp1*exp(svp2*(t-svpt0)/(T-svp3)))
+    rh = 1.E2 * (p*q/(q*(1.-eps) + eps))/(svp1*np.exp(svp2*(t-svpt0)/(t-svp3)))
+    #print 'rh', rh
+
     return hur
 
@@ -2320 +2332 @@
             tdps = c*gamma/(b-gamma)
          hur = 100*expr[b(tdps/(c+tdps)-tas/(c+tas))]
-
+    >>> var_hur_tas_tds(10.+273.15, 8.+273.15, None, None)
+    0.8737387866043359
     """
     fname = 'compute_hur_tas_tds'
@@ -2336 +2349 @@
     hur = np.exp(b*(tdsC/(c+tdsC)-tasC/(c+tasC)))
                 
+    Xlv = 2.5e6
+    Rv = 461.5
+
+    rh = np.exp(Xlv/Rv*(1./tas - 1./tds)) * 1.E2
+    print 'rh', rh
+
     return hur, dnamesvar, dvnamesvar