modd_chs_aerosol.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!!     ######################
       MODULE modd_chs_aerosol
!!     ######################
!!
!!     PURPOSE
!!     -------
!!
!!     declaration of variables and types for the aerosol system in the
!!     surface meteo-france scheme
!!
!!     METHOD
!!     ------
!!
!!
!!     REFERENCE
!!     ---------
!!     none
!!
!!
!!     AUTHOR
!!     ------
!!     P. Tulet (CNRM)
!!
!!
!!     MODIFICATIONS
!!     -------------
!!
!!--------------------------------------------------------------------
!!     DECLARATIONS
!!     ------------
IMPLICIT NONE
!
! aerosol mode parameters
LOGICAL      :: lch_aero_flux     = .false. ! switch to active pronostic aerosols
!
LOGICAL      :: lco2pm            = .false. ! switch to active primary emission derived from CO 

LOGICAL      :: lvarsigi  = .false.   ! switch to active pronostic dispersion for I mode
LOGICAL      :: lvarsigj  = .false.   ! switch to active pronostic dispersion for J mode

!
INTEGER, PARAMETER         :: jpmode=2      ! number of modes
INTEGER, PARAMETER         :: jpin=jpmode*3 ! number of differential equations
INTEGER, SAVE, DIMENSION(JPMODE) :: nm0,nm3,nm6   ! index of the moments in arrays
!
!* indices of Aerosol chemical parameters
!
INTEGER, PARAMETER :: nsp=4        ! number of chemical species
                                   ! for ARES or isorropia NSP=4 these are
INTEGER, PARAMETER :: jp_aer_so4 = 1
INTEGER, PARAMETER :: jp_aer_no3 = 2
INTEGER, PARAMETER :: jp_aer_nh3 = 3
INTEGER, PARAMETER :: jp_aer_h2o = 4
!
INTEGER, PARAMETER :: jp_aer_so4g = jp_aer_so4
INTEGER, PARAMETER :: jp_aer_no3g = jp_aer_no3
INTEGER, PARAMETER :: jp_aer_nh3g = jp_aer_nh3
!
INTEGER, PARAMETER :: ncarb=3     ! number of chemically inert species
                                  ! (like black carbon)
INTEGER, PARAMETER :: jp_aer_oc = 5
INTEGER, PARAMETER :: jp_aer_bc = 6
INTEGER, PARAMETER :: jp_aer_dst = 7

INTEGER            :: nsoa = 10    ! number of condensable species that may form
                                   ! secondary aerosols
INTEGER            :: nm6_aer = 2  ! number of condensable species that may form
                                   ! secondary aerosols
INTEGER            :: jp_aer_soa1 = 8 
INTEGER            :: jp_aer_soa2 = 9
INTEGER            :: jp_aer_soa3 = 10
INTEGER            :: jp_aer_soa4 = 11
INTEGER            :: jp_aer_soa5 = 12
INTEGER            :: jp_aer_soa6 = 13
INTEGER            :: jp_aer_soa7 = 14
INTEGER            :: jp_aer_soa8 = 15
INTEGER            :: jp_aer_soa9 = 16 
INTEGER            :: jp_aer_soa10 = 17

 CHARACTER(LEN=32),DIMENSION(:), ALLOCATABLE :: caeronames

INTEGER            :: jp_ch_so4i = 1  
INTEGER            :: jp_ch_so4j = 2  
INTEGER            :: jp_ch_no3i = 3  
INTEGER            :: jp_ch_no3j = 4  
INTEGER            :: jp_ch_nh3i = 5  
INTEGER            :: jp_ch_nh3j = 6  
INTEGER            :: jp_ch_h2oi = 7  
INTEGER            :: jp_ch_h2oj = 8  
INTEGER            :: jp_ch_oci  = 9  
INTEGER            :: jp_ch_ocj  = 10  
INTEGER            :: jp_ch_bci  = 11  
INTEGER            :: jp_ch_bcj  = 12 
INTEGER            :: jp_ch_dsti  = 13 
INTEGER            :: jp_ch_dstj  = 14 
INTEGER            :: jp_ch_soa1i  = 15  
INTEGER            :: jp_ch_soa1j  = 16
INTEGER            :: jp_ch_soa2i  = 17
INTEGER            :: jp_ch_soa2j  = 18
INTEGER            :: jp_ch_soa3i  = 19
INTEGER            :: jp_ch_soa3j  = 20  
INTEGER            :: jp_ch_soa4i  = 21  
INTEGER            :: jp_ch_soa4j  = 22 
INTEGER            :: jp_ch_soa5i  = 23  
INTEGER            :: jp_ch_soa5j  = 24  
INTEGER            :: jp_ch_soa6i  = 25  
INTEGER            :: jp_ch_soa6j  = 26  
INTEGER            :: jp_ch_soa7i  = 27  
INTEGER            :: jp_ch_soa7j  = 28  
INTEGER            :: jp_ch_soa8i  = 29  
INTEGER            :: jp_ch_soa8j  = 30  
INTEGER            :: jp_ch_soa9i  = 31  
INTEGER            :: jp_ch_soa9j  = 32  
INTEGER            :: jp_ch_soa10i = 33  
INTEGER            :: jp_ch_soa10j = 34  
INTEGER            :: jp_ch_m0i = 35  
INTEGER            :: jp_ch_m0j = 36  
INTEGER            :: jp_ch_m6i = 37  
INTEGER            :: jp_ch_m6j = 38  

!INTEGER, PARAMETER :: JPNN=NSP+NSOA+NCARB
                                   
! Index for gas species which interact with aerosols
INTEGER :: jp_ch_hno3,  jp_ch_h2so4, jp_ch_nh3, jp_ch_o3, jp_ch_co,  &
             jp_ch_urg1, jp_ch_urg2, jp_ch_rpg2, jp_ch_rp18, jp_ch_ur26,&
             jp_ch_rpg3, jp_ch_urg4, jp_ch_ur8, jp_ch_ur17, jp_ch_ur7, jp_ch_urg6, &
             jp_ch_arac, jp_ch_urg7, jp_ch_rpg7, jp_ch_rpr7, jp_ch_urg8, jp_ch_ur19, &
             jp_ch_urg9, jp_ch_urg10, jp_ch_pan8, jp_ch_ur22, jp_ch_rpr4, jp_ch_ap7, &
             jp_ch_rpr3, jp_ch_ur21, jp_ch_ur28, jp_ch_ur29,  jp_ch_ur30, &
             jp_ch_rpr9, jp_ch_rp12, jp_ch_ur3, jp_ch_ur23, jp_ch_ur31, jp_ch_ap1, &
             jp_ch_ap6, jp_ch_adac, jp_ch_ur2, jp_ch_ur14, jp_ch_ur27, jp_ch_rp14, &
             jp_ch_rp19, jp_ch_ur11, jp_ch_ur15, jp_ch_ap10, jp_ch_ur20, jp_ch_ur34, &
             jp_ch_ap11, jp_ch_ap12, jp_ch_ur5, jp_ch_ur6, jp_ch_ap8  

INTEGER :: jp_ch_h2o2,  jp_ch_so2, jp_ch_so42m
!



! Molar mass of each aerosols parents (in kg/mol)
REAL, PARAMETER :: xhno3=63.01287
REAL, PARAMETER :: xnh3 =17.03061
REAL, PARAMETER :: xurg1=88.
REAL, PARAMETER :: xurg2=1.76981e+02
REAL, PARAMETER :: xrpg2=1.68000e+02
REAL, PARAMETER :: xrp18=1.84000e+02
REAL, PARAMETER :: xrpg3=1.53772e+02
REAL, PARAMETER :: xurg4=1.95867e+02
REAL, PARAMETER :: xur17=1.72000e+02
REAL, PARAMETER :: xrpr3=1.86000e+02
REAL, PARAMETER :: xap7 =2.33000e+02
REAL, PARAMETER :: xurg6=1.89153e+02
REAL, PARAMETER :: xur22=2.12000e+02
REAL, PARAMETER :: xurg7=1.56781e+02
REAL, PARAMETER :: xrpr4=1.67000e+02
REAL, PARAMETER :: xrpr7=1.50000e+02
REAL, PARAMETER :: xrpg7=1.96059e+02
REAL, PARAMETER :: xurg8=1.73777e+02
REAL, PARAMETER :: xurg9=2.61676e+02
REAL, PARAMETER :: xur26=1.68000e+02
REAL, PARAMETER :: xurg10=2.14834e+02
REAL, PARAMETER :: xur7=1.68000e+02
REAL, PARAMETER :: xur8=1.84000e+02
REAL, PARAMETER :: xpan8=2.63000e+02
REAL, PARAMETER :: xarac=1.32000e+02
REAL, PARAMETER :: xur19=1.70000e+02



!
REAL         :: xemisradiusi  = 0.036   ! mean radius of primary aerosol
                                       ! emission for I mode
REAL         :: xemisradiusj  = 0.385    ! mean radius of primary aerosol
                                       ! emission for J mode
REAL         :: xemissigi     = 1.86   ! dispersion of primary aerosol
                                       ! emission for I mode
REAL         :: xemissigj     = 1.29   ! dispersion of primary aerosol
                                       ! emission for J mode
 CHARACTER*4  :: crgunit   = 'NUMB'    ! type of log-normal geometric mean radius given
!                                     ! in nameliste (mass on number)



!----------------------------------------------------------------------------
!
!*  constants
!
REAL, PARAMETER :: xpboltz=1.380658e-23    ! Boltzmann constant (J/K)
REAL, PARAMETER :: xavogadro=6.0221367e+23 ! Avogadro constant
REAL, PARAMETER :: xmd    = 28.9644e-3     ! Air mass molarity

!
END MODULE modd_chs_aerosol