SURFEX v7.3
General documentation of Surfex
|
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