SURFEX v7.3
General documentation of Surfex
 All Classes Files Functions Variables Typedefs
/home/dasprezs/EXPORT_v7_3/src/SURFEX/modd_chs_aerosol.F90
Go to the documentation of this file.
00001 !!     ######################
00002        MODULE MODD_CHS_AEROSOL
00003 !!     ######################
00004 !!
00005 !!     PURPOSE
00006 !!     -------
00007 !!
00008 !!     declaration of variables and types for the aerosol system in the
00009 !!     surface meteo-france scheme
00010 !!
00011 !!     METHOD
00012 !!     ------
00013 !!
00014 !!
00015 !!     REFERENCE
00016 !!     ---------
00017 !!     none
00018 !!
00019 !!
00020 !!     AUTHOR
00021 !!     ------
00022 !!     P. Tulet (CNRM)
00023 !!
00024 !!
00025 !!     MODIFICATIONS
00026 !!     -------------
00027 !!
00028 !!--------------------------------------------------------------------
00029 !!     DECLARATIONS
00030 !!     ------------
00031 IMPLICIT NONE
00032 !
00033 ! aerosol mode parameters
00034 LOGICAL      :: LCH_AERO_FLUX     = .FALSE. ! switch to active pronostic aerosols
00035 !
00036 LOGICAL      :: LCO2PM            = .FALSE. ! switch to active primary emission derived from CO 
00037 
00038 LOGICAL      :: LVARSIGI  = .FALSE.   ! switch to active pronostic dispersion for I mode
00039 LOGICAL      :: LVARSIGJ  = .FALSE.   ! switch to active pronostic dispersion for J mode
00040 
00041 !
00042 INTEGER, PARAMETER         :: JPMODE=2      ! number of modes
00043 INTEGER, PARAMETER         :: JPIN=JPMODE*3 ! number of differential equations
00044 INTEGER, SAVE, DIMENSION(JPMODE) :: NM0,NM3,NM6   ! index of the moments in arrays
00045 !
00046 !* indices of Aerosol chemical parameters
00047 !
00048 INTEGER, PARAMETER :: NSP=4        ! number of chemical species
00049                                    ! for ARES or isorropia NSP=4 these are
00050 INTEGER, PARAMETER :: JP_AER_SO4 = 1
00051 INTEGER, PARAMETER :: JP_AER_NO3 = 2
00052 INTEGER, PARAMETER :: JP_AER_NH3 = 3
00053 INTEGER, PARAMETER :: JP_AER_H2O = 4
00054 !
00055 INTEGER, PARAMETER :: JP_AER_SO4g = JP_AER_SO4
00056 INTEGER, PARAMETER :: JP_AER_NO3g = JP_AER_NO3
00057 INTEGER, PARAMETER :: JP_AER_NH3g = JP_AER_NH3
00058 !
00059 INTEGER, PARAMETER :: NCARB=2     ! number of chemically inert species
00060                                   ! (like black carbon)
00061 INTEGER, PARAMETER :: JP_AER_OC = 5
00062 INTEGER, PARAMETER :: JP_AER_BC = 6
00063 
00064 INTEGER            :: NSOA = 10    ! number of condensable species that may form
00065                                    ! secondary aerosols
00066 INTEGER            :: NM6_AER = 2  ! number of condensable species that may form
00067                                    ! secondary aerosols
00068 INTEGER            :: JP_AER_SOA1 = 7 
00069 INTEGER            :: JP_AER_SOA2 = 8
00070 INTEGER            :: JP_AER_SOA3 = 9
00071 INTEGER            :: JP_AER_SOA4 = 10
00072 INTEGER            :: JP_AER_SOA5 = 11
00073 INTEGER            :: JP_AER_SOA6 = 12
00074 INTEGER            :: JP_AER_SOA7 = 13
00075 INTEGER            :: JP_AER_SOA8 = 14
00076 INTEGER            :: JP_AER_SOA9 = 15 
00077 INTEGER            :: JP_AER_SOA10 = 16
00078 
00079  CHARACTER(LEN=32),DIMENSION(:), ALLOCATABLE :: CAERONAMES
00080 
00081 INTEGER            :: JP_CH_SO4I = 1  
00082 INTEGER            :: JP_CH_SO4J = 2  
00083 INTEGER            :: JP_CH_NO3I = 3  
00084 INTEGER            :: JP_CH_NO3J = 4  
00085 INTEGER            :: JP_CH_NH3I = 5  
00086 INTEGER            :: JP_CH_NH3J = 6  
00087 INTEGER            :: JP_CH_H2OI = 7  
00088 INTEGER            :: JP_CH_H2OJ = 8  
00089 INTEGER            :: JP_CH_OCI  = 9  
00090 INTEGER            :: JP_CH_OCJ  = 10  
00091 INTEGER            :: JP_CH_BCI  = 11  
00092 INTEGER            :: JP_CH_BCJ  = 12 
00093 INTEGER            :: JP_CH_SOA1I  = 13  
00094 INTEGER            :: JP_CH_SOA1J  = 14
00095 INTEGER            :: JP_CH_SOA2I  = 15
00096 INTEGER            :: JP_CH_SOA2J  = 16
00097 INTEGER            :: JP_CH_SOA3I  = 17
00098 INTEGER            :: JP_CH_SOA3J  = 18  
00099 INTEGER            :: JP_CH_SOA4I  = 19  
00100 INTEGER            :: JP_CH_SOA4J  = 20 
00101 INTEGER            :: JP_CH_SOA5I  = 21  
00102 INTEGER            :: JP_CH_SOA5J  = 22  
00103 INTEGER            :: JP_CH_SOA6I  = 23  
00104 INTEGER            :: JP_CH_SOA6J  = 24  
00105 INTEGER            :: JP_CH_SOA7I  = 25  
00106 INTEGER            :: JP_CH_SOA7J  = 26  
00107 INTEGER            :: JP_CH_SOA8I  = 27  
00108 INTEGER            :: JP_CH_SOA8J  = 28  
00109 INTEGER            :: JP_CH_SOA9I  = 29  
00110 INTEGER            :: JP_CH_SOA9J  = 30  
00111 INTEGER            :: JP_CH_SOA10I = 31  
00112 INTEGER            :: JP_CH_SOA10J = 32  
00113 INTEGER            :: JP_CH_M0I = 33  
00114 INTEGER            :: JP_CH_M0J = 34  
00115 INTEGER            :: JP_CH_M6I = 35  
00116 INTEGER            :: JP_CH_M6J = 36  
00117 
00118 !INTEGER, PARAMETER :: JPNN=NSP+NSOA+NCARB
00119                                    
00120 ! Index for gas species which interact with aerosols
00121 INTEGER :: JP_CH_HNO3,  JP_CH_H2SO4, JP_CH_NH3, JP_CH_O3, JP_CH_CO,  
00122              JP_CH_URG1, JP_CH_URG2, JP_CH_RPG2, JP_CH_RP18, JP_CH_UR26,
00123              JP_CH_RPG3, JP_CH_URG4, JP_CH_UR8, JP_CH_UR17, JP_CH_UR7, JP_CH_URG6, 
00124              JP_CH_ARAC, JP_CH_URG7, JP_CH_RPG7, JP_CH_RPR7, JP_CH_URG8, JP_CH_UR19, 
00125              JP_CH_URG9, JP_CH_URG10, JP_CH_PAN8, JP_CH_UR22, JP_CH_RPR4, JP_CH_AP7, 
00126              JP_CH_RPR3, JP_CH_UR21, JP_CH_UR28, JP_CH_UR29,  JP_CH_UR30, 
00127              JP_CH_RPR9, JP_CH_RP12, JP_CH_UR3, JP_CH_UR23, JP_CH_UR31, JP_CH_AP1, 
00128              JP_CH_AP6, JP_CH_ADAC, JP_CH_UR2, JP_CH_UR14, JP_CH_UR27, JP_CH_RP14, 
00129              JP_CH_RP19, JP_CH_UR11, JP_CH_UR15, JP_CH_AP10, JP_CH_UR20, JP_CH_UR34, 
00130              JP_CH_AP11, JP_CH_AP12, JP_CH_UR5, JP_CH_UR6, JP_CH_AP8  
00131 
00132 INTEGER :: JP_CH_H2O2,  JP_CH_SO2, JP_CH_SO42M
00133 !
00134 
00135 
00136 
00137 ! Molar mass of each aerosols parents (in kg/mol)
00138 REAL, PARAMETER :: XHNO3=63.01287
00139 REAL, PARAMETER :: XNH3 =17.03061
00140 REAL, PARAMETER :: XURG1=88.
00141 REAL, PARAMETER :: XURG2=1.76981E+02
00142 REAL, PARAMETER :: XRPG2=1.68000E+02
00143 REAL, PARAMETER :: XRP18=1.84000E+02
00144 REAL, PARAMETER :: XRPG3=1.53772E+02
00145 REAL, PARAMETER :: XURG4=1.95867E+02
00146 REAL, PARAMETER :: XUR17=1.72000E+02
00147 REAL, PARAMETER :: XRPR3=1.86000E+02
00148 REAL, PARAMETER :: XAP7 =2.33000E+02
00149 REAL, PARAMETER :: XURG6=1.89153E+02
00150 REAL, PARAMETER :: XUR22=2.12000E+02
00151 REAL, PARAMETER :: XURG7=1.56781E+02
00152 REAL, PARAMETER :: XRPR4=1.67000E+02
00153 REAL, PARAMETER :: XRPR7=1.50000E+02
00154 REAL, PARAMETER :: XRPG7=1.96059E+02
00155 REAL, PARAMETER :: XURG8=1.73777E+02
00156 REAL, PARAMETER :: XURG9=2.61676E+02
00157 REAL, PARAMETER :: XUR26=1.68000E+02
00158 REAL, PARAMETER :: XURG10=2.14834E+02
00159 REAL, PARAMETER :: XUR7=1.68000E+02
00160 REAL, PARAMETER :: XUR8=1.84000E+02
00161 REAL, PARAMETER :: XPAN8=2.63000E+02
00162 REAL, PARAMETER :: XARAC=1.32000E+02
00163 REAL, PARAMETER :: XUR19=1.70000E+02
00164 
00165 
00166 
00167 !
00168 REAL         :: XEMISRADIUSI  = 0.036   ! mean radius of primary aerosol
00169                                        ! emission for I mode
00170 REAL         :: XEMISRADIUSJ  = 0.385    ! mean radius of primary aerosol
00171                                        ! emission for J mode
00172 REAL         :: XEMISSIGI     = 1.86   ! dispersion of primary aerosol
00173                                        ! emission for I mode
00174 REAL         :: XEMISSIGJ     = 1.29   ! dispersion of primary aerosol
00175                                        ! emission for J mode
00176  CHARACTER*4  :: CRGUNIT   = 'NUMB'    ! type of log-normal geometric mean radius given
00177 !                                     ! in nameliste (mass on number)
00178 
00179 
00180 
00181 !----------------------------------------------------------------------------
00182 !
00183 !*  constants
00184 !
00185 REAL, PARAMETER :: XPBOLTZ=1.380658e-23    ! Boltzmann constant (J/K)
00186 REAL, PARAMETER :: XAVOGADRO=6.0221367E+23 ! Avogadro constant
00187 REAL, PARAMETER :: XMD    = 28.9644E-3     ! Air mass molarity
00188 
00189 !
00190 END MODULE MODD_CHS_AEROSOL