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

Timestamp:

Aug 29, 2017 2:10:28 PM (7 years ago)

Author:

schwenkel

Message:

improved aerosol initialization for bulk microphysics

File:

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

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Initialization of chemical aerosol composition
+! 
+! 2346 2017-08-09 16:39:17Z suehring
 ! Bugfix, correct determination of topography top index
 ! 
@@ -970 +973 @@
 
     USE cloud_parameters,                                                      &
-        ONLY: l_d_rv, rho_l, r_v
+        ONLY: l_d_rv, molecular_weight_of_solute,                              &
+              molecular_weight_of_water, rho_l, r_v, rho_s, vanthoff
 
     USE constants,                                                             &
@@ -978 +982 @@
 
     USE particle_attributes,                                                   &
-        ONLY: aero_type, aero_weight, log_sigma, molecular_weight_of_solute,   &
-              molecular_weight_of_water, na, rho_s, rm, vanthoff
+        ONLY: aero_species, aero_type, aero_weight, log_sigma, na, rm
 
     IMPLICIT NONE
@@ -1003 +1006 @@
     INTEGER(iwp)  :: kp             !<
 
+!
+!-- Set constants for different aerosol species
+    IF ( TRIM(aero_species) .EQ. 'nacl' ) THEN
+       molecular_weight_of_solute = 0.05844_wp 
+       rho_s                      = 2165.0_wp
+       vanthoff                   = 2.0_wp
+    ELSEIF ( TRIM(aero_species) .EQ. 'c3h4o4' ) THEN
+       molecular_weight_of_solute = 0.10406_wp 
+       rho_s                      = 1600.0_wp
+       vanthoff                   = 1.37_wp
+    ELSEIF ( TRIM(aero_species) .EQ. 'nh4o3' ) THEN
+       molecular_weight_of_solute = 0.08004_wp 
+       rho_s                      = 1720.0_wp
+       vanthoff                   = 2.31_wp
+    ELSE
+       WRITE( message_string, * ) 'unknown aerosol species ',   &
+                                'aero_species = "', TRIM( aero_species ), '"'
+       CALL message( 'lpm_init', 'PA0470', 1, 2, 0, 6, 0 )
+    ENDIF
 !
 !-- The following typical aerosol spectra are taken from Jaenicke (1993):