Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

optaer2_5wv.F90 @ 6610

Last change on this file since 6610 was 6610, checked in by acosce, 10 months ago
INCA used for ICOLMDZORINCA_CO2_Transport_GMD_2023
File size: 48.1 KB

Line
1	!! =========================================================================
2	!! INCA - INteraction with Chemistry and Aerosols
3	!!
4	!! Copyright Laboratoire des Sciences du Climat et de l'Environnement (LSCE)
5	!! Unite mixte CEA-CNRS-UVSQ
6	!!
7	!! Contributors to this INCA subroutine:
8	!!
9	!! Yves Balkanski, LSCE,
10	!! Celine Deandreis, LSCE
11	!! Johannes Quaas
12	!! Nicolas Yan, LSCE
13	!!
14	!!
15	!! Anne Cozic, LSCE, anne.cozic@cea.fr
16	!! Yann Meurdesoif, LSCE, yann.meurdesoif@cea.fr
17	!!
18	!!
19	!! This software is a computer program whose purpose is to simulate the
20	!! atmospheric gas phase and aerosol composition. The model is designed to be
21	!! used within a transport model or a general circulation model. This version
22	!! of INCA was designed to be coupled to the LMDz GCM. LMDz-INCA accounts
23	!! for emissions, transport (resolved and sub-grid scale), photochemical
24	!! transformations, and scavenging (dry deposition and washout) of chemical
25	!! species and aerosols interactively in the GCM. Several versions of the INCA
26	!! model are currently used depending on the envisaged applications with the
27	!! chemistry-climate model.
28	!!
29	!! This software is governed by the CeCILL license under French law and
30	!! abiding by the rules of distribution of free software. You can use,
31	!! modify and/ or redistribute the software under the terms of the CeCILL
32	!! license as circulated by CEA, CNRS and INRIA at the following URL
33	!! "http://www.cecill.info".
34	!!
35	!! As a counterpart to the access to the source code and rights to copy,
36	!! modify and redistribute granted by the license, users are provided only
37	!! with a limited warranty and the software's author, the holder of the
38	!! economic rights, and the successive licensors have only limited
39	!! liability.
40	!!
41	!! In this respect, the user's attention is drawn to the risks associated
42	!! with loading, using, modifying and/or developing or reproducing the
43	!! software by the user in light of its specific status of free software,
44	!! that may mean that it is complicated to manipulate, and that also
45	!! therefore means that it is reserved for developers and experienced
46	!! professionals having in-depth computer knowledge. Users are therefore
47	!! encouraged to load and test the software's suitability as regards their
48	!! requirements in conditions enabling the security of their systems and/or
49	!! data to be ensured and, more generally, to use and operate it in the
50	!! same conditions as regards security.
51	!!
52	!! The fact that you are presently reading this means that you have had
53	!! knowledge of the CeCILL license and that you accept its terms.
54	!! =========================================================================
55	#include <inca_define.h>
56
57	#ifdef AER
58
59	SUBROUTINE OPTAER_5WV(pdel, tr_seri, delt, RHcl, ai)
60
61	USE INCA_DIM
62	USE SPECIES_NAMES
63	USE CHEM_TRACNM, ONLY: solsym
64	USE AEROSOL_MOD, ONLY: nmodes,mass1index,massNindex, &
65	numberindex,ntr,trmx,trnx,aimode,asmode, &
66	cimode,csmode,ssmode
67	USE AEROSOL_DIAG , ONLY : tausum,angst,lambda,las,cla,tau
68	USE AEROSOL_METEO, ONLY : zdp1
69	USE CHEM_CONS , ONLY : gravit
70	!
71	! (avril 2006 - mai 2009)
72	! a partir d'une sous-routine de Johannes Quaas pour les sulfates
73	!
74	!
75	! Refractive indices for seasalt come from Shettle and Fenn (1979)
76	!
77	! Refractive indices from water come from Hale and Querry (1973)
78	!
79	! Refractive indices from Ammonium Sulfate Toon and Pollack (1976)
80	!
81	! Refractive indices for Dust, internal mixture of minerals coated with 1.5% hematite
82	! by Volume (Balkanski et al., 2006)
83	!
84	! Refractive indices for POM: Kinne (pers. Communication
85	!
86	! Refractive index for BC from Shettle and Fenn (1979)
87	!
88	! Shettle, E. P., & Fenn, R. W. (1979), Models for the aerosols of the lower atmosphere and
89	! the effects of humidity variations on their optical properties, U.S. Air Force Geophysics
90	! Laboratory Rept. AFGL-TR-79-0214, Hanscomb Air Force Base, MA.
91	!
92	! Hale, G. M. and M. R. Querry, Optical constants of water in the 200-nm to 200-m
93	! wavelength region, Appl. Opt., 12, 555-563, 1973.
94	!
95	! Toon, O. B. and J. B. Pollack, The optical constants of several atmospheric aerosol species:
96	! Ammonium sulfate, aluminum oxide, and sodium chloride, J. Geohys. Res., 81, 5733-5748,
97	! 1976.
98	!
99	! Balkanski, Y., M. Schulz, T. Claquin And O. Boucher, Reevaluation of mineral aerosol
100	! radiative forcings suggests a better agreement with satellite and AERONET data, Atmospheric
101	! Chemistry and Physics Discussions., 6, pp 8383-8419, 2006.
102	!
103	IMPLICIT NONE
104	!
105	! Arguments:
106	!
107	REAL, DIMENSION(PLON,PLEV), INTENT(in) :: pdel
108	REAL, INTENT(in) :: delt
109	REAL, DIMENSION(PLON,PLEV,PCNST), INTENT(in) :: tr_seri ! Mass Mixing ratio kg/kg
110	REAL, DIMENSION(PLON,PLEV), INTENT(in) :: RHcl ! humidite relative ciel clair
111	REAL, DIMENSION(PLON), INTENT(out) :: ai ! POLDER aerosol index
112
113	!
114	! Local
115	!
116	LOGICAL :: soluble
117	REAL :: cg_temp, piz_temp, cg_temp2, piz_temp2
118
119	INTEGER :: i, k, m, nx
120	INTEGER :: spsol, spinsol, la, na
121	INTEGER :: spss
122	INTEGER :: RH_num
123	INTEGER, PARAMETER :: nbre_RH = 12
124	INTEGER, PARAMETER :: naero_soluble = 13 ! 1- BC soluble; 2- POM soluble; 3- SO4 accumulation
125	! 4- SO4 Coarse 5- Seasalt SS: 6- Seasalt CS; 7- Seasalt AS
126	! 8- Nitrate accumulation Coarse 9- Nitrate coarse 10 SOAA1 11 SOAA2 12 SOAB1 13 SOAB2
127	INTEGER, PARAMETER :: naero_insoluble = 4 ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
128
129	REAL :: A1_ASSSM(PLEV), A2_ASSSM(PLEV), A3_ASSSM(PLEV),&
130	B1_ASSSM(PLEV), B2_ASSSM(PLEV), C1_ASSSM(PLEV), C2_ASSSM(PLEV),&
131	A1_CSSSM(PLEV), A2_CSSSM(PLEV), A3_CSSSM(PLEV),&
132	B1_CSSSM(PLEV), B2_CSSSM(PLEV), C1_CSSSM(PLEV), C2_CSSSM(PLEV), &
133	A1_SSSSM(PLEV), A2_SSSSM(PLEV), A3_SSSSM(PLEV),&
134	B1_SSSSM(PLEV), B2_SSSSM(PLEV), C1_SSSSM(PLEV), C2_SSSSM(PLEV)
135
136
137	REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
138	REAL :: DELTA(PLON,PLEV), rh(PLON,PLEV), H
139	REAL, PARAMETER :: RH_MAX=95.
140	REAL :: tau_ae5wv_int ! Intermediate computation of epaisseur optique aerosol
141	REAL :: piz_ae5wv_int ! Intermediate single scattering albedo aerosol
142	REAL :: cg_ae5wv_int ! Intermediate asymmetry parameter aerosol
143	REAL :: fac
144
145	REAL :: tau3d(PLON,PLEV), piz3d(PLON,PLEV), cg3d(PLON,PLEV)
146	REAL :: abs3d(PLON,PLEV) ! epaisseur optique d'absorption
147
148
149	REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble) ! ext. coeff. Soluble comp. units *** m2/g
150	REAL :: alpha_aeri_5wv(las,naero_insoluble) ! ext. coeff. Insoluble comp.
151	REAL :: cg_aers_5wv(nbre_RH,las,naero_soluble) ! Asym. param. soluble comp.
152	REAL :: cg_aeri_5wv(las,naero_insoluble) ! Asym. param. insoluble comp
153	REAL :: piz_aers_5wv(nbre_RH,las,naero_soluble)
154	REAL :: piz_aeri_5wv(las,naero_insoluble) ! Insoluble comp.
155
156	REAL :: fact_RH(nbre_RH)
157	!
158	! Proprietes optiques
159	!
160	REAL :: radry = 287.054 ! dry air mass constant
161	!
162	!
163	!
164	! From here on we look at the optical parameters at 5 wavelengths:
165	! 443nm, 550, 670, 765 and 865 nm
166	! le 12 AVRIL 2006
167	!
168
169	!!ACCUMULATION MODE
170
171	DATA A1_ASSSM/ 4.373E+00, 4.373E+00, 4.361E+00, &
172	4.346E+00, 4.331E+00, 4.278E+00, &
173	4.249E+00, 4.193E+00, 4.162E+00, &
174	4.103E+00, 4.069E+00, 4.035E+00, &
175	3.962E+00, 3.933E+00, 3.904E+00, &
176	3.893E+00, 3.860E+00, 3.847E+00, 3.806E+00, &
177	3.824E+00, 3.780E+00, 3.713E+00, 3.646E+00, &
178	3.547E+00, 3.448E+00, 3.314E+00, 3.179E+00, &
179	3.017E+00, 2.855E+00, 2.855E+00, 2.743E+00, &
180	2.630E+00, 2.630E+00, 2.630E+00, 2.630E+00, &
181	2.630E+00, 2.630E+00, 2.630E+00, 2.630E+00/
182	DATA A2_ASSSM/ 2.496E+00, 2.496E+00, 2.489E+00, &
183	2.481E+00, 2.472E+00, 2.442E+00, &
184	2.427E+00, 2.394E+00, 2.376E+00, &
185	2.342E+00, 2.322E+00, 2.303E+00, &
186	2.261E+00, 2.245E+00, 2.228E+00, &
187	2.222E+00, 2.205E+00, 2.196E+00, 2.189E+00, &
188	2.183E+00, 2.158E+00, 2.120E+00, 2.081E+00, &
189	2.025E+00, 1.968E+00, 1.880E+00, 1.814E+00, &
190	1.720E+00, 1.680E+00, 1.630E+00, 1.501E+00, &
191	1.501E+00, 1.501E+00, 1.501E+00, 1.501E+00, &
192	1.501E+00, 1.501E+00, 1.501E+00, 1.501E+00/
193	DATA A3_ASSSM/-4.688E-02, -4.688E-02, -4.676E-02, &
194	-4.660E-02, -4.644E-02, -4.587E-02, &
195	-4.558E-02, -4.496E-02, -4.463E-02, &
196	-4.399E-02, -4.363E-02, -4.326E-02, &
197	-4.248E-02, -4.217E-02, -4.186E-02, &
198	-4.174E-02, -4.138E-02, -4.125E-02, -4.109E-02, &
199	-4.100E-02, -4.053E-02, -3.982E-02, -3.910E-02, &
200	-3.804E-02, -3.697E-02, -3.553E-02, -3.408E-02, &
201	-3.235E-02, -3.061E-02, -3.061E-02, -2.940E-02, &
202	-2.819E-02, -2.819E-02, -2.819E-02, -2.819E-02, &
203	-2.819E-02, -2.819E-02, -2.819E-02, -2.819E-02/
204	DATA B1_ASSSM/ 1.165E-08, 1.165E-08, 1.145E-08, &
205	1.121E-08, 1.097E-08, 1.012E-08, &
206	9.696E-09, 8.747E-09, 8.261E-09, &
207	7.297E-09, 6.749E-09, 6.201E-09, &
208	5.026E-09, 4.562E-09, 4.098E-09, &
209	3.921E-09, 3.377E-09, 3.187E-09, 2.948E-09, &
210	2.807E-09, 2.291E-09, 2.183E-09, 2.075E-09, &
211	1.916E-09, 1.756E-09, 1.539E-09, 1.322E-09, &
212	1.062E-09, 8.011E-10, 8.011E-10, 6.195E-10, &
213	4.379E-10, 4.379E-10, 4.379E-10, 4.379E-10, &
214	4.379E-10, 4.379E-10, 4.379E-10, 4.379E-10/
215	DATA B2_ASSSM/ 2.193E-08, 2.193E-08, 2.192E-08, &
216	2.190E-08, 2.187E-08, 2.179E-08, &
217	2.175E-08, 2.166E-08, 2.162E-08, &
218	2.153E-08, 2.148E-08, 2.143E-08, &
219	2.132E-08, 2.128E-08, 2.124E-08, &
220	2.122E-08, 2.117E-08, 2.115E-08, 2.113E-08, &
221	2.112E-08, 2.106E-08, 2.103E-08, 2.100E-08, &
222	2.095E-08, 2.090E-08, 2.083E-08, 2.077E-08, &
223	2.069E-08, 2.061E-08, 2.061E-08, 2.055E-08, &
224	2.049E-08, 2.049E-08, 2.049E-08, 2.049E-08, &
225	2.049E-08, 2.049E-08, 2.049E-08, 2.049E-08/
226	DATA C1_ASSSM/ 7.365E-01, 7.365E-01, 7.365E-01, &
227	7.365E-01, 7.365E-01, 7.364E-01, &
228	7.363E-01, 7.363E-01, 7.362E-01, &
229	7.361E-01, 7.360E-01, 7.359E-01, &
230	7.358E-01, 7.357E-01, 7.357E-01, &
231	7.357E-01, 7.357E-01, 7.356E-01, 7.356E-01, &
232	7.356E-01, 7.355E-01, 7.354E-01, 7.354E-01, &
233	7.353E-01, 7.352E-01, 7.351E-01, 7.350E-01, &
234	7.348E-01, 7.347E-01, 7.347E-01, 7.346E-01, &
235	7.345E-01, 7.345E-01, 7.345E-01, 7.345E-01, &
236	7.345E-01, 7.345E-01, 7.345E-01, 7.345E-01/
237	DATA C2_ASSSM/ 5.833E-02, 5.833E-02, 5.835E-02, &
238	5.838E-02, 5.841E-02, 5.850E-02, &
239	5.854E-02, 5.864E-02, 5.870E-02, &
240	5.880E-02, 5.886E-02, 5.891E-02, &
241	5.904E-02, 5.909E-02, 5.914E-02, &
242	5.916E-02, 5.922E-02, 5.924E-02, 5.926E-02, &
243	5.928E-02, 5.934E-02, 5.939E-02, 5.944E-02, &
244	5.952E-02, 5.959E-02, 5.969E-02, 5.979E-02, &
245	5.991E-02, 6.003E-02, 6.003E-02, 6.012E-02, &
246	6.020E-02, 6.020E-02, 6.020E-02, 6.020E-02, &
247	6.020E-02, 6.020E-02, 6.020E-02, 6.020E-02/
248	!COARSE MODE
249	DATA A1_CSSSM/ 7.403E-01, 7.403E-01, 7.422E-01, &
250	7.524E-01, 7.626E-01, 8.019E-01, &
251	8.144E-01, 8.398E-01, 8.527E-01, &
252	8.702E-01, 8.754E-01, 8.806E-01, &
253	8.937E-01, 9.213E-01, 9.489E-01, &
254	9.759E-01, 1.068E+00, 1.105E+00, 1.164E+00, &
255	1.199E+00, 1.359E+00, 1.510E+00, 1.660E+00, &
256	2.100E+00, 2.540E+00, 3.480E+00, 4.421E+00, &
257	3.286E+00, 2.151E+00, 2.151E+00, 1.551E+00, &
258	9.518E-01, 9.518E-01, 9.518E-01, 9.518E-01, &
259	9.518E-01, 9.518E-01, 9.518E-01, 9.518E-01/
260	DATA A2_CSSSM/ 4.522E-01, 4.522E-01, 4.532E-01, &
261	4.588E-01, 4.644E-01, 4.859E-01, &
262	4.928E-01, 5.066E-01, 5.137E-01, &
263	5.233E-01, 5.262E-01, 5.290E-01, &
264	5.361E-01, 5.508E-01, 5.655E-01, &
265	5.798E-01, 6.284E-01, 6.483E-01, 6.795E-01, &
266	6.979E-01, 7.819E-01, 8.654E-01, 9.488E-01, &
267	1.082E+00, 1.450E+00, 1.986E+00, 2.523E+00, &
268	1.876E+00, 1.228E+00, 1.228E+00, 8.856E-01, &
269	5.433E-01, 5.433E-01, 5.433E-01, 5.433E-01, &
270	5.433E-01, 5.433E-01, 5.433E-01, 5.433E-01/
271	DATA A3_CSSSM/-8.516E-03, -8.516E-03, -8.535E-03, &
272	-8.640E-03, -8.744E-03, -9.148E-03, &
273	-9.277E-03, -9.537E-03, -9.668E-03, &
274	-9.848E-03, -9.902E-03, -9.955E-03, &
275	-1.009E-02, -1.036E-02, -1.064E-02, &
276	-1.091E-02, -1.182E-02, -1.219E-02, -1.278E-02, &
277	-1.312E-02, -1.470E-02, -1.626E-02, -1.783E-02, &
278	-2.254E-02, -2.724E-02, -3.732E-02, -4.740E-02, &
279	-3.523E-02, -2.306E-02, -2.306E-02, -1.664E-02, &
280	-1.021E-02, -1.021E-02, -1.021E-02, -1.021E-02, &
281	-1.021E-02, -1.021E-02, -1.021E-02, -1.021E-02/
282	DATA B1_CSSSM/ 2.535E-07, 2.535E-07, 2.530E-07, &
283	2.504E-07, 2.479E-07, 2.380E-07, &
284	2.348E-07, 2.284E-07, 2.252E-07, &
285	2.208E-07, 2.195E-07, 2.182E-07, &
286	2.149E-07, 2.100E-07, 2.051E-07, &
287	2.005E-07, 1.848E-07, 1.784E-07, 1.683E-07, &
288	1.624E-07, 1.353E-07, 1.182E-07, 1.012E-07, &
289	8.068E-08, 6.016E-08, 4.059E-08, 2.102E-08, &
290	1.051E-08, 0.000E+00, 0.000E+00, 0.000E+00, &
291	0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, &
292	0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00/
293	DATA B2_CSSSM/ 1.221E-07, 1.221E-07, 1.217E-07, &
294	1.198E-07, 1.179E-07, 1.104E-07, &
295	1.080E-07, 1.032E-07, 1.008E-07, &
296	9.744E-08, 9.645E-08, 9.546E-08, &
297	9.299E-08, 9.053E-08, 8.807E-08, &
298	8.588E-08, 7.847E-08, 7.544E-08, 7.066E-08, &
299	6.786E-08, 5.504E-08, 4.792E-08, 4.080E-08, &
300	3.520E-08, 2.960E-08, 2.630E-08, 2.300E-08, &
301	2.165E-08, 2.030E-08, 2.030E-08, 2.014E-08, &
302	1.997E-08, 1.997E-08, 1.997E-08, 1.997E-08, &
303	1.997E-08, 1.997E-08, 1.997E-08, 1.997E-08/
304	DATA C1_CSSSM/ 7.659E-01, 7.659E-01, 7.658E-01, &
305	7.655E-01, 7.652E-01, 7.639E-01, &
306	7.635E-01, 7.627E-01, 7.623E-01, &
307	7.618E-01, 7.616E-01, 7.614E-01, &
308	7.610E-01, 7.604E-01, 7.598E-01, &
309	7.592E-01, 7.574E-01, 7.566E-01, 7.553E-01, &
310	7.546E-01, 7.513E-01, 7.492E-01, 7.472E-01, &
311	7.448E-01, 7.423E-01, 7.400E-01, 7.376E-01, &
312	7.359E-01, 7.342E-01, 7.342E-01, 7.338E-01, &
313	7.334E-01, 7.334E-01, 7.334E-01, 7.334E-01, &
314	7.334E-01, 7.334E-01, 7.334E-01, 7.334E-01/
315	DATA C2_CSSSM/ 3.691E-02, 3.691E-02, 3.694E-02, &
316	3.712E-02, 3.729E-02, 3.796E-02, &
317	3.818E-02, 3.861E-02, 3.883E-02, &
318	3.913E-02, 3.922E-02, 3.931E-02, &
319	3.953E-02, 3.994E-02, 4.035E-02, &
320	4.074E-02, 4.208E-02, 4.263E-02, 4.349E-02, &
321	4.400E-02, 4.631E-02, 4.782E-02, 4.933E-02, &
322	5.132E-02, 5.331E-02, 5.332E-02, 5.734E-02, &
323	5.894E-02, 6.053E-02, 6.053E-02, 6.090E-02, &
324	6.128E-02, 6.128E-02, 6.128E-02, 6.128E-02, &
325	6.128E-02, 6.128E-02, 6.128E-02, 6.128E-02/
326	!SUPER COARSE MODE
327	DATA A1_SSSSM/ 2.836E-01, 2.836E-01, 2.876E-01, &
328	2.720E-01, 2.563E-01, 2.414E-01, &
329	2.478E-01, 2.544E-01, 2.546E-01, &
330	2.572E-01, 2.605E-01, 2.638E-01, &
331	2.781E-01, 2.974E-01, 3.167E-01, &
332	3.154E-01, 4.748E-01, 5.286E-01, 6.340E-01, &
333	6.959E-01, 9.233E-01, 1.103E+00, 1.282E+00, &
334	1.559E+00, 1.836E+00, 2.408E+00, 2.981E+00, &
335	3.668E+00, 4.355E+00, 4.355E+00, 4.207E+00, &
336	4.059E+00, 4.059E+00, 4.059E+00, 4.059E+00, &
337	4.059E+00, 4.059E+00, 4.059E+00, 4.059E+00/
338	DATA A2_SSSSM/ 1.608E-01, 1.608E-01, 1.651E-01, &
339	1.614E-01, 1.577E-01, 1.587E-01, &
340	1.636E-01, 1.688E-01, 1.690E-01, &
341	1.711E-01, 1.736E-01, 1.762E-01, &
342	1.874E-01, 2.006E-01, 2.138E-01, &
343	2.342E-01, 3.057E-01, 3.363E-01, 3.940E-01, &
344	4.279E-01, 5.519E-01, 6.470E-01, 7.421E-01, &
345	8.950E-01, 1.048E+00, 1.375E+00, 1.702E+00, &
346	2.094E+00, 2.485E+00, 2.485E+00, 2.401E+00, &
347	2.317E+00, 2.317E+00, 2.317E+00, 2.317E+00, &
348	2.317E+00, 2.317E+00, 2.317E+00, 2.317E+00/
349	DATA A3_SSSSM/-3.025E-03, -3.025E-03, -3.111E-03, &
350	-3.046E-03, -2.981E-03, -3.005E-03, &
351	-3.099E-03, -3.196E-03, -3.200E-03, &
352	-3.239E-03, -3.288E-03, -3.336E-03, &
353	-3.548E-03, -3.798E-03, -4.047E-03, &
354	-4.430E-03, -5.770E-03, -6.345E-03, -7.426E-03, &
355	-8.061E-03, -1.038E-02, -1.216E-02, -1.395E-02, &
356	-1.682E-02, -1.970E-02, -2.584E-02, -3.197E-02, &
357	-3.933E-02, -4.669E-02, -4.669E-02, -4.510E-02, &
358	-4.352E-02, -4.352E-02, -4.352E-02, -4.352E-02, &
359	-4.352E-02, -4.352E-02, -4.352E-02, -4.352E-02/
360	DATA B1_SSSSM/ 6.759E-07, 6.759E-07, 6.246E-07, &
361	5.894E-07, 5.542E-07, 4.953E-07, &
362	4.850E-07, 4.742E-07, 4.738E-07, &
363	4.695E-07, 4.642E-07, 4.588E-07, &
364	4.354E-07, 4.150E-07, 3.947E-07, &
365	3.785E-07, 3.264E-07, 3.067E-07, 2.802E-07, &
366	2.646E-07, 2.095E-07, 1.788E-07, 1.481E-07, &
367	1.192E-07, 9.024E-08, 7.386E-08, 5.747E-08, &
368	4.066E-08, 2.384E-08, 2.384E-08, 1.521E-08, &
369	6.599E-09, 6.599E-09, 6.599E-09, 6.599E-09, &
370	6.599E-09, 6.599E-09, 6.599E-09, 6.599E-09/
371	DATA B2_SSSSM/ 5.977E-07, 5.977E-07, 5.390E-07, &
372	4.929E-07, 4.468E-07, 3.696E-07, &
373	3.570E-07, 3.438E-07, 3.433E-07, &
374	3.380E-07, 3.314E-07, 3.249E-07, &
375	2.962E-07, 2.722E-07, 2.483E-07, &
376	2.318E-07, 1.806E-07, 1.623E-07, 1.423E-07, &
377	1.305E-07, 9.015E-08, 7.563E-08, 6.111E-08, &
378	4.936E-08, 3.761E-08, 3.332E-08, 2.903E-08, &
379	2.620E-08, 2.337E-08, 2.337E-08, 2.242E-08, &
380	2.147E-08, 2.147E-08, 2.147E-08, 2.147E-08, &
381	2.147E-08, 2.147E-08, 2.147E-08, 2.147E-08/
382	DATA C1_SSSSM/ 8.120E-01, 8.120E-01, 8.084E-01, &
383	8.050E-01, 8.016E-01, 7.953E-01, &
384	7.941E-01, 7.928E-01, 7.928E-01, &
385	7.923E-01, 7.916E-01, 7.910E-01, &
386	7.882E-01, 7.858E-01, 7.834E-01, &
387	7.814E-01, 7.750E-01, 7.725E-01, 7.692E-01, &
388	7.673E-01, 7.604E-01, 7.566E-01, 7.529E-01, &
389	7.494E-01, 7.458E-01, 7.438E-01, 7.419E-01, &
390	7.399E-01, 7.379E-01, 7.379E-01, 7.370E-01, &
391	7.360E-01, 7.360E-01, 7.360E-01, 7.360E-01, &
392	7.360E-01, 7.360E-01, 7.360E-01, 7.360E-01/
393	DATA C2_SSSSM/ 2.388E-02, 2.388E-02, 2.392E-02, &
394	2.424E-02, 2.457E-02, 2.552E-02, &
395	2.584E-02, 2.616E-02, 2.618E-02, &
396	2.631E-02, 2.647E-02, 2.663E-02, &
397	2.735E-02, 2.805E-02, 2.875E-02, &
398	2.954E-02, 3.222E-02, 3.330E-02, 3.510E-02, &
399	3.615E-02, 3.997E-02, 4.259E-02, 4.521E-02, &
400	4.780E-02, 5.038E-02, 5.198E-02, 5.358E-02, &
401	5.532E-02, 5.705E-02, 5.705E-02, 5.796E-02, &
402	5.887E-02, 5.358E-02, 5.705E-02, 5.887E-02, &
403	5.038E-02, 5.358E-02, 5.705E-02, 5.887E-02/
404	!! ENDIF
405	!*********************************************************************
406	!
407	DATA alpha_aers_5wv/ &
408	! bc soluble
409	7.930,7.930,7.930,7.930,7.930,7.930, &
410	7.930,7.930,10.893,12.618,14.550,16.613, &
411	7.658,7.658,7.658,7.658,7.658,7.658, &
412	7.658,7.658,10.351,11.879,13.642,15.510, &
413	7.195,7.195,7.195,7.195,7.195,7.195, &
414	7.195,7.195,9.551,10.847,12.381,13.994, &
415	6.736,6.736,6.736,6.736,6.736,6.736, &
416	6.736,6.736,8.818,9.938,11.283,12.687, &
417	6.277,6.277,6.277,6.277,6.277,6.277, &
418	6.277,6.277,8.123,9.094,10.275,11.501, &
419	! pom soluble
420	6.676,6.676,6.676,6.676,6.710,6.934, &
421	7.141,7.569,8.034,8.529,9.456,10.511, &
422	5.109,5.109,5.109,5.109,5.189,5.535, &
423	5.960,6.852,8.008,9.712,12.897,19.676, &
424	3.718,3.718,3.718,3.718,3.779,4.042, &
425	4.364,5.052,5.956,7.314,9.896,15.688, &
426	2.849,2.849,2.849,2.849,2.897,3.107, &
427	3.365,3.916,4.649,5.760,7.900,12.863, &
428	2.229,2.229,2.229,2.229,2.268,2.437, &
429	2.645,3.095,3.692,4.608,6.391,10.633, &
430	! Sulfate Accumulation
431	5.751,6.215,6.690,7.024,7.599,8.195, &
432	9.156,10.355,12.660,14.823,18.908,24.508, &
433	4.320,4.675,5.052,5.375,5.787,6.274, &
434	7.066,8.083,10.088,12.003,15.697,21.133, &
435	3.079,3.351,3.639,3.886,4.205,4.584, &
436	5.206,6.019,7.648,9.234,12.391,17.220, &
437	2.336,2.552,2.781,2.979,3.236,3.540, &
438	4.046,4.711,6.056,7.388,10.093,14.313, &
439	1.777,1.949,2.134,2.292,2.503,2.751, &
440	3.166,3.712,4.828,5.949,8.264,11.922, &
441	! Sulfate Coarse
442	5.751,6.215,6.690,7.024,7.599,8.195, &
443	9.156,10.355,12.660,14.823,18.908,24.508, &
444	4.320,4.675,5.052,5.375,5.787,6.274, &
445	7.066,8.083,10.088,12.003,15.697,21.133, &
446	3.079,3.351,3.639,3.886,4.205,4.584, &
447	5.206,6.019,7.648,9.234,12.391,17.220, &
448	2.336,2.552,2.781,2.979,3.236,3.540, &
449	4.046,4.711,6.056,7.388,10.093,14.313, &
450	1.777,1.949,2.134,2.292,2.503,2.751, &
451	3.166,3.712,4.828,5.949,8.264,11.922, &
452	! seasalt Super Coarse soluble SS
453	0.50,0.90,1.05,1.21,1.40,2.41, &
454	2.66,3.11,3.88,4.52,5.69,8.84, &
455	0.00,0.00,0.00,0.00,0.00,0.00, &
456	0.00,0.00,0.00,0.00,0.00,0.00, &
457	0.52,0.93,1.08,1.24,1.43,2.47, &
458	2.73,3.20,3.99,4.64,5.84,9.04, &
459	0.52,0.93,1.09,1.25,1.44,2.50, &
460	2.76,3.23,4.03,4.68,5.89,9.14, &
461	0.52,0.94,1.09,1.26,1.45,2.51, &
462	2.78,3.25,4.06,4.72,5.94,9.22, &
463	! seasalt soluble CS
464	0.50,0.90,1.05,1.21,1.40,2.41, &
465	2.66,3.11,3.88,4.52,5.69,8.84, &
466	0.00,0.00,0.00,0.00,0.00,0.00, &
467	0.00,0.00,0.00,0.00,0.00,0.00, &
468	0.52,0.93,1.08,1.24,1.43,2.47, &
469	2.73,3.20,3.99,4.64,5.84,9.04, &
470	0.52,0.93,1.09,1.25,1.44,2.50, &
471	2.76,3.23,4.03,4.68,5.89,9.14, &
472	0.52,0.94,1.09,1.26,1.45,2.51, &
473	2.78,3.25,4.06,4.72,5.94,9.22, &
474	! seasalt soluble AS
475	4.28, 7.17, 8.44, 9.85,11.60,22.44, &
476	25.34,30.54,39.38,46.52,59.33,91.77, &
477	0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
478	0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
479	2.48, 4.22, 5.02, 5.94, 7.11,15.29, &
480	17.70,22.31,30.73,38.06,52.15,90.59, &
481	1.90, 3.29, 3.94, 4.69, 5.65, 12.58, &
482	14.68,18.77,26.41,33.25,46.77,85.50, &
483	1.47, 2.59, 3.12, 3.74, 4.54, 10.42, &
484	12.24,15.82,22.66,28.91,41.54,79.33, &
485	! Nitrate Accumulation
486	4.752,4.752,4.752,5.341,5.973,6.648, &
487	7.508,9.390,12.436,15.501,20.153,34.333, &
488	3.232,3.232,3.232,3.636,4.063,4.522, &
489	5.107,6.386,8.458,10.543,13.707,23.351, &
490	2.217,2.217,2.217,2.494,2.787,3.102, &
491	3.503,4.381,5.802,7.232,9.402,16.018, &
492	1.630,1.630,1.630,1.834,2.049,2.280, &
493	2.575,3.221,4.266,5.317,6.913,11.777, &
494	1.234,1.234,1.234,1.388,1.551,1.726, &
495	1.950,2.438,3.229,4.025,5.233, 8.916, &
496	! Nitrate Coarse
497	0.720,0.720,0.720,0.720,0.720,0.720, &
498	0.720,0.720,0.720,0.720,0.720,0.720, &
499	0.729,0.729,0.729,0.729,0.729,0.729, &
500	0.729,0.729,0.729,0.729,0.729,0.729, &
501	0.738,0.738,0.738,0.738,0.738,0.738, &
502	0.738,0.738,0.738,0.738,0.738,0.738, &
503	0.743,0.743,0.743,0.743,0.743,0.743, &
504	0.743,0.743,0.743,0.743,0.743,0.743, &
505	0.746,0.746,0.746,0.746,0.746,0.746, &
506	0.746,0.746,0.746,0.746,0.746,0.746, &
507	! SOA
508	6.676,6.676,6.676,6.676,6.710,6.934, &
509	7.141,7.569,8.034,8.529,9.456,10.511, &
510	5.109,5.109,5.109,5.109,5.189,5.535, &
511	5.960,6.852,8.008,9.712,12.897,19.676, &
512	3.718,3.718,3.718,3.718,3.779,4.042, &
513	4.364,5.052,5.956,7.314,9.896,15.688, &
514	2.849,2.849,2.849,2.849,2.897,3.107, &
515	3.365,3.916,4.649,5.760,7.900,12.863, &
516	2.229,2.229,2.229,2.229,2.268,2.437, &
517	2.645,3.095,3.692,4.608,6.391,10.633, &
518	! SOA
519	6.676,6.676,6.676,6.676,6.710,6.934, &
520	7.141,7.569,8.034,8.529,9.456,10.511, &
521	5.109,5.109,5.109,5.109,5.189,5.535, &
522	5.960,6.852,8.008,9.712,12.897,19.676, &
523	3.718,3.718,3.718,3.718,3.779,4.042, &
524	4.364,5.052,5.956,7.314,9.896,15.688, &
525	2.849,2.849,2.849,2.849,2.897,3.107, &
526	3.365,3.916,4.649,5.760,7.900,12.863, &
527	2.229,2.229,2.229,2.229,2.268,2.437, &
528	2.645,3.095,3.692,4.608,6.391,10.633, &
529	! SOA
530	6.676,6.676,6.676,6.676,6.710,6.934, &
531	7.141,7.569,8.034,8.529,9.456,10.511, &
532	5.109,5.109,5.109,5.109,5.189,5.535, &
533	5.960,6.852,8.008,9.712,12.897,19.676, &
534	3.718,3.718,3.718,3.718,3.779,4.042, &
535	4.364,5.052,5.956,7.314,9.896,15.688, &
536	2.849,2.849,2.849,2.849,2.897,3.107, &
537	3.365,3.916,4.649,5.760,7.900,12.863, &
538	2.229,2.229,2.229,2.229,2.268,2.437, &
539	2.645,3.095,3.692,4.608,6.391,10.633, &
540	! SOA
541	6.676,6.676,6.676,6.676,6.710,6.934, &
542	7.141,7.569,8.034,8.529,9.456,10.511, &
543	5.109,5.109,5.109,5.109,5.189,5.535, &
544	5.960,6.852,8.008,9.712,12.897,19.676, &
545	3.718,3.718,3.718,3.718,3.779,4.042, &
546	4.364,5.052,5.956,7.314,9.896,15.688, &
547	2.849,2.849,2.849,2.849,2.897,3.107, &
548	3.365,3.916,4.649,5.760,7.900,12.863, &
549	2.229,2.229,2.229,2.229,2.268,2.437, &
550	2.645,3.095,3.692,4.608,6.391,10.633 /
551
552	DATA alpha_aeri_5wv/ &
553	! dust insoluble
554	0.759, 0.770, 0.775, 0.775, 0.772, &
555	!!jb bc insoluble
556	11.536,10.033, 8.422, 7.234, 6.270, &
557	! pom insoluble
558	5.042, 3.101, 1.890, 1.294, 0.934, &
559	! nitrate insoluble
560	0.720, 0.729, 0.738, 0.743, 0.746 /
561
562	DATA cg_aers_5wv/ &
563	! bc soluble
564	.651, .651, .651, .651, .651, .651, &
565	.651, .651, .738, .764, .785, .800, &
566	.597, .597, .597, .597, .597, .597, &
567	.597, .597, .695, .725, .751, .770, &
568	.543, .543, .543, .543, .543, .543, &
569	.543, .543, .650, .684, .714, .736, &
570	.504, .504, .504, .504, .504, .504, &
571	.504, .504, .614, .651, .683, .708, &
572	.469, .469, .469, .469, .469, .469, &
573	.469, .469, .582, .620, .655, .681, &
574	! pom soluble
575	.679, .679, .679, .679, .683, .691, &
576	.703, .720, .736, .751, .766, .784, &
577	.656, .656, .656, .656, .659, .669, &
578	.681, .699, .717, .735, .750, .779, &
579	.623, .623, .623, .623, .627, .637, &
580	.649, .668, .688, .709, .734, .762, &
581	.592, .592, .592, .592, .595, .605, &
582	.618, .639, .660, .682, .711, .743, &
583	.561, .561, .561, .561, .565, .575, &
584	.588, .609, .632, .656, .688, .724, &
585	! sulfate Accumulation
586	.671, .684, .697, .704, .714, .723, &
587	.734, .746, .762, .771, .781, .789, &
588	.653, .666, .678, .687, .697, .707, &
589	.719, .732, .751, .762, .775, .789, &
590	.622, .635, .648, .657, .667, .678, &
591	.691, .705, .728, .741, .758, .777, &
592	.591, .604, .617, .627, .638, .650, &
593	.664, .679, .704, .719, .739, .761, &
594	.560, .574, .587, .597, .609, .621, &
595	.637, .653, .680, .697, .719, .745, &
596	! sulfate Coarse
597	.671, .684, .697, .704, .714, .723, &
598	.734, .746, .762, .771, .781, .789, &
599	.653, .666, .678, .687, .697, .707, &
600	.719, .732, .751, .762, .775, .789, &
601	.622, .635, .648, .657, .667, .678, &
602	.691, .705, .728, .741, .758, .777, &
603	.591, .604, .617, .627, .638, .650, &
604	.664, .679, .704, .719, .739, .761, &
605	.560, .574, .587, .597, .609, .621, &
606	.637, .653, .680, .697, .719, .745, &
607	! seasalt Super Coarse soluble (SS)
608	0.730,0.753,0.760,0.766,0.772,0.793, &
609	0.797,0.802,0.809,0.813,0.820,0.830, &
610	0.719,0.744,0.751,0.757,0.764,0.786, &
611	0.791,0.796,0.803,0.808,0.815,0.826, &
612	0.721,0.744,0.750,0.756,0.762,0.784, &
613	0.787,0.793,0.800,0.804,0.811,0.822, &
614	0.717,0.741,0.747,0.753,0.759,0.780, &
615	0.784,0.789,0.795,0.800,0.806,0.817, &
616	0.715,0.739,0.745,0.751,0.757,0.777, &
617	0.781,0.786,0.793,0.797,0.803,0.814, &
618	! seasalt soluble (CS)
619	0.730,0.753,0.760,0.766,0.772,0.793, &
620	0.797,0.802,0.809,0.813,0.820,0.830, &
621	0.719,0.744,0.751,0.757,0.764,0.786, &
622	0.791,0.796,0.803,0.808,0.815,0.826, &
623	0.721,0.744,0.750,0.756,0.762,0.784, &
624	0.787,0.793,0.800,0.804,0.811,0.822, &
625	0.717,0.741,0.747,0.753,0.759,0.780, &
626	0.784,0.789,0.795,0.800,0.806,0.817, &
627	0.715,0.739,0.745,0.751,0.757,0.777, &
628	0.781,0.786,0.793,0.797,0.803,0.814, &
629	! seasalt soluble (AS)
630	0.698,0.722,0.729,0.736,0.743,0.765, &
631	0.768,0.773,0.777,0.779,0.781,0.779, &
632	0.682,0.710,0.719,0.727,0.735,0.764, &
633	0.769,0.776,0.783,0.787,0.791,0.792, &
634	0.658,0.691,0.701,0.710,0.720,0.756, &
635	0.763,0.771,0.782,0.788,0.795,0.801, &
636	0.632,0.668,0.679,0.690,0.701,0.743, &
637	0.750,0.762,0.775,0.783,0.792,0.804, &
638	0.605,0.644,0.656,0.669,0.681,0.729, &
639	0.737,0.750,0.765,0.775,0.787,0.803, &
640	! Nitrate Accumulation
641	.712, .712, .712, .724, .730, .732, &
642	.736, .741, .744, .743, .740, .732, &
643	.685, .685, .685, .695, .704, .708, &
644	.715, .728, .738, .740, .744, .743, &
645	.642, .642, .642, .655, .668, .672, &
646	.684, .703, .721, .728, .735, .747, &
647	.604, .604, .604, .619, .634, .640, &
648	.654, .677, .701, .711, .728, .744, &
649	.566, .566, .566, .584, .601, .606, &
650	.624, .652, .680, .695, .718, .738, &
651	! Nitrate Coarse
652	!!! Nitrate does not lead to hygroscopic growth on this mode
653	!!! it is either with seasalt that is doing the growht
654	!!! or associated with dust and not growing
655	.754, .754, .754, .754, .754, .754, &
656	.754, .754, .754, .754, .754, .754, &
657	.748, .748, .748, .748, .748, .748, &
658	.748, .748, .748, .748, .748, .748, &
659	.741, .741, .741, .741, .741, .741, &
660	.741, .741, .741, .741, .741, .741, &
661	.737, .737, .737, .737, .737, .737, &
662	.737, .737, .737, .737, .737, .737, &
663	.735, .735, .735, .735, .735, .735, &
664	.735, .735, .735, .735, .735, .735, &
665	! SOA
666	.679, .679, .679, .679, .683, .691, &
667	.703, .720, .736, .751, .766, .784, &
668	.656, .656, .656, .656, .659, .669, &
669	.681, .699, .717, .735, .750, .779, &
670	.623, .623, .623, .623, .627, .637, &
671	.649, .668, .688, .709, .734, .762, &
672	.592, .592, .592, .592, .595, .605, &
673	.618, .639, .660, .682, .711, .743, &
674	.561, .561, .561, .561, .565, .575, &
675	.588, .609, .632, .656, .688, .724, &
676	! SOA
677	.679, .679, .679, .679, .683, .691, &
678	.703, .720, .736, .751, .766, .784, &
679	.656, .656, .656, .656, .659, .669, &
680	.681, .699, .717, .735, .750, .779, &
681	.623, .623, .623, .623, .627, .637, &
682	.649, .668, .688, .709, .734, .762, &
683	.592, .592, .592, .592, .595, .605, &
684	.618, .639, .660, .682, .711, .743, &
685	.561, .561, .561, .561, .565, .575, &
686	.588, .609, .632, .656, .688, .724, &
687	! SOA
688	.679, .679, .679, .679, .683, .691, &
689	.703, .720, .736, .751, .766, .784, &
690	.656, .656, .656, .656, .659, .669, &
691	.681, .699, .717, .735, .750, .779, &
692	.623, .623, .623, .623, .627, .637, &
693	.649, .668, .688, .709, .734, .762, &
694	.592, .592, .592, .592, .595, .605, &
695	.618, .639, .660, .682, .711, .743, &
696	.561, .561, .561, .561, .565, .575, &
697	.588, .609, .632, .656, .688, .724, &
698	! SOA
699	.679, .679, .679, .679, .683, .691, &
700	.703, .720, .736, .751, .766, .784, &
701	.656, .656, .656, .656, .659, .669, &
702	.681, .699, .717, .735, .750, .779, &
703	.623, .623, .623, .623, .627, .637, &
704	.649, .668, .688, .709, .734, .762, &
705	.592, .592, .592, .592, .595, .605, &
706	.618, .639, .660, .682, .711, .743, &
707	.561, .561, .561, .561, .565, .575, &
708	.588, .609, .632, .656, .688, .724 /
709	!
710
711	DATA cg_aeri_5wv/&
712	! dust insoluble
713	0.714, 0.697, 0.688, 0.683, 0.679, &
714	! bc insoluble
715	0.511, 0.445, 0.384, 0.342, 0.307, &
716	! pom insoluble
717	0.596, 0.536, 0.466, 0.409, 0.359, &
718	! nitrate insoluble
719	0.754, 0.748, 0.741, 0.737, 0.735 /
720	!
721	DATA piz_aers_5wv/&
722	! bc soluble
723	.445, .445, .445, .445, .445, .445, &
724	.445, .445, .470, .487, .508, .531, &
725	.442, .442, .442, .442, .442, .442, &
726	.442, .442, .462, .481, .506, .533, &
727	.427, .427, .427, .427, .427, .427, &
728	.427, .427, .449, .470, .497, .526, &
729	.413, .413, .413, .413, .413, .413, &
730	.413, .413, .437, .458, .486, .516, &
731	.399, .399, .399, .399, .399, .399, &
732	.399, .399, .423, .445, .473, .506, &
733	! pom soluble
734	.975, .975, .975, .975, .975, .977, &
735	.979, .982, .984, .987, .990, .994, &
736	.972, .972, .972, .972, .973, .974, &
737	.977, .980, .983, .986, .989, .993, &
738	.963, .963, .963, .963, .964, .966, &
739	.969, .974, .977, .982, .986, .991, &
740	.955, .955, .955, .955, .955, .958, &
741	.962, .967, .972, .977, .983, .989, &
742	.944, .944, .944, .944, .944, .948, &
743	.952, .959, .962, .972, .979, .987, &
744	! sulfate Coarse
745	1.000,1.000,1.000,1.000,1.000,1.000, &
746	1.000,1.000,1.000,1.000,1.000,1.000, &
747	1.000,1.000,1.000,1.000,1.000,1.000, &
748	1.000,1.000,1.000,1.000,1.000,1.000, &
749	1.000,1.000,1.000,1.000,1.000,1.000, &
750	1.000,1.000,1.000,1.000,1.000,1.000, &
751	1.000,1.000,1.000,1.000,1.000,1.000, &
752	1.000,1.000,1.000,1.000,1.000,1.000, &
753	1.000,1.000,1.000,1.000,1.000,1.000, &
754	1.000,1.000,1.000,1.000,1.000,1.000, &
755	! sulfate Accumulation
756	1.000,1.000,1.000,1.000,1.000,1.000, &
757	1.000,1.000,1.000,1.000,1.000,1.000, &
758	1.000,1.000,1.000,1.000,1.000,1.000, &
759	1.000,1.000,1.000,1.000,1.000,1.000, &
760	1.000,1.000,1.000,1.000,1.000,1.000, &
761	1.000,1.000,1.000,1.000,1.000,1.000, &
762	1.000,1.000,1.000,1.000,1.000,1.000, &
763	1.000,1.000,1.000,1.000,1.000,1.000, &
764	1.000,1.000,1.000,1.000,1.000,1.000, &
765	1.000,1.000,1.000,1.000,1.000,1.000, &
766	! seasalt soluble (SS)
767	1.000,1.000,1.000,1.000,1.000,1.000, &
768	1.000,1.000,1.000,1.000,1.000,1.000, &
769	1.000,1.000,1.000,1.000,1.000,1.000, &
770	1.000,1.000,1.000,1.000,1.000,1.000, &
771	1.000,1.000,1.000,1.000,1.000,1.000, &
772	1.000,1.000,1.000,1.000,1.000,1.000, &
773	1.000,1.000,1.000,1.000,1.000,1.000, &
774	1.000,1.000,1.000,1.000,1.000,1.000, &
775	1.000,1.000,1.000,1.000,1.000,1.000, &
776	1.000,1.000,1.000,1.000,1.000,1.000, &
777	! seasalt soluble (CS)
778	1.000,1.000,1.000,1.000,1.000,1.000, &
779	1.000,1.000,1.000,1.000,1.000,1.000, &
780	1.000,1.000,1.000,1.000,1.000,1.000, &
781	1.000,1.000,1.000,1.000,1.000,1.000, &
782	1.000,1.000,1.000,1.000,1.000,1.000, &
783	1.000,1.000,1.000,1.000,1.000,1.000, &
784	1.000,1.000,1.000,1.000,1.000,1.000, &
785	1.000,1.000,1.000,1.000,1.000,1.000, &
786	1.000,1.000,1.000,1.000,1.000,1.000, &
787	1.000,1.000,1.000,1.000,1.000,1.000, &
788	! seasalt soluble (AS)
789	1.000,1.000,1.000,1.000,1.000,1.000, &
790	1.000,1.000,1.000,1.000,1.000,1.000, &
791	1.000,1.000,1.000,1.000,1.000,1.000, &
792	1.000,1.000,1.000,1.000,1.000,1.000, &
793	1.000,1.000,1.000,1.000,1.000,1.000, &
794	1.000,1.000,1.000,1.000,1.000,1.000, &
795	1.000,1.000,1.000,1.000,1.000,1.000, &
796	1.000,1.000,1.000,1.000,1.000,1.000, &
797	1.000,1.000,1.000,1.000,1.000,1.000, &
798	1.000,1.000,1.000,1.000,1.000,1.000, &
799	! Nitrate Accumulation
800	1.000,1.000,1.000,1.000,1.000,1.000, &
801	1.000,1.000,1.000,1.000,1.000,1.000, &
802	1.000,1.000,1.000,1.000,1.000,1.000, &
803	1.000,1.000,1.000,1.000,1.000,1.000, &
804	1.000,1.000,1.000,1.000,1.000,1.000, &
805	1.000,1.000,1.000,1.000,1.000,1.000, &
806	1.000,1.000,1.000,1.000,1.000,1.000, &
807	1.000,1.000,1.000,1.000,1.000,1.000, &
808	1.000,1.000,1.000,1.000,1.000,1.000, &
809	1.000,1.000,1.000,1.000,1.000,1.000, &
810	! Nitrate soluble (CS)
811	1.000,1.000,1.000,1.000,1.000,1.000, &
812	1.000,1.000,1.000,1.000,1.000,1.000, &
813	1.000,1.000,1.000,1.000,1.000,1.000, &
814	1.000,1.000,1.000,1.000,1.000,1.000, &
815	1.000,1.000,1.000,1.000,1.000,1.000, &
816	1.000,1.000,1.000,1.000,1.000,1.000, &
817	1.000,1.000,1.000,1.000,1.000,1.000, &
818	1.000,1.000,1.000,1.000,1.000,1.000, &
819	1.000,1.000,1.000,1.000,1.000,1.000, &
820	1.000,1.000,1.000,1.000,1.000,1.000, &
821	! SOA
822	.975, .975, .975, .975, .975, .977, &
823	.979, .982, .984, .987, .990, .994, &
824	.972, .972, .972, .972, .973, .974, &
825	.977, .980, .983, .986, .989, .993, &
826	.963, .963, .963, .963, .964, .966, &
827	.969, .974, .977, .982, .986, .991, &
828	.955, .955, .955, .955, .955, .958, &
829	.962, .967, .972, .977, .983, .989, &
830	.944, .944, .944, .944, .944, .948, &
831	.952, .959, .962, .972, .979, .987, &
832	! SOA
833	.975, .975, .975, .975, .975, .977, &
834	.979, .982, .984, .987, .990, .994, &
835	.972, .972, .972, .972, .973, .974, &
836	.977, .980, .983, .986, .989, .993, &
837	.963, .963, .963, .963, .964, .966, &
838	.969, .974, .977, .982, .986, .991, &
839	.955, .955, .955, .955, .955, .958, &
840	.962, .967, .972, .977, .983, .989, &
841	.944, .944, .944, .944, .944, .948, &
842	.952, .959, .962, .972, .979, .987, &
843	! SOA
844	.975, .975, .975, .975, .975, .977, &
845	.979, .982, .984, .987, .990, .994, &
846	.972, .972, .972, .972, .973, .974, &
847	.977, .980, .983, .986, .989, .993, &
848	.963, .963, .963, .963, .964, .966, &
849	.969, .974, .977, .982, .986, .991, &
850	.955, .955, .955, .955, .955, .958, &
851	.962, .967, .972, .977, .983, .989, &
852	.944, .944, .944, .944, .944, .948, &
853	.952, .959, .962, .972, .979, .987, &
854	! SOA
855	.975, .975, .975, .975, .975, .977, &
856	.979, .982, .984, .987, .990, .994, &
857	.972, .972, .972, .972, .973, .974, &
858	.977, .980, .983, .986, .989, .993, &
859	.963, .963, .963, .963, .964, .966, &
860	.969, .974, .977, .982, .986, .991, &
861	.955, .955, .955, .955, .955, .958, &
862	.962, .967, .972, .977, .983, .989, &
863	.944, .944, .944, .944, .944, .948, &
864	.952, .959, .962, .972, .979, .987 /
865	!
866	DATA piz_aeri_5wv/&
867	! dust insoluble
868	0.944, 0.970, 0.977, 0.982, 0.987, &
869	! bc insoluble
870	0.415, 0.387, 0.355, 0.328, 0.301, &
871	! pom insoluble
872	0.972, 0.963, 0.943, 0.923, 0.897, &
873	! Nitrate insoluble
874	1.000, 1.000, 1.000, 1.000, 1.000 /
875	!
876	! loop over modes, use of precalculated nmd and corresponding sigma
877	! loop over wavelengths
878	! for each mass species in mode
879	! interpolate from Sext to retrieve Sext_at_gridpoint_per_species
880	! compute optical_thickness_at_gridpoint_per_species
881
882	ai(:)=0.
883	tau_ae5wv_int=0.
884	piz_ae5wv_int=0.
885	cg_ae5wv_int=0.
886
887	tausum(:,:,:)=0.
888	tau(:,:,:,:)=0.
889
890	!
891	! Calculations that need to be done since we are not in the subroutines INCA
892	!
893	! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)]
894	zdp1=pdel/(gravit*delt)
895
896	DO RH_Num=1,nbre_RH-1
897	fact_RH(RH_num)=1./(RH_tab(RH_num+1)-RH_tab(RH_num))
898	ENDDO
899
900	DO k=1, PLEV
901	DO i=1, PLON
902	rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
903	RH_num = INT( rh(i,k)/10. + 1.)
904	IF (rh(i,k).GT.85.) RH_num=10
905	IF (rh(i,k).GT.90.) RH_num=11
906	DELTA(i,k)=(rh(i,k)-RH_tab(RH_num))*fact_RH(RH_num)
907	ENDDO
908	ENDDO
909
910
911	DO na=1,nmodes
912	IF (numberindex(na) .EQ. 0) CYCLE
913	nx=numberindex(na)
914	DO m=mass1index(na),massNindex(na) ! tau is only computed for each mass
915	! sum of tau in mode is in (tausum(:,la,nx))
916
917	fac=1.0
918	#ifndef DUSS
919	IF (m.EQ.id_ASBCM) THEN
920	soluble=.TRUE.
921	spsol=1
922	spss=0
923	ELSEIF (m.EQ.id_ASPOMM) THEN
924	soluble=.TRUE.
925	spsol=2
926	spss=0
927	#ifdef NMHC
928	ELSEIF (m.EQ.id_ASARp1a) THEN
929	soluble=.TRUE.
930	spsol=10
931	spss=0
932	ELSEIF (m.EQ.id_ASARp2a) THEN
933	soluble=.TRUE.
934	spsol=11
935	spss=0
936	ELSEIF (m.EQ.id_ASAPp1a) THEN
937	soluble=.TRUE.
938	spsol=12
939	spss=0
940	ELSEIF (m.EQ.id_ASAPp2a) THEN
941	soluble=.TRUE.
942	spsol=13
943	spss=0
944	#endif
945	ELSEIF (m.EQ.id_ASSO4M) THEN
946	soluble=.TRUE.
947	spsol=3
948	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
949	spss=0
950	ELSEIF (m.EQ.id_CSSO4M) THEN
951	soluble=.TRUE.
952	spsol=4
953	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
954	spss=0
955	ELSEIF (m.EQ.id_SSSSM) THEN
956	soluble=.TRUE.
957	spsol=5
958	spss=3
959	ELSEIF (m.EQ.id_CSSSM) THEN
960	soluble=.TRUE.
961	spsol=6
962	spss=2
963	ELSEIF (m.EQ.id_ASSSM) THEN
964	soluble=.TRUE.
965	spsol=7
966	spss=1
967	ELSEIF (m.EQ.id_ASNO3M) THEN
968	soluble=.TRUE.
969	spsol=8
970	spss=0
971	fac=1.29 ! (NH4)NO3/NO3 80/62 mass conversion factor for OD
972	ELSEIF (m.EQ.id_CSNO3M) THEN
973	soluble=.TRUE.
974	spsol=9
975	spss=0
976	ELSEIF (m.EQ.id_CIDUSTM) THEN
977	soluble=.FALSE.
978	spinsol=1
979	spss=0
980	ELSEIF (m.EQ.id_AIBCM) THEN
981	soluble=.FALSE.
982	spinsol=2
983	spss=0
984	ELSEIF (m.EQ.id_AIPOMM) THEN
985	soluble=.FALSE.
986	spinsol=3
987	spss=0
988	ELSEIF (m.EQ.id_CINO3M) THEN
989	soluble=.FALSE.
990	spinsol=4
991	spss=0
992	ELSE
993	CYCLE
994	ENDIF
995	#else
996	IF (m.EQ.id_CIDUSTM) THEN
997	soluble=.FALSE.
998	spinsol=1
999	spss=0
1000	ELSEIF (m.EQ.id_SSSSM) THEN
1001	soluble=.TRUE.
1002	spsol=5
1003	spss=3
1004	ELSEIF (m.EQ.id_CSSSM) THEN
1005	soluble=.TRUE.
1006	spsol=6
1007	spss=2
1008	ELSEIF (m.EQ.id_ASSSM) THEN
1009	soluble=.TRUE.
1010	spsol=7
1011	spss=1
1012	ELSE
1013	CYCLE
1014	ENDIF
1015	#endif
1016
1017
1018	DO la=1,las
1019
1020	IF (soluble) THEN
1021
1022	IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm
1023	IF (spss.EQ.1) THEN !accumulation mode
1024	DO k=1,PLEV
1025	DO i=1,PLON
1026	H=rh(i,k)/100
1027
1028	tau_ae5wv_int = A1_ASSSM(k) + A2_ASSSM(k)*H + A3_ASSSM(k)/(H-1.05)
1029	piz_ae5wv_int = 1 - B1_ASSSM(k) - B2_ASSSM(k)*H
1030	cg_ae5wv_int = C1_ASSSM(k) + C2_ASSSM(k)*H
1031
1032	tau(i,k,la,m) = tr_seri(i,k,m)1000.zdp1(i,k)tau_ae5wv_intdelt*fac
1033	tau(i,k,la,nx) = tau(i,k,la,nx) + tau(i,k,la,m)
1034	tausum(i,la,m)=tausum(i,la,m)+tau(i,k,la,m)
1035	tausum(i,la,nx)=tausum(i,la,nx)+tau(i,k,la,m)
1036	ENDDO
1037	ENDDO
1038	ENDIF
1039
1040	IF (spss.EQ.2) THEN !coarse mode
1041	DO k=1,PLEV
1042	DO i=1,PLON
1043	H=rh(i,k)/100
1044
1045	tau_ae5wv_int = A1_CSSSM(k) + A2_CSSSM(k)*H + A3_CSSSM(k)/(H-1.05)
1046	piz_ae5wv_int = 1 - B1_CSSSM(k) - B2_CSSSM(k)*H
1047	cg_ae5wv_int = C1_CSSSM(k) + C2_CSSSM(k)*H
1048
1049	tau(i,k,la,m) = tr_seri(i,k,m)1000.zdp1(i,k)tau_ae5wv_intdelt*fac
1050	tau(i,k,la,nx) = tau(i,k,la,nx) + tau(i,k,la,m)
1051	tausum(i,la,m)=tausum(i,la,m)+tau(i,k,la,m)
1052	tausum(i,la,nx)=tausum(i,la,nx)+tau(i,k,la,m)
1053	ENDDO
1054	ENDDO
1055	ENDIF
1056	IF (spss.EQ.3) THEN !super coarse mode
1057	DO k=1,PLEV
1058	DO i=1,PLON
1059	H=rh(i,k)/100
1060
1061	tau_ae5wv_int = A1_SSSSM(k) + A2_SSSSM(k)*H + A3_SSSSM(k)/(H-1.05)
1062	piz_ae5wv_int = 1 - B1_SSSSM(k) - B2_SSSSM(k)*H
1063	cg_ae5wv_int = C1_SSSSM(k) + C2_SSSSM(k)*H
1064
1065	tau(i,k,la,m) = tr_seri(i,k,m)1000.zdp1(i,k)tau_ae5wv_intdelt*fac
1066	tau(i,k,la,nx) = tau(i,k,la,nx) + tau(i,k,la,m)
1067	tausum(i,la,m)=tausum(i,la,m)+tau(i,k,la,m)
1068	tausum(i,la,nx)=tausum(i,la,nx)+tau(i,k,la,m)
1069	ENDDO
1070	ENDDO
1071	ENDIF
1072	ELSE ! if la .eq 2 & spss .ne. 0
1073	DO k=1,PLEV
1074	DO i=1,PLON
1075	rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
1076	RH_num = INT( rh(i,k)/10. + 1.)
1077	IF (rh(i,k).GT.85.) RH_num=10
1078	IF (rh(i,k).GT.90.) RH_num=11
1079
1080	tau_ae5wv_int = &
1081	alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
1082	(alpha_aers_5wv(RH_num+1,la,spsol) - &
1083	alpha_aers_5wv(RH_num,la,spsol))
1084
1085	piz_ae5wv_int = &
1086	piz_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
1087	(piz_aers_5wv(RH_num+1,la,spsol) - &
1088	piz_aers_5wv(RH_num,la,spsol))
1089
1090	cg_ae5wv_int = &
1091	cg_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
1092	(cg_aers_5wv(RH_num+1,la,spsol) - &
1093	cg_aers_5wv(RH_num,la,spsol))
1094
1095	tau(i,k,la,m) = tr_seri(i,k,m)1000.zdp1(i,k)tau_ae5wv_intdelt*fac
1096	tau(i,k,la,nx) = tau(i,k,la,nx) + tau(i,k,la,m)
1097	tausum(i,la,m)=tausum(i,la,m)+tau(i,k,la,m)
1098	tausum(i,la,nx)=tausum(i,la,nx)+tau(i,k,la,m)
1099	ENDDO
1100	ENDDO
1101	ENDIF ! if la .eq 2 & spss .ne. 0
1102	ELSE ! if soluble
1103	DO k=1,PLEV
1104	DO i=1,PLON
1105	tau_ae5wv_int = alpha_aeri_5wv(la,spinsol)
1106	piz_ae5wv_int = piz_aeri_5wv(la,spinsol)
1107	cg_ae5wv_int = cg_aeri_5wv(la,spinsol)
1108
1109	tau(i,k,la,m) = tr_seri(i,k,m)1000.zdp1(i,k)tau_ae5wv_intdelt*fac
1110	tau(i,k,la,nx) = tau(i,k,la,nx) + tau(i,k,la,m)
1111	tausum(i,la,m)=tausum(i,la,m)+tau(i,k,la,m)
1112	tausum(i,la,nx)=tausum(i,la,nx)+tau(i,k,la,m)
1113	ENDDO
1114	ENDDO
1115	ENDIF ! if soluble
1116
1117
1118	ENDDO ! boucle sur les longueurs d'onde
1119
1120	ENDDO ! Boucle sur les masses de traceurs
1121	ENDDO ! Boucle sur les modes
1122
1123
1124
1125	! DO k=1, PLEV
1126	! DO i=1, PLON
1127	! !
1128	! ! 5 prochaines lignes à revoir
1129	! ! piz_temp2=piz3d(i,k)
1130	! ! cg_temp2=cg3d(i,k)
1131	! !
1132	! ! piz3d(i,k)=piz_temp2/MAX(tau3d(i,k),1.e-20)
1133	! ! cg3d(i,k)=cg_temp2/MAX(piz_temp,1.e-20)
1134	! !
1135	! ENDDO
1136	! ENDDO
1137
1138	RETURN
1139	END SUBROUTINE OPTAER_5WV
1140	#endif

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: CONFIG/publications/ICOLMDZORINCA_CO2_Transport_GMD_2023/INCA/src/INCA_SRC/optaer2_5wv.F90 @ 6610

Download in other formats: