Changeset 1849 for palm/trunk/SOURCE/init_cloud_physics.f90

Timestamp:

Apr 8, 2016 11:33:18 AM (8 years ago)

Author:

hoffmann

Message:

lpm_droplet_condensation improved, microphysics partially modularized

File:

• palm/trunk/SOURCE/init_cloud_physics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_cloud_physics.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-!
+! bfactor removed
 ! 
 ! Former revisions:
@@ -85 +85 @@
         
     USE cloud_parameters,                                                      &
-        ONLY:  bfactor, cp, c_sedimentation, dpirho_l, dt_precipitation,       &
-               hyrho, k_st, l_d_cp, l_d_r, l_d_rv, l_v, mass_of_solute,  &
-               molecular_weight_of_solute, molecular_weight_of_water, pirho_l, &
-               pt_d_t, rho_l, r_d, r_v, sed_qc_const, schmidt, schmidt_p_1d3,  &
-               sigma_gc, t_d_pt, vanthoff, w_precipitation
-        
-    USE constants,                                                             &
-        ONLY:  pi
-        
+        ONLY:  cp, hyrho, l_d_cp, l_d_r, l_d_rv, l_v, pt_d_t, rho_l, r_d, r_v, &
+               t_d_pt
+                
     USE control_parameters,                                                    &
-        ONLY:  g, message_string, microphysics_seifert, pt_surface,            &
-               rho_surface, surface_pressure
+        ONLY:  g, message_string, pt_surface, rho_surface, surface_pressure
     
     USE indices,                                                               &
@@ -118 +111 @@
     l_d_rv = l_v / r_v
 
-    schmidt_p_1d3 = schmidt**( 1.0_wp / 3.0_wp )
-
-    pirho_l  = pi * rho_l / 6.0_wp
-    dpirho_l = 1.0_wp / pirho_l
-!
-!-- constant for the sedimentation of cloud water (2-moment cloud physics)
-    sed_qc_const = k_st * ( 3.0_wp / ( 4.0_wp * pi * rho_l )                   &
-                          )**( 2.0_wp / 3.0_wp ) *                             &
-                   EXP( 5.0_wp * LOG( sigma_gc )**2 )
-!
-!-- Calculate timestep according to precipitation
-    IF ( microphysics_seifert )  THEN
-       dt_precipitation = c_sedimentation * MINVAL( dzu(nzb+2:nzt) ) /         &
-                          w_precipitation
-    ENDIF
-!
-!-- Calculate factor used in equation for droplet growth by condensation
-    bfactor = 3.0_wp * vanthoff * mass_of_solute * molecular_weight_of_water   &
-              / ( 4.0_wp * pi * rho_l * molecular_weight_of_solute )
 !
 !-- Calculate: