Changeset 2522 for palm/trunk/SOURCE/microphysics_mod.f90

Timestamp:

Oct 5, 2017 2:20:37 PM (7 years ago)

Author:

schwenkel

Message:

minor bugfix

File:

• palm/trunk/SOURCE/microphysics_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/microphysics_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor bugfix
+! 
+! 2375 2017-08-29 14:10:28Z schwenkel
 ! Improved aerosol initilization and some minor bugfixes 
 ! for droplet sedimenation
@@ -609 +612 @@
 !--             (see: Bott + Trautmann, 2002, Atm. Res., 64)
                 k_act  = 0.7_wp
-
-                IF ( sat >= 0.0 .AND. .NOT. curvature_solution_effects_bulk ) THEN
+                activ  = 0.0_wp
+
+
+                IF ( sat > 0.0 .AND. .NOT. curvature_solution_effects_bulk ) THEN
 !
 !--                Compute the number of activated Aerosols
@@ -627 +632 @@
                       ( 1.0_wp, ( sat / sat_max )**k_act) - nc(k,j,i) ) ) /    &
                        dt_micro
-                ELSEIF ( sat >= 0.0 .AND. curvature_solution_effects_bulk ) THEN
+                ELSEIF ( sat > 0.0 .AND. curvature_solution_effects_bulk ) THEN
 !
 !--                Curvature effect (afactor) with surface tension
@@ -850 +855 @@
                 nc_auto = MERGE( nc(k,j,i), nc_const, microphysics_morrison )
 
-                IF ( qc(k,j,i) > eps_sb )  THEN
+                IF ( qc(k,j,i) > eps_sb  .AND.  nc_auto > eps_mr )  THEN
 
                    k_au = k_cc / ( 20.0_wp * x0 )
@@ -856 +861 @@
 !--                Intern time scale of coagulation (Seifert and Beheng, 2006):
 !--                (1.0_wp - qc(k,j,i) / ( qc(k,j,i) + qr(k,j,i) ))
-                   tau_cloud = 1.0_wp - qc(k,j,i) / ( qr(k,j,i) + qc(k,j,i) )
+                   tau_cloud = MAX( 1.0_wp - qc(k,j,i) / ( qr(k,j,i) +         &
+                                    qc(k,j,i) ), 0.0_wp )
 !
 !--                Universal function for autoconversion process
@@ -868 +874 @@
 !
 !--                Mean weight of cloud droplets:
-                   xc = hyrho(k) * qc(k,j,i) / nc_auto
+                   xc = MAX( hyrho(k) * qc(k,j,i) / nc_auto, xcmin) 
 !
 !--                Parameterized turbulence effects on autoconversion (Seifert,
@@ -2041 +2047 @@
 !--       (see: Bott + Trautmann, 2002, Atm. Res., 64)
           k_act  = 0.7_wp
+          activ  = 0.0_wp
 
           IF ( sat >= 0.0 .AND. .NOT. curvature_solution_effects_bulk )  THEN
@@ -2262 +2269 @@
           nc_auto = MERGE ( nc_1d(k), nc_const, microphysics_morrison )
 
-          IF ( qc_1d(k) > eps_sb )  THEN
+          IF ( qc_1d(k) > eps_sb  .AND.  nc_auto > eps_mr )  THEN
 
              k_au = k_cc / ( 20.0_wp * x0 )
@@ -2268 +2275 @@
 !--          Intern time scale of coagulation (Seifert and Beheng, 2006):
 !--          (1.0_wp - qc(k,j,i) / ( qc(k,j,i) + qr_1d(k) ))
-             tau_cloud = 1.0_wp - qc_1d(k) / ( qr_1d(k) + qc_1d(k) )
+             tau_cloud = MAX( 1.0_wp - qc_1d(k) / ( qr_1d(k) + qc_1d(k) ),     & 
+                              0.0_wp )
 !
 !--          Universal function for autoconversion process