Changeset 4506 for palm/trunk/SOURCE

Timestamp:

Apr 21, 2020 10:57:45 AM (5 years ago)

Author:

schwenkel

Message:

Use magnus formula for liquid water temperature

File:

• palm/trunk/SOURCE/bulk_cloud_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Use correct magnus formula for liquid water temperature
+! 
+! 4502 2020-04-17 16:14:16Z schwenkel
 ! Implementation of ice microphysics
 ! 
@@ -134 +137 @@
 
     USE basic_constants_and_equations_mod,                                     &
-        ONLY:  c_p, g, lv_d_cp, lv_d_rd, l_v, magnus, magnus_ice,              &
+        ONLY:  c_p, g, lv_d_cp, lv_d_rd, l_v, magnus, magnus_ice, magnus_tl,   &
                molecular_weight_of_solute,                                     &
                molecular_weight_of_water, pi, rho_l, rho_s, r_d, r_v, vanthoff,&
@@ -1280 +1283 @@
              d_exner = exner_function_invers(hyp)
              exner = 1.0_wp / exner_function_invers(hyp)
-             hyrho  = ideal_gas_law_rho_pt(hyp, pt(:, nys, nxl) )
+             hyrho  = ideal_gas_law_rho_pt(hyp, pt_init)
 !
 !--          Compute reference density
@@ -1396 +1399 @@
              d_exner = exner_function_invers(hyp)
              exner = 1.0_wp / exner_function_invers(hyp)
-             hyrho  = ideal_gas_law_rho_pt(hyp, pt(:, nys, nxl) )
+             hyrho  = ideal_gas_law_rho_pt(hyp, pt_init)
 !
 !--          Compute reference density
@@ -6434 +6437 @@
 !
 !--    Calculate water vapor saturation pressure
-       e_s = magnus( t_l )
+       e_s = magnus_tl( t_l )
 !
 !--    Computation of saturation mixing ratio: