1 |
|
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 |
✓✗ |
1 |
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 |
|
1 |
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 |
|
1 |
R_R = 0.001_JPRB |
125 |
|
1 |
R_S = 0.001_JPRB |
126 |
|
|
|
127 |
|
1 |
RFRAER = 0.5_JPRB |
128 |
|
1 |
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 |
|
1 |
NBRH=12 |
135 |
|
1 |
RRHMAX = 95._JPRB |
136 |
|
|
RRHTAB = (/ 0._JPRB, 10._JPRB, 20._JPRB, 30._JPRB, 40._JPRB, 50._JPRB & |
137 |
|
1 |
& , 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 |
|
1 |
& , 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 |
& , 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 |
|
1 |
RMMD_SS = (/ 0.30_JPRB, 3.00_JPRB, 10.00_JPRB /) |
152 |
|
1 |
RFRSS = (/ 0.7_JPRB, 0.7_JPRB, 0.7_JPRB /) |
153 |
|
1 |
RHO_WAT = 1000._JPRB |
154 |
|
1 |
RHO_ICE = 500._JPRB |
155 |
|
|
|
156 |
|
|
!- computed off-line by gems_ss.f (m s-1) |
157 |
|
|
|
158 |
✓✓ |
4 |
RVSEDOCE(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
159 |
✓✓ |
4 |
RVSEDSIC(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
160 |
✓✓ |
4 |
RVSEDLND(1:3) = (/ 0.24E-04_JPRB, 0.20E-02_JPRB, 0.20E-01_JPRB /) |
161 |
✓✓ |
4 |
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 |
✓✓ |
4 |
RVDPOCE(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
166 |
✓✓ |
4 |
RVDPSIC(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
167 |
✓✓ |
4 |
RVDPLND(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB, 0.145E-01_JPRB /) |
168 |
✓✓ |
4 |
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 |
|
1 |
RMMD_DD = (/ 0.32_JPRB, 0.75_JPRB, 9.0_JPRB /) |
177 |
|
1 |
RRHO_DD = (/ 2600._JPRB, 2600._JPRB, 2600._JPRB /) |
178 |
|
1 |
RFRDD = (/ 0.7_JPRB, 0.7_JPRB, 0.7_JPRB /) |
179 |
|
|
|
180 |
|
|
!- computed off-line by gems_dust.f |
181 |
|
|
|
182 |
✓✓ |
4 |
RVSEDOCE(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
183 |
✓✓ |
4 |
RVSEDSIC(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
184 |
✓✓ |
4 |
RVSEDLND(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /) |
185 |
✓✓ |
4 |
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 |
✓✓ |
4 |
RVDPOCE(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
190 |
✓✓ |
4 |
RVDPSIC(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
191 |
✓✓ |
4 |
RVDPLND(4:6) = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /) |
192 |
✓✓ |
4 |
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 |
|
1 |
RFROM = (/ 0.0_JPRB, 0.7_JPRB /) |
200 |
|
|
|
201 |
✓✓ |
3 |
RVSEDOCE(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
202 |
✓✓ |
3 |
RVSEDSIC(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
203 |
✓✓ |
3 |
RVSEDLND(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
204 |
✓✓ |
3 |
RVSEDLIC(7:8) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
205 |
|
|
|
206 |
|
|
! adapted from LMDZ (m s-1) |
207 |
|
|
|
208 |
✓✓ |
3 |
RVDPOCE(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
209 |
✓✓ |
3 |
RVDPSIC(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
210 |
✓✓ |
3 |
RVDPLND(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
211 |
✓✓ |
3 |
RVDPLIC(7:8) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
212 |
|
|
|
213 |
|
|
|
214 |
|
|
!- black carbon |
215 |
|
1 |
RFRBC = (/ 0.0_JPRB, 0.7_JPRB /) |
216 |
|
|
|
217 |
✓✓ |
3 |
RVSEDOCE(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
218 |
✓✓ |
3 |
RVSEDSIC(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
219 |
✓✓ |
3 |
RVSEDLND(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
220 |
✓✓ |
3 |
RVSEDLIC(9:10) = (/ 0.10E+00_JPRB, 0.10E+00_JPRB /) |
221 |
|
|
|
222 |
|
|
! adapted from LMDZ (m s-1) |
223 |
|
|
|
224 |
✓✓ |
3 |
RVDPOCE(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
225 |
✓✓ |
3 |
RVDPSIC(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
226 |
✓✓ |
3 |
RVDPLND(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
227 |
✓✓ |
3 |
RVDPLIC(9:10) = (/ 0.10E-02_JPRB, 0.10E-02_JPRB /) |
228 |
|
|
|
229 |
|
|
!- sulfate |
230 |
|
1 |
RFRSO4 = 0.7_JPRB |
231 |
|
|
|
232 |
|
1 |
RVSEDOCE(11) = 0.05_JPRB |
233 |
|
1 |
RVSEDSIC(11) = 0.25_JPRB |
234 |
|
1 |
RVSEDLND(11) = 0.25_JPRB |
235 |
|
1 |
RVSEDLIC(11) = 0.25_JPRB |
236 |
|
|
|
237 |
|
|
! adapted from LMDZ (m s-1) |
238 |
|
|
|
239 |
|
1 |
RVDPOCE(11) = 0.05E-02_JPRB |
240 |
|
1 |
RVDPSIC(11) = 0.25E-02_JPRB |
241 |
|
1 |
RVDPLND(11) = 0.25E-02_JPRB |
242 |
|
1 |
RVDPLIC(11) = 0.25E-02_JPRB |
243 |
|
|
|
244 |
|
|
!- fly ash |
245 |
|
1 |
RFRIF = 0.7_JPRB |
246 |
|
|
|
247 |
|
1 |
RVSEDOCE(12) = 0.20E+00_JPRB |
248 |
|
1 |
RVSEDSIC(12) = 0.20E+00_JPRB |
249 |
|
1 |
RVSEDLND(12) = 0.20E+00_JPRB |
250 |
|
1 |
RVSEDLIC(12) = 0.20E+00_JPRB |
251 |
|
|
|
252 |
|
|
! adapted from LMDZ (m s-1) |
253 |
|
|
|
254 |
|
1 |
RVDPOCE(12) = 0.20E-02_JPRB |
255 |
|
1 |
RVDPSIC(12) = 0.20E-02_JPRB |
256 |
|
1 |
RVDPLND(12) = 0.20E-02_JPRB |
257 |
|
1 |
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 |
|
1 |
& 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 |
|
1 |
& 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 |
|
1 |
NDD1=4 |
275 |
|
1 |
NSS1=1 |
276 |
|
|
|
277 |
|
1 |
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 |
|
1 |
&, 6.02E+23_JPRB, 6.02E+23_JPRB, 6.02E+23_JPRB /) |
283 |
|
|
|
284 |
|
|
|
285 |
|
|
! ---------------------------------------------------------------- |
286 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SU_AERP',1,ZHOOK_HANDLE) |
287 |
|
1 |
END SUBROUTINE SU_AERP |
288 |
|
|
|