Changeset 1871 for palm/trunk/SOURCE/lpm_droplet_condensation.f90

Timestamp:

Apr 15, 2016 11:46:09 AM (8 years ago)

Author:

hoffmann

Message:

initialization of aerosol spectra added

File:

• palm/trunk/SOURCE/lpm_droplet_condensation.f90 (modified) (15 diffs)

palm/trunk/SOURCE/lpm_droplet_condensation.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!
+! Initialization of aerosols added.
 ! 
 ! Former revisions:
 ! -----------------
 ! $Id$
-!
 ! 
 ! 1849 2016-04-08 11:33:18Z hoffmann
@@ -121 +120 @@
 
     USE particle_attributes,                                                   &
-        ONLY:  curvature_solution_effects, hall_kernel, mass_of_solute,        &
+        ONLY:  curvature_solution_effects, hall_kernel,                        &
                molecular_weight_of_solute, molecular_weight_of_water,          &
-               number_of_particles, particles, radius_classes,                 &
+               number_of_particles, particles, radius_classes, rho_s,          &
                use_kernel_tables, vanthoff, wang_kernel
 
@@ -140 +139 @@
     INTEGER(iwp) :: ros_count                  !<
  
-    INTEGER(iwp), PARAMETER ::  maxtry = 40    !<
+    INTEGER(iwp), PARAMETER ::  maxtry = 100   !<
 
     LOGICAL ::  repeat                         !<
@@ -195 +194 @@
     REAL(wp), PARAMETER ::  e3 = 0.0_wp                !<
     REAL(wp), PARAMETER ::  e4 = 125.0_wp / 108.0_wp   !<
+    REAL(wp), PARAMETER ::  eps_ros = 1.0E-3_wp        !< accuracy of Rosenbrock method
     REAL(wp), PARAMETER ::  gam = 0.5_wp               !<
     REAL(wp), PARAMETER ::  grow = 1.5_wp              !<
@@ -200 +200 @@
     REAL(wp), PARAMETER ::  pshrnk = -1.0_wp /3.0_wp   !<
     REAL(wp), PARAMETER ::  shrnk = 0.5_wp             !<
-    REAL(wp), PARAMETER ::  eps_ros = 1.0E-4_wp        !< accuracy of Rosenbrock method
 
 !
@@ -311 +310 @@
              afactor = 2.0_wp * sigma / ( rho_l * r_v * t_int )
 !
-!--          Solute effect (bfactor), mass of solute to be replaced by variable
-!--          aerosol radius
-             bfactor = 3.0_wp * vanthoff * mass_of_solute *                     &
-                       molecular_weight_of_water / ( 4.0_wp * pi * rho_l *      &
-                                                     molecular_weight_of_solute &
-                                                   )
+!--          Solute effect (bfactor)
+             bfactor = vanthoff * rho_s * particles(n)%rvar2**3 *              &
+                       molecular_weight_of_water /                             &
+                       ( rho_l * molecular_weight_of_solute )
 
              r_ros = particles(n)%radius
@@ -326 +323 @@
 
 !
-!--          Internal time step should not be > 1.0E-2 in case of evaporation
-!--          because larger values may lead to secondary solutions which are
-!--          physically unlikely
-             IF ( dt_ros_next > 1.0E-2_wp  .AND.  e_a / e_s < 1.0_wp )  THEN
+!--          Internal time step should not be > 1.0E-2 and < 1.0E-6 
+             IF ( dt_ros_next > 1.0E-2_wp )  THEN
                 dt_ros_next = 1.0E-3_wp
+             ELSEIF ( dt_ros_next < 1.0E-6_wp )  THEN
+                dt_ros_next = 1.0E-6_wp
              ENDIF
+
 !
 !--          If calculation of Rosenbrock method is repeated due to unreasonalble
@@ -337 +335 @@
 !--          reduced
              IF ( ros_count > 1 )  THEN
-                dt_ros_next = dt_ros_next - ( 0.2_wp * dt_ros_next )
+                dt_ros_next = dt_ros_next * 0.1_wp
              ELSEIF ( ros_count > 5 )  THEN
 !
@@ -349 +347 @@
 !--          Internal time step must not be larger than PALM time step
              dt_ros = MIN( dt_ros_next, dt_3d )
+
 !
 !--          Integrate growth equation in time unless PALM time step is reached
@@ -364 +363 @@
                 drdt_ini       = drdt
                 dt_ros_sum_ini = dt_ros_sum
-                r_ros_ini      = r_ros
+                r_ros_ini      = MAX( r_ros, particles(n)%rvar2 )
 
 !
@@ -378 +377 @@
 
                    IF ( jtry == maxtry+1 )  THEN
-                      message_string = 'maxtry > 40 in Rosenbrock method'
+                      message_string = 'maxtry > 100 in Rosenbrock method'
                       CALL message( 'lpm_droplet_condensation', 'PA0347', 2,   &
                                     2, 0, 6, 0 )
@@ -385 +384 @@
                    aa    = 1.0_wp / ( gam * dt_ros ) - d2rdtdr
                    g1    = drdt_ini / aa
-                   r_ros = r_ros_ini + a21 * g1
+                   r_ros = MAX( r_ros_ini + a21 * g1, particles(n)%rvar2 )
                    drdt  = ddenom * ventilation_effect(n) * ( e_a / e_s - 1.0_wp - &
                                                               afactor / r_ros +    &
@@ -393 +392 @@
                    g2    = ( drdt + c21 * g1 / dt_ros )&
                            / aa
-                   r_ros = r_ros_ini + a31 * g1 + a32 * g2
+                   r_ros = MAX( r_ros_ini + a31 * g1 + a32 * g2, particles(n)%rvar2 )
                    drdt  = ddenom * ventilation_effect(n) * ( e_a / e_s - 1.0_wp - &
                                                               afactor / r_ros +    &
@@ -403 +402 @@
                    g4    = ( drdt +  &
                              ( c41 * g1 + c42 * g2 + c43 * g3 ) / dt_ros ) / aa
-                   r_ros = r_ros_ini + b1 * g1 + b2 * g2 + b3 * g3 + b4 * g4
+                   r_ros = MAX( r_ros_ini + b1 * g1 + b2 * g2 + b3 * g3 +      &
+                                b4 * g4, particles(n)%rvar2 )
 
                    dt_ros_sum = dt_ros_sum_ini + dt_ros
@@ -451 +451 @@
              particles(n)%rvar1 = dt_ros_next
 
-             new_r(n) = r_ros
 !
 !--          Radius should not fall below 1E-8 because Rosenbrock method may
 !--          lead to errors otherwise
-             new_r(n) = MAX( new_r(n), 1.0E-8_wp )
+             new_r(n) = MAX( r_ros, particles(n)%rvar2 )
 !
 !--          Check if calculated droplet radius change is reasonable since in