Changeset 3361 for palm/trunk/SOURCE/bulk_cloud_model_mod.f90

Timestamp:

Oct 16, 2018 8:39:37 PM (6 years ago)

Author:

knoop

Message:

Introduced global constant rd_d_rv=0.622

File:

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

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -186 +186 @@
                molecular_weight_of_water, pi, rho_l, rho_s, r_d, r_v, vanthoff,&
                exner_function, exner_function_invers, ideal_gas_law_rho,       &
-               ideal_gas_law_rho_pt, barometric_formula
+               ideal_gas_law_rho_pt, barometric_formula, rd_d_rv
 
     USE control_parameters,                                                    &
@@ -3308 +3308 @@
 !
 !--       Computation of saturation mixing ratio:
-          q_s   = 0.622_wp * e_s / ( hyp(k) -  e_s )
-          alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
+          q_s   = rd_d_rv * e_s / ( hyp(k) -  e_s )
+          alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
           q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) /                        &
                         ( 1.0_wp + alpha * q_s )
@@ -3423 +3423 @@
 !
 !--       Computation of saturation mixing ratio:
-          q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
-          alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
+          q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
+          alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
           q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) /               &
                         ( 1.0_wp + alpha * q_s )
@@ -3803 +3803 @@
 !
 !--          Computation of saturation mixing ratio:
-             q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
-             alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
+             q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
+             alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
              q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
 !
@@ -4247 +4247 @@
 !
 !--    Computation of saturation mixing ratio:
-       q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
+       q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
 !
 !--    Correction factor
-       alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
+       alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
 !
 !--    Correction of the approximated value