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

palm/trunk/SOURCE/modules.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Moved mass_of_solute, molecular_weight_of_solute, molecular_weight_of_water,
+! vanthoff back from particle attributes because they can now also be used in
+! bulk microphysics. 
+! Added aerosol_bulk, aerosol_nacl, aerosol_c3h4o4, aerosol_nh4no3
+! 
+! 2372 2017-08-25 12:37:32Z sward
 ! y_shift namelist parameter added
 ! 
@@ -878 +884 @@
     USE kinds
 
-    REAL(wp) ::  cp = 1005.0_wp         !< heat capacity of dry air (J kg-1 K-1)
-    REAL(wp) ::  l_v = 2.5E+06_wp       !< latent heat of vaporization (J kg-1)
-    REAL(wp) ::  l_d_cp                 !< l_v / cp
-    REAL(wp) ::  l_d_r                  !< l_v / r_d
-    REAL(wp) ::  l_d_rv                 !< l_v / r_v    
-    REAL(wp) ::  rho_l = 1.0E3_wp       !< density of water (kg m-3)
-    REAL(wp) ::  r_d = 287.0_wp         !< sp. gas const. dry air (J kg-1 K-1)
-    REAL(wp) ::  r_v = 461.51_wp        !< sp. gas const. water vapor (J kg-1 K-1)
+    REAL(wp) ::  cp = 1005.0_wp                            !< heat capacity of dry air (J kg-1 K-1)
+    REAL(wp) ::  l_v = 2.5E+06_wp                          !< latent heat of vaporization (J kg-1)
+    REAL(wp) ::  l_d_cp                                    !< l_v / cp
+    REAL(wp) ::  l_d_r                                     !< l_v / r_d
+    REAL(wp) ::  l_d_rv                                    !< l_v / r_v    
+    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) ::  rho_l = 1.0E3_wp                          !< density of water (kg m-3)
+    REAL(wp) ::  rho_s = 2165.0_wp                         !< density of NaCl (kg m-3)
+    REAL(wp) ::  r_d = 287.0_wp                            !< sp. gas const. dry air (J kg-1 K-1)
+    REAL(wp) ::  r_v = 461.51_wp                           !< sp. gas const. water vapor (J kg-1 K-1)
+    REAL(wp) ::  vanthoff = 2.0_wp                         !< van't Hoff factor for NaCl
+
 
 
@@ -974 +985 @@
     CHARACTER (LEN=20)   ::  bc_uv_b = 'dirichlet'                        !< namelist parameter
     CHARACTER (LEN=20)   ::  bc_uv_t = 'dirichlet'                        !< namelist parameter
+    CHARACTER (LEN=20)   ::  aerosol_bulk = 'nacl'                        !< namelist parameter
     CHARACTER (LEN=20)   ::  cloud_scheme = 'saturation_adjust'           !< namelist parameter
     CHARACTER (LEN=20)   ::  coupling_mode = 'uncoupled'                  !< coupling mode for atmosphere-ocean coupling  
@@ -1119 +1131 @@
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  mask_k_global  !< global grid index of masked output point on z-dimension
 
+    LOGICAL ::  aerosol_nacl =.TRUE.                         !< nacl aerosol for bulk scheme
+    LOGICAL ::  aerosol_c3h4o4 =.FALSE.                      !< malonic acid aerosol for bulk scheme
+    LOGICAL ::  aerosol_nh4no3 =.FALSE.                      !< malonic acid aerosol for bulk scheme
     LOGICAL ::  bc_lr_cyc =.TRUE.                            !< left-right boundary condition cyclic?
     LOGICAL ::  bc_lr_dirrad = .FALSE.                       !< left-right boundary condition dirichlet/radiation?