Ignore:
Timestamp:
Aug 29, 2017 2:10:28 PM (4 years ago)
Author:
schwenkel
Message:

improved aerosol initialization for bulk microphysics

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/modules.f90

    r2372 r2375  
    2525! -----------------
    2626! $Id$
     27! Moved mass_of_solute, molecular_weight_of_solute, molecular_weight_of_water,
     28! vanthoff back from particle attributes because they can now also be used in
     29! bulk microphysics.
     30! Added aerosol_bulk, aerosol_nacl, aerosol_c3h4o4, aerosol_nh4no3
     31!
     32! 2372 2017-08-25 12:37:32Z sward
    2733! y_shift namelist parameter added
    2834!
     
    878884    USE kinds
    879885
    880     REAL(wp) ::  cp = 1005.0_wp         !< heat capacity of dry air (J kg-1 K-1)
    881     REAL(wp) ::  l_v = 2.5E+06_wp       !< latent heat of vaporization (J kg-1)
    882     REAL(wp) ::  l_d_cp                 !< l_v / cp
    883     REAL(wp) ::  l_d_r                  !< l_v / r_d
    884     REAL(wp) ::  l_d_rv                 !< l_v / r_v   
    885     REAL(wp) ::  rho_l = 1.0E3_wp       !< density of water (kg m-3)
    886     REAL(wp) ::  r_d = 287.0_wp         !< sp. gas const. dry air (J kg-1 K-1)
    887     REAL(wp) ::  r_v = 461.51_wp        !< sp. gas const. water vapor (J kg-1 K-1)
     886    REAL(wp) ::  cp = 1005.0_wp                            !< heat capacity of dry air (J kg-1 K-1)
     887    REAL(wp) ::  l_v = 2.5E+06_wp                          !< latent heat of vaporization (J kg-1)
     888    REAL(wp) ::  l_d_cp                                    !< l_v / cp
     889    REAL(wp) ::  l_d_r                                     !< l_v / r_d
     890    REAL(wp) ::  l_d_rv                                    !< l_v / r_v   
     891    REAL(wp) ::  molecular_weight_of_solute = 0.05844_wp   !< mol. m. NaCl (kg mol-1)
     892    REAL(wp) ::  molecular_weight_of_water = 0.01801528_wp !< mol. m. H2O (kg mol-1)
     893    REAL(wp) ::  rho_l = 1.0E3_wp                          !< density of water (kg m-3)
     894    REAL(wp) ::  rho_s = 2165.0_wp                         !< density of NaCl (kg m-3)
     895    REAL(wp) ::  r_d = 287.0_wp                            !< sp. gas const. dry air (J kg-1 K-1)
     896    REAL(wp) ::  r_v = 461.51_wp                           !< sp. gas const. water vapor (J kg-1 K-1)
     897    REAL(wp) ::  vanthoff = 2.0_wp                         !< van't Hoff factor for NaCl
     898
    888899
    889900
     
    974985    CHARACTER (LEN=20)   ::  bc_uv_b = 'dirichlet'                        !< namelist parameter
    975986    CHARACTER (LEN=20)   ::  bc_uv_t = 'dirichlet'                        !< namelist parameter
     987    CHARACTER (LEN=20)   ::  aerosol_bulk = 'nacl'                        !< namelist parameter
    976988    CHARACTER (LEN=20)   ::  cloud_scheme = 'saturation_adjust'           !< namelist parameter
    977989    CHARACTER (LEN=20)   ::  coupling_mode = 'uncoupled'                  !< coupling mode for atmosphere-ocean coupling 
     
    11191131    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  mask_k_global  !< global grid index of masked output point on z-dimension
    11201132
     1133    LOGICAL ::  aerosol_nacl =.TRUE.                         !< nacl aerosol for bulk scheme
     1134    LOGICAL ::  aerosol_c3h4o4 =.FALSE.                      !< malonic acid aerosol for bulk scheme
     1135    LOGICAL ::  aerosol_nh4no3 =.FALSE.                      !< malonic acid aerosol for bulk scheme
    11211136    LOGICAL ::  bc_lr_cyc =.TRUE.                            !< left-right boundary condition cyclic?
    11221137    LOGICAL ::  bc_lr_dirrad = .FALSE.                       !< left-right boundary condition dirichlet/radiation?
Note: See TracChangeset for help on using the changeset viewer.