Context Navigation

aeropt_2bands.F90 @ 3814

Last change on this file since 3814 was 3814, checked in by ymipsl, 10 years ago

remove all dynamic dependency in LMDZ physics except for the include "dimensions.h"

File size: 54.2 KB

Line
1	!
2	! $Id: aeropt_2bands.F90 2148 2014-11-18 17:55:52Z fhourdin $
3	!
4	SUBROUTINE AEROPT_2BANDS( &
5	pdel, m_allaer, delt, RHcl, &
6	tau_allaer, piz_allaer, &
7	cg_allaer, m_allaer_pi, &
8	flag_aerosol, pplay, t_seri, presnivs)
9
10	USE dimphy
11	USE aero_mod
12	USE phys_local_var_mod, only: absvisaer
13	USE pres2lev_mod
14
15
16	! Yves Balkanski le 12 avril 2006
17	! Celine Deandreis
18	! Anne Cozic Avril 2009
19	! a partir d'une sous-routine de Johannes Quaas pour les sulfates
20	!
21	IMPLICIT NONE
22
23	INCLUDE "YOMCST.h"
24	INCLUDE "iniprint.h"
25
26	!
27	! Input arguments:
28	!
29	REAL, DIMENSION(klon,klev), INTENT(in) :: pdel
30	REAL, INTENT(in) :: delt
31	REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer
32	!RAF
33	REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer_pi
34	REAL, DIMENSION(klon,klev), INTENT(in) :: RHcl ! humidite relative ciel clair
35	!RAF REAL, DIMENSION(klon,naero_tot),INTENT(in) :: fractnat_allaer
36	INTEGER, INTENT(in) :: flag_aerosol
37	REAL, DIMENSION(klon,klev), INTENT(in) :: pplay
38	REAL, DIMENSION(klon,klev), INTENT(in) :: t_seri
39	REAL, DIMENSION(klev), INTENT(in) :: presnivs
40	!
41	! Output arguments:
42	!
43	REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(out) :: tau_allaer ! epaisseur optique aerosol
44	REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(out) :: piz_allaer ! single scattering albedo aerosol
45	REAL, DIMENSION(klon,klev,naero_grp,nbands), INTENT(out) :: cg_allaer ! asymmetry parameter aerosol
46
47	!
48	! Local
49	!
50	REAL, DIMENSION(klon,klev,naero_tot,nbands) :: tau_ae
51	!RAF
52	REAL, DIMENSION(klon,klev,naero_tot,nbands) :: tau_ae_pi
53	REAL, DIMENSION(klon,klev,naero_tot,nbands) :: piz_ae
54	REAL, DIMENSION(klon,klev,naero_tot,nbands) :: cg_ae
55	LOGICAL :: soluble
56	INTEGER :: i, k,n, ierr, inu, m, mrfspecies
57	INTEGER :: spsol, spinsol, spss
58	INTEGER :: RH_num(klon,klev)
59	INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
60
61	INTEGER, PARAMETER :: nbre_RH=12
62	INTEGER, PARAMETER :: naero_soluble=7 ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
63	! 5- seasalt super coarse 6- seasalt coarse 7- seasalt acc.
64	INTEGER, PARAMETER :: naero_insoluble=3 ! 1- Dust; 2- BC insoluble; 3- POM insoluble
65	LOGICAL, SAVE :: firstcall=.TRUE.
66	!$OMP THREADPRIVATE(firstcall)
67
68	! Coefficient optiques sur 19 niveaux
69	REAL, SAVE, DIMENSION(nb_level) :: presnivs_19 ! Pression milieux couche pour 19 niveaux (nb_level)
70	!$OMP THREADPRIVATE(presnivs_19)
71
72	REAL, SAVE, DIMENSION(nb_level) :: A1_ASSSM_b1_19, A2_ASSSM_b1_19, A3_ASSSM_b1_19,&
73	B1_ASSSM_b1_19, B2_ASSSM_b1_19, C1_ASSSM_b1_19, C2_ASSSM_b1_19,&
74	A1_CSSSM_b1_19, A2_CSSSM_b1_19, A3_CSSSM_b1_19,&
75	B1_CSSSM_b1_19, B2_CSSSM_b1_19, C1_CSSSM_b1_19, C2_CSSSM_b1_19,&
76	A1_SSSSM_b1_19, A2_SSSSM_b1_19, A3_SSSSM_b1_19,&
77	B1_SSSSM_b1_19, B2_SSSSM_b1_19, C1_SSSSM_b1_19, C2_SSSSM_b1_19,&
78	A1_ASSSM_b2_19, A2_ASSSM_b2_19, A3_ASSSM_b2_19,&
79	B1_ASSSM_b2_19, B2_ASSSM_b2_19, C1_ASSSM_b2_19, C2_ASSSM_b2_19,&
80	A1_CSSSM_b2_19, A2_CSSSM_b2_19, A3_CSSSM_b2_19,&
81	B1_CSSSM_b2_19, B2_CSSSM_b2_19, C1_CSSSM_b2_19, C2_CSSSM_b2_19,&
82	A1_SSSSM_b2_19, A2_SSSSM_b2_19, A3_SSSSM_b2_19,&
83	B1_SSSSM_b2_19, B2_SSSSM_b2_19, C1_SSSSM_b2_19, C2_SSSSM_b2_19
84	!$OMP THREADPRIVATE(A1_ASSSM_b1_19, A2_ASSSM_b1_19, A3_ASSSM_b1_19)
85	!$OMP THREADPRIVATE(B1_ASSSM_b1_19, B2_ASSSM_b1_19, C1_ASSSM_b1_19, C2_ASSSM_b1_19)
86	!$OMP THREADPRIVATE(A1_CSSSM_b1_19, A2_CSSSM_b1_19, A3_CSSSM_b1_19)
87	!$OMP THREADPRIVATE(B1_CSSSM_b1_19, B2_CSSSM_b1_19, C1_CSSSM_b1_19, C2_CSSSM_b1_19)
88	!$OMP THREADPRIVATE(A1_SSSSM_b1_19, A2_SSSSM_b1_19, A3_SSSSM_b1_19)
89	!$OMP THREADPRIVATE(B1_SSSSM_b1_19, B2_SSSSM_b1_19, C1_SSSSM_b1_19, C2_SSSSM_b1_19)
90	!$OMP THREADPRIVATE(A1_ASSSM_b2_19, A2_ASSSM_b2_19, A3_ASSSM_b2_19)
91	!$OMP THREADPRIVATE(B1_ASSSM_b2_19, B2_ASSSM_b2_19, C1_ASSSM_b2_19, C2_ASSSM_b2_19)
92	!$OMP THREADPRIVATE(A1_CSSSM_b2_19, A2_CSSSM_b2_19, A3_CSSSM_b2_19)
93	!$OMP THREADPRIVATE(B1_CSSSM_b2_19, B2_CSSSM_b2_19, C1_CSSSM_b2_19, C2_CSSSM_b2_19)
94	!$OMP THREADPRIVATE(A1_SSSSM_b2_19, A2_SSSSM_b2_19, A3_SSSSM_b2_19)
95	!$OMP THREADPRIVATE(B1_SSSSM_b2_19, B2_SSSSM_b2_19, C1_SSSSM_b2_19, C2_SSSSM_b2_19)
96
97
98	! Coefficient optiques interpole sur le nombre de niveau du modele
99	REAL, ALLOCATABLE, DIMENSION(:), SAVE :: &
100	A1_ASSSM_b1, A2_ASSSM_b1, A3_ASSSM_b1,&
101	B1_ASSSM_b1, B2_ASSSM_b1, C1_ASSSM_b1, C2_ASSSM_b1,&
102	A1_CSSSM_b1, A2_CSSSM_b1, A3_CSSSM_b1,&
103	B1_CSSSM_b1, B2_CSSSM_b1, C1_CSSSM_b1, C2_CSSSM_b1,&
104	A1_SSSSM_b1, A2_SSSSM_b1, A3_SSSSM_b1,&
105	B1_SSSSM_b1, B2_SSSSM_b1, C1_SSSSM_b1, C2_SSSSM_b1,&
106	A1_ASSSM_b2, A2_ASSSM_b2, A3_ASSSM_b2,&
107	B1_ASSSM_b2, B2_ASSSM_b2, C1_ASSSM_b2, C2_ASSSM_b2,&
108	A1_CSSSM_b2, A2_CSSSM_b2, A3_CSSSM_b2,&
109	B1_CSSSM_b2, B2_CSSSM_b2, C1_CSSSM_b2, C2_CSSSM_b2,&
110	A1_SSSSM_b2, A2_SSSSM_b2, A3_SSSSM_b2,&
111	B1_SSSSM_b2, B2_SSSSM_b2, C1_SSSSM_b2, C2_SSSSM_b2
112	!$OMP THREADPRIVATE(A1_ASSSM_b1, A2_ASSSM_b1, A3_ASSSM_b1)
113	!$OMP THREADPRIVATE(B1_ASSSM_b1, B2_ASSSM_b1, C1_ASSSM_b1, C2_ASSSM_b1)
114	!$OMP THREADPRIVATE(A1_CSSSM_b1, A2_CSSSM_b1, A3_CSSSM_b1)
115	!$OMP THREADPRIVATE(B1_CSSSM_b1, B2_CSSSM_b1, C1_CSSSM_b1, C2_CSSSM_b1)
116	!$OMP THREADPRIVATE(A1_SSSSM_b1, A2_SSSSM_b1, A3_SSSSM_b1)
117	!$OMP THREADPRIVATE(B1_SSSSM_b1, B2_SSSSM_b1, C1_SSSSM_b1, C2_SSSSM_b1)
118	!$OMP THREADPRIVATE(A1_ASSSM_b2, A2_ASSSM_b2, A3_ASSSM_b2)
119	!$OMP THREADPRIVATE(B1_ASSSM_b2, B2_ASSSM_b2, C1_ASSSM_b2, C2_ASSSM_b2)
120	!$OMP THREADPRIVATE(A1_CSSSM_b2, A2_CSSSM_b2, A3_CSSSM_b2)
121	!$OMP THREADPRIVATE(B1_CSSSM_b2, B2_CSSSM_b2, C1_CSSSM_b2, C2_CSSSM_b2)
122	!$OMP THREADPRIVATE(A1_SSSSM_b2, A2_SSSSM_b2, A3_SSSSM_b2)
123	!$OMP THREADPRIVATE(B1_SSSSM_b2, B2_SSSSM_b2, C1_SSSSM_b2, C2_SSSSM_b2)
124
125	REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
126	REAL, PARAMETER :: RH_MAX=95.
127	REAL:: DELTA(klon,klev), rh(klon,klev), H
128	REAL:: tau_ae2b_int ! Intermediate computation of epaisseur optique aerosol
129	REAL:: piz_ae2b_int ! Intermediate computation of Single scattering albedo
130	REAL:: cg_ae2b_int ! Intermediate computation of Assymetry parameter
131	REAL :: Fact_RH(nbre_RH)
132	REAL :: zrho
133	REAL :: fac
134	REAL :: zdp1(klon,klev)
135	REAL, PARAMETER :: gravit = 9.80616 ! m2/s
136	INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name
137	INTEGER :: nb_aer
138	REAL, DIMENSION(klon,klev,naero_tot) :: mass_temp
139	!RAF
140	REAL, DIMENSION(klon,klev,naero_tot) :: mass_temp_pi
141
142	!
143	! Proprietes optiques
144	!
145	REAL:: alpha_aers_2bands(nbre_RH,nbands,naero_soluble) !--unit m2/g SO4
146	REAL:: alpha_aeri_2bands(nbands,naero_insoluble)
147	REAL:: cg_aers_2bands(nbre_RH,nbands,naero_soluble) !--unit
148	REAL:: cg_aeri_2bands(nbands,naero_insoluble)
149	REAL:: piz_aers_2bands(nbre_RH,nbands,naero_soluble) !-- unit
150	REAL:: piz_aeri_2bands(nbands,naero_insoluble) !-- unit
151
152	INTEGER :: id
153	LOGICAL :: used_aer(naero_tot)
154	REAL :: tmp_var, tmp_var_pi
155
156	DATA presnivs_19/&
157	100426.5, 98327.6, 95346.5, 90966.8, 84776.9, &
158	76536.5, 66292.2, 54559.3, 42501.8, 31806, &
159	23787.5, 18252.7, 13996, 10320.8, 7191.1, &
160	4661.7, 2732.9, 1345.6, 388.2/
161
162
163	!*********************BAND 1*********************************
164	!ACCUMULATION MODE
165	DATA A1_ASSSM_b1_19/ 4.373E+00, 4.361E+00, 4.331E+00, &
166	4.278E+00, 4.223E+00, 4.162E+00, &
167	4.103E+00, 4.035E+00, 3.962E+00, &
168	3.904E+00, 3.871E+00, 3.847E+00, &
169	3.824E+00, 3.780E+00, 3.646E+00, &
170	3.448E+00, 3.179E+00, 2.855E+00, 2.630E+00/
171	DATA A2_ASSSM_b1_19/ 2.496E+00, 2.489E+00, 2.472E+00, &
172	2.442E+00, 2.411E+00, 2.376E+00, &
173	2.342E+00, 2.303E+00, 2.261E+00, &
174	2.228E+00, 2.210E+00, 2.196E+00, &
175	2.183E+00, 2.158E+00, 2.081E+00, &
176	1.968E+00, 1.814E+00, 1.630E+00, 1.501E+00/
177	DATA A3_ASSSM_b1_19/-4.688E-02, -4.676E-02, -4.644E-02, &
178	-4.587E-02, -4.528E-02, -4.463E-02, &
179	-4.399E-02, -4.326E-02, -4.248E-02, &
180	-4.186E-02, -4.151E-02, -4.125E-02, &
181	-4.100E-02, -4.053E-02, -3.910E-02, &
182	-3.697E-02, -3.408E-02, -3.061E-02, -2.819E-02/
183	DATA B1_ASSSM_b1_19/ 1.165E-08, 1.145E-08, 1.097E-08, &
184	1.012E-08, 9.233E-09, 8.261E-09, &
185	7.297E-09, 6.201E-09, 5.026E-09, &
186	4.098E-09, 3.567E-09, 3.187E-09, &
187	2.807E-09, 2.291E-09, 2.075E-09, &
188	1.756E-09, 1.322E-09, 8.011E-10, 4.379E-10/
189	DATA B2_ASSSM_b1_19/ 2.193E-08, 2.192E-08, 2.187E-08, &
190	2.179E-08, 2.171E-08, 2.162E-08, &
191	2.153E-08, 2.143E-08, 2.132E-08, &
192	2.124E-08, 2.119E-08, 2.115E-08, &
193	2.112E-08, 2.106E-08, 2.100E-08, &
194	2.090E-08, 2.077E-08, 2.061E-08, 2.049E-08/
195	DATA C1_ASSSM_b1_19/ 7.365E-01, 7.365E-01, 7.365E-01, &
196	7.364E-01, 7.363E-01, 7.362E-01, &
197	7.361E-01, 7.359E-01, 7.358E-01, &
198	7.357E-01, 7.356E-01, 7.356E-01, &
199	7.356E-01, 7.355E-01, 7.354E-01, &
200	7.352E-01, 7.350E-01, 7.347E-01, 7.345E-01/
201	DATA C2_ASSSM_b1_19/ 5.833E-02, 5.835E-02, 5.841E-02, &
202	5.850E-02, 5.859E-02, 5.870E-02, &
203	5.880E-02, 5.891E-02, 5.904E-02, &
204	5.914E-02, 5.920E-02, 5.924E-02, &
205	5.928E-02, 5.934E-02, 5.944E-02, &
206	5.959E-02, 5.979E-02, 6.003E-02, 6.020E-02/
207	!COARSE MODE
208	DATA A1_CSSSM_b1_19/ 7.403E-01, 7.422E-01, 7.626E-01, &
209	8.019E-01, 8.270E-01, 8.527E-01, &
210	8.702E-01, 8.806E-01, 8.937E-01, &
211	9.489E-01, 1.030E+00, 1.105E+00, &
212	1.199E+00, 1.357E+00, 1.660E+00, &
213	2.540E+00, 4.421E+00, 2.151E+00, 9.518E-01/
214	DATA A2_CSSSM_b1_19/ 4.522E-01, 4.532E-01, 4.644E-01, &
215	4.859E-01, 4.996E-01, 5.137E-01, &
216	5.233E-01, 5.290E-01, 5.361E-01, &
217	5.655E-01, 6.085E-01, 6.483E-01, &
218	6.979E-01, 7.819E-01, 9.488E-01, &
219	1.450E+00, 2.523E+00, 1.228E+00, 5.433E-01/
220	DATA A3_CSSSM_b1_19/-8.516E-03, -8.535E-03, -8.744E-03, &
221	-9.148E-03, -9.406E-03, -9.668E-03, &
222	-9.848E-03, -9.955E-03, -1.009E-02, &
223	-1.064E-02, -1.145E-02, -1.219E-02, &
224	-1.312E-02, -1.470E-02, -1.783E-02, &
225	-2.724E-02, -4.740E-02, -2.306E-02, -1.021E-02/
226	DATA B1_CSSSM_b1_19/ 2.535E-07, 2.530E-07, 2.479E-07, &
227	2.380E-07, 2.317E-07, 2.252E-07, &
228	2.208E-07, 2.182E-07, 2.149E-07, &
229	2.051E-07, 1.912E-07, 1.784E-07, &
230	1.624E-07, 1.353E-07, 1.012E-07, &
231	6.016E-08, 2.102E-08, 0.000E+00, 0.000E+00/
232	DATA B2_CSSSM_b1_19/ 1.221E-07, 1.217E-07, 1.179E-07, &
233	1.104E-07, 1.056E-07, 1.008E-07, &
234	9.744E-08, 9.546E-08, 9.299E-08, &
235	8.807E-08, 8.150E-08, 7.544E-08, &
236	6.786E-08, 5.504E-08, 4.080E-08, &
237	2.960E-08, 2.300E-08, 2.030E-08, 1.997E-08/
238	DATA C1_CSSSM_b1_19/ 7.659E-01, 7.658E-01, 7.652E-01, &
239	7.639E-01, 7.631E-01, 7.623E-01, &
240	7.618E-01, 7.614E-01, 7.610E-01, &
241	7.598E-01, 7.581E-01, 7.566E-01, &
242	7.546E-01, 7.513E-01, 7.472E-01, &
243	7.423E-01, 7.376E-01, 7.342E-01, 7.334E-01/
244	DATA C2_CSSSM_b1_19/ 3.691E-02, 3.694E-02, 3.729E-02, &
245	3.796E-02, 3.839E-02, 3.883E-02, &
246	3.913E-02, 3.931E-02, 3.953E-02, &
247	4.035E-02, 4.153E-02, 4.263E-02, &
248	4.400E-02, 4.631E-02, 4.933E-02, &
249	5.331E-02, 5.734E-02, 6.053E-02, 6.128E-02/
250	!SUPER COARSE MODE
251	DATA A1_SSSSM_b1_19/ 2.836E-01, 2.876E-01, 2.563E-01, &
252	2.414E-01, 2.541E-01, 2.546E-01, &
253	2.572E-01, 2.638E-01, 2.781E-01, &
254	3.167E-01, 4.209E-01, 5.286E-01, &
255	6.959E-01, 9.233E-01, 1.282E+00, &
256	1.836E+00, 2.981E+00, 4.355E+00, 4.059E+00/
257	DATA A2_SSSSM_b1_19/ 1.608E-01, 1.651E-01, 1.577E-01, &
258	1.587E-01, 1.686E-01, 1.690E-01, &
259	1.711E-01, 1.762E-01, 1.874E-01, &
260	2.138E-01, 2.751E-01, 3.363E-01, &
261	4.279E-01, 5.519E-01, 7.421E-01, &
262	1.048E+00, 1.702E+00, 2.485E+00, 2.317E+00/
263	DATA A3_SSSSM_b1_19/-3.025E-03, -3.111E-03, -2.981E-03, &
264	-3.005E-03, -3.193E-03, -3.200E-03, &
265	-3.239E-03, -3.336E-03, -3.548E-03, &
266	-4.047E-03, -5.196E-03, -6.345E-03, &
267	-8.061E-03, -1.038E-02, -1.395E-02, &
268	-1.970E-02, -3.197E-02, -4.669E-02, -4.352E-02/
269	DATA B1_SSSSM_b1_19/ 6.759E-07, 6.246E-07, 5.542E-07, &
270	4.953E-07, 4.746E-07, 4.738E-07, &
271	4.695E-07, 4.588E-07, 4.354E-07, &
272	3.947E-07, 3.461E-07, 3.067E-07, &
273	2.646E-07, 2.095E-07, 1.481E-07, &
274	9.024E-08, 5.747E-08, 2.384E-08, 6.599E-09/
275	DATA B2_SSSSM_b1_19/ 5.977E-07, 5.390E-07, 4.468E-07, &
276	3.696E-07, 3.443E-07, 3.433E-07, &
277	3.380E-07, 3.249E-07, 2.962E-07, &
278	2.483E-07, 1.989E-07, 1.623E-07, &
279	1.305E-07, 9.015E-08, 6.111E-08, &
280	3.761E-08, 2.903E-08, 2.337E-08, 2.147E-08/
281	DATA C1_SSSSM_b1_19/ 8.120E-01, 8.084E-01, 8.016E-01, &
282	7.953E-01, 7.929E-01, 7.928E-01, &
283	7.923E-01, 7.910E-01, 7.882E-01, &
284	7.834E-01, 7.774E-01, 7.725E-01, &
285	7.673E-01, 7.604E-01, 7.529E-01, &
286	7.458E-01, 7.419E-01, 7.379E-01, 7.360E-01/
287	DATA C2_SSSSM_b1_19/ 2.388E-02, 2.392E-02, 2.457E-02, 2.552E-02, &
288	2.615E-02, 2.618E-02, 2.631E-02, 2.663E-02, &
289	2.735E-02, 2.875E-02, 3.113E-02, 3.330E-02, &
290	3.615E-02, 3.997E-02, 4.521E-02, 5.038E-02, &
291	5.358E-02, 5.705E-02, 5.887E-02/
292	!*******************BAND 2**********************************************
293	!ACCUMULATION MODE
294	DATA A1_ASSSM_b2_19/1.256E+00, 1.246E+00, 1.226E+00, 1.187E+00, 1.148E+00, &
295	1.105E+00, 1.062E+00, 1.014E+00, 9.616E-01, 9.205E-01, &
296	8.970E-01, 8.800E-01, 8.632E-01, 8.371E-01, 7.943E-01, &
297	7.308E-01, 6.448E-01, 5.414E-01, 4.693E-01/
298	DATA A2_ASSSM_b2_19/5.321E-01, 5.284E-01, 5.196E-01, 5.036E-01, 4.872E-01, &
299	4.691E-01, 4.512E-01, 4.308E-01, 4.089E-01, 3.917E-01, &
300	3.818E-01, 3.747E-01, 3.676E-01, 3.567E-01, 3.385E-01, &
301	3.116E-01, 2.751E-01, 2.312E-01, 2.006E-01/
302	DATA A3_ASSSM_b2_19/-1.053E-02, -1.046E-02, -1.028E-02, -9.964E-03, -9.637E-03, &
303	-9.279E-03, -8.923E-03, -8.518E-03, -8.084E-03, -7.741E-03, &
304	-7.545E-03, -7.405E-03, -7.265E-03, -7.048E-03, -6.687E-03, &
305	-6.156E-03, -5.433E-03, -4.565E-03, -3.961E-03/
306	DATA B1_ASSSM_b2_19/1.560E-02, 1.560E-02, 1.561E-02, 1.565E-02, 1.568E-02, &
307	1.572E-02, 1.576E-02, 1.580E-02, 1.584E-02, 1.588E-02, &
308	1.590E-02, 1.592E-02, 1.593E-02, 1.595E-02, 1.599E-02, &
309	1.605E-02, 1.612E-02, 1.621E-02, 1.627E-02/
310	DATA B2_ASSSM_b2_19/1.073E-02, 1.074E-02, 1.076E-02, 1.079E-02, 1.082E-02, &
311	1.085E-02, 1.089E-02, 1.093E-02, 1.097E-02, 1.100E-02, &
312	1.102E-02, 1.103E-02, 1.105E-02, 1.107E-02, 1.110E-02, &
313	1.115E-02, 1.122E-02, 1.130E-02, 1.136E-02/
314	DATA C1_ASSSM_b2_19/7.429E-01, 7.429E-01, 7.429E-01, 7.427E-01, 7.427E-01, &
315	7.424E-01, 7.423E-01, 7.422E-01, 7.421E-01, 7.420E-01, &
316	7.419E-01, 7.419E-01, 7.418E-01, 7.417E-01, 7.416E-01, &
317	7.415E-01, 7.413E-01, 7.409E-01, 7.408E-01/
318	DATA C2_ASSSM_b2_19/3.031E-02, 3.028E-02, 3.022E-02, 3.011E-02, 2.999E-02, &
319	2.986E-02, 2.973E-02, 2.959E-02, 2.943E-02, 2.931E-02, &
320	2.924E-02, 2.919E-02, 2.913E-02, 2.905E-02, 2.893E-02, &
321	2.874E-02, 2.847E-02, 2.817E-02, 2.795E-02/
322	!COARSE MODE
323	DATA A1_CSSSM_b2_19/7.061E-01, 7.074E-01, 7.211E-01, 7.476E-01, 7.647E-01, &
324	7.817E-01, 7.937E-01, 8.007E-01, 8.095E-01, 8.436E-01, &
325	8.932E-01, 9.390E-01, 9.963E-01, 1.093E+00, 1.256E+00, &
326	1.668E+00, 1.581E+00, 3.457E-01, 1.331E-01/
327	DATA A2_CSSSM_b2_19/3.617E-01, 3.621E-01, 3.662E-01, 3.739E-01, 3.789E-01, &
328	3.840E-01, 3.874E-01, 3.895E-01, 3.921E-01, 4.001E-01, &
329	4.117E-01, 4.223E-01, 4.356E-01, 4.581E-01, 5.099E-01, &
330	6.831E-01, 6.663E-01, 1.481E-01, 5.703E-02/
331	DATA A3_CSSSM_b2_19/-6.953E-03, -6.961E-03, -7.048E-03, -7.216E-03, -7.322E-03, &
332	-7.431E-03, -7.506E-03, -7.551E-03, -7.606E-03, -7.791E-03, &
333	-8.059E-03, -8.305E-03, -8.613E-03, -9.134E-03, -1.023E-02, &
334	-1.365E-02, -1.320E-02, -2.922E-03, -1.125E-03/
335	DATA B1_CSSSM_b2_19/1.007E-02, 1.008E-02, 1.012E-02, 1.019E-02, 1.024E-02, &
336	1.029E-02, 1.033E-02, 1.035E-02, 1.038E-02, 1.056E-02, &
337	1.083E-02, 1.109E-02, 1.140E-02, 1.194E-02, 1.270E-02, &
338	1.390E-02, 1.524E-02, 1.639E-02, 1.667E-02/
339	DATA B2_CSSSM_b2_19/4.675E-03, 4.682E-03, 4.760E-03, 4.908E-03, 5.004E-03, &
340	5.102E-03, 5.168E-03, 5.207E-03, 5.256E-03, 5.474E-03, &
341	5.793E-03, 6.089E-03, 6.457E-03, 7.081E-03, 7.923E-03, &
342	9.127E-03, 1.041E-02, 1.147E-02, 1.173E-02/
343	DATA C1_CSSSM_b2_19/7.571E-01, 7.571E-01, 7.570E-01, 7.568E-01, 7.565E-01, &
344	7.564E-01, 7.563E-01, 7.562E-01, 7.562E-01, 7.557E-01, &
345	7.552E-01, 7.545E-01, 7.539E-01, 7.527E-01, 7.509E-01, &
346	7.478E-01, 7.440E-01, 7.404E-01, 7.394E-01/
347	DATA C2_CSSSM_b2_19/4.464E-02, 4.465E-02, 4.468E-02, 4.474E-02, 4.477E-02, &
348	4.480E-02, 4.482E-02, 4.484E-02, 4.486E-02, 4.448E-02, &
349	4.389E-02, 4.334E-02, 4.264E-02, 4.148E-02, 3.957E-02, &
350	3.588E-02, 3.149E-02, 2.751E-02, 2.650E-02/
351	!SUPER COARSE MODE
352	DATA A1_SSSSM_b2_19/2.357E-01, 2.490E-01, 2.666E-01, 2.920E-01, 3.120E-01, &
353	3.128E-01, 3.169E-01, 3.272E-01, 3.498E-01, 3.960E-01, &
354	4.822E-01, 5.634E-01, 6.763E-01, 8.278E-01, 1.047E+00, &
355	1.340E+00, 1.927E+00, 1.648E+00, 1.031E+00/
356	DATA A2_SSSSM_b2_19/1.219E-01, 1.337E-01, 1.633E-01, 1.929E-01, 2.057E-01, &
357	2.062E-01, 2.089E-01, 2.155E-01, 2.300E-01, 2.560E-01, &
358	2.908E-01, 3.199E-01, 3.530E-01, 3.965E-01, 4.475E-01, &
359	5.443E-01, 7.943E-01, 6.928E-01, 4.381E-01/
360	DATA A3_SSSSM_b2_19/-2.387E-03, -2.599E-03, -3.092E-03, -3.599E-03, -3.832E-03, &
361	-3.842E-03, -3.890E-03, -4.012E-03, -4.276E-03, -4.763E-03, &
362	-5.455E-03, -6.051E-03, -6.763E-03, -7.708E-03, -8.887E-03, &
363	-1.091E-02, -1.585E-02, -1.373E-02, -8.665E-03/
364	DATA B1_SSSSM_b2_19/1.260E-02, 1.211E-02, 1.126E-02, 1.056E-02, 1.038E-02, &
365	1.037E-02, 1.033E-02, 1.023E-02, 1.002E-02, 9.717E-03, &
366	9.613E-03, 9.652E-03, 9.983E-03, 1.047E-02, 1.168E-02, &
367	1.301E-02, 1.399E-02, 1.514E-02, 1.578E-02/
368	DATA B2_SSSSM_b2_19/2.336E-03, 2.419E-03, 2.506E-03, 2.610E-03, 2.690E-03, &
369	2.694E-03, 2.711E-03, 2.752E-03, 2.844E-03, 3.043E-03, &
370	3.455E-03, 3.871E-03, 4.507E-03, 5.373E-03, 6.786E-03, &
371	8.238E-03, 9.208E-03, 1.032E-02, 1.091E-02/
372	DATA C1_SSSSM_b2_19/7.832E-01, 7.787E-01, 7.721E-01, 7.670E-01, 7.657E-01, &
373	7.657E-01, 7.654E-01, 7.648E-01, 7.634E-01, 7.613E-01, &
374	7.596E-01, 7.585E-01, 7.574E-01, 7.560E-01, 7.533E-01, &
375	7.502E-01, 7.476E-01, 7.443E-01, 7.423E-01/
376	DATA C2_SSSSM_b2_19/3.144E-02, 3.268E-02, 3.515E-02, 3.748E-02, 3.837E-02, &
377	3.840E-02, 3.860E-02, 3.906E-02, 4.006E-02, 4.173E-02, &
378	4.338E-02, 4.435E-02, 4.459E-02, 4.467E-02, 4.202E-02, &
379	3.864E-02, 3.559E-02, 3.183E-02, 2.964E-02/
380	!***************************************************************************
381
382	spsol = 0
383	spinsol = 0
384	spss = 0
385
386	DATA alpha_aers_2bands/ &
387	! bc soluble
388	7.675,7.675,7.675,7.675,7.675,7.675, &
389	7.675,7.675,10.433,11.984,13.767,15.567,&
390	4.720,4.720,4.720,4.720,4.720,4.720, &
391	4.720,4.720,6.081,6.793,7.567,9.344, &
392	! pom soluble
393	5.503,5.503,5.503,5.503,5.588,5.957, &
394	6.404,7.340,8.545,10.319,13.595,20.398, &
395	1.402,1.402,1.402,1.402,1.431,1.562, &
396	1.715,2.032,2.425,2.991,4.193,7.133, &
397	! sulfate
398	4.681,5.062,5.460,5.798,6.224,6.733, &
399	7.556,8.613,10.687,12.265,16.32,21.692, &
400	1.107,1.239,1.381,1.490,1.635,1.8030, &
401	2.071,2.407,3.126,3.940,5.539,7.921, &
402	! sulfate coarse
403	4.681,5.062,5.460,5.798,6.224,6.733, &
404	7.556,8.613,10.687,12.265,16.32,21.692, &
405	1.107,1.239,1.381,1.490,1.635,1.8030, &
406	2.071,2.407,3.126,3.940,5.539,7.921, &
407	! seasalt Super Coarse Soluble (SS)
408	0.5090,0.6554,0.7129,0.7767,0.8529,1.2728, &
409	1.3820,1.5792,1.9173,2.2002,2.7173,4.1487, &
410	0.5167,0.6613,0.7221,0.7868,0.8622,1.3027, &
411	1.4227,1.6317,1.9887,2.2883,2.8356,4.3453, &
412	! seasalt Coarse Soluble (CS)
413	0.5090,0.6554,0.7129,0.7767,0.8529,1.2728, &
414	1.3820,1.5792,1.9173,2.2002,2.7173,4.1487, &
415	0.5167,0.6613,0.7221,0.7868,0.8622,1.3027, &
416	1.4227,1.6317,1.9887,2.2883,2.8356,4.3453, &
417	! seasalt Accumulation Soluble (AS)
418	4.125, 4.674, 5.005, 5.434, 5.985, 10.006, &
419	11.175,13.376,17.264,20.540,26.604, 42.349,&
420	4.187, 3.939, 3.919, 3.937, 3.995, 5.078, &
421	5.511, 6.434, 8.317,10.152,14.024, 26.537/
422
423	DATA alpha_aeri_2bands/ &
424	! dust insoluble
425	0.7661,0.7123,&
426	! bc insoluble
427	10.360,4.437, &
428	! pom insoluble
429	3.741,0.606/
430
431	DATA cg_aers_2bands/ &
432	! bc soluble
433	.612, .612, .612, .612, .612, .612, &
434	.612, .612, .702, .734, .760, .796, &
435	.433, .433, .433, .433, .433, .433, &
436	.433, .433, .534, .575, .613, .669, &
437	! pom soluble
438	.663, .663, .663, .663, .666, .674, &
439	.685, .702, .718, .737, .757, .777, &
440	.544, .544, .544, .544, .547, .554, &
441	.565, .583, .604, .631, .661, .698, &
442	! sulfate
443	.658, .669, .680, .688, .698, .707, &
444	.719, .733, .752, .760, .773, .786, &
445	.544, .555, .565, .573, .583, .593, &
446	.610, .628, .655, .666, .692, .719, &
447	! sulfate coarse
448	.658, .669, .680, .688, .698, .707, &
449	.719, .733, .752, .760, .773, .786, &
450	.544, .555, .565, .573, .583, .593, &
451	.610, .628, .655, .666, .692, .719, &
452	! seasalt Super Coarse soluble (SS)
453	.727, .747, .755, .761, .770, .788, &
454	.792, .799, .805, .809, .815, .826, &
455	.717, .738, .745, .752, .761, .779, &
456	.781, .786, .793, .797, .803, .813, &
457	! seasalt Coarse soluble (CS)
458	.727, .747, .755, .761, .770, .788, &
459	.792, .799, .805, .809, .815, .826, &
460	.717, .738, .745, .752, .761, .779, &
461	.781, .786, .793, .797, .803, .813, &
462	! Sesalt Accumulation Soluble (AS)
463	.727, .741, .748, .754, .761, .782, &
464	.787, .792, .797, .799, .801, .799, &
465	.606, .645, .658, .669, .681, .726, &
466	.734, .746, .761, .770, .782, .798/
467
468	DATA cg_aeri_2bands/ &
469	! dust insoluble
470	.701, .670, &
471	! bc insoluble
472	.471, .297, &
473	! pom insoluble
474	.568, .365/
475
476	DATA piz_aers_2bands/&
477	! bc soluble
478	.445, .445, .445, .445, .445, .445, &
479	.445, .445, .461, .480, .505, .528, &
480	.362, .362, .362, .362, .362, .362, &
481	.362, .362, .381, .405, .437, .483, &
482	! pom soluble
483	.972, .972, .972, .972, .972, .974, &
484	.976, .979, .982, .986, .989, .992, &
485	.924, .924, .924, .924, .925, .927, &
486	.932, .938, .945, .952, .961, .970, &
487	! sulfate
488	1.000,1.000,1.000,1.000,1.000,1.000, &
489	1.000,1.000,1.000,1.000,1.000,1.000, &
490	.992, .988, .988, .987, .986, .985, &
491	.985, .985, .984, .984, .984, .984, &
492	! sulfate coarse
493	1.000,1.000,1.000,1.000,1.000,1.000, &
494	1.000,1.000,1.000,1.000,1.000,1.000, &
495	.992, .988, .988, .987, .986, .985, &
496	.985, .985, .984, .984, .984, .984, &
497	! seasalt Super Coarse Soluble (SS)
498	1.000,1.000,1.000,1.000,1.000,1.000, &
499	1.000,1.000,1.000,1.000,1.000,1.000, &
500	0.992,0.989,0.987,0.986,0.986,0.980, &
501	0.980,0.978,0.976,0.976,0.974,0.971, &
502	! seasalt Coarse soluble (CS)
503	1.000,1.000,1.000,1.000,1.000,1.000, &
504	1.000,1.000,1.000,1.000,1.000,1.000, &
505	0.992,0.989,0.987,0.986,0.986,0.980, &
506	0.980,0.978,0.976,0.976,0.974,0.971, &
507	! seasalt Accumulation Soluble (AS)
508	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
509	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
510	0.970, 0.975, 0.976, 0.977, 0.978, 0.982, &
511	0.982, 0.983, 0.984, 0.984, 0.985, 0.985/
512
513	DATA piz_aeri_2bands/ &
514	! dust insoluble
515	.963, .987, &
516	! bc insoluble
517	.395, .264, &
518	! pom insoluble
519	.966, .859/
520
521	! Interpolation des coefficients optiques de 19 niveaux vers le nombre des niveaux du model
522	IF (firstcall) THEN
523	firstcall=.FALSE.
524
525	IF (.NOT. ALLOCATED(A1_ASSSM_b1)) THEN
526	ALLOCATE(A1_ASSSM_b1(klev),A2_ASSSM_b1(klev), A3_ASSSM_b1(klev),&
527	B1_ASSSM_b1(klev), B2_ASSSM_b1(klev), C1_ASSSM_b1(klev), C2_ASSSM_b1(klev),&
528	A1_CSSSM_b1(klev), A2_CSSSM_b1(klev), A3_CSSSM_b1(klev),&
529	B1_CSSSM_b1(klev), B2_CSSSM_b1(klev), C1_CSSSM_b1(klev), C2_CSSSM_b1(klev),&
530	A1_SSSSM_b1(klev), A2_SSSSM_b1(klev), A3_SSSSM_b1(klev),&
531	B1_SSSSM_b1(klev), B2_SSSSM_b1(klev), C1_SSSSM_b1(klev), C2_SSSSM_b1(klev),&
532	A1_ASSSM_b2(klev), A2_ASSSM_b2(klev), A3_ASSSM_b2(klev),&
533	B1_ASSSM_b2(klev), B2_ASSSM_b2(klev), C1_ASSSM_b2(klev), C2_ASSSM_b2(klev),&
534	A1_CSSSM_b2(klev), A2_CSSSM_b2(klev), A3_CSSSM_b2(klev),&
535	B1_CSSSM_b2(klev), B2_CSSSM_b2(klev), C1_CSSSM_b2(klev), C2_CSSSM_b2(klev),&
536	A1_SSSSM_b2(klev), A2_SSSSM_b2(klev), A3_SSSSM_b2(klev),&
537	B1_SSSSM_b2(klev), B2_SSSSM_b2(klev), C1_SSSSM_b2(klev), C2_SSSSM_b2(klev), stat=ierr)
538	IF (ierr /= 0) CALL abort_physic('aeropt_2bands', 'pb in allocation 1',1)
539	END IF
540
541	! bande 1
542	CALL pres2lev(A1_ASSSM_b1_19, A1_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
543	CALL pres2lev(A2_ASSSM_b1_19, A2_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
544	CALL pres2lev(A3_ASSSM_b1_19, A3_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
545	CALL pres2lev(B1_ASSSM_b1_19, B1_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
546	CALL pres2lev(B2_ASSSM_b1_19, B2_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
547	CALL pres2lev(C1_ASSSM_b1_19, C1_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
548	CALL pres2lev(C2_ASSSM_b1_19, C2_ASSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
549
550	CALL pres2lev(A1_CSSSM_b1_19, A1_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
551	CALL pres2lev(A2_CSSSM_b1_19, A2_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
552	CALL pres2lev(A3_CSSSM_b1_19, A3_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
553	CALL pres2lev(B1_CSSSM_b1_19, B1_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
554	CALL pres2lev(B2_CSSSM_b1_19, B2_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
555	CALL pres2lev(C1_CSSSM_b1_19, C1_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
556	CALL pres2lev(C2_CSSSM_b1_19, C2_CSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
557
558	CALL pres2lev(A1_SSSSM_b1_19, A1_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
559	CALL pres2lev(A2_SSSSM_b1_19, A2_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
560	CALL pres2lev(A3_SSSSM_b1_19, A3_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
561	CALL pres2lev(B1_SSSSM_b1_19, B1_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
562	CALL pres2lev(B2_SSSSM_b1_19, B2_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
563	CALL pres2lev(C1_SSSSM_b1_19, C1_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
564	CALL pres2lev(C2_SSSSM_b1_19, C2_SSSSM_b1, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
565
566	! bande 2
567	CALL pres2lev(A1_ASSSM_b2_19, A1_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
568	CALL pres2lev(A2_ASSSM_b2_19, A2_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
569	CALL pres2lev(A3_ASSSM_b2_19, A3_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
570	CALL pres2lev(B1_ASSSM_b2_19, B1_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
571	CALL pres2lev(B2_ASSSM_b2_19, B2_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
572	CALL pres2lev(C1_ASSSM_b2_19, C1_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
573	CALL pres2lev(C2_ASSSM_b2_19, C2_ASSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
574
575	CALL pres2lev(A1_CSSSM_b2_19, A1_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
576	CALL pres2lev(A2_CSSSM_b2_19, A2_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
577	CALL pres2lev(A3_CSSSM_b2_19, A3_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
578	CALL pres2lev(B1_CSSSM_b2_19, B1_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
579	CALL pres2lev(B2_CSSSM_b2_19, B2_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
580	CALL pres2lev(C1_CSSSM_b2_19, C1_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
581	CALL pres2lev(C2_CSSSM_b2_19, C2_CSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
582
583	CALL pres2lev(A1_SSSSM_b2_19, A1_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
584	CALL pres2lev(A2_SSSSM_b2_19, A2_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
585	CALL pres2lev(A3_SSSSM_b2_19, A3_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
586	CALL pres2lev(B1_SSSSM_b2_19, B1_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
587	CALL pres2lev(B2_SSSSM_b2_19, B2_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
588	CALL pres2lev(C1_SSSSM_b2_19, C1_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
589	CALL pres2lev(C2_SSSSM_b2_19, C2_SSSSM_b2, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
590
591	END IF ! firstcall
592
593
594	DO k=1, klev
595	DO i=1, klon
596	zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3
597	!CDIR UNROLL=naero_spc
598	mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
599	!RAF zrho
600	!CDIR UNROLL=naero_spc
601	mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9
602	zdp1(i,k)=pdel(i,k)/(gravitdelt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 s)]
603	ENDDO
604	ENDDO
605
606	IF (flag_aerosol .EQ. 1) THEN
607	nb_aer = 2
608	ALLOCATE (aerosol_name(nb_aer))
609	aerosol_name(1) = id_ASSO4M_phy
610	aerosol_name(2) = id_CSSO4M_phy
611	ELSEIF (flag_aerosol .EQ. 2) THEN
612	nb_aer = 2
613	ALLOCATE (aerosol_name(nb_aer))
614	aerosol_name(1) = id_ASBCM_phy
615	aerosol_name(2) = id_AIBCM_phy
616	ELSEIF (flag_aerosol .EQ. 3) THEN
617	nb_aer = 2
618	ALLOCATE (aerosol_name(nb_aer))
619	aerosol_name(1) = id_ASPOMM_phy
620	aerosol_name(2) = id_AIPOMM_phy
621	ELSEIF (flag_aerosol .EQ. 4) THEN
622	nb_aer = 3
623	ALLOCATE (aerosol_name(nb_aer))
624	aerosol_name(1) = id_CSSSM_phy
625	aerosol_name(2) = id_SSSSM_phy
626	aerosol_name(3) = id_ASSSM_phy
627	ELSEIF (flag_aerosol .EQ. 5) THEN
628	nb_aer = 1
629	ALLOCATE (aerosol_name(nb_aer))
630	aerosol_name(1) = id_CIDUSTM_phy
631	ELSEIF (flag_aerosol .EQ. 6) THEN
632	nb_aer = 10
633	ALLOCATE (aerosol_name(nb_aer))
634	aerosol_name(1) = id_ASSO4M_phy
635	aerosol_name(2) = id_ASBCM_phy
636	aerosol_name(3) = id_AIBCM_phy
637	aerosol_name(4) = id_ASPOMM_phy
638	aerosol_name(5) = id_AIPOMM_phy
639	aerosol_name(6) = id_CSSSM_phy
640	aerosol_name(7) = id_SSSSM_phy
641	aerosol_name(8) = id_ASSSM_phy
642	aerosol_name(9) = id_CIDUSTM_phy
643	aerosol_name(10)= id_CSSO4M_phy
644	ENDIF
645
646
647	!
648	! loop over modes, use of precalculated nmd and corresponding sigma
649	! loop over wavelengths
650	! for each mass species in mode
651	! interpolate from Sext to retrieve Sext_at_gridpoint_per_species
652	! compute optical_thickness_at_gridpoint_per_species
653
654
655
656	!!CDIR ON_ADB(RH_tab)
657	!CDIR ON_ADB(fact_RH)
658	!CDIR SHORTLOOP
659	DO n=1,nbre_RH-1
660	fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n))
661	ENDDO
662
663	DO k=1, KLEV
664	!!CDIR ON_ADB(RH_tab)
665	!CDIR ON_ADB(fact_RH)
666	DO i=1, KLON
667	rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
668	RH_num(i,k) = INT( rh(i,k)/10. + 1.)
669	IF (rh(i,k).GT.85.) RH_num(i,k)=10
670	IF (rh(i,k).GT.90.) RH_num(i,k)=11
671
672	DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k))
673	ENDDO
674	ENDDO
675
676	used_aer(:)=.FALSE.
677
678	DO m=1,nb_aer ! tau is only computed for each mass
679	fac=1.0
680	IF (aerosol_name(m).EQ.id_ASBCM_phy) THEN
681	soluble=.TRUE.
682	spsol=1
683	spss=0
684	ELSEIF (aerosol_name(m).EQ.id_ASPOMM_phy) THEN
685	soluble=.TRUE.
686	spsol=2
687	spss=0
688	ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
689	soluble=.TRUE.
690	spsol=3
691	spss=0
692	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
693	ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
694	soluble=.TRUE.
695	spsol=4
696	spss=0
697	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
698	ELSEIF (aerosol_name(m).EQ.id_SSSSM_phy) THEN
699	soluble=.TRUE.
700	spsol=5
701	spss=3
702	ELSEIF (aerosol_name(m).EQ.id_CSSSM_phy) THEN
703	soluble=.TRUE.
704	spsol=6
705	spss=2
706	ELSEIF (aerosol_name(m).EQ.id_ASSSM_phy) THEN
707	soluble=.TRUE.
708	spsol=7
709	spss=1
710	ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
711	soluble=.FALSE.
712	spinsol=1
713	spss=0
714	ELSEIF (aerosol_name(m).EQ.id_AIBCM_phy) THEN
715	soluble=.FALSE.
716	spinsol=2
717	spss=0
718	ELSEIF (aerosol_name(m).EQ.id_AIPOMM_phy) THEN
719	soluble=.FALSE.
720	spinsol=3
721	spss=0
722	ELSE
723	CYCLE
724	ENDIF
725
726	id=aerosol_name(m)
727	used_aer(id)=.TRUE.
728
729
730	IF (soluble) THEN
731
732	IF (spss.NE.0) THEN
733
734	IF (spss.EQ.1) THEN !accumulation mode
735	DO k=1, KLEV
736	!CDIR ON_ADB(A1_ASSSM_b1)
737	!CDIR ON_ADB(A2_ASSSM_b1)
738	!CDIR ON_ADB(A3_ASSSM_b1)
739	!CDIR ON_ADB(B1_ASSSM_b1)
740	!CDIR ON_ADB(B2_ASSSM_b1)
741	!CDIR ON_ADB(C1_ASSSM_b1)
742	!CDIR ON_ADB(C2_ASSSM_b2)
743	!CDIR ON_ADB(A1_ASSSM_b2)
744	!CDIR ON_ADB(A2_ASSSM_b2)
745	!CDIR ON_ADB(A3_ASSSM_b2)
746	!CDIR ON_ADB(B1_ASSSM_b2)
747	!CDIR ON_ADB(B2_ASSSM_b2)
748	!CDIR ON_ADB(C1_ASSSM_b2)
749	!CDIR ON_ADB(C2_ASSSM_b2)
750	DO i=1, KLON
751	H=rh(i,k)/100
752	tmp_var=mass_temp(i,k,spsol)1000.zdp1(i,k)deltfac
753	tmp_var_pi=mass_temp_pi(i,k,spsol)1000.zdp1(i,k)deltfac
754
755	! band 1
756	tau_ae2b_int=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05)
757	piz_ae2b_int=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H
758	cg_ae2b_int=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H
759
760	tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
761	tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int
762	piz_ae(i,k,id,1) = piz_ae2b_int
763	cg_ae(i,k,id,1)= cg_ae2b_int
764
765	!band 2
766	tau_ae2b_int=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05)
767	piz_ae2b_int=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H
768	cg_ae2b_int=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H
769
770	tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
771	tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int
772	piz_ae(i,k,id,2) = piz_ae2b_int
773	cg_ae(i,k,id,2)= cg_ae2b_int
774
775	ENDDO
776	ENDDO
777	ENDIF
778
779	IF (spss.EQ.2) THEN !coarse mode
780	DO k=1, KLEV
781	!CDIR ON_ADB(A1_CSSSM_b1)
782	!CDIR ON_ADB(A2_CSSSM_b1)
783	!CDIR ON_ADB(A3_CSSSM_b1)
784	!CDIR ON_ADB(B1_CSSSM_b1)
785	!CDIR ON_ADB(B2_CSSSM_b1)
786	!CDIR ON_ADB(C1_CSSSM_b1)
787	!CDIR ON_ADB(C2_CSSSM_b2)
788	!CDIR ON_ADB(A1_CSSSM_b2)
789	!CDIR ON_ADB(A2_CSSSM_b2)
790	!CDIR ON_ADB(A3_CSSSM_b2)
791	!CDIR ON_ADB(B1_CSSSM_b2)
792	!CDIR ON_ADB(B2_CSSSM_b2)
793	!CDIR ON_ADB(C1_CSSSM_b2)
794	!CDIR ON_ADB(C2_CSSSM_b2)
795	DO i=1, KLON
796	H=rh(i,k)/100
797	tmp_var=mass_temp(i,k,spsol)1000.zdp1(i,k)deltfac
798	tmp_var_pi=mass_temp_pi(i,k,spsol)1000.zdp1(i,k)deltfac
799	! band 1
800	tau_ae2b_int=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05)
801	piz_ae2b_int=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H
802	cg_ae2b_int=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H
803
804	tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
805	tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int
806	piz_ae(i,k,id,1) = piz_ae2b_int
807	cg_ae(i,k,id,1)= cg_ae2b_int
808
809	! band 2
810	tau_ae2b_int=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05)
811	piz_ae2b_int=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H
812	cg_ae2b_int=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H
813
814	tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
815	tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int
816	piz_ae(i,k,id,2) = piz_ae2b_int
817	cg_ae(i,k,id,2)= cg_ae2b_int
818
819	ENDDO
820	ENDDO
821	ENDIF
822
823	IF (spss.EQ.3) THEN !super coarse mode
824	DO k=1, KLEV
825	!CDIR ON_ADB(A1_SSSSM_b1)
826	!CDIR ON_ADB(A2_SSSSM_b1)
827	!CDIR ON_ADB(A3_SSSSM_b1)
828	!CDIR ON_ADB(B1_SSSSM_b1)
829	!CDIR ON_ADB(B2_SSSSM_b1)
830	!CDIR ON_ADB(C1_SSSSM_b1)
831	!CDIR ON_ADB(C2_SSSSM_b2)
832	!CDIR ON_ADB(A1_SSSSM_b2)
833	!CDIR ON_ADB(A2_SSSSM_b2)
834	!CDIR ON_ADB(A3_SSSSM_b2)
835	!CDIR ON_ADB(B1_SSSSM_b2)
836	!CDIR ON_ADB(B2_SSSSM_b2)
837	!CDIR ON_ADB(C1_SSSSM_b2)
838	!CDIR ON_ADB(C2_SSSSM_b2)
839	DO i=1, KLON
840	H=rh(i,k)/100
841	tmp_var=mass_temp(i,k,spsol)1000.zdp1(i,k)deltfac
842	tmp_var_pi=mass_temp_pi(i,k,spsol)1000.zdp1(i,k)deltfac
843
844	! band 1
845	tau_ae2b_int=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05)
846	piz_ae2b_int=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H
847	cg_ae2b_int=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H
848
849	tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
850	tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int
851	piz_ae(i,k,id,1) = piz_ae2b_int
852	cg_ae(i,k,id,1)= cg_ae2b_int
853
854	! band 2
855	tau_ae2b_int=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05)
856	piz_ae2b_int=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H
857	cg_ae2b_int=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H
858
859	tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
860	tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int
861	piz_ae(i,k,id,2) = piz_ae2b_int
862	cg_ae(i,k,id,2)= cg_ae2b_int
863
864	ENDDO
865	ENDDO
866	ENDIF
867
868	ELSE
869
870	!CDIR ON_ADB(alpha_aers_2bands)
871	!CDIR ON_ADB(piz_aers_2bands)
872	!CDIR ON_ADB(cg_aers_2bands)
873	DO k=1, KLEV
874	DO i=1, KLON
875	tmp_var=mass_temp(i,k,spsol)1000.zdp1(i,k)deltfac
876	tmp_var_pi=mass_temp_pi(i,k,spsol)1000.zdp1(i,k)deltfac
877	!CDIR UNROLL=nbands
878	DO inu=1,nbands
879
880	tau_ae2b_int= alpha_aers_2bands(RH_num(i,k),inu,spsol)+ &
881	DELTA(i,k)* (alpha_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
882	alpha_aers_2bands(RH_num(i,k),inu,spsol))
883
884	piz_ae2b_int = piz_aers_2bands(RH_num(i,k),inu,spsol) + &
885	DELTA(i,k)* (piz_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
886	piz_aers_2bands(RH_num(i,k),inu,spsol))
887
888	cg_ae2b_int = cg_aers_2bands(RH_num(i,k),inu,spsol) + &
889	DELTA(i,k)* (cg_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
890	cg_aers_2bands(RH_num(i,k),inu,spsol))
891
892	tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
893	tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
894	piz_ae(i,k,id,inu) = piz_ae2b_int
895	cg_ae(i,k,id,inu)= cg_ae2b_int
896
897	ENDDO
898	ENDDO
899	ENDDO
900
901	ENDIF
902
903	ELSE ! For all aerosol insoluble components
904
905	!CDIR ON_ADB(alpha_aers_2bands)
906	!CDIR ON_ADB(piz_aers_2bands)
907	!CDIR ON_ADB(cg_aers_2bands)
908	DO k=1, KLEV
909	DO i=1, KLON
910	tmp_var=mass_temp(i,k,naero_soluble+ spinsol)1000.zdp1(i,k)deltfac
911	tmp_var_pi=mass_temp_pi(i,k,naero_soluble+spinsol)1000.zdp1(i,k)deltfac
912	!CDIR UNROLL=nbands
913	DO inu=1,nbands
914	tau_ae2b_int = alpha_aeri_2bands(inu,spinsol)
915	piz_ae2b_int = piz_aeri_2bands(inu,spinsol)
916	cg_ae2b_int = cg_aeri_2bands(inu,spinsol)
917
918	tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
919	tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int
920	piz_ae(i,k,id,inu) = piz_ae2b_int
921	cg_ae(i,k,id,inu)= cg_ae2b_int
922	ENDDO
923	ENDDO
924	ENDDO
925
926	ENDIF ! soluble
927
928	ENDDO ! nb_aer
929
930	!correction bug OB
931	! DO m=1,nb_aer
932	DO m=1,naero_tot
933	IF (.NOT. used_aer(m)) THEN
934	tau_ae(:,:,m,:)=0.
935	tau_ae_pi(:,:,m,:)=0.
936	piz_ae(:,:,m,:)=0.
937	cg_ae(:,:,m,:)=0.
938	ENDIF
939	ENDDO
940
941	DO inu=1, nbands
942	DO mrfspecies=1,naero_grp
943	IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER
944	DO k=1, KLEV
945	DO i=1, KLON
946	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)+ &
947	tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+ &
948	tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ &
949	tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+ &
950	tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu)
951	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
952
953	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
954	tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ &
955	tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ &
956	tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ &
957	tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ &
958	tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ &
959	tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ &
960	tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ &
961	tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ &
962	tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
963	/tau_allaer(i,k,mrfspecies,inu)
964	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
965
966	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)* &
967	piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ &
968	tau_ae(i,k,id_CSSO4M_phy,inu)piz_ae(i,k,id_CSSO4M_phy,inu) &
969	cg_ae(i,k,id_CSSO4M_phy,inu)+ &
970	tau_ae(i,k,id_ASBCM_phy,inu)piz_ae(i,k,id_ASBCM_phy,inu) &
971	cg_ae(i,k,id_ASBCM_phy,inu)+ &
972	tau_ae(i,k,id_AIBCM_phy,inu)piz_ae(i,k,id_AIBCM_phy,inu) &
973	cg_ae(i,k,id_AIBCM_phy,inu)+ &
974	tau_ae(i,k,id_ASPOMM_phy,inu)piz_ae(i,k,id_ASPOMM_phy,inu) &
975	cg_ae(i,k,id_ASPOMM_phy,inu)+ &
976	tau_ae(i,k,id_AIPOMM_phy,inu)piz_ae(i,k,id_AIPOMM_phy,inu) &
977	cg_ae(i,k,id_AIPOMM_phy,inu)+ &
978	tau_ae(i,k,id_ASSSM_phy,inu)piz_ae(i,k,id_ASSSM_phy,inu) &
979	cg_ae(i,k,id_ASSSM_phy,inu)+ &
980	tau_ae(i,k,id_CSSSM_phy,inu)piz_ae(i,k,id_CSSSM_phy,inu) &
981	cg_ae(i,k,id_CSSSM_phy,inu)+ &
982	tau_ae(i,k,id_SSSSM_phy,inu)piz_ae(i,k,id_SSSSM_phy,inu) &
983	cg_ae(i,k,id_SSSSM_phy,inu)+ &
984	tau_ae(i,k,id_CIDUSTM_phy,inu)piz_ae(i,k,id_CIDUSTM_phy,inu) &
985	cg_ae(i,k,id_CIDUSTM_phy,inu))/ &
986	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
987	ENDDO
988	ENDDO
989
990	ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT
991
992	DO k=1, KLEV
993	DO i=1, KLON
994	!RAF
995	tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+ &
996	tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ &
997	tau_ae_pi(i,k,id_ASBCM_phy,inu)+ &
998	tau_ae_pi(i,k,id_AIBCM_phy,inu)+ &
999	tau_ae_pi(i,k,id_ASPOMM_phy,inu)+ &
1000	tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ &
1001	tau_ae_pi(i,k,id_ASSSM_phy,inu)+ &
1002	tau_ae_pi(i,k,id_CSSSM_phy,inu)+ &
1003	tau_ae_pi(i,k,id_SSSSM_phy,inu)+ &
1004	tau_ae_pi(i,k,id_CIDUSTM_phy,inu)
1005	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1006
1007	piz_allaer(i,k,mrfspecies,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
1008	tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ &
1009	tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ &
1010	tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ &
1011	tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ &
1012	tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ &
1013	tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ &
1014	tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ &
1015	tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ &
1016	tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
1017	/tau_allaer(i,k,mrfspecies,inu)
1018	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1019
1020	cg_allaer(i,k,mrfspecies,inu)=(&
1021	tau_ae_pi(i,k,id_ASSO4M_phy,inu)piz_ae(i,k,id_ASSO4M_phy,inu)cg_ae(i,k,id_ASSO4M_phy,inu)+ &
1022	tau_ae_pi(i,k,id_CSSO4M_phy,inu)piz_ae(i,k,id_CSSO4M_phy,inu)cg_ae(i,k,id_CSSO4M_phy,inu)+ &
1023	tau_ae_pi(i,k,id_ASBCM_phy,inu)piz_ae(i,k,id_ASBCM_phy,inu)cg_ae(i,k,id_ASBCM_phy,inu)+ &
1024	tau_ae_pi(i,k,id_AIBCM_phy,inu)piz_ae(i,k,id_AIBCM_phy,inu)cg_ae(i,k,id_AIBCM_phy,inu)+ &
1025	tau_ae_pi(i,k,id_ASPOMM_phy,inu)piz_ae(i,k,id_ASPOMM_phy,inu)cg_ae(i,k,id_ASPOMM_phy,inu)+ &
1026	tau_ae_pi(i,k,id_AIPOMM_phy,inu)piz_ae(i,k,id_AIPOMM_phy,inu)cg_ae(i,k,id_AIPOMM_phy,inu)+ &
1027	tau_ae_pi(i,k,id_ASSSM_phy,inu)piz_ae(i,k,id_ASSSM_phy,inu)cg_ae(i,k,id_ASSSM_phy,inu)+ &
1028	tau_ae_pi(i,k,id_CSSSM_phy,inu)piz_ae(i,k,id_CSSSM_phy,inu)cg_ae(i,k,id_CSSSM_phy,inu)+ &
1029	tau_ae_pi(i,k,id_SSSSM_phy,inu)piz_ae(i,k,id_SSSSM_phy,inu)cg_ae(i,k,id_SSSSM_phy,inu)+ &
1030	tau_ae_pi(i,k,id_CIDUSTM_phy,inu)piz_ae(i,k,id_CIDUSTM_phy,inu)&
1031	cg_ae(i,k,id_CIDUSTM_phy,inu))/ &
1032	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1033	ENDDO
1034	ENDDO
1035
1036	ELSEIF (mrfspecies .EQ. 4) THEN ! = BC
1037	DO k=1, KLEV
1038	DO i=1, KLON
1039	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)
1040	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1041	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu) &
1042	+tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu))/ &
1043	tau_allaer(i,k,mrfspecies,inu)
1044	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1045	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM_phy,inu)piz_ae(i,k,id_ASBCM_phy,inu) cg_ae(i,k,id_ASBCM_phy,inu)&
1046	+tau_ae(i,k,id_AIBCM_phy,inu)piz_ae(i,k,id_AIBCM_phy,inu)cg_ae(i,k,id_AIBCM_phy,inu))/ &
1047	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1048	ENDDO
1049	ENDDO
1050
1051	ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4
1052
1053	DO k=1, KLEV
1054	DO i=1, KLON
1055	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)
1056	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1057	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu) &
1058	+tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu))/ &
1059	tau_allaer(i,k,mrfspecies,inu)
1060	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1061	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M_phy,inu)* &
1062	piz_ae(i,k,id_CSSO4M_phy,inu) *cg_ae(i,k,id_CSSO4M_phy,inu)&
1063	+tau_ae(i,k,id_ASSO4M_phy,inu)piz_ae(i,k,id_ASSO4M_phy,inu) &
1064	cg_ae(i,k,id_ASSO4M_phy,inu))/ &
1065	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1066	ENDDO
1067	ENDDO
1068
1069	ELSEIF (mrfspecies .EQ. 6) THEN ! = POM
1070
1071	DO k=1, KLEV
1072	DO i=1, KLON
1073	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)
1074	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1075	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu) &
1076	+tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu))/ &
1077	tau_allaer(i,k,mrfspecies,inu)
1078	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1079	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM_phy,inu)piz_ae(i,k,id_ASPOMM_phy,inu) cg_ae(i,k,id_ASPOMM_phy,inu)&
1080	+tau_ae(i,k,id_AIPOMM_phy,inu)piz_ae(i,k,id_AIPOMM_phy,inu)cg_ae(i,k,id_AIPOMM_phy,inu))/ &
1081	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1082	ENDDO
1083	ENDDO
1084
1085	ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST
1086
1087	DO k=1, KLEV
1088	DO i=1, KLON
1089	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM_phy,inu)
1090	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1091	piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM_phy,inu)
1092	cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM_phy,inu)
1093	ENDDO
1094	ENDDO
1095
1096	ELSEIF (mrfspecies .EQ. 8) THEN ! = SS
1097
1098	DO k=1, KLEV
1099	DO i=1, KLON
1100	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+tau_ae(i,k,id_SSSSM_phy,inu)
1101	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1102	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu) &
1103	+tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu) &
1104	+tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu))/ &
1105	tau_allaer(i,k,mrfspecies,inu)
1106	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1107	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM_phy,inu)piz_ae(i,k,id_ASSSM_phy,inu) cg_ae(i,k,id_ASSSM_phy,inu)&
1108	+tau_ae(i,k,id_CSSSM_phy,inu)piz_ae(i,k,id_CSSSM_phy,inu)cg_ae(i,k,id_CSSSM_phy,inu) &
1109	+tau_ae(i,k,id_SSSSM_phy,inu)piz_ae(i,k,id_SSSSM_phy,inu)cg_ae(i,k,id_SSSSM_phy,inu))/ &
1110	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1111	ENDDO
1112	ENDDO
1113
1114	ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3
1115
1116	DO k=1, KLEV
1117	DO i=1, KLON
1118	tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary
1119	piz_allaer(i,k,mrfspecies,inu)=0.
1120	cg_allaer(i,k,mrfspecies,inu)=0.
1121	ENDDO
1122	ENDDO
1123
1124	ELSE
1125
1126	DO k=1, KLEV
1127	DO i=1, KLON
1128	tau_allaer(i,k,mrfspecies,inu)=0.
1129	piz_allaer(i,k,mrfspecies,inu)=0.
1130	cg_allaer(i,k,mrfspecies,inu)=0.
1131	ENDDO
1132	ENDDO
1133
1134	ENDIF
1135
1136	ENDDO
1137	ENDDO
1138
1139
1140	inu=1
1141	DO i=1, KLON
1142	absvisaer(i)=SUM((1-piz_allaer(i,:,:,inu))*tau_allaer(i,:,:,inu))
1143	END DO
1144
1145	DEALLOCATE(aerosol_name)
1146
1147	END SUBROUTINE AEROPT_2BANDS

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

Download in other formats:

Original Format