Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

optaer2_2bands.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: 84.5 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
58
59
60	#ifdef AER
61	SUBROUTINE OPTAER_2BANDS(pdel, tr_seri, delt, RHcl, ai)
62
63	USE INCA_DIM
64	USE SPECIES_NAMES
65	USE CHEM_TRACNM, ONLY: solsym
66	USE AEROSOL_MOD, ONLY: nmodes,mass1index,massNindex,numberindex
67	USE AEROSOL_DIAG , ONLY : lambda, tau_ae, piz_ae, cg_ae,tau_allaer, piz_allaer, cg_allaer,fractnat_allaer, naero_grp
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	IMPLICIT NONE
75
76	!
77	! Arguments:
78	!
79	REAL, DIMENSION(PLON,PLEV), INTENT(in) :: pdel
80	REAL, INTENT(in) :: delt
81	REAL, DIMENSION(PLON,PLEV,PCNST), INTENT(in) :: tr_seri ! Mass Mixing ratio kg/kg
82	REAL, DIMENSION(PLON,PLEV), INTENT(in) :: RHcl ! humidite relative ciel clair
83	REAL, DIMENSION(PLON), INTENT(out) :: ai ! POLDER aerosol index
84
85	!
86	! Local
87	!
88	REAL :: TAUSUM_b1(naero_grp)
89	LOGICAL :: soluble
90	INTEGER :: i, k, inu, m, nx, mrfspecies, na, n
91	INTEGER :: spsol, spinsol, sp, spss
92	INTEGER, PARAMETER :: nbre_RH=12
93	INTEGER, PARAMETER :: naero_soluble=13 ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
94	! 5- seasalt super coarse 6- seasalt coarse 7- seasalt acc.
95	! 8- Nitrate acc. 9- nitrate coarse 10- SOAA1 11-SOAA2 12-SOAB1 13-SOAB2
96	INTEGER, PARAMETER :: naero_insoluble=4 ! 1- Dust; 2- BC insoluble; 3- POM insoluble 4- NO3 insoluble
97	INTEGER :: RH_num(PLON,PLEV)
98	INTEGER :: la443, la550, la670, la765, la865
99	REAL :: A1_ASSSM_b1(PLEV), A2_ASSSM_b1(PLEV), A3_ASSSM_b1(PLEV),&
100	B1_ASSSM_b1(PLEV), B2_ASSSM_b1(PLEV), C1_ASSSM_b1(PLEV), C2_ASSSM_b1(PLEV),&
101	A1_CSSSM_b1(PLEV), A2_CSSSM_b1(PLEV), A3_CSSSM_b1(PLEV),&
102	B1_CSSSM_b1(PLEV), B2_CSSSM_b1(PLEV), C1_CSSSM_b1(PLEV), C2_CSSSM_b1(PLEV),&
103	A1_SSSSM_b1(PLEV), A2_SSSSM_b1(PLEV), A3_SSSSM_b1(PLEV),&
104	B1_SSSSM_b1(PLEV), B2_SSSSM_b1(PLEV), C1_SSSSM_b1(PLEV), C2_SSSSM_b1(PLEV),&
105	A1_ASSSM_b2(PLEV), A2_ASSSM_b2(PLEV), A3_ASSSM_b2(PLEV),&
106	B1_ASSSM_b2(PLEV), B2_ASSSM_b2(PLEV), C1_ASSSM_b2(PLEV), C2_ASSSM_b2(PLEV),&
107	A1_CSSSM_b2(PLEV), A2_CSSSM_b2(PLEV), A3_CSSSM_b2(PLEV),&
108	B1_CSSSM_b2(PLEV), B2_CSSSM_b2(PLEV), C1_CSSSM_b2(PLEV), C2_CSSSM_b2(PLEV),&
109	A1_SSSSM_b2(PLEV), A2_SSSSM_b2(PLEV), A3_SSSSM_b2(PLEV),&
110	B1_SSSSM_b2(PLEV), B2_SSSSM_b2(PLEV), C1_SSSSM_b2(PLEV), C2_SSSSM_b2(PLEV)
111	REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
112	REAL, PARAMETER :: RH_MAX=95.
113	REAL:: DELTA(PLON,PLEV), rh(PLON,PLEV), H
114	REAL:: tau_ae2b_int ! Intermediate computation of epaisseur optique aerosol
115	REAL:: piz_ae2b_int ! Intermediate computation of Single scattering albedo
116	REAL:: cg_ae2b_int ! Intermediate computation of Assymetry parameter
117	REAL:: Fact_RH(nbre_RH)
118	REAL:: taue670(PLON) ! epaisseur optique aerosol absorption 550 nm
119	REAL:: taue865(PLON) ! epaisseur optique aerosol extinction 865 nm
120	REAL:: fac
121	REAL :: tmp_var
122	!
123	! Proprietes optiques
124	!
125	REAL:: alpha_aers_2bands(nbre_RH,2,naero_soluble) !--unit m2/g SO4
126	REAL:: alpha_aeri_2bands(2,naero_insoluble)
127	REAL:: cg_aers_2bands(nbre_RH,2,naero_soluble) !--unit
128	REAL:: cg_aeri_2bands(2,naero_insoluble)
129	REAL:: piz_aers_2bands(nbre_RH,2,naero_soluble) !-- unit
130	REAL:: piz_aeri_2bands(2,naero_insoluble) !-- unit
131	!
132	REAL, PARAMETER :: mintau = 1.e-8
133	REAL, PARAMETER :: minpiz = 1.e-3
134	!
135	!*********************BAND 1*********************************
136	!ACCUMULATION MODE
137	DATA A1_ASSSM_b1/ 4.373E+00, 4.373E+00, 4.361E+00, &
138	4.346E+00, 4.331E+00, 4.278E+00, &
139	4.250E+00, 4.177E+00, 3.132E+00, &
140	4.103E+00, 4.069E+00, 4.035E+00, &
141	3.962E+00, 3.933E+00, 3.904E+00, &
142	3.893E+00, 3.857E+00, 3.847E+00, 3.832E+00, &
143	3.824E+00, 3.780E+00, 3.713E+00, 3.646E+00, &
144	3.547E+00, 3.448E+00, 3.312E+00, 3.179E+00, &
145	3.017E+00, 2.855E+00, 2.855E+00, 2.742E+00, &
146	2.630E+00, 2.630E+00, 2.630E+00, 2.630E+00, &
147	2.630E+00, 2.630E+00, 2.630E+00, 2.630E+00/
148
149	DATA A2_ASSSM_b1/ 2.496E+00, 2.496E+00, 2.489E+00, &
150	2.480E+00, 2.472E+00, 2.442E+00, &
151	2.426E+00, 2.392E+00, 2.372E+00, &
152	2.342E+00, 2.322E+00, 2.303E+00, &
153	2.261E+00, 2.244E+00, 2.228E+00, &
154	2.222E+00, 2.203E+00, 2.196E+00, 2.188E+00, &
155	2.183E+00, 2.158E+00, 2.120E+00, 2.081E+00, &
156	2.024E+00, 1.968E+00, 1.891E+00, 1.814E+00, &
157	1.722E+00, 1.630E+00, 1.630E+00, 1.566E+00, &
158	1.501E+00, 1.501E+00, 1.501E+00, 1.501E+00, &
159	1.501E+00, 1.501E+00, 1.501E+00, 1.501E+00/
160
161	DATA A3_ASSSM_b1/-4.688E-02, -4.688E-02, -4.676E-02, &
162	-4.660E-02, -4.644E-02, -4.587E-02, &
163	-4.558E-02, -4.496E-02, -4.463E-02, &
164	-4.399E-02, -4.362E-02, -4.326E-02, &
165	-4.248E-02, -4.217E-02, -4.186E-02, &
166	-4.174E-02, -4.138E-02, -4.125E-02, -4.112E-02, &
167	-4.100E-02, -4.053E-02, -3.982E-02, -3.910E-02, &
168	-3.804E-02, -3.697E-02, -3.552E-02, -3.408E-02, &
169	-3.234E-02, -3.061E-02, -3.061E-02, -2.940E-02, &
170	-2.819E-02, -2.819E-02, -2.819E-02, -2.819E-02, &
171	-2.819E-02, -2.819E-02, -2.819E-02, -2.819E-02/
172
173	DATA B1_ASSSM_b1/ 1.165E-08, 1.165E-08, 1.145E-08, &
174	1.121E-08, 1.097E-08, 1.012E-08, &
175	9.676E-09, 8.747E-09, 8.261E-09, &
176	7.297E-09, 6.749E-09, 6.201E-09, &
177	5.026E-09, 4.562E-09, 4.098E-09, &
178	3.921E-09, 3.377E-09, 3.187E-09, 2.948E-09, &
179	2.807E-09, 2.291E-09, 2.183E-09, 2.075E-09, &
180	1.916E-09, 1.756E-09, 1.539E-09, 1.322E-09, &
181	1.062E-09, 8.011E-10, 8.011E-10, 6.195E-10, &
182	4.379E-10, 4.379E-10, 4.379E-10, 4.379E-10, &
183	4.379E-10, 4.379E-10, 4.379E-10, 4.379E-10/
184
185	DATA B2_ASSSM_b1/ 2.193E-08, 2.193E-08, 2.192E-08, &
186	2.190E-08, 2.187E-08, 2.179E-08, &
187	2.175E-08, 2.166E-08, 2.162E-08, &
188	2.153E-08, 2.148E-08, 2.143E-08, &
189	2.132E-08, 2.128E-08, 2.124E-08, &
190	2.122E-08, 2.117E-08, 2.115E-08, 2.113E-08, &
191	2.112E-08, 2.106E-08, 2.103E-08, 2.100E-08, &
192	2.095E-08, 2.090E-08, 2.083E-08, 2.077E-08, &
193	2.069E-08, 2.061E-08, 2.061E-08, 2.055E-08, &
194	2.049E-08, 2.049E-08, 2.049E-08, 2.049E-08, &
195	2.049E-08, 2.049E-08, 2.049E-08, 2.049E-08/
196
197	DATA C1_ASSSM_b1/ 7.365E-01, 7.365E-01, 7.365E-01, &
198	7.365E-01, 7.365E-01, 7.364E-01, &
199	7.363E-01, 7.363E-01, 7.362E-01, &
200	7.361E-01, 7.360E-01, 7.359E-01, &
201	7.358E-01, 7.357E-01, 7.357E-01, &
202	7.357E-01, 7.357E-01, 7.356E-01, 7.356E-01, &
203	7.356E-01, 7.355E-01, 7.354E-01, 7.354E-01, &
204	7.353E-01, 7.352E-01, 7.351E-01, 7.350E-01, &
205	7.348E-01, 7.347E-01, 7.347E-01, 7.346E-01, &
206	7.345E-01, 7.345E-01, 7.345E-01, 7.345E-01, &
207	7.345E-01, 7.345E-01, 7.345E-01, 7.345E-01/
208
209	DATA C2_ASSSM_b1/ 5.833E-02, 5.833E-02, 5.835E-02, &
210	5.838E-02, 5.841E-02, 5.850E-02, &
211	5.854E-02, 5.864E-02, 5.870E-02, &
212	5.880E-02, 5.886E-02, 5.891E-02, &
213	5.904E-02, 5.909E-02, 5.914E-02, &
214	5.916E-02, 5.922E-02, 5.924E-02, 5.926E-02, &
215	5.928E-02, 5.934E-02, 5.939E-02, 5.944E-02, &
216	5.952E-02, 5.959E-02, 5.969E-02, 5.979E-02, &
217	5.991E-02, 6.003E-02, 6.003E-02, 6.012E-02, &
218	6.020E-02, 6.020E-02, 6.020E-02, 6.020E-02, &
219	6.020E-02, 6.020E-02, 6.020E-02, 6.020E-02/
220
221	!COARSE MODE
222	DATA A1_CSSSM_b1/ 7.403E-01, 7.403E-01, 7.422E-01, &
223	7.524E-01, 7.626E-01, 8.019E-01, &
224	8.144E-01, 8.398E-01, 8.527E-01, &
225	8.702E-01, 8.754E-01, 8.806E-01, &
226	8.937E-01, 9.213E-01, 9.489E-01, &
227	9.759E-01, 1.068E+00, 1.105E+00, 1.164E+00, &
228	1.199E+00, 1.359E+00, 1.510E+00, 1.660E+00, &
229	2.100E+00, 2.540E+00, 3.480E+00, 4.421E+00, &
230	3.286E+00, 2.151E+00, 2.151E+00, 1.551E+00, &
231	9.518E-01, 9.518E-01, 9.518E-01, 9.518E-01, &
232	9.518E-01, 9.518E-01, 9.518E-01, 9.518E-01/
233
234	DATA A2_CSSSM_b1/ 4.522E-01, 4.522E-01, 4.532E-01, &
235	4.588E-01, 4.644E-01, 4.859E-01, &
236	4.928E-01, 5.066E-01, 5.137E-01, &
237	5.233E-01, 5.262E-01, 5.290E-01, &
238	5.361E-01, 5.508E-01, 5.655E-01, &
239	5.798E-01, 6.284E-01, 6.483E-01, 6.795E-01, &
240	6.979E-01, 7.819E-01, 8.654E-01, 9.488E-01, &
241	1.082E+00, 1.450E+00, 1.986E+00, 2.523E+00, &
242	1.876E+00, 1.228E+00, 1.228E+00, 8.856E-01, &
243	5.433E-01, 5.433E-01, 5.433E-01, 5.433E-01, &
244	5.433E-01, 5.433E-01, 5.433E-01, 5.433E-01/
245
246	DATA A3_CSSSM_b1/-8.516E-03, -8.516E-03, -8.535E-03, &
247	-8.640E-03, -8.744E-03, -9.148E-03, &
248	-9.277E-03, -9.537E-03, -9.668E-03, &
249	-9.848E-03, -9.902E-03, -9.955E-03, &
250	-1.009E-02, -1.036E-02, -1.064E-02, &
251	-1.091E-02, -1.182E-02, -1.219E-02, -1.278E-02, &
252	-1.312E-02, -1.470E-02, -1.626E-02, -1.783E-02, &
253	-2.254E-02, -2.724E-02, -3.732E-02, -4.740E-02, &
254	-3.523E-02, -2.306E-02, -2.306E-02, -1.664E-02, &
255	-1.021E-02, -1.021E-02, -1.021E-02, -1.021E-02, &
256	-1.021E-02, -1.021E-02, -1.021E-02, -1.021E-02/
257
258	DATA B1_CSSSM_b1/ 2.535E-07, 2.535E-07, 2.530E-07, &
259	2.504E-07, 2.479E-07, 2.380E-07, &
260	2.348E-07, 2.284E-07, 2.252E-07, &
261	2.208E-07, 2.195E-07, 2.182E-07, &
262	2.149E-07, 2.100E-07, 2.051E-07, &
263	2.005E-07, 1.848E-07, 1.784E-07, 1.683E-07, &
264	1.624E-07, 1.353E-07, 1.182E-07, 1.012E-07, &
265	8.068E-08, 6.016E-08, 4.059E-08, 2.102E-08, &
266	1.051E-08, 0.000E+00, 0.000E+00, 0.000E+00, &
267	0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, &
268	0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00/
269
270	DATA B2_CSSSM_b1/ 1.221E-07, 1.221E-07, 1.217E-07, &
271	1.198E-07, 1.179E-07, 1.104E-07, &
272	1.080E-07, 1.032E-07, 1.008E-07, &
273	9.744E-08, 9.645E-08, 9.546E-08, &
274	9.299E-08, 9.053E-08, 8.807E-08, &
275	8.588E-08, 7.847E-08, 7.544E-08, 7.066E-08, &
276	6.786E-08, 5.504E-08, 4.792E-08, 4.080E-08, &
277	3.520E-08, 2.960E-08, 2.630E-08, 2.300E-08, &
278	2.165E-08, 2.030E-08, 2.030E-08, 2.014E-08, &
279	1.997E-08, 1.997E-08, 1.997E-08, 1.997E-08, &
280	1.997E-08, 1.997E-08, 1.997E-08, 1.997E-08/
281
282	DATA C1_CSSSM_b1/ 7.659E-01, 7.659E-01, 7.658E-01, &
283	7.655E-01, 7.652E-01, 7.639E-01, &
284	7.635E-01, 7.627E-01, 7.623E-01, &
285	7.618E-01, 7.616E-01, 7.614E-01, &
286	7.610E-01, 7.604E-01, 7.598E-01, &
287	7.592E-01, 7.574E-01, 7.566E-01, 7.553E-01, &
288	7.546E-01, 7.513E-01, 7.492E-01, 7.472E-01, &
289	7.448E-01, 7.423E-01, 7.400E-01, 7.376E-01, &
290	7.359E-01, 7.342E-01, 7.342E-01, 7.338E-01, &
291	7.334E-01, 7.334E-01, 7.334E-01, 7.334E-01, &
292	7.334E-01, 7.334E-01, 7.334E-01, 7.334E-01/
293
294	DATA C2_CSSSM_b1/ 3.691E-02, 3.691E-02, 3.694E-02, &
295	3.712E-02, 3.729E-02, 3.796E-02, &
296	3.818E-02, 3.861E-02, 3.883E-02, &
297	3.913E-02, 3.922E-02, 3.931E-02, &
298	3.953E-02, 3.994E-02, 4.035E-02, &
299	4.074E-02, 4.208E-02, 4.263E-02, 4.349E-02, &
300	4.400E-02, 4.631E-02, 4.782E-02, 4.933E-02, &
301	5.132E-02, 5.331E-02, 5.332E-02, 5.734E-02, &
302	5.894E-02, 6.053E-02, 6.053E-02, 6.090E-02, &
303	6.128E-02, 6.128E-02, 6.128E-02, 6.128E-02, &
304	6.128E-02, 6.128E-02, 6.128E-02, 6.128E-02/
305
306	!SUPER COARSE MODE
307	DATA A1_SSSSM_b1/ 2.836E-01, 2.836E-01, 2.876E-01, &
308	2.720E-01, 2.563E-01, 2.414E-01, &
309	2.478E-01, 2.544E-01, 2.546E-01, &
310	2.572E-01, 2.605E-01, 2.638E-01, &
311	2.781E-01, 2.974E-01, 3.167E-01, &
312	3.154E-01, 4.748E-01, 5.286E-01, 6.340E-01, &
313	6.959E-01, 9.233E-01, 1.103E+00, 1.282E+00, &
314	1.559E+00, 1.836E+00, 2.408E+00, 2.981E+00, &
315	3.668E+00, 4.355E+00, 4.355E+00, 4.207E+00, &
316	4.059E+00, 4.059E+00, 4.059E+00, 4.059E+00, &
317	4.059E+00, 4.059E+00, 4.059E+00, 4.059E+00/
318
319	DATA A2_SSSSM_b1/ 1.608E-01, 1.608E-01, 1.651E-01, &
320	1.614E-01, 1.577E-01, 1.587E-01, &
321	1.636E-01, 1.688E-01, 1.690E-01, &
322	1.711E-01, 1.736E-01, 1.762E-01, &
323	1.874E-01, 2.006E-01, 2.138E-01, &
324	2.342E-01, 3.057E-01, 3.363E-01, 3.940E-01, &
325	4.279E-01, 5.519E-01, 6.470E-01, 7.421E-01, &
326	8.950E-01, 1.048E+00, 1.375E+00, 1.702E+00, &
327	2.094E+00, 2.485E+00, 2.485E+00, 2.401E+00, &
328	2.317E+00, 2.317E+00, 2.317E+00, 2.317E+00, &
329	2.317E+00, 2.317E+00, 2.317E+00, 2.317E+00/
330
331	DATA A3_SSSSM_b1/-3.025E-03, -3.025E-03, -3.111E-03, &
332	-3.046E-03, -2.981E-03, -3.005E-03, &
333	-3.099E-03, -3.196E-03, -3.200E-03, &
334	-3.239E-03, -3.288E-03, -3.336E-03, &
335	-3.548E-03, -3.798E-03, -4.047E-03, &
336	-4.430E-03, -5.770E-03, -6.345E-03, -7.426E-03, &
337	-8.061E-03, -1.038E-02, -1.216E-02, -1.395E-02, &
338	-1.682E-02, -1.970E-02, -2.584E-02, -3.197E-02, &
339	-3.933E-02, -4.669E-02, -4.669E-02, -4.510E-02, &
340	-4.352E-02, -4.352E-02, -4.352E-02, -4.352E-02, &
341	-4.352E-02, -4.352E-02, -4.352E-02, -4.352E-02/
342
343	DATA B1_SSSSM_b1/ 6.759E-07, 6.759E-07, 6.246E-07, &
344	5.894E-07, 5.542E-07, 4.953E-07, &
345	4.850E-07, 4.742E-07, 4.738E-07, &
346	4.695E-07, 4.642E-07, 4.588E-07, &
347	4.354E-07, 4.150E-07, 3.947E-07, &
348	3.785E-07, 3.264E-07, 3.067E-07, 2.802E-07, &
349	2.646E-07, 2.095E-07, 1.788E-07, 1.481E-07, &
350	1.192E-07, 9.024E-08, 7.386E-08, 5.747E-08, &
351	4.066E-08, 2.384E-08, 2.384E-08, 1.521E-08, &
352	6.599E-09, 6.599E-09, 6.599E-09, 6.599E-09, &
353	6.599E-09, 6.599E-09, 6.599E-09, 6.599E-09/
354
355	DATA B2_SSSSM_b1/ 5.977E-07, 5.977E-07, 5.390E-07, &
356	4.929E-07, 4.468E-07, 3.696E-07, &
357	3.570E-07, 3.438E-07, 3.433E-07, &
358	3.380E-07, 3.314E-07, 3.249E-07, &
359	2.962E-07, 2.722E-07, 2.483E-07, &
360	2.318E-07, 1.806E-07, 1.623E-07, 1.423E-07, &
361	1.305E-07, 9.015E-08, 7.563E-08, 6.111E-08, &
362	4.936E-08, 3.761E-08, 3.332E-08, 2.903E-08, &
363	2.620E-08, 2.337E-08, 2.337E-08, 2.242E-08, &
364	2.147E-08, 2.147E-08, 2.147E-08, 2.147E-08, &
365	2.147E-08, 2.147E-08, 2.147E-08, 2.147E-08/
366
367	DATA C1_SSSSM_b1/ 8.120E-01, 8.120E-01, 8.084E-01, &
368	8.050E-01, 8.016E-01, 7.953E-01, &
369	7.941E-01, 7.928E-01, 7.928E-01, &
370	7.923E-01, 7.916E-01, 7.910E-01, &
371	7.882E-01, 7.858E-01, 7.834E-01, &
372	7.814E-01, 7.750E-01, 7.725E-01, 7.692E-01, &
373	7.673E-01, 7.604E-01, 7.566E-01, 7.529E-01, &
374	7.494E-01, 7.458E-01, 7.438E-01, 7.419E-01, &
375	7.399E-01, 7.379E-01, 7.379E-01, 7.370E-01, &
376	7.360E-01, 7.360E-01, 7.360E-01, 7.360E-01, &
377	7.360E-01, 7.360E-01, 7.360E-01, 7.360E-01/
378
379	DATA C2_SSSSM_b1/ 2.388E-02, 2.388E-02, 2.392E-02, &
380	2.424E-02, 2.457E-02, 2.552E-02, &
381	2.584E-02, 2.616E-02, 2.618E-02, &
382	2.631E-02, 2.647E-02, 2.663E-02, &
383	2.735E-02, 2.805E-02, 2.875E-02, &
384	2.954E-02, 3.222E-02, 3.330E-02, 3.510E-02, &
385	3.615E-02, 3.997E-02, 4.259E-02, 4.521E-02, &
386	4.780E-02, 5.038E-02, 5.198E-02, 5.358E-02, &
387	5.532E-02, 5.705E-02, 5.705E-02, 5.796E-02, &
388	5.887E-02, 5.358E-02, 5.705E-02, 5.887E-02, &
389	5.038E-02, 5.358E-02, 5.705E-02, 5.887E-02/
390
391	!*******************BAND 2**********************************************
392	!ACCUMULATION MODE
393	DATA A1_ASSSM_b2/1.256E+00, 1.256E+00, 1.246E+00, &
394	1.236E+00, 1.226E+00, 1.187E+00, &
395	1.168E+00, 1.126E+00, 1.105E+00, &
396	1.062E+00, 1.038E+00, 1.014E+00, &
397	9.616E-01, 9.410E-01, 9.205E-01, &
398	9.127E-01, 8.885E-01, 8.800E-01, 8.694E-01, &
399	8.632E-01, 8.371E-01, 8.157E-01, 7.943E-01, &
400	7.626E-01, 7.308E-01, 6.878E-01, 6.448E-01, &
401	5.931E-01, 5.414E-01, 5.414E-01, 5.054E-01, &
402	4.693E-01, 4.693E-01, 4.693E-01, 4.693E-01, &
403	4.693E-01, 4.693E-01, 4.693E-01, 4.693E-01/
404
405	DATA A2_ASSSM_b2/5.321E-01, 5.321E-01, 5.284E-01, &
406	5.240E-01, 5.196E-01, 5.036E-01, &
407	4.954E-01, 4.782E-01, 4.691E-01, &
408	4.512E-01, 4.410E-01, 4.308E-01, &
409	4.089E-01, 4.003E-01, 3.917E-01, &
410	3.884E-01, 3.782E-01, 3.747E-01, 3.702E-01, &
411	3.676E-01, 3.567E-01, 3.476E-01, 3.385E-01, &
412	3.250E-01, 3.116E-01, 2.934E-01, 2.751E-01, &
413	2.532E-01, 2.312E-01, 2.312E-01, 2.159E-01, &
414	2.006E-01, 2.006E-01, 2.006E-01, 2.006E-01, &
415	2.006E-01, 2.006E-01, 2.006E-01, 2.006E-01/
416
417	DATA A3_ASSSM_b2/-1.053E-02, -1.053E-02, -1.046E-02, &
418	-1.037E-02, -1.028E-02, -9.964E-03, &
419	-9.800E-03, -9.458E-03, -9.279E-03, &
420	-8.923E-03, -8.720E-03, -8.518E-03, &
421	-8.084E-03, -7.912E-03, -7.741E-03, &
422	-7.676E-03, -7.475E-03, -7.405E-03, -7.317E-03, &
423	-7.265E-03, -7.048E-03, -6.868E-03, -6.687E-03, &
424	-6.422E-03, -6.156E-03, -5.794E-03, -5.433E-03, &
425	-4.999E-03, -4.565E-03, -4.565E-03, -4.263E-03, &
426	-3.961E-03, -3.961E-03, -3.961E-03, -3.961E-03, &
427	-3.961E-03, -3.961E-03, -3.961E-03, -3.961E-03/
428
429	DATA B1_ASSSM_b2/1.560E-02, 1.560E-02, 1.560E-02, &
430	1.560E-02, 1.561E-02, 1.565E-02, &
431	1.566E-02, 1.570E-02, 1.572E-02, &
432	1.576E-02, 1.578E-02, 1.580E-02, &
433	1.584E-02, 1.586E-02, 1.588E-02, &
434	1.589E-02, 1.591E-02, 1.592E-02, 1.593E-02, &
435	1.593E-02, 1.595E-02, 1.597E-02, 1.599E-02, &
436	1.602E-02, 1.605E-02, 1.608E-02, 1.612E-02, &
437	1.616E-02, 1.621E-02, 1.621E-02, 1.624E-02, &
438	1.627E-02, 1.627E-02, 1.627E-02, 1.627E-02, &
439	1.627E-02, 1.627E-02, 1.627E-02, 1.627E-02/
440
441	DATA B2_ASSSM_b2/1.073E-02, 1.073E-02, 1.074E-02, &
442	1.075E-02, 1.076E-02, 1.079E-02, &
443	1.080E-02, 1.084E-02, 1.085E-02, &
444	1.089E-02, 1.091E-02, 1.093E-02, &
445	1.097E-02, 1.098E-02, 1.100E-02, &
446	1.101E-02, 1.102E-02, 1.103E-02, 1.104E-02, &
447	1.105E-02, 1.107E-02, 1.108E-02, 1.110E-02, &
448	1.112E-02, 1.115E-02, 1.118E-02, 1.122E-02, &
449	1.126E-02, 1.130E-02, 1.130E-02, 1.133E-02, &
450	1.136E-02, 1.136E-02, 1.136E-02, 1.136E-02, &
451	1.136E-02, 1.136E-02, 1.136E-02, 1.136E-02/
452
453	DATA C1_ASSSM_b2/7.429E-01, 7.429E-01, 7.429E-01, &
454	7.429E-01, 7.429E-01, 7.427E-01, &
455	7.427E-01, 7.426E-01, 7.424E-01, &
456	7.423E-01, 7.422E-01, 7.422E-01, &
457	7.421E-01, 7.420E-01, 7.420E-01, &
458	7.420E-01, 7.419E-01, 7.419E-01, 7.418E-01, &
459	7.418E-01, 7.417E-01, 7.416E-01, 7.416E-01, &
460	7.416E-01, 7.415E-01, 7.414E-01, 7.413E-01, &
461	7.411E-01, 7.409E-01, 7.409E-01, 7.408E-01, &
462	7.408E-01, 7.408E-01, 7.408E-01, 7.408E-01, &
463	7.408E-01, 7.408E-01, 7.408E-01, 7.408E-01/
464
465
466	DATA C2_ASSSM_b2/3.031E-02, 3.031E-02, 3.028E-02, &
467	3.025E-02, 3.022E-02, 3.011E-02, &
468	3.005E-02, 2.992E-02, 2.986E-02, &
469	2.973E-02, 2.966E-02, 2.959E-02, &
470	2.943E-02, 2.937E-02, 2.931E-02, &
471	2.929E-02, 2.922E-02, 2.919E-02, 2.915E-02, &
472	2.913E-02, 2.905E-02, 2.899E-02, 2.893E-02, &
473	2.884E-02, 2.874E-02, 2.860E-02, 2.847E-02, &
474	2.832E-02, 2.817E-02, 2.817E-02, 2.806E-02, &
475	2.795E-02, 2.795E-02, 2.795E-02, 2.795E-02, &
476	2.795E-02, 2.795E-02, 2.795E-02, 2.795E-02/
477
478	!COARSE MODE
479	DATA A1_CSSSM_b2/7.061E-01, 7.061E-01, 7.074E-01, &
480	7.142E-01, 7.211E-01, 7.476E-01, &
481	7.562E-01, 7.732E-01, 7.817E-01, &
482	7.937E-01, 7.972E-01, 8.007E-01, &
483	8.095E-01, 8.266E-01, 8.436E-01, &
484	8.601E-01, 9.161E-01, 9.390E-01, 9.751E-01, &
485	9.963E-01, 1.093E+00, 1.174E+00, 1.256E+00, &
486	1.462E+00, 1.668E+00, 1.624E+00, 1.581E+00, &
487	9.634E-01, 3.457E-01, 3.457E-01, 2.394E-01, &
488	1.331E-01, 1.331E-01, 1.331E-01, 1.331E-01, &
489	1.331E-01, 1.331E-01, 1.331E-01, 1.331E-01/
490
491	DATA A2_CSSSM_b2/3.617E-01, 3.617E-01, 3.621E-01, &
492	3.642E-01, 3.662E-01, 3.739E-01, &
493	3.764E-01, 3.814E-01, 3.840E-01, &
494	3.874E-01, 3.884E-01, 3.895E-01, &
495	3.921E-01, 3.961E-01, 4.001E-01, &
496	4.040E-01, 4.170E-01, 4.223E-01, 4.307E-01, &
497	4.356E-01, 4.581E-01, 4.840E-01, 5.099E-01, &
498	5.965E-01, 6.831E-01, 6.747E-01, 6.663E-01, &
499	4.072E-01, 1.481E-01, 1.481E-01, 1.026E-01, &
500	5.703E-02, 5.703E-02, 5.703E-02, 5.703E-02, &
501	5.703E-02, 5.703E-02, 5.703E-02, 5.703E-02/
502
503	DATA A3_CSSSM_b2/-6.953E-03, -6.953E-03, -6.961E-03, &
504	-7.004E-03, -7.048E-03, -7.216E-03, &
505	-7.269E-03, -7.376E-03, -7.431E-03, &
506	-7.506E-03, -7.528E-03, -7.551E-03, &
507	-7.606E-03, -7.698E-03, -7.791E-03, &
508	-7.880E-03, -8.182E-03, -8.305E-03, -8.499E-03, &
509	-8.613E-03, -9.134E-03, -9.682E-03, -1.023E-02, &
510	-1.194E-02, -1.365E-02, -1.342E-02, -1.320E-02, &
511	-8.061E-03, -2.922E-03, -2.922E-03, -2.023E-03, &
512	-1.125E-03, -1.125E-03, -1.125E-03, -1.125E-03, &
513	-1.125E-03, -1.125E-03, -1.125E-03, -1.125E-03/
514
515	DATA B1_CSSSM_b2/1.007E-02, 1.007E-02, 1.008E-02, &
516	1.010E-02, 1.012E-02, 1.019E-02, &
517	1.022E-02, 1.026E-02, 1.029E-02, &
518	1.033E-02, 1.034E-02, 1.035E-02, &
519	1.038E-02, 1.047E-02, 1.056E-02, &
520	1.065E-02, 1.096E-02, 1.109E-02, 1.128E-02, &
521	1.140E-02, 1.194E-02, 1.232E-02, 1.270E-02, &
522	1.330E-02, 1.390E-02, 1.457E-02, 1.524E-02, &
523	1.582E-02, 1.639E-02, 1.639E-02, 1.653E-02, &
524	1.667E-02, 1.667E-02, 1.667E-02, 1.667E-02, &
525	1.667E-02, 1.667E-02, 1.667E-02, 1.667E-02/
526
527	DATA B2_CSSSM_b2/4.675E-03, 4.675E-03, 4.782E-03, &
528	4.721E-03, 4.760E-03, 4.908E-03, &
529	4.956E-03, 5.053E-03, 5.102E-03, &
530	5.168E-03, 5.188E-03, 5.207E-03, &
531	5.256E-03, 5.365E-03, 5.474E-03, &
532	5.580E-03, 5.941E-03, 6.089E-03, 6.321E-03, &
533	6.457E-03, 7.081E-03, 7.472E-03, 7.923E-03, &
534	8.525E-03, 9.127E-03, 9.768E-03, 1.041E-02, &
535	1.094E-02, 1.147E-02, 1.147E-02, 1.160E-02, &
536	1.173E-02, 1.173E-02, 1.173E-02, 1.173E-02, &
537	1.173E-02, 1.173E-02, 1.173E-02, 1.173E-02/
538
539	DATA C1_CSSSM_b2/7.571E-01, 7.571E-01, 7.571E-01, &
540	7.570E-01, 7.570E-01, 7.568E-01, &
541	7.566E-01, 7.564E-01, 7.564E-01, &
542	7.563E-01, 7.562E-01, 7.562E-01, &
543	7.562E-01, 7.560E-01, 7.557E-01, &
544	7.555E-01, 7.548E-01, 7.545E-01, 7.541E-01, &
545	7.539E-01, 7.527E-01, 7.518E-01, 7.509E-01, &
546	7.494E-01, 7.478E-01, 7.459E-01, 7.440E-01, &
547	7.422E-01, 7.404E-01, 7.404E-01, 7.399E-01, &
548	7.394E-01, 7.394E-01, 7.394E-01, 7.394E-01, &
549	7.394E-01, 7.394E-01, 7.394E-01, 7.394E-01/
550
551	DATA C2_CSSSM_b2/4.464E-02, 4.464E-02, 4.465E-02, &
552	4.466E-02, 4.468E-02, 4.474E-02, &
553	4.476E-02, 4.478E-02, 4.480E-02, &
554	4.482E-02, 4.483E-02, 4.484E-02, &
555	4.486E-02, 4.467E-02, 4.448E-02, &
556	4.428E-02, 4.362E-02, 4.334E-02, 4.290E-02, &
557	4.264E-02, 4.148E-02, 4.052E-02, 3.957E-02, &
558	3.772E-02, 3.588E-02, 3.368E-02, 3.149E-02, &
559	2.950E-02, 2.751E-02, 2.751E-02, 2.700E-02, &
560	2.650E-02, 2.650E-02, 2.650E-02, 2.650E-02, &
561	2.650E-02, 2.650E-02, 2.650E-02, 2.650E-02/
562
563	!SUPER COARSE MODE
564	DATA A1_SSSSM_b2/2.357E-01, 2.357E-01, 2.490E-01, &
565	2.578E-01, 2.666E-01, 2.920E-01, &
566	3.020E-01, 3.124E-01, 3.128E-01, &
567	3.169E-01, 3.220E-01, 3.272E-01, &
568	3.498E-01, 3.729E-01, 3.960E-01, &
569	4.247E-01, 5.228E-01, 5.634E-01, 6.345E-01, &
570	6.763E-01, 8.278E-01, 9.374E-01, 1.047E+00, &
571	1.194E+00, 1.340E+00, 1.634E+00, 1.927E+00, &
572	1.788E+00, 1.648E+00, 1.648E+00, 1.340E+00, &
573	1.031E+00, 1.031E+00, 1.031E+00, 1.031E+00, &
574	1.031E+00, 1.031E+00, 1.031E+00, 1.031E+00/
575
576	DATA A2_SSSSM_b2/1.219E-01, 1.219E-01, 1.337E-01, &
577	1.485E-01, 1.633E-01, 1.929E-01, &
578	1.993E-01, 2.060E-01, 2.062E-01, &
579	2.089E-01, 2.122E-01, 2.155E-01, &
580	2.300E-01, 2.430E-01, 2.560E-01, &
581	2.676E-01, 3.054E-01, 3.199E-01, 3.408E-01, &
582	3.530E-01, 3.965E-01, 4.220E-01, 4.475E-01, &
583	4.959E-01, 4.443E-01, 6.693E-01, 7.943E-01, &
584	7.436E-01, 6.928E-01, 6.928E-01, 5.654E-01, &
585	4.381E-01, 4.381E-01, 4.381E-01, 4.381E-01, &
586	4.381E-01, 4.381E-01, 4.381E-01, 4.381E-01/
587
588	DATA A3_SSSSM_b2/-2.387E-03, -2.387E-03, -2.599E-03, &
589	-2.846E-03, -3.092E-03, -3.599E-03, &
590	-3.716E-03, -3.837E-03, -3.842E-03, &
591	-3.890E-03, -3.951E-03, -4.012E-03, &
592	-4.276E-03, -4.520E-03, -4.763E-03, &
593	-4.994E-03, -5.753E-03, -6.051E-03, -6.500E-03, &
594	-6.763E-03, -7.708E-03, -8.298E-03, -8.887E-03, &
595	-9.898E-03, -1.091E-02, -1.338E-02, -1.585E-02, &
596	-1.479E-02, -1.373E-02, -1.373E-02, -1.120E-02, &
597	-8.665E-03, -8.665E-03, -8.665E-03, -8.665E-03, &
598	-8.665E-03, -8.665E-03, -8.665E-03, -8.665E-03/
599
600	DATA B1_SSSSM_b2/1.260E-02, 1.260E-02, 1.211E-02, &
601	1.168E-02, 1.126E-02, 1.056E-02, &
602	1.047E-02, 1.038E-02, 1.037E-02, &
603	1.033E-02, 1.028E-02, 1.023E-02, &
604	1.002E-02, 9.868E-03, 9.717E-03, &
605	9.682E-03, 9.632E-03, 9.652E-03, 9.860E-03, &
606	9.983E-03, 1.047E-02, 1.108E-02, 1.168E-02, &
607	1.234E-02, 1.301E-02, 1.350E-02, 1.399E-02, &
608	1.456E-02, 1.514E-02, 1.514E-02, 1.546E-02, &
609	1.578E-02, 1.578E-02, 1.578E-02, 1.578E-02, &
610	1.578E-02, 1.578E-02, 1.578E-02, 1.578E-02/
611
612	DATA B2_SSSSM_b2/2.336E-03, 2.336E-03, 2.419E-03, &
613	2.462E-03, 2.506E-03, 2.610E-03, &
614	2.650E-03, 2.692E-03, 2.694E-03, &
615	2.711E-03, 2.732E-03, 2.752E-03, &
616	2.844E-03, 2.944E-03, 3.043E-03, &
617	3.180E-03, 3.663E-03, 3.871E-03, 4.272E-03, &
618	4.507E-03, 5.373E-03, 6.080E-03, 6.786E-03, &
619	7.512E-03, 8.238E-03, 8.723E-03, 9.208E-03, &
620	9.764E-03, 1.032E-02, 1.032E-02, 1.062E-02, &
621	1.091E-02, 1.091E-02, 1.091E-02, 1.091E-02, &
622	1.091E-02, 1.091E-02, 1.091E-02, 1.091E-02/
623
624	DATA C1_SSSSM_b2/7.832E-01, 7.832E-01, 7.787E-01, &
625	7.754E-01, 7.721E-01, 7.670E-01, &
626	7.664E-01, 7.657E-01, 7.657E-01, &
627	7.654E-01, 7.651E-01, 7.648E-01, &
628	7.634E-01, 7.624E-01, 7.613E-01, &
629	7.607E-01, 7.590E-01, 7.585E-01, 7.578E-01, &
630	7.574E-01, 7.560E-01, 7.546E-01, 7.533E-01, &
631	7.518E-01, 7.502E-01, 7.489E-01, 7.476E-01, &
632	7.460E-01, 7.443E-01, 7.443E-01, 7.433E-01, &
633	7.423E-01, 7.423E-01, 7.423E-01, 7.423E-01, &
634	7.423E-01, 7.423E-01, 7.423E-01, 7.423E-01/
635
636	DATA C2_SSSSM_b2/3.144E-02, 3.144E-02, 3.268E-02, &
637	3.392E-02, 3.515E-02, 3.748E-02, &
638	3.792E-02, 3.838E-02, 3.840E-02, &
639	3.860E-02, 3.883E-02, 3.906E-02, &
640	4.006E-02, 4.090E-02, 4.173E-02, &
641	4.228E-02, 4.386E-02, 4.435E-02, 4.450E-02, &
642	4.459E-02, 4.467E-02, 4.334E-02, 4.202E-02, &
643	4.033E-02, 3.864E-02, 3.712E-02, 3.559E-02, &
644	3.371E-02, 3.183E-02, 3.183E-02, 3.074E-02, &
645	2.964E-02, 2.964E-02, 2.964E-02, 2.964E-02, &
646	2.964E-02, 2.964E-02, 2.964E-02, 2.964E-02/
647	!
648	DATA alpha_aers_2bands/ &
649	! bc soluble
650	7.675,7.675,7.675,7.675,7.675,7.675, &
651	7.675,7.675,10.433,11.984,13.767,15.567,&
652	4.720,4.720,4.720,4.720,4.720,4.720, &
653	4.720,4.720,6.081,6.793,7.567,9.344, &
654	! pom soluble
655	5.503,5.503,5.503,5.503,5.588,5.957, &
656	6.404,7.340,8.545,10.319,13.595,20.398, &
657	1.402,1.402,1.402,1.402,1.431,1.562, &
658	1.715,2.032,2.425,2.991,4.193,7.133, &
659	! sulfate accumulation
660	4.681,5.062,5.460,5.798,6.224,6.733, &
661	7.556,8.613,10.687,12.265,16.32,21.692, &
662	1.107,1.239,1.381,1.490,1.635,1.8030, &
663	2.071,2.407,3.126,3.940,5.539,7.921, &
664	! sulfate Coarse
665	4.681,5.062,5.460,5.798,6.224,6.733, &
666	7.556,8.613,10.687,12.265,16.32,21.692, &
667	1.107,1.239,1.381,1.490,1.635,1.8030, &
668	2.071,2.407,3.126,3.940,5.539,7.921, &
669	! seasalt Super Coarse Soluble (SS)
670	0.5090,0.6554,0.7129,0.7767,0.8529,1.2728, &
671	1.3820,1.5792,1.9173,2.2002,2.7173,4.1487, &
672	0.5167,0.6613,0.7221,0.7868,0.8622,1.3027, &
673	1.4227,1.6317,1.9887,2.2883,2.8356,4.3453, &
674	! seasalt Coarse Soluble (CS)
675	0.5090,0.6554,0.7129,0.7767,0.8529,1.2728, &
676	1.3820,1.5792,1.9173,2.2002,2.7173,4.1487, &
677	0.5167,0.6613,0.7221,0.7868,0.8622,1.3027, &
678	1.4227,1.6317,1.9887,2.2883,2.8356,4.3453, &
679	! seasalt Accumulation Soluble (AS)
680	4.125, 4.674, 5.005, 5.434, 5.985, 10.006, &
681	11.175,13.376,17.264,20.540,26.604, 42.349,&
682	4.187, 3.939, 3.919, 3.937, 3.995, 5.078, &
683	5.511, 6.434, 8.317,10.152,14.024, 26.537, &
684	! nitrate accumulation
685	3.672,3.672,3.672,4.131,4.616,5.137, &
686	5.802,7.256,9.610,11.978,15.573,26.530, &
687	0.807,0.807,0.807,0.908,1.014,1.129, &
688	1.275,1.595,2.112,2.632,3.422, 5.831, &
689	! nitrate Coarse
690	0.725,0.725,0.725,0.816,0.911,1.014, &
691	1.146,1.433,1.897,2.365,3.075,5.238, &
692	0.722,0.722,0.722,0.812,0.908,1.010, &
693	1.141,1.427,1.889,2.355,3.062,5.216, &
694	! SOA
695	5.503,5.503,5.503,5.503,5.588,5.957, &
696	6.404,7.340,8.545,10.319,13.595,20.398, &
697	1.402,1.402,1.402,1.402,1.431,1.562, &
698	1.715,2.032,2.425,2.991,4.193,7.133, &
699	! SOA
700	5.503,5.503,5.503,5.503,5.588,5.957, &
701	6.404,7.340,8.545,10.319,13.595,20.398, &
702	1.402,1.402,1.402,1.402,1.431,1.562, &
703	1.715,2.032,2.425,2.991,4.193,7.133, &
704	! SOA
705	5.503,5.503,5.503,5.503,5.588,5.957, &
706	6.404,7.340,8.545,10.319,13.595,20.398, &
707	1.402,1.402,1.402,1.402,1.431,1.562, &
708	1.715,2.032,2.425,2.991,4.193,7.133, &
709	! SOA
710	5.503,5.503,5.503,5.503,5.588,5.957, &
711	6.404,7.340,8.545,10.319,13.595,20.398, &
712	1.402,1.402,1.402,1.402,1.431,1.562, &
713	1.715,2.032,2.425,2.991,4.193,7.133 /
714	!
715	DATA alpha_aeri_2bands/ &
716	! dust insoluble
717	0.7661,0.7123,&
718	! bc insoluble
719	10.360,4.437, &
720	! pom insoluble
721	3.741,0.606, &
722	! nitrate insoluble
723	0.725,0.722/
724
725	DATA cg_aers_2bands/ &
726	! bc soluble
727	.612, .612, .612, .612, .612, .612, &
728	.612, .612, .702, .734, .760, .796, &
729	.433, .433, .433, .433, .433, .433, &
730	.433, .433, .534, .575, .613, .669, &
731	! pom soluble
732	.663, .663, .663, .663, .666, .674, &
733	.685, .702, .718, .737, .757, .777, &
734	.544, .544, .544, .544, .547, .554, &
735	.565, .583, .604, .631, .661, .698, &
736	! sulfate accumulation
737	.658, .669, .680, .688, .698, .707, &
738	.719, .733, .752, .760, .773, .786, &
739	.544, .555, .565, .573, .583, .593, &
740	.610, .628, .655, .666, .692, .719, &
741	! sulfate coarse
742	.658, .669, .680, .688, .698, .707, &
743	.719, .733, .752, .760, .773, .786, &
744	.544, .555, .565, .573, .583, .593, &
745	.610, .628, .655, .666, .692, .719, &
746	! seasalt Super Coarse soluble (SS)
747	0.727, 0.747, 0.755, 0.761, 0.770, 0.788, &
748	0.792, 0.799, 0.805, 0.809, 0.815, 0.826, &
749	0.717, 0.738, 0.745, 0.752, 0.761, 0.779, &
750	0.781, 0.786, 0.793, 0.797, 0.803, 0.813, &
751	! seasalt Coarse soluble (CS)
752	0.727, 0.747, 0.755, 0.761, 0.770, 0.788, &
753	0.792, 0.799, 0.805, 0.809, 0.815, 0.826, &
754	0.717, 0.738, 0.745, 0.752, 0.761, 0.779, &
755	0.781, 0.786, 0.793, 0.797, 0.803, 0.813, &
756	! Sesalt Accumulation Soluble (AS)
757	0.727, 0.741, 0.748, 0.754, 0.761, 0.782, &
758	0.787, 0.792, 0.797, 0.799, 0.801, 0.799, &
759	0.606, 0.645, 0.658, 0.669, 0.681, 0.726, &
760	0.734, 0.746, 0.761, 0.770, 0.782, 0.798, &
761	! nitrate accumulation
762	.683, .683, .683, .698, .712, .716, &
763	.721, .730, .737, .738, .740, .737, &
764	.518, .518, .518, .550, .581, .586, &
765	.603, .625, .655, .670, .684, .714, &
766	! nitrate coarse
767	.749, .749, .749, .749, .749, .749, &
768	.749, .749, .749, .749, .749, .749, &
769	.731, .731, .731, .731, .731, .731, &
770	.731, .731, .731, .731, .731, .731, &
771	! SOA
772	.663, .663, .663, .663, .666, .674, &
773	.685, .702, .718, .737, .757, .777, &
774	.544, .544, .544, .544, .547, .554, &
775	.565, .583, .604, .631, .661, .698, &
776	! SOA
777	.663, .663, .663, .663, .666, .674, &
778	.685, .702, .718, .737, .757, .777, &
779	.544, .544, .544, .544, .547, .554, &
780	.565, .583, .604, .631, .661, .698, &
781	! SOA
782	.663, .663, .663, .663, .666, .674, &
783	.685, .702, .718, .737, .757, .777, &
784	.544, .544, .544, .544, .547, .554, &
785	.565, .583, .604, .631, .661, .698, &
786	! SOA
787	.663, .663, .663, .663, .666, .674, &
788	.685, .702, .718, .737, .757, .777, &
789	.544, .544, .544, .544, .547, .554, &
790	.565, .583, .604, .631, .661, .698 /
791
792	DATA cg_aeri_2bands/ &
793	! dust insoluble
794	.701, .670, &
795	! bc insoluble
796	.471, .297, &
797	! pom insoluble
798	.568, .365, &
799	!nitrate insoluble
800	.749, .731/
801
802	DATA piz_aers_2bands/&
803	! bc soluble
804	.445, .445, .445, .445, .445, .445, &
805	.445, .445, .461, .480, .505, .528, &
806	.362, .362, .362, .362, .362, .362, &
807	.362, .362, .381, .405, .437, .483, &
808	! pom soluble
809	.972, .972, .972, .972, .972, .974, &
810	.976, .979, .982, .986, .989, .992, &
811	.924, .924, .924, .924, .925, .927, &
812	.932, .938, .945, .952, .961, .970, &
813	! sulfate accumulation
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	.992, .988, .988, .987, .986, .985, &
817	.985, .985, .984, .984, .984, .984, &
818	! sulfate Coarse
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	.992, .988, .988, .987, .986, .985, &
822	.985, .985, .984, .984, .984, .984, &
823	! seasalt Super Coarse soluble (SS)
824	1.000,1.000,1.000,1.000,1.000,1.000, &
825	1.000,1.000,1.000,1.000,1.000,1.000, &
826	0.992,0.989,0.987,0.986,0.986,0.980, &
827	0.980,0.978,0.976,0.976,0.974,0.971, &
828	! seasalt Coarse soluble (CS)
829	1.000,1.000,1.000,1.000,1.000,1.000, &
830	1.000,1.000,1.000,1.000,1.000,1.000, &
831	0.992,0.989,0.987,0.986,0.986,0.980, &
832	0.980,0.978,0.976,0.976,0.974,0.971, &
833	! seasalt Accumulation Soluble (AS)
834	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
835	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
836	0.970, 0.975, 0.976, 0.977, 0.978, 0.982, &
837	0.982, 0.983, 0.984, 0.984, 0.985, 0.985, &
838	! Nitrate accumulation
839	1.000,1.000,1.000,1.000,1.000,1.000, &
840	1.000,1.000,1.000,1.000,1.000,1.000, &
841	.984, .984, .984, .985, .986, .986, &
842	.988, .989, .990, .991, .991, .992, &
843	! Nitrate Coarse
844	1.000,1.000,1.000,1.000,1.000,1.000, &
845	1.000,1.000,1.000,1.000,1.000,1.000, &
846	.988, .988, .988, .988, .988, .988, &
847	.988, .988, .988, .988, .988, .988, &
848	! SOA
849	.972, .972, .972, .972, .972, .974, &
850	.976, .979, .982, .986, .989, .992, &
851	.924, .924, .924, .924, .925, .927, &
852	.932, .938, .945, .952, .961, .970, &
853	! SOA
854	.972, .972, .972, .972, .972, .974, &
855	.976, .979, .982, .986, .989, .992, &
856	.924, .924, .924, .924, .925, .927, &
857	.932, .938, .945, .952, .961, .970, &
858	! SOA
859	.972, .972, .972, .972, .972, .974, &
860	.976, .979, .982, .986, .989, .992, &
861	.924, .924, .924, .924, .925, .927, &
862	.932, .938, .945, .952, .961, .970, &
863	! SOA
864	.972, .972, .972, .972, .972, .974, &
865	.976, .979, .982, .986, .989, .992, &
866	.924, .924, .924, .924, .925, .927, &
867	.932, .938, .945, .952, .961, .970 /
868
869	DATA piz_aeri_2bands/ &
870	! dust insoluble
871	.963, .987, &
872	! bc insoluble
873	.395, .264, &
874	! pom insoluble
875	.966, .859, &
876	! Nitrate Coarse insoluble
877	1.000, .988/
878
879	DO i=1, PLON
880	taue670(i)=0.0
881	taue865(i)=0.0
882	ENDDO
883
884	!
885	! loop over modes, use of precalculated nmd and corresponding sigma
886	! loop over wavelengths
887	! for each mass species in mode
888	! interpolate from Sext to retrieve Sext_at_gridpoint_per_species
889	! compute optical_thickness_at_gridpoint_per_species
890	ai(:)=0.
891	tau_ae(:,:,:,:)=0.
892	piz_ae(:,:,:,:)=1.
893	cg_ae(:,:,:,:)=0.5
894	tau_allaer(:,:,:,:)=0.
895	piz_allaer(:,:,:,:)=1.
896	cg_allaer(:,:,:,:)=0.5
897
898	!
899	! Calculations that need to be done since we are not in the subroutines INCA
900	!
901	! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)]
902	zdp1=pdel/(gravit*delt)
903
904
905	DO n=1,nbre_RH-1
906	fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n))
907	ENDDO
908
909	DO k=1, PLEV
910	DO i=1, PLON
911	rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
912	RH_num(i,k) = INT( rh(i,k)/10. + 1.)
913	IF (rh(i,k).GT.85.) RH_num(i,k)=10
914	IF (rh(i,k).GT.90.) RH_num(i,k)=11
915
916	DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k))
917	ENDDO
918	ENDDO
919
920
921
922	DO na=1,nmodes ! Loop on modes
923	IF (numberindex(na) .EQ. 0) CYCLE
924	nx=numberindex(na)
925	DO m=mass1index(na),massNindex(na) ! tau is only computed for each mass
926	! sum of tau in mode is in (tausum(:,la,nx))
927
928	fac=1.0
929	#ifndef DUSS
930	IF (m.EQ.id_ASBCM) THEN
931	soluble=.TRUE.
932	spsol=1
933	spss=0
934	ELSEIF (m.EQ.id_ASPOMM) THEN
935	soluble=.TRUE.
936	spsol=2
937	spss=0
938	#ifdef NMHC
939	ELSEIF (m.EQ.id_ASARp1a) THEN
940	soluble=.TRUE.
941	spsol=10
942	spss=0
943	ELSEIF (m.EQ.id_ASARp2a) THEN
944	soluble=.TRUE.
945	spsol=11
946	spss=0
947	ELSEIF (m.EQ.id_ASAPp1a) THEN
948	soluble=.TRUE.
949	spsol=12
950	spss=0
951	ELSEIF (m.EQ.id_ASAPp2a) THEN
952	soluble=.TRUE.
953	spsol=13
954	spss=0
955	#endif
956	ELSEIF (m.EQ.id_ASSO4M) THEN
957	soluble=.TRUE.
958	spsol=3
959	spss=0
960	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
961	ELSEIF (m.EQ.id_CSSO4M) THEN
962	soluble=.TRUE.
963	spsol=4
964	spss=0
965	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
966	ELSEIF (m.EQ.id_SSSSM) THEN
967	soluble=.TRUE.
968	spsol=5
969	spss=3
970	ELSEIF (m.EQ.id_CSSSM) THEN
971	soluble=.TRUE.
972	spsol=6
973	spss=2
974	ELSEIF (m.EQ.id_ASSSM) THEN
975	soluble=.TRUE.
976	spsol=7
977	spss=1
978	ELSEIF (m.EQ.id_ASNO3M) THEN
979	soluble=.TRUE.
980	spsol=8
981	spss=0
982	fac=1.29 ! (NH4)NO3/NO3 80/62 mass conversion factor for OD
983	ELSEIF (m.EQ.id_CSNO3M) THEN
984	soluble=.TRUE.
985	spsol=9
986	spss=0
987	ELSEIF (m.EQ.id_CIDUSTM) THEN
988	soluble=.FALSE.
989	spinsol=1
990	spss=0
991	ELSEIF (m.EQ.id_AIBCM) THEN
992	soluble=.FALSE.
993	spinsol=2
994	spss=0
995	ELSEIF (m.EQ.id_AIPOMM) THEN
996	soluble=.FALSE.
997	spinsol=3
998	spss=0
999	ELSEIF (m.EQ.id_CINO3M) THEN
1000	soluble=.FALSE.
1001	spinsol=4
1002	spss=0
1003	ELSE
1004	CYCLE
1005	ENDIF
1006	#else
1007	IF (m.EQ.id_CIDUSTM) THEN
1008	soluble=.FALSE.
1009	spinsol=1
1010	spss=0
1011	ELSEIF (m.EQ.id_SSSSM) THEN
1012	soluble=.TRUE.
1013	spsol=5
1014	spss=3
1015	ELSEIF (m.EQ.id_CSSSM) THEN
1016	soluble=.TRUE.
1017	spsol=6
1018	spss=2
1019	ELSEIF (m.EQ.id_ASSSM) THEN
1020	soluble=.TRUE.
1021	spsol=7
1022	spss=1
1023	ELSE
1024	CYCLE
1025	ENDIF
1026	#endif
1027	!--------------------------------------------------
1028	IF (soluble) THEN
1029	IF (spss .NE. 0) THEN
1030	IF (spss .EQ. 1) THEN ! Accumulation Mode
1031	DO k= 1,PLEV
1032	DO i=1,PLON
1033	H=rh(i,k)/100
1034	tmp_var = tr_seri(i,k,m)1000.zdp1(i,k)deltfac
1035
1036	! band 1
1037	tau_ae2b_int = A1_ASSSM_b1(k) + A2_ASSSM_b1(k)*H + A3_ASSSM_b1(k)/(H-1.05)
1038	piz_ae2b_int = 1 - B1_ASSSM_b1(k) - B2_ASSSM_b1(k)*H
1039	cg_ae2b_int = C1_ASSSM_b1(k) + C2_ASSSM_b1(k)*H
1040
1041	tau_ae(i,k,m,1) = tau_ae(i,k,m,1) + tmp_var*tau_ae2b_int
1042	piz_ae(i,k,m,1) = piz_ae2b_int
1043	cg_ae(i,k,m,1) = cg_ae2b_int
1044
1045	!band 2
1046	tau_ae2b_int = A1_ASSSM_b2(k) + A2_ASSSM_b2(k)*H + A3_ASSSM_b2(k)/(H-1.05)
1047	piz_ae2b_int = 1 - B1_ASSSM_b2(k) - B2_ASSSM_b2(k)*H
1048	cg_ae2b_int = C1_ASSSM_b2(k) + C2_ASSSM_b2(k)*H
1049
1050	tau_ae(i,k,m,2) = tau_ae(i,k,m,2) + tmp_var*tau_ae2b_int
1051	piz_ae(i,k,m,2) = piz_ae2b_int
1052	cg_ae(i,k,m,2) = cg_ae2b_int
1053	ENDDO
1054	ENDDO
1055	ENDIF ! if spss .eq. 1
1056
1057	IF (spss.EQ.2) THEN ! Coarse mode
1058	DO k= 1,PLEV
1059	DO i=1,PLON
1060	H=rh(i,k)/100
1061	tmp_var = tr_seri(i,k,m)1000.zdp1(i,k)deltfac
1062
1063	! band 1
1064	tau_ae2b_int = A1_CSSSM_b1(k) + A2_CSSSM_b1(k)*H + A3_CSSSM_b1(k)/(H-1.05)
1065	piz_ae2b_int = 1 - B1_CSSSM_b1(k) - B2_CSSSM_b1(k)*H
1066	cg_ae2b_int = C1_CSSSM_b1(k) + C2_CSSSM_b1(k)*H
1067
1068	tau_ae(i,k,m,1) = tau_ae(i,k,m,1) + tmp_var*tau_ae2b_int
1069	piz_ae(i,k,m,1) = piz_ae2b_int
1070	cg_ae(i,k,m,1) = cg_ae2b_int
1071
1072
1073	! band 2
1074	tau_ae2b_int = A1_CSSSM_b2(k) + A2_CSSSM_b2(k)*H + A3_CSSSM_b2(k)/(H-1.05)
1075	piz_ae2b_int = 1 - B1_CSSSM_b2(k) - B2_CSSSM_b2(k)*H
1076	cg_ae2b_int = C1_CSSSM_b2(k) + C2_CSSSM_b2(k)*H
1077
1078	tau_ae(i,k,m,2) = tau_ae(i,k,m,2) + tmp_var*tau_ae2b_int
1079	piz_ae(i,k,m,2) = piz_ae2b_int
1080	cg_ae(i,k,m,2) = cg_ae2b_int
1081	ENDDO
1082	ENDDO
1083	ENDIF !if spss .eq. 2
1084
1085	IF (spss.EQ.3) THEN !Super Coarse mode
1086	DO k= 1,PLEV
1087	DO i=1,PLON
1088	H=rh(i,k)/100
1089	tmp_var = tr_seri(i,k,m)1000.zdp1(i,k)deltfac
1090
1091	! band 1
1092	tau_ae2b_int = A1_SSSSM_b1(k) + A2_SSSSM_b1(k)*H + A3_SSSSM_b1(k)/(H-1.05)
1093	piz_ae2b_int = 1 - B1_SSSSM_b1(k) - B2_SSSSM_b1(k)*H
1094	cg_ae2b_int = C1_SSSSM_b1(k) + C2_SSSSM_b1(k)*H
1095
1096	tau_ae(i,k,m,1) = tau_ae(i,k,m,1) + tmp_var*tau_ae2b_int
1097	piz_ae(i,k,m,1) = piz_ae2b_int
1098	cg_ae(i,k,m,1) = cg_ae2b_int
1099
1100	! band 2
1101	tau_ae2b_int = A1_SSSSM_b2(k) + A2_SSSSM_b2(k)*H + A3_SSSSM_b2(k)/(H-1.05)
1102	piz_ae2b_int = 1 - B1_SSSSM_b2(k) - B2_SSSSM_b2(k)*H
1103	cg_ae2b_int = C1_SSSSM_b2(k) + C2_SSSSM_b2(k)*H
1104
1105	tau_ae(i,k,m,2) = tau_ae(i,k,m,2) + tmp_var*tau_ae2b_int
1106	piz_ae(i,k,m,2) = piz_ae2b_int
1107	cg_ae(i,k,m,2) = cg_ae2b_int
1108	ENDDO
1109	ENDDO
1110	ENDIF
1111
1112	ELSE
1113	DO k= 1,PLEV
1114	DO i=1,PLON
1115	tmp_var = tr_seri(i,k,m)1000.zdp1(i,k)deltfac
1116
1117	DO inu=1,2
1118	tau_ae2b_int = alpha_aers_2bands(RH_num(i,k),inu,spsol)+ &
1119	DELTA(i,k)* (alpha_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
1120	alpha_aers_2bands(RH_num(i,k),inu,spsol))
1121
1122	piz_ae2b_int = piz_aers_2bands(RH_num(i,k),inu,spsol) + &
1123	DELTA(i,k)* (piz_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
1124	piz_aers_2bands(RH_num(i,k),inu,spsol))
1125
1126	cg_ae2b_int = cg_aers_2bands(RH_num(i,k),inu,spsol) + &
1127	DELTA(i,k)* (cg_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
1128	cg_aers_2bands(RH_num(i,k),inu,spsol))
1129
1130	tau_ae(i,k,m,inu) = tau_ae(i,k,m,inu) + tmp_var*tau_ae2b_int
1131	piz_ae(i,k,m,inu) = piz_ae2b_int
1132	cg_ae(i,k,m,inu) = cg_ae2b_int
1133	ENDDO
1134
1135	ENDDO
1136	ENDDO
1137	ENDIF ! if spss .ne. 0
1138
1139	ELSE ! if soluble
1140	DO k= 1,PLEV
1141	DO i=1,PLON
1142	tmp_var = tr_seri(i,k,m)1000.zdp1(i,k)deltfac
1143
1144	DO inu=1,2
1145	tau_ae2b_int = alpha_aeri_2bands(inu,spinsol)
1146	piz_ae2b_int = piz_aeri_2bands(inu,spinsol)
1147	cg_ae2b_int = cg_aeri_2bands(inu,spinsol)
1148
1149	tau_ae(i,k,m,inu) = tau_ae(i,k,m,inu) + tmp_var*tau_ae2b_int
1150	piz_ae(i,k,m,inu) = piz_ae2b_int
1151	cg_ae(i,k,m,inu) = cg_ae2b_int
1152	ENDDO
1153
1154	ENDDO
1155	ENDDO
1156	ENDIF ! if soluble
1157
1158	!--------------------------------------------------
1159
1160
1161	ENDDO ! Boucle sur les masses de traceurs
1162	ENDDO ! Boucle sur les modes
1163
1164
1165	#ifndef DUSS
1166	DO inu=1, 2
1167	DO mrfspecies=1,naero_grp
1168	IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER
1169	DO k=1, PLEV
1170	DO i=1, PLON
1171	#ifdef NMHC
1172	tau_allaer(i,k,mrfspecies,inu) = &
1173	tau_ae(i,k,id_ASSO4M,inu)+ tau_ae(i,k,id_CSSO4M,inu)+ &
1174	tau_ae(i,k,id_ASBCM,inu) + tau_ae(i,k,id_AIBCM,inu)+ &
1175	tau_ae(i,k,id_ASPOMM,inu) + tau_ae(i,k,id_AIPOMM,inu)+ &
1176	tau_ae(i,k,id_ASARp1a,inu) + tau_ae(i,k,id_ASARp2a,inu)+ &
1177	tau_ae(i,k,id_ASAPp1a,inu) + tau_ae(i,k,id_ASAPp2a,inu)+ &
1178	tau_ae(i,k,id_ASNO3M,inu) + tau_ae(i,k,id_CSNO3M,inu)+ &
1179	tau_ae(i,k,id_ASSSM,inu) + tau_ae(i,k,id_CSSSM,inu)+ &
1180	tau_ae(i,k,id_SSSSM,inu) + tau_ae(i,k,id_CIDUSTM,inu)+ &
1181	tau_ae(i,k,id_CINO3M,inu)
1182	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1183
1184	piz_allaer(i,k,mrfspecies,inu)= &
1185	(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
1186	tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
1187	tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
1188	tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
1189	tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
1190	tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &
1191	tau_ae(i,k,id_ASARp1a,inu)*piz_ae(i,k,id_ASARp1a,inu)+ &
1192	tau_ae(i,k,id_ASARp2a,inu)*piz_ae(i,k,id_ASARp2a,inu)+ &
1193	tau_ae(i,k,id_ASAPp1a,inu)*piz_ae(i,k,id_ASAPp1a,inu)+ &
1194	tau_ae(i,k,id_ASAPp2a,inu)*piz_ae(i,k,id_ASAPp2a,inu)+ &
1195	tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
1196	tau_ae(i,k,id_ASNO3M,inu)*piz_ae(i,k,id_ASNO3M,inu)+ &
1197	tau_ae(i,k,id_CSNO3M,inu)*piz_ae(i,k,id_CSNO3M,inu)+ &
1198	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
1199	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
1200	tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)+ &
1201	tau_ae(i,k,id_CINO3M,inu)*piz_ae(i,k,id_CINO3M,inu)) &
1202	/tau_allaer(i,k,mrfspecies,inu)
1203	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1204
1205	cg_allaer(i,k,mrfspecies,inu)= &
1206	(tau_ae(i,k,id_ASSO4M,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu)+ &
1207	tau_ae(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu)cg_ae(i,k,id_CSSO4M,inu)+ &
1208	tau_ae(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu)cg_ae(i,k,id_ASBCM,inu)+ &
1209	tau_ae(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu)+ &
1210	tau_ae(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu)cg_ae(i,k,id_ASPOMM,inu)+ &
1211	tau_ae(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu)+ &
1212	tau_ae(i,k,id_ASARp1a,inu)piz_ae(i,k,id_ASARp1a,inu)cg_ae(i,k,id_ASARp1a,inu)+ &
1213	tau_ae(i,k,id_ASARp2a,inu)piz_ae(i,k,id_ASARp2a,inu)cg_ae(i,k,id_ASARp2a,inu)+ &
1214	tau_ae(i,k,id_ASAPp1a,inu)piz_ae(i,k,id_ASAPp1a,inu)cg_ae(i,k,id_ASAPp1a,inu)+ &
1215	tau_ae(i,k,id_ASAPp2a,inu)piz_ae(i,k,id_ASAPp2a,inu)cg_ae(i,k,id_ASAPp2a,inu)+ &
1216	tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu)+ &
1217	tau_ae(i,k,id_ASNO3M,inu)piz_ae(i,k,id_ASNO3M,inu)cg_ae(i,k,id_ASNO3M,inu)+ &
1218	tau_ae(i,k,id_CSNO3M,inu)piz_ae(i,k,id_CSNO3M,inu)cg_ae(i,k,id_CSNO3M,inu)+ &
1219	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu)+ &
1220	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu)+ &
1221	tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu)cg_ae(i,k,id_CIDUSTM,inu)+ &
1222	tau_ae(i,k,id_CINO3M,inu)piz_ae(i,k,id_CINO3M,inu)cg_ae(i,k,id_CINO3M,inu))/ &
1223	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1224	#else
1225	tau_allaer(i,k,mrfspecies,inu) = &
1226	tau_ae(i,k,id_ASSO4M,inu)+ tau_ae(i,k,id_CSSO4M,inu)+ &
1227	tau_ae(i,k,id_ASBCM,inu) + tau_ae(i,k,id_AIBCM,inu)+ &
1228	tau_ae(i,k,id_ASPOMM,inu) + tau_ae(i,k,id_AIPOMM,inu)+ &
1229	tau_ae(i,k,id_ASNO3M,inu) + tau_ae(i,k,id_CSNO3M,inu)+ &
1230	tau_ae(i,k,id_ASSSM,inu) + tau_ae(i,k,id_CSSSM,inu)+ &
1231	tau_ae(i,k,id_SSSSM,inu) + tau_ae(i,k,id_CIDUSTM,inu)+ &
1232	tau_ae(i,k,id_CINO3M,inu)
1233	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1234
1235	piz_allaer(i,k,mrfspecies,inu)= &
1236	(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
1237	tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
1238	tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
1239	tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
1240	tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
1241	tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &
1242	tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
1243	tau_ae(i,k,id_ASNO3M,inu)*piz_ae(i,k,id_ASNO3M,inu)+ &
1244	tau_ae(i,k,id_CSNO3M,inu)*piz_ae(i,k,id_CSNO3M,inu)+ &
1245	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
1246	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
1247	tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)+ &
1248	tau_ae(i,k,id_CINO3M,inu)*piz_ae(i,k,id_CINO3M,inu)) &
1249	/tau_allaer(i,k,mrfspecies,inu)
1250	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1251
1252	cg_allaer(i,k,mrfspecies,inu)= &
1253	(tau_ae(i,k,id_ASSO4M,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu)+ &
1254	tau_ae(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu)cg_ae(i,k,id_CSSO4M,inu)+ &
1255	tau_ae(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu)cg_ae(i,k,id_ASBCM,inu)+ &
1256	tau_ae(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu)+ &
1257	tau_ae(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu)cg_ae(i,k,id_ASPOMM,inu)+ &
1258	tau_ae(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu)+ &
1259	tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu)+ &
1260	tau_ae(i,k,id_ASNO3M,inu)piz_ae(i,k,id_ASNO3M,inu)cg_ae(i,k,id_ASNO3M,inu)+ &
1261	tau_ae(i,k,id_CSNO3M,inu)piz_ae(i,k,id_CSNO3M,inu)cg_ae(i,k,id_CSNO3M,inu)+ &
1262	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu)+ &
1263	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu)+ &
1264	tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu)cg_ae(i,k,id_CIDUSTM,inu)+ &
1265	tau_ae(i,k,id_CINO3M,inu)piz_ae(i,k,id_CINO3M,inu)cg_ae(i,k,id_CINO3M,inu))/ &
1266	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1267
1268	#endif
1269	ENDDO
1270	ENDDO
1271
1272	ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT
1273	DO k=1, PLEV
1274	DO i=1, PLON
1275	#ifdef NMHC
1276	tau_allaer(i,k,mrfspecies,inu)= &
1277	tau_ae(i,k,id_ASSO4M,inu)*fractnat_allaer(i,id_ASSO4M)+ &
1278	tau_ae(i,k,id_CSSO4M,inu)*fractnat_allaer(i,id_CSSO4M)+ &
1279	tau_ae(i,k,id_ASBCM,inu)*fractnat_allaer(i,id_ASBCM)+ &
1280	tau_ae(i,k,id_AIBCM,inu)*fractnat_allaer(i,id_AIBCM)+ &
1281	tau_ae(i,k,id_ASPOMM,inu)*fractnat_allaer(i,id_ASPOMM)+ &
1282	tau_ae(i,k,id_AIPOMM,inu)*fractnat_allaer(i,id_AIPOMM)+ &
1283	tau_ae(i,k,id_ASAPp1a,inu)*fractnat_allaer(i,id_ASAPp1a)+ &
1284	tau_ae(i,k,id_ASAPp2a,inu)*fractnat_allaer(i,id_ASAPp2a)+ &
1285	tau_ae(i,k,id_ASARp1a,inu)*fractnat_allaer(i,id_ASARp1a)+ &
1286	tau_ae(i,k,id_ASARp2a,inu)*fractnat_allaer(i,id_ASARp2a)+ &
1287	tau_ae(i,k,id_ASSSM,inu)*fractnat_allaer(i,id_ASSSM)+ &
1288	tau_ae(i,k,id_CSSSM,inu)*fractnat_allaer(i,id_CSSSM)+ &
1289	tau_ae(i,k,id_SSSSM,inu)*fractnat_allaer(i,id_SSSSM)+ &
1290	tau_ae(i,k,id_CIDUSTM,inu)*fractnat_allaer(i,id_CIDUSTM)
1291	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1292
1293	piz_allaer(i,k,mrfspecies,inu)= &
1294	(tau_ae(i,k,id_ASSO4M,inu)fractnat_allaer(i,id_ASSO4M)piz_ae(i,k,id_ASSO4M,inu)+ &
1295	tau_ae(i,k,id_CSSO4M,inu)fractnat_allaer(i,id_CSSO4M)piz_ae(i,k,id_CSSO4M,inu)+ &
1296	tau_ae(i,k,id_ASBCM,inu)fractnat_allaer(i,id_ASBCM)piz_ae(i,k,id_ASBCM,inu)+ &
1297	tau_ae(i,k,id_AIBCM,inu)fractnat_allaer(i,id_AIBCM)piz_ae(i,k,id_AIBCM,inu)+ &
1298	tau_ae(i,k,id_ASPOMM,inu)fractnat_allaer(i,id_ASPOMM)piz_ae(i,k,id_ASPOMM,inu)+ &
1299	tau_ae(i,k,id_AIPOMM,inu)fractnat_allaer(i,id_AIPOMM)piz_ae(i,k,id_AIPOMM,inu)+ &
1300	tau_ae(i,k,id_ASAPp1a,inu)fractnat_allaer(i,id_ASAPp1a)piz_ae(i,k,id_ASAPp1a,inu)+ &
1301	tau_ae(i,k,id_ASAPp2a,inu)fractnat_allaer(i,id_ASAPp2a)piz_ae(i,k,id_ASAPp2a,inu)+ &
1302	tau_ae(i,k,id_ASARp1a,inu)fractnat_allaer(i,id_ASARp1a)piz_ae(i,k,id_ASARp1a,inu)+ &
1303	tau_ae(i,k,id_ASARp2a,inu)fractnat_allaer(i,id_ASARp2a)piz_ae(i,k,id_ASARp2a,inu)+ &
1304	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_ASSSM)piz_ae(i,k,id_ASSSM,inu)+ &
1305	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CSSSM)piz_ae(i,k,id_CSSSM,inu)+ &
1306	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_SSSSM)piz_ae(i,k,id_SSSSM,inu)+ &
1307	tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CIDUSTM,inu)) &
1308	/tau_allaer(i,k,mrfspecies,inu)
1309	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1310
1311	cg_allaer(i,k,mrfspecies,inu)= &
1312	(tau_ae(i,k,id_ASSO4M,inu)fractnat_allaer(i,id_ASSO4M)piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &
1313	tau_ae(i,k,id_CSSO4M,inu)fractnat_allaer(i,id_CSSO4M)piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ &
1314	tau_ae(i,k,id_ASBCM,inu)fractnat_allaer(i,id_ASBCM)piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ &
1315	tau_ae(i,k,id_AIBCM,inu)fractnat_allaer(i,id_AIBCM)piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu)+ &
1316	tau_ae(i,k,id_ASPOMM,inu)fractnat_allaer(i,id_ASPOMM)piz_ae(i,k,id_ASPOMM,inu)*cg_ae(i,k,id_ASPOMM,inu)+ &
1317	tau_ae(i,k,id_AIPOMM,inu)fractnat_allaer(i,id_AIPOMM)piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu)+ &
1318	tau_ae(i,k,id_ASAPp1a,inu)fractnat_allaer(i,id_ASAPp1a)piz_ae(i,k,id_ASAPp1a,inu)*cg_ae(i,k,id_ASAPp1a,inu)+ &
1319	tau_ae(i,k,id_ASAPp2a,inu)fractnat_allaer(i,id_ASAPp2a)piz_ae(i,k,id_ASAPp2a,inu)*cg_ae(i,k,id_ASAPp2a,inu)+ &
1320	tau_ae(i,k,id_ASARp1a,inu)fractnat_allaer(i,id_ASARp1a)piz_ae(i,k,id_ASARp1a,inu)*cg_ae(i,k,id_ASARp1a,inu)+ &
1321	tau_ae(i,k,id_ASARp2a,inu)fractnat_allaer(i,id_ASARp2a)piz_ae(i,k,id_ASARp2a,inu)*cg_ae(i,k,id_ASARp2a,inu)+ &
1322	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_ASSSM)piz_ae(i,k,id_ASSSM,inu)*cg_ae(i,k,id_ASSSM,inu)+ &
1323	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CSSSM)piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu)+ &
1324	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_SSSSM)piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu)+ &
1325	tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ &
1326	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1327
1328	#else
1329	tau_allaer(i,k,mrfspecies,inu)= &
1330	tau_ae(i,k,id_ASSO4M,inu)*fractnat_allaer(i,id_ASSO4M)+ &
1331	tau_ae(i,k,id_CSSO4M,inu)*fractnat_allaer(i,id_CSSO4M)+ &
1332	tau_ae(i,k,id_ASBCM,inu)*fractnat_allaer(i,id_ASBCM)+ &
1333	tau_ae(i,k,id_AIBCM,inu)*fractnat_allaer(i,id_AIBCM)+ &
1334	tau_ae(i,k,id_ASPOMM,inu)*fractnat_allaer(i,id_ASPOMM)+ &
1335	tau_ae(i,k,id_AIPOMM,inu)*fractnat_allaer(i,id_AIPOMM)+ &
1336	tau_ae(i,k,id_ASSSM,inu)*fractnat_allaer(i,id_ASSSM)+ &
1337	tau_ae(i,k,id_CSSSM,inu)*fractnat_allaer(i,id_CSSSM)+ &
1338	tau_ae(i,k,id_SSSSM,inu)*fractnat_allaer(i,id_SSSSM)+ &
1339	tau_ae(i,k,id_CIDUSTM,inu)*fractnat_allaer(i,id_CIDUSTM)
1340	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1341
1342	piz_allaer(i,k,mrfspecies,inu)= &
1343	(tau_ae(i,k,id_ASSO4M,inu)fractnat_allaer(i,id_ASSO4M)piz_ae(i,k,id_ASSO4M,inu)+ &
1344	tau_ae(i,k,id_CSSO4M,inu)fractnat_allaer(i,id_CSSO4M)piz_ae(i,k,id_CSSO4M,inu)+ &
1345	tau_ae(i,k,id_ASBCM,inu)fractnat_allaer(i,id_ASBCM)piz_ae(i,k,id_ASBCM,inu)+ &
1346	tau_ae(i,k,id_AIBCM,inu)fractnat_allaer(i,id_AIBCM)piz_ae(i,k,id_AIBCM,inu)+ &
1347	tau_ae(i,k,id_ASPOMM,inu)fractnat_allaer(i,id_ASPOMM)piz_ae(i,k,id_ASPOMM,inu)+ &
1348	tau_ae(i,k,id_AIPOMM,inu)fractnat_allaer(i,id_AIPOMM)piz_ae(i,k,id_AIPOMM,inu)+ &
1349	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_ASSSM)piz_ae(i,k,id_ASSSM,inu)+ &
1350	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CSSSM)piz_ae(i,k,id_CSSSM,inu)+ &
1351	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_SSSSM)piz_ae(i,k,id_SSSSM,inu)+ &
1352	tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CIDUSTM,inu)) &
1353	/tau_allaer(i,k,mrfspecies,inu)
1354	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1355
1356	cg_allaer(i,k,mrfspecies,inu)= &
1357	(tau_ae(i,k,id_ASSO4M,inu)fractnat_allaer(i,id_ASSO4M)piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &
1358	tau_ae(i,k,id_CSSO4M,inu)fractnat_allaer(i,id_CSSO4M)piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ &
1359	tau_ae(i,k,id_ASBCM,inu)fractnat_allaer(i,id_ASBCM)piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ &
1360	tau_ae(i,k,id_AIBCM,inu)fractnat_allaer(i,id_AIBCM)piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu)+ &
1361	tau_ae(i,k,id_ASPOMM,inu)fractnat_allaer(i,id_ASPOMM)piz_ae(i,k,id_ASPOMM,inu)*cg_ae(i,k,id_ASPOMM,inu)+ &
1362	tau_ae(i,k,id_AIPOMM,inu)fractnat_allaer(i,id_AIPOMM)piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu)+ &
1363	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_ASSSM)piz_ae(i,k,id_ASSSM,inu)*cg_ae(i,k,id_ASSSM,inu)+ &
1364	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CSSSM)piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu)+ &
1365	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_SSSSM)piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu)+ &
1366	tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ &
1367	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1368
1369	#endif
1370	ENDDO
1371	ENDDO
1372
1373	ELSEIF (mrfspecies .EQ. 4) THEN ! = BC
1374	DO k=1, PLEV
1375	DO i=1, PLON
1376	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)
1377	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1378
1379	piz_allaer(i,k,mrfspecies,inu) = &
1380	(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) + &
1381	tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ &
1382	tau_allaer(i,k,mrfspecies,inu)
1383	piz_allaer(i,k,mrfspecies,inu) = MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1384
1385	cg_allaer(i,k,mrfspecies,inu) = &
1386	(tau_ae(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu)cg_ae(i,k,id_ASBCM,inu) + &
1387	tau_ae(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu))/ &
1388	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1389
1390	ENDDO
1391	ENDDO
1392	ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4
1393	DO k=1, PLEV
1394	DO i=1, PLON
1395	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)
1396	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1397
1398	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) &
1399	+tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ &
1400	tau_allaer(i,k,mrfspecies,inu)
1401	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1402
1403	cg_allaer(i,k,mrfspecies,inu)=&
1404	(tau_ae(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu) cg_ae(i,k,id_CSSO4M,inu)&
1405	+tau_ae(i,k,id_ASSO4M,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu))/ &
1406	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1407	ENDDO
1408	ENDDO
1409
1410	ELSEIF (mrfspecies .EQ. 6) THEN ! = POM
1411	DO k=1, PLEV
1412	DO i=1, PLON
1413	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)
1414	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1415
1416	piz_allaer(i,k,mrfspecies,inu) = &
1417	(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) &
1418	+tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ &
1419	tau_allaer(i,k,mrfspecies,inu)
1420	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1421
1422	cg_allaer(i,k,mrfspecies,inu) = &
1423	(tau_ae(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu) cg_ae(i,k,id_ASPOMM,inu)&
1424	+tau_ae(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu))/ &
1425	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1426
1427	ENDDO
1428	ENDDO
1429
1430	#ifdef NMHC
1431	ELSEIF (mrfspecies .EQ. 12) THEN ! = SOAA
1432	DO k=1, PLEV
1433	DO i=1, PLON
1434	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_ASARp1a,inu)+tau_ae(i,k,id_ASARp2a,inu)
1435	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1436
1437	piz_allaer(i,k,mrfspecies,inu) = &
1438	(tau_ae(i,k,id_ASARp1a,inu)*piz_ae(i,k,id_ASARp1a,inu) &
1439	+tau_ae(i,k,id_ASARp2a,inu)*piz_ae(i,k,id_ASARp2a,inu))/ &
1440	tau_allaer(i,k,mrfspecies,inu)
1441	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1442
1443	cg_allaer(i,k,mrfspecies,inu) = &
1444	(tau_ae(i,k,id_ASARp1a,inu)piz_ae(i,k,id_ASARp1a,inu) cg_ae(i,k,id_ASARp1a,inu)&
1445	+tau_ae(i,k,id_ASARp2a,inu)piz_ae(i,k,id_ASARp2a,inu)cg_ae(i,k,id_ASARp2a,inu))/ &
1446	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1447
1448	ENDDO
1449	ENDDO
1450	ELSEIF (mrfspecies .EQ. 13) THEN ! = SOAB
1451	DO k=1, PLEV
1452	DO i=1, PLON
1453	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_ASAPp1a,inu)+tau_ae(i,k,id_ASAPp2a,inu)
1454	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1455
1456	piz_allaer(i,k,mrfspecies,inu) = &
1457	(tau_ae(i,k,id_ASAPp1a,inu)*piz_ae(i,k,id_ASAPp1a,inu) &
1458	+tau_ae(i,k,id_ASAPp2a,inu)*piz_ae(i,k,id_ASAPp2a,inu))/ &
1459	tau_allaer(i,k,mrfspecies,inu)
1460	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1461
1462	cg_allaer(i,k,mrfspecies,inu) = &
1463	(tau_ae(i,k,id_ASAPp1a,inu)piz_ae(i,k,id_ASAPp1a,inu)cg_ae(i,k,id_ASAPp1a,inu)&
1464	+tau_ae(i,k,id_ASAPp2a,inu)piz_ae(i,k,id_ASAPp2a,inu)cg_ae(i,k,id_ASAPp2a,inu))/ &
1465	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1466
1467	ENDDO
1468	ENDDO
1469	#endif
1470	ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST
1471	DO k=1, PLEV
1472	DO i=1, PLON
1473	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_CIDUSTM,inu)
1474	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1475
1476	piz_allaer(i,k,mrfspecies,inu)= &
1477	tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu) / &
1478	tau_allaer(i,k,mrfspecies,inu)
1479	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1480
1481	cg_allaer(i,k,mrfspecies,inu)= &
1482	tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu) cg_ae(i,k,id_CIDUSTM,inu)/ &
1483	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1484
1485
1486	ENDDO
1487	ENDDO
1488	ELSEIF (mrfspecies .EQ. 8) THEN ! = SS
1489	DO k=1, PLEV
1490	DO i=1, PLON
1491	tau_allaer(i,k,mrfspecies,inu) = &
1492	tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)
1493	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1494
1495	piz_allaer(i,k,mrfspecies,inu)= &
1496	(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) + &
1497	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) + &
1498	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
1499	tau_allaer(i,k,mrfspecies,inu)
1500	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1501
1502	cg_allaer(i,k,mrfspecies,inu)= &
1503	(tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu) cg_ae(i,k,id_ASSSM,inu) + &
1504	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu) + &
1505	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu))/ &
1506	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1507
1508	ENDDO
1509	ENDDO
1510
1511	ELSEIF (mrfspecies .EQ. 9) THEN ! = FNO3
1512	DO k=1, PLEV
1513	DO i=1, PLON
1514
1515	tau_allaer(i,k,mrfspecies,inu)= tau_ae(i,k,id_ASNO3M,inu)
1516	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1517
1518	piz_allaer(i,k,mrfspecies,inu)= piz_ae(i,k,id_ASNO3M,inu)
1519	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1520
1521	cg_allaer(i,k,mrfspecies,inu)= cg_ae(i,k,id_ASNO3M,inu)
1522
1523	ENDDO
1524	ENDDO
1525	ELSEIF (mrfspecies .EQ. 10) THEN ! = DNO3
1526	DO k=1, PLEV
1527	DO i=1, PLON
1528	tau_allaer(i,k,mrfspecies,inu)= tau_ae(i,k,id_CINO3M,inu)
1529	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1530
1531	piz_allaer(i,k,mrfspecies,inu)= piz_ae(i,k,id_CINO3M,inu)
1532	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1533
1534	cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CINO3M,inu)
1535
1536	ENDDO
1537	ENDDO
1538	ELSEIF (mrfspecies .EQ. 11) THEN ! = SNO3
1539	DO k=1, PLEV
1540	DO i=1, PLON
1541
1542	tau_allaer(i,k,mrfspecies,inu)= tau_ae(i,k,id_CSNO3M,inu)
1543	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1544
1545	piz_allaer(i,k,mrfspecies,inu)= piz_ae(i,k,id_CSNO3M,inu)
1546	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1547
1548	cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CSNO3M,inu)
1549
1550	ENDDO
1551	ENDDO
1552	ENDIF
1553	ENDDO
1554	ENDDO
1555	#else
1556	DO inu=1, 2
1557	DO mrfspecies=1,naero_grp
1558	IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER
1559	DO k=1, PLEV
1560	DO i=1, PLON
1561
1562	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_CIDUSTM,inu)+ tau_ae(i,k,id_ASSSM,inu)+ &
1563	tau_ae(i,k,id_CSSSM,inu)+ tau_ae(i,k,id_SSSSM,inu)
1564	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1565
1566	piz_allaer(i,k,mrfspecies,inu)= &
1567	(tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)+ &
1568	tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
1569	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
1570	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)) &
1571	/tau_allaer(i,k,mrfspecies,inu)
1572	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1573
1574	cg_allaer(i,k,mrfspecies,inu)= &
1575	(tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu)cg_ae(i,k,id_CIDUSTM,inu)+ &
1576	tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu)+ &
1577	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu)+ &
1578	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu)) &
1579	/(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1580
1581	ENDDO
1582	ENDDO
1583
1584	ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT (DUST + SS)
1585	DO k=1, PLEV
1586	DO i=1, PLON
1587	tau_allaer(i,k,mrfspecies,inu)= &
1588	(tau_ae(i,k,id_CIDUSTM,inu)*fractnat_allaer(i,id_CIDUSTM)+ &
1589	tau_ae(i,k,id_ASSSM,inu)*fractnat_allaer(i,id_ASSSM)+ &
1590	tau_ae(i,k,id_CSSSM,inu)*fractnat_allaer(i,id_CSSSM)+ &
1591	tau_ae(i,k,id_SSSSM,inu)*fractnat_allaer(i,id_SSSSM))
1592	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1593
1594	piz_allaer(i,k,mrfspecies,inu)= &
1595	(tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CIDUSTM,inu)+ &
1596	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_ASSSM,inu)+ &
1597	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_CSSSM,inu)+ &
1598	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_CIDUSTM)piz_ae(i,k,id_SSSSM,inu)) &
1599	/tau_allaer(i,k,mrfspecies,inu)
1600	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1601
1602	cg_allaer(i,k,mrfspecies,inu)= &
1603	(tau_ae(i,k,id_CIDUSTM,inu)fractnat_allaer(i,id_CIDUSTM) &
1604	piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu)+ &
1605	tau_ae(i,k,id_ASSSM,inu)fractnat_allaer(i,id_ASSSM) &
1606	piz_ae(i,k,id_ASSSM,inu)*cg_ae(i,k,id_ASSSM,inu)+ &
1607	tau_ae(i,k,id_CSSSM,inu)fractnat_allaer(i,id_CSSSM) &
1608	piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu)+ &
1609	tau_ae(i,k,id_SSSSM,inu)fractnat_allaer(i,id_SSSSM) &
1610	piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ &
1611	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1612
1613	ENDDO
1614	ENDDO
1615
1616	ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST ONLY
1617	DO k=1, PLEV
1618	DO i=1, PLON
1619	tau_allaer(i,k,mrfspecies,inu) = tau_ae(i,k,id_CIDUSTM,inu)
1620	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1621
1622	piz_allaer(i,k,mrfspecies,inu)= &
1623	(tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu))/ tau_allaer(i,k,mrfspecies,inu)
1624	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1625
1626	cg_allaer(i,k,mrfspecies,inu)= &
1627	(tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu) cg_ae(i,k,id_CIDUSTM,inu))/ &
1628	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1629
1630	ENDDO
1631	ENDDO
1632	ELSEIF (mrfspecies .EQ. 8) THEN ! = SS
1633	DO k=1, PLEV
1634	DO i=1, PLON
1635
1636	tau_allaer(i,k,mrfspecies,inu) = &
1637	tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)
1638	tau_allaer(i,k,mrfspecies,inu) = MAX(tau_allaer(i,k,mrfspecies,inu),mintau)
1639
1640	piz_allaer(i,k,mrfspecies,inu)= &
1641	(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) + &
1642	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) + &
1643	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
1644	tau_allaer(i,k,mrfspecies,inu)
1645	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),minpiz)
1646
1647	cg_allaer(i,k,mrfspecies,inu)= &
1648	(tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu) + &
1649	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu) + &
1650	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu))/ &
1651	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1652
1653	ENDDO
1654	ENDDO
1655	ENDIF
1656	ENDDO
1657	ENDDO
1658	#endif
1659	!
1660	! DO i=1, PLON
1661	! ai(i)=(-log(MAX(taue670(i),0.0001)/
1662	! . MAX(taue865(i),0.0001))/log(670./865.)) *
1663	! . taue865(i)
1664	! ENDDO
1665	!
1666	RETURN
1667	END SUBROUTINE OPTAER_2BANDS
1668	#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_2bands.F90 @ 6610

Download in other formats: