Changeset 2608 for palm/trunk/SOURCE/calc_liquid_water_content.f90

Timestamp:

Nov 13, 2017 2:04:26 PM (6 years ago)

Author:

schwenkel

Message:

Inital revision of diagnostic_quantities_mod allows unified calculation of magnus equation and saturion mixing ratio

File:

• palm/trunk/SOURCE/calc_liquid_water_content.f90 (modified) (4 diffs)

palm/trunk/SOURCE/calc_liquid_water_content.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Calculation of supersaturation in external module (diagnostic_quantities_mod).
+! Change: correct calculation of saturation specific humidity to saturation
+! mixing ratio (the factor of 0.378 vanishes).
+! 
+! 2292 2017-06-20 09:51:42Z schwenkel
 ! Implementation of new microphysic scheme: cloud_scheme = 'morrison' 
 ! includes two more prognostic equations for cloud drop concentration (nc)  
@@ -82 +87 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  hyp, pt, q, qc, ql, qr
-
-    USE cloud_parameters,                                                      &
-        ONLY:  l_d_cp, l_d_r, t_d_pt
+        ONLY:  q, qc, ql, qr
 
     USE control_parameters,                                                    &
         ONLY:  microphysics_morrison, microphysics_seifert
+
+    USE diagnostic_quantities_mod,                                             &
+        ONLY:  e_s, magnus, q_s, sat, supersaturation, t_l
 
     USE indices,                                                               &
@@ -104 +109 @@
     INTEGER(iwp) ::  k !<
 
-    REAL(wp) ::  alpha !<
-    REAL(wp) ::  e_s   !<
-    REAL(wp) ::  q_s   !<
-    REAL(wp) ::  t_l   !<
 
     DO  i = nxlg, nxrg
@@ -114 +115 @@
 
 !
-!--          Compute the liquid water temperature
-             t_l = t_d_pt(k) * pt(k,j,i)
-
-!
-!--          Compute saturation water vapor pressure at t_l
-             e_s = 610.78_wp * EXP( 17.269_wp * ( t_l - 273.16_wp ) /          &
-                                                ( t_l - 35.86_wp ) )
-
-!
-!--          Compute approximation of saturation humidity
-             q_s = 0.622_wp * e_s / ( hyp(k) - 0.378_wp * e_s )
-
-!
-!--          Correction factor
-             alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l )
-
-!
-!--          Correction of the approximated value
-!--          (see: Cuijpers + Duynkerke, 1993, JAS, 23)
-             q_s = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
+!--          Call calculation of supersaturation located 
+!--          in diagnostic_quantities_mod
+             CALL supersaturation ( i, j, k )
 
 !