| 1 | SUBROUTINE SU_AERP |
|---|
| 2 | |
|---|
| 3 | !**** *SU_AERP* - INITIALIZE MODULES YOEAERSRC, YOEAERSNK |
|---|
| 4 | |
|---|
| 5 | ! PURPOSE. |
|---|
| 6 | ! -------- |
|---|
| 7 | ! INITIALIZE YOEAERSRC AND YOEAERSNK, THE MODULES THAT CONTAINS |
|---|
| 8 | ! COEFFICIENTS NEEDED TO RUN THE PROGNOSTIC AEROSOLS |
|---|
| 9 | |
|---|
| 10 | !** INTERFACE. |
|---|
| 11 | ! ---------- |
|---|
| 12 | ! *CALL* *SU_AERP |
|---|
| 13 | |
|---|
| 14 | ! EXPLICIT ARGUMENTS : |
|---|
| 15 | ! -------------------- |
|---|
| 16 | ! NONE |
|---|
| 17 | |
|---|
| 18 | ! IMPLICIT ARGUMENTS : |
|---|
| 19 | ! -------------------- |
|---|
| 20 | ! YOEAERSRC, YOEAERSNK, YOEAERATM |
|---|
| 21 | |
|---|
| 22 | ! METHOD. |
|---|
| 23 | ! ------- |
|---|
| 24 | ! SEE DOCUMENTATION |
|---|
| 25 | |
|---|
| 26 | ! EXTERNALS. |
|---|
| 27 | ! ---------- |
|---|
| 28 | |
|---|
| 29 | ! REFERENCE. |
|---|
| 30 | ! ---------- |
|---|
| 31 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS |
|---|
| 32 | |
|---|
| 33 | ! AUTHOR. |
|---|
| 34 | ! ------- |
|---|
| 35 | ! JEAN-JACQUES MORCRETTE *ECMWF* |
|---|
| 36 | ! from O.BOUCHER (LOA, 1998-03) |
|---|
| 37 | |
|---|
| 38 | ! MODIFICATIONS. |
|---|
| 39 | ! -------------- |
|---|
| 40 | ! ORIGINAL : 2004-05-10 |
|---|
| 41 | |
|---|
| 42 | ! ------------------------------------------------------------------ |
|---|
| 43 | |
|---|
| 44 | USE PARKIND1 ,ONLY : JPRB |
|---|
| 45 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
|---|
| 46 | |
|---|
| 47 | USE YOEAERSRC , ONLY : RSSFLX |
|---|
| 48 | |
|---|
| 49 | USE YOEAERSNK , ONLY : R_R, R_S, RALPHAR, RALPHAS, RFRAER, RFRGAS, & |
|---|
| 50 | & RRHMAX, RRHTAB, RRHO_SS, RSSGROW, RMMD_SS, RMMD_DD, RRHO_DD, & |
|---|
| 51 | & RFRBC , RFRIF , RFROM , RFRSO4 , RFRDD , RFRSS , RHO_WAT, RHO_ICE, & |
|---|
| 52 | & RVDPOCE, RVDPSIC, RVDPLND, RVDPLIC, RVSEDOCE, RVSEDSIC, RVSEDLND, RVSEDLIC, & |
|---|
| 53 | & NBRH |
|---|
| 54 | |
|---|
| 55 | USE YOEAERATM , ONLY : RMASSE, RMFMIN, NDD1, NSS1 |
|---|
| 56 | |
|---|
| 57 | IMPLICIT NONE |
|---|
| 58 | |
|---|
| 59 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
|---|
| 60 | ! ---------------------------------------------------------------- |
|---|
| 61 | IF (LHOOK) CALL DR_HOOK('SU_AERP',0,ZHOOK_HANDLE) |
|---|
| 62 | |
|---|
| 63 | !-- For the ECMWF model, the following tables when dimensioned to 12 |
|---|
| 64 | ! can refer to 12 values of RH |
|---|
| 65 | ! (RHTAB, RHHO_SS, RSSGROW) |
|---|
| 66 | ! or to 12 types/bins of aerosols with the following mapping: NTYPAER |
|---|
| 67 | ! 1- 3 sea-salt 0.03 - 0.5 - 5 - 20 microns 1 |
|---|
| 68 | ! 4- 6 dust 0.03 - 0.5 - 0.9 - 20 microns 2 |
|---|
| 69 | ! 7- 8 POM hydrophilic, hydrophobic 3 |
|---|
| 70 | ! 9-10 BC hydrophilic, hydrophobic 4 |
|---|
| 71 | ! 11 sulfate 5 |
|---|
| 72 | ! 12 fly ash 6 |
|---|
| 73 | ! 13 stratospheric aerosols 7 |
|---|
| 74 | ! 14 volcanic aerosols 8 |
|---|
| 75 | ! 15 9 |
|---|
| 76 | ! (RVDPOCE, RVDSIC, RVDPLND, RVDPLIC) |
|---|
| 77 | ! (RVSEDOCE,RVSEDSIC,RVSEDLND,RVSEDLIC) |
|---|
| 78 | |
|---|
| 79 | |
|---|
| 80 | |
|---|
| 81 | !* 1. PARAMETERS RELATED TO SOURCES |
|---|
| 82 | ! ----------------------------- |
|---|
| 83 | |
|---|
| 84 | !-- OB's original 12 types and 24 different values had the following mapping |
|---|
| 85 | ! DMS SO2 SO4 H2S DMSO MSA H2O2 |
|---|
| 86 | ! BC(2) OM(2) FlyAsh DU(2) SS(10) |
|---|
| 87 | |
|---|
| 88 | !-- OB's original SS 10 bins |
|---|
| 89 | ! bin sizes: 0.03-0.06-0.13-0.25-0.5-1.0-2.0-5.0-10.-15.-20 |
|---|
| 90 | !RSSFLX = (/ & |
|---|
| 91 | ! & 0.20526E-09_JPRB, 0.49292E-09_JPRB, 0.97079E-09_JPRB, 0.31938E-08_JPRB & |
|---|
| 92 | ! &, 0.16245E-07_JPRB, 0.86292E-07_JPRB, 0.31326E-06_JPRB, 0.24671E-06_JPRB & |
|---|
| 93 | ! &, 0.14109E-06_JPRB, 0.11784E-06_JPRB /) |
|---|
| 94 | |
|---|
| 95 | ! maximum possible number of aerosol types |
|---|
| 96 | !NMAXTAER=9 already defined in SU_AERW |
|---|
| 97 | |
|---|
| 98 | !N.B. Fluxes of sea salt for each size bin are given in mg m-2 s-1 at wind |
|---|
| 99 | ! speed of 1 m s-1 at 10m height (at 80% RH) in OB's seasalt.F |
|---|
| 100 | ! RSSFLX also in mg m-2 s-1 |
|---|
| 101 | !-- OB's ECMWF 3 bins of sea salt: 0.03, 0.5, 5, 20 microns |
|---|
| 102 | RSSFLX = (/ 4.85963536E-09_JPRB, 4.15358556E-07_JPRB, 5.04905813E-07_JPRB /) |
|---|
| 103 | |
|---|
| 104 | ! OB's original vdep_oce, vdep_sic, vdep_ter, vdep_lic were given over 24 values |
|---|
| 105 | |
|---|
| 106 | ! following 12 values for 10 SS and 2 DU in cm s-1 |
|---|
| 107 | !RVDPOCE = (/ 0.1_JPRB, 1.2_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB & |
|---|
| 108 | ! &, 0.1_JPRB, 1.2_JPRB, 1.2_JPRB, 1.2_JPRB, 1.5_JPRB, 1.5_JPRB /) |
|---|
| 109 | ! |
|---|
| 110 | !RVDPSIC = (/ 0.1_JPRB, 1.2_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB & |
|---|
| 111 | ! &, 0.1_JPRB, 1.2_JPRB, 1.2_JPRB, 1.2_JPRB, 1.5_JPRB, 1.5_JPRB /) |
|---|
| 112 | ! |
|---|
| 113 | !RVDPLND = (/ 0.1_JPRB, 1.2_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB & |
|---|
| 114 | ! &, 0.1_JPRB, 1.2_JPRB, 1.2_JPRB, 1.2_JPRB, 1.5_JPRB, 1.5_JPRB /) |
|---|
| 115 | ! |
|---|
| 116 | !RVDPLIC = (/ 0.1_JPRB, 1.2_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB, 0.1_JPRB & |
|---|
| 117 | ! &, 0.1_JPRB, 1.2_JPRB, 1.2_JPRB, 1.2_JPRB, 1.5_JPRB, 1.5_JPRB /) |
|---|
| 118 | |
|---|
| 119 | |
|---|
| 120 | |
|---|
| 121 | !* 2. PARAMETERS RELATED TO SINKS |
|---|
| 122 | ! --------------------------- |
|---|
| 123 | |
|---|
| 124 | R_R = 0.001_JPRB |
|---|
| 125 | R_S = 0.001_JPRB |
|---|
| 126 | |
|---|
| 127 | RFRAER = 0.5_JPRB |
|---|
| 128 | RFRGAS = 1.0_JPRB |
|---|
| 129 | |
|---|
| 130 | !* 2.1 SEA SALT |
|---|
| 131 | ! -------- |
|---|
| 132 | !-- parameters related to SEA SALT: 12 relates to 12 values of relative humidity |
|---|
| 133 | |
|---|
| 134 | NBRH=12 |
|---|
| 135 | RRHMAX = 95._JPRB |
|---|
| 136 | RRHTAB = (/ 0._JPRB, 10._JPRB, 20._JPRB, 30._JPRB, 40._JPRB, 50._JPRB & |
|---|
| 137 | & , 60._JPRB, 70._JPRB, 80._JPRB, 85._JPRB, 90._JPRB, 95._JPRB /) |
|---|
| 138 | RRHO_SS = (/ 2160._JPRB, 2160._JPRB, 2160._JPRB, 2160._JPRB, 1451.6_JPRB & |
|---|
| 139 | & , 1367.9_JPRB, 1302.9_JPRB, 1243.2_JPRB, 1182.7_JPRB, 1149.5_JPRB & |
|---|
| 140 | & , 1111.6_JPRB, 1063.1_JPRB /) |
|---|
| 141 | RSSGROW = (/ 0.503_JPRB, 0.503_JPRB, 0.503_JPRB, 0.503_JPRB, 0.724_JPRB & |
|---|
| 142 | & , 0.782_JPRB, 0.838_JPRB, 0.905_JPRB, 1.000_JPRB, 1.072_JPRB & |
|---|
| 143 | & , 1.188_JPRB, 1.447_JPRB /) |
|---|
| 144 | |
|---|
| 145 | !-- OB's original 10 bins !RMMD_SS = (/ 0.09_JPRB, 0.19_JPRB, 0.38_JPRB, 0.75_JPRB, 1.50_JPRB & |
|---|
| 146 | ! & , 3.00_JPRB, 7.00_JPRB, 15.0_JPRB, 25.0_JPRB, 35.0_JPRB /) |
|---|
| 147 | |
|---|
| 148 | !-- OB's ECMWF 3 bins of sea salt |
|---|
| 149 | ! bins are 0.03 - 0.5 - 5.0 - 20 microns |
|---|
| 150 | |
|---|
| 151 | RMMD_SS = (/ 0.30_JPRB, 3.00_JPRB, 10.00_JPRB /) |
|---|
| 152 | RFRSS = (/ 0.7_JPRB, 0.7_JPRB, 0.7_JPRB /) |
|---|
| 153 | RHO_WAT = 1000._JPRB |
|---|
| 154 | RHO_ICE = 500._JPRB |
|---|
| 155 | |
|---|
| 156 | !- computed off-line by gems_ss.f (m s-1) |
|---|
| 157 | |
|---|
| 158 | RVSEDOCE(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
|---|
| 159 | RVSEDSIC(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
|---|
| 160 | RVSEDLND(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
|---|
| 161 | RVSEDLIC(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
|---|
| 162 | |
|---|
| 163 | ! adapted from LMDZ (m s-1) |
|---|
| 164 | |
|---|
| 165 | RVDPOCE(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
|---|
| 166 | RVDPSIC(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
|---|
| 167 | RVDPLND(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
|---|
| 168 | RVDPLIC(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
|---|
| 169 | |
|---|
| 170 | |
|---|
| 171 | !* 2.2 DESERT DUST |
|---|
| 172 | ! ----------- |
|---|
| 173 | !- parameters related to DESERT DUST (OB's ECMWF 3 bins) |
|---|
| 174 | ! bins are 0.03 - 0.55 - 0.9 - 20 microns |
|---|
| 175 | |
|---|
| 176 | RMMD_DD = (/ 0.32_JPRB, 0.75_JPRB, 9.0_JPRB /) |
|---|
| 177 | RRHO_DD = (/ 2600._JPRB, 2600._JPRB, 2600._JPRB /) |
|---|
| 178 | RFRDD = (/ 0.7_JPRB, 0.7_JPRB, 0.7_JPRB /) |
|---|
| 179 | |
|---|
| 180 | !- computed off-line by gems_dust.f |
|---|
| 181 | |
|---|
| 182 | RVSEDOCE(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
|---|
| 183 | RVSEDSIC(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
|---|
| 184 | RVSEDLND(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
|---|
| 185 | RVSEDLIC(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
|---|
| 186 | |
|---|
| 187 | ! adapted from LMDZ (m s-1) |
|---|
| 188 | |
|---|
| 189 | RVDPOCE(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
|---|
| 190 | RVDPSIC(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
|---|
| 191 | RVDPLND(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
|---|
| 192 | RVDPLIC(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
|---|
| 193 | |
|---|
| 194 | |
|---|
| 195 | !* 2.3 OTHER AEROSOLS (to be improved later!) |
|---|
| 196 | ! -------------- |
|---|
| 197 | !- parameters related to other aerosol types |
|---|
| 198 | !- particulate organic matter POM |
|---|
| 199 | RFROM = (/ 0.0_JPRB, 0.7_JPRB /) |
|---|
| 200 | |
|---|
| 201 | RVSEDOCE(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 202 | RVSEDSIC(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 203 | RVSEDLND(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 204 | RVSEDLIC(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 205 | |
|---|
| 206 | ! adapted from LMDZ (m s-1) |
|---|
| 207 | |
|---|
| 208 | RVDPOCE(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 209 | RVDPSIC(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 210 | RVDPLND(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 211 | RVDPLIC(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 212 | |
|---|
| 213 | |
|---|
| 214 | !- black carbon |
|---|
| 215 | RFRBC = (/ 0.0_JPRB, 0.7_JPRB /) |
|---|
| 216 | |
|---|
| 217 | RVSEDOCE(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 218 | RVSEDSIC(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 219 | RVSEDLND(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 220 | RVSEDLIC(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
|---|
| 221 | |
|---|
| 222 | ! adapted from LMDZ (m s-1) |
|---|
| 223 | |
|---|
| 224 | RVDPOCE(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 225 | RVDPSIC(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 226 | RVDPLND(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 227 | RVDPLIC(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
|---|
| 228 | |
|---|
| 229 | !- sulfate |
|---|
| 230 | RFRSO4 = 0.7_JPRB |
|---|
| 231 | |
|---|
| 232 | RVSEDOCE(11) = 0.05_JPRB |
|---|
| 233 | RVSEDSIC(11) = 0.25_JPRB |
|---|
| 234 | RVSEDLND(11) = 0.25_JPRB |
|---|
| 235 | RVSEDLIC(11) = 0.25_JPRB |
|---|
| 236 | |
|---|
| 237 | ! adapted from LMDZ (m s-1) |
|---|
| 238 | |
|---|
| 239 | RVDPOCE(11) = 0.05E-02_JPRB |
|---|
| 240 | RVDPSIC(11) = 0.25E-02_JPRB |
|---|
| 241 | RVDPLND(11) = 0.25E-02_JPRB |
|---|
| 242 | RVDPLIC(11) = 0.25E-02_JPRB |
|---|
| 243 | |
|---|
| 244 | !- fly ash |
|---|
| 245 | RFRIF = 0.7_JPRB |
|---|
| 246 | |
|---|
| 247 | RVSEDOCE(12) = 0.20E+00_JPRB |
|---|
| 248 | RVSEDSIC(12) = 0.20E+00_JPRB |
|---|
| 249 | RVSEDLND(12) = 0.20E+00_JPRB |
|---|
| 250 | RVSEDLIC(12) = 0.20E+00_JPRB |
|---|
| 251 | |
|---|
| 252 | ! adapted from LMDZ (m s-1) |
|---|
| 253 | |
|---|
| 254 | RVDPOCE(12) = 0.20E-02_JPRB |
|---|
| 255 | RVDPSIC(12) = 0.20E-02_JPRB |
|---|
| 256 | RVDPLND(12) = 0.20E-02_JPRB |
|---|
| 257 | RVDPLIC(12) = 0.20E-02_JPRB |
|---|
| 258 | |
|---|
| 259 | |
|---|
| 260 | |
|---|
| 261 | |
|---|
| 262 | !- NB: 15 values for all possible types of ECMWF aerosols |
|---|
| 263 | RALPHAR = (/ & |
|---|
| 264 | & 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, & |
|---|
| 265 | & 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB, 0.001_JPRB /) |
|---|
| 266 | RALPHAS = (/ & |
|---|
| 267 | & 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, & |
|---|
| 268 | & 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB, 0.01_JPRB /) |
|---|
| 269 | |
|---|
| 270 | |
|---|
| 271 | !* 3. PARAMETERS RELATED TO TRANSPORT WITHIN THE FREE ATMOSPHERE |
|---|
| 272 | ! ---------------------------------------------------------- |
|---|
| 273 | |
|---|
| 274 | NDD1=4 |
|---|
| 275 | NSS1=1 |
|---|
| 276 | |
|---|
| 277 | RMFMIN = 1.E-10_JPRB |
|---|
| 278 | |
|---|
| 279 | RMASSE = (/ & |
|---|
| 280 | & 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB & |
|---|
| 281 | &, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB & |
|---|
| 282 | &, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB /) |
|---|
| 283 | |
|---|
| 284 | |
|---|
| 285 | ! ---------------------------------------------------------------- |
|---|
| 286 | IF (LHOOK) CALL DR_HOOK('SU_AERP',1,ZHOOK_HANDLE) |
|---|
| 287 | END SUBROUTINE SU_AERP |
|---|
| 288 | |
|---|
