source: LMDZ5/branches/testing/libf/phylmd/rrtm/su_aerp.F90 @ 2056

Last change on this file since 2056 was 1999, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1920:1997 into testing branch
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 9.4 KB
Line 
1SUBROUTINE 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
44USE PARKIND1  ,ONLY : JPRB
45USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
46
47USE YOEAERSRC , ONLY : RSSFLX
48
49USE 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
55USE YOEAERATM , ONLY : RMASSE, RMFMIN, NDD1, NSS1
56
57IMPLICIT NONE
58
59REAL(KIND=JPRB) :: ZHOOK_HANDLE
60!     ----------------------------------------------------------------
61IF (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 
102RSSFLX = (/ 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
124R_R = 0.001_JPRB
125R_S = 0.001_JPRB
126
127RFRAER = 0.5_JPRB
128RFRGAS = 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
134NBRH=12
135RRHMAX = 95._JPRB
136RRHTAB = (/ 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 /)
138RRHO_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 /)
141RSSGROW = (/ 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
151RMMD_SS = (/ 0.30_JPRB, 3.00_JPRB, 10.00_JPRB /)
152RFRSS   = (/  0.7_JPRB,  0.7_JPRB,   0.7_JPRB /)
153RHO_WAT = 1000._JPRB
154RHO_ICE = 500._JPRB
155
156!-  computed off-line by gems_ss.f  (m s-1)
157
158RVSEDOCE(1:3) = (/   0.24E-04_JPRB,   0.20E-02_JPRB,   0.20E-01_JPRB /)
159RVSEDSIC(1:3) = (/   0.24E-04_JPRB,   0.20E-02_JPRB,   0.20E-01_JPRB /)
160RVSEDLND(1:3) = (/   0.24E-04_JPRB,   0.20E-02_JPRB,   0.20E-01_JPRB /)
161RVSEDLIC(1:3) = (/   0.24E-04_JPRB,   0.20E-02_JPRB,   0.20E-01_JPRB /)
162
163! adapted from LMDZ   (m s-1)
164
165RVDPOCE(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB,   0.145E-01_JPRB /)
166RVDPSIC(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB,   0.145E-01_JPRB /)
167RVDPLND(1:3) = (/ 0.100E-02_JPRB, 0.110E-01_JPRB,   0.145E-01_JPRB /)
168RVDPLIC(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
176RMMD_DD = (/  0.32_JPRB,  0.75_JPRB,   9.0_JPRB /)
177RRHO_DD = (/ 2600._JPRB, 2600._JPRB, 2600._JPRB /)
178RFRDD   = (/   0.7_JPRB,   0.7_JPRB,   0.7_JPRB /)
179
180!-  computed off-line by gems_dust.f
181
182RVSEDOCE(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /)
183RVSEDSIC(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /)
184RVSEDLND(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /)
185RVSEDLIC(4:6) = (/ 0.70E-04_JPRB, 0.20E-03_JPRB, 0.20E-02_JPRB /)
186
187! adapted from LMDZ   (m s-1)
188
189RVDPOCE(4:6)  = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /)
190RVDPSIC(4:6)  = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /)
191RVDPLND(4:6)  = (/ 0.100E-02_JPRB, 0.500E-02_JPRB, 0.110E-01_JPRB /)
192RVDPLIC(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
199RFROM = (/  0.0_JPRB,  0.7_JPRB /)
200
201RVSEDOCE(7:8) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
202RVSEDSIC(7:8) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
203RVSEDLND(7:8) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
204RVSEDLIC(7:8) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
205
206! adapted from LMDZ   (m s-1)
207
208RVDPOCE(7:8) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
209RVDPSIC(7:8) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
210RVDPLND(7:8) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
211RVDPLIC(7:8) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
212
213
214!- black carbon
215RFRBC = (/  0.0_JPRB,  0.7_JPRB /)
216
217RVSEDOCE(9:10) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
218RVSEDSIC(9:10) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
219RVSEDLND(9:10) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
220RVSEDLIC(9:10) = (/   0.10E+00_JPRB,   0.10E+00_JPRB /)
221
222! adapted from LMDZ   (m s-1)
223
224RVDPOCE(9:10) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
225RVDPSIC(9:10) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
226RVDPLND(9:10) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
227RVDPLIC(9:10) = (/   0.10E-02_JPRB,   0.10E-02_JPRB /)
228
229!- sulfate
230RFRSO4 = 0.7_JPRB
231
232RVSEDOCE(11) = 0.05_JPRB
233RVSEDSIC(11) = 0.25_JPRB
234RVSEDLND(11) = 0.25_JPRB
235RVSEDLIC(11) = 0.25_JPRB
236
237! adapted from LMDZ   (m s-1)
238
239RVDPOCE(11) = 0.05E-02_JPRB
240RVDPSIC(11) = 0.25E-02_JPRB
241RVDPLND(11) = 0.25E-02_JPRB
242RVDPLIC(11) = 0.25E-02_JPRB
243
244!- fly ash
245RFRIF  = 0.7_JPRB
246
247RVSEDOCE(12) = 0.20E+00_JPRB
248RVSEDSIC(12) = 0.20E+00_JPRB
249RVSEDLND(12) = 0.20E+00_JPRB
250RVSEDLIC(12) = 0.20E+00_JPRB
251
252! adapted from LMDZ   (m s-1)
253
254RVDPOCE(12) = 0.20E-02_JPRB
255RVDPSIC(12) = 0.20E-02_JPRB
256RVDPLND(12) = 0.20E-02_JPRB
257RVDPLIC(12) = 0.20E-02_JPRB
258
259
260
261
262!- NB: 15 values for all possible types of ECMWF aerosols
263RALPHAR = (/ &
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 /)
266RALPHAS = (/ &
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
274NDD1=4
275NSS1=1
276
277RMFMIN = 1.E-10_JPRB
278
279RMASSE = (/ &
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!     ----------------------------------------------------------------
286IF (LHOOK) CALL DR_HOOK('SU_AERP',1,ZHOOK_HANDLE)
287END SUBROUTINE SU_AERP
288
Note: See TracBrowser for help on using the repository browser.