Changeset 1890 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Apr 22, 2016 8:52:11 AM (8 years ago)

Author:

hoffmann

Message:

improvements for the consideration of aerosols in the LCM

File:

• palm/trunk/SOURCE/lpm_init.f90 (modified) (4 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Initialization of aerosol equilibrium radius not possible in supersaturated 
+! environments. Therefore, a maximum supersaturation of -1 % is assumed during
+! initialization.
 ! 
 ! Former revisions:
@@ -811 +813 @@
     REAL(wp)  :: weight_sum         !< sum of all weighting factors
 
-    INTEGER(iwp), DIMENSION(:,:,:), INTENT(IN) ::  local_start !<
+    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg), INTENT(IN) ::  local_start !<
 
     INTEGER(iwp)  :: n              !<
@@ -980 +982 @@
 
 !
-!--          The formula is only valid for subsaturated environments. In (super-)
-!--          saturated air, the inital radius is used.
-             IF ( e_a / e_s < 1.0_wp )  THEN
+!--          The formula is only valid for subsaturated environments. For 
+!--          supersaturations higher than -1 %, the supersaturation is set to -1%.
+             IF ( e_a / e_s < 0.99_wp )  THEN
 
                 DO  n = local_start(kp,jp,ip), number_of_particles  !only new particles 
@@ -995 +997 @@
                 ENDDO
 
+             ELSE
+
+                DO  n = local_start(kp,jp,ip), number_of_particles  !only new particles 
+
+                   bfactor             = vanthoff * molecular_weight_of_water *    &
+                                         rho_s * particles(n)%rvar2**3 /           &
+                                         ( molecular_weight_of_solute * rho_l )
+                   particles(n)%radius = particles(n)%rvar2 * ( bfactor /          &
+                                         particles(n)%rvar2**3 )**(1.0_wp/3.0_wp) *&
+                                         0.01_wp**(-1.0_wp/3.0_wp)
+
+                ENDDO
+
              ENDIF