Changeset 2375 for palm/trunk/SOURCE/mod_particle_attributes.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/mod_particle_attributes.f90 (modified) (4 diffs)

palm/trunk/SOURCE/mod_particle_attributes.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
-!
+! 
+! 
 ! Former revisions:
 ! -----------------
 ! $Id$
+! molecular_weight_of_solute, molecular_weight_of_water, vanthoff removed and 
+! added in modules. Parameters are also used in bulk-microphysics.
+!
+! 2312 2017-07-14 20:26:51Z hoffmann
 ! Aerosol initialization improved.
 !
@@ -95 +99 @@
     USE kinds
 
-    CHARACTER(LEN=15) ::  aero_type = 'maritime'                   !< aerosol type
-    CHARACTER(LEN=15) ::  bc_par_lr = 'cyclic'                     !< left/right boundary condition
-    CHARACTER(LEN=15) ::  bc_par_ns = 'cyclic'                     !< north/south boundary condition
-    CHARACTER(LEN=15) ::  bc_par_b  = 'reflect'                    !< bottom boundary condition
-    CHARACTER(LEN=15) ::  bc_par_t  = 'absorb'                     !< top boundary condition
-    CHARACTER(LEN=15) ::  collision_kernel = 'none'                !< collision kernel
+    CHARACTER(LEN=15) ::  aero_species = 'nacl'                    !< aerosol species
+    CHARACTER(LEN=15) ::  aero_type    = 'maritime'                !< aerosol type
+    CHARACTER(LEN=15) ::  bc_par_lr    = 'cyclic'                  !< left/right boundary condition
+    CHARACTER(LEN=15) ::  bc_par_ns    = 'cyclic'                  !< north/south boundary condition
+    CHARACTER(LEN=15) ::  bc_par_b     = 'reflect'                 !< bottom boundary condition
+    CHARACTER(LEN=15) ::  bc_par_t     = 'absorb'                  !< top boundary condition
+    CHARACTER(LEN=15) ::  collision_kernel   = 'none'              !< collision kernel
     CHARACTER(LEN=5)  ::  splitting_function = 'gamma'             !< function for calculation critical weighting factor
-    CHARACTER(LEN=5)  ::  splitting_mode = 'const'                 !< splitting mode
+    CHARACTER(LEN=5)  ::  splitting_mode     = 'const'             !< splitting mode
 
     INTEGER(iwp) ::  deleted_particles = 0                        !< number of deleted particles per time step
@@ -172 +177 @@
     REAL(wp) ::  initial_weighting_factor = 1.0_wp         !< namelist parameter (see documentation)
     REAL(wp) ::  log_sigma(3) = 1.0_wp                     !< namelist parameter (see documentation)
-    REAL(wp) ::  molecular_weight_of_solute = 0.05844_wp   !< mol. m. NaCl (kg mol-1)
-    REAL(wp) ::  molecular_weight_of_water = 0.01801528_wp !< mol. m. H2O (kg mol-1)
     REAL(wp) ::  na(3) = 0.0_wp                            !< namelist parameter (see documentation)
     REAL(wp) ::  number_concentration = -1.0_wp            !< namelist parameter (see documentation)
@@ -179 +182 @@
     REAL(wp) ::  radius_merge = 1.0E-7_wp                  !< namelist parameter (see documentation)
     REAL(wp) ::  radius_split = 40.0E-6_wp                 !< namelist parameter (see documentation)
-    REAL(wp) ::  rho_s = 2165.0_wp                         !< density of NaCl (kg m-3)
     REAL(wp) ::  rm(3) = 1.0E-6_wp                         !< namelist parameter (see documentation)
     REAL(wp) ::  sgs_wf_part                               !< parameter for sgs
     REAL(wp) ::  time_prel = 0.0_wp                        !< time for particle release
     REAL(wp) ::  time_write_particle_data = 0.0_wp         !< write particle data at current time on file
-    REAL(wp) ::  vanthoff = 2.0_wp                         !< van't Hoff factor for NaCl
     REAL(wp) ::  weight_factor_merge = -1.0_wp             !< namelist parameter (see documentation)
     REAL(wp) ::  weight_factor_split = -1.0_wp             !< namelist parameter (see documentation)