Context Navigation

aeropt_2bands.F90 @ 5085

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

Nouvelle version qui inclut les effets des aérosols et propose les mêmes diagnostics des effets
directs et indirects que l'ancienne version du rayonnement.
OB

New RRTM version that includes the effects of aerosols and outputs the same direct and indirect effects
diagnostics as the old version
OB

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory
Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 53.3 KB

Line
1	!
2	! $Id: aeropt_2bands.F90 2003 2014-04-04 12:51:02Z fairhead $
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_spc), INTENT(in) :: m_allaer
32	!RAF
33	REAL, DIMENSION(klon,klev,naero_spc), 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_spc) :: mass_temp
139	!RAF
140	REAL, DIMENSION(klon,klev,naero_spc) :: 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_gcm('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
610	aerosol_name(2) = id_CSSO4M
611	ELSEIF (flag_aerosol .EQ. 2) THEN
612	nb_aer = 2
613	ALLOCATE (aerosol_name(nb_aer))
614	aerosol_name(1) = id_ASBCM
615	aerosol_name(2) = id_AIBCM
616	ELSEIF (flag_aerosol .EQ. 3) THEN
617	nb_aer = 2
618	ALLOCATE (aerosol_name(nb_aer))
619	aerosol_name(1) = id_ASPOMM
620	aerosol_name(2) = id_AIPOMM
621	ELSEIF (flag_aerosol .EQ. 4) THEN
622	nb_aer = 3
623	ALLOCATE (aerosol_name(nb_aer))
624	aerosol_name(1) = id_CSSSM
625	aerosol_name(2) = id_SSSSM
626	aerosol_name(3) = id_ASSSM
627	ELSEIF (flag_aerosol .EQ. 5) THEN
628	nb_aer = 1
629	ALLOCATE (aerosol_name(nb_aer))
630	aerosol_name(1) = id_CIDUSTM
631	ELSEIF (flag_aerosol .EQ. 6) THEN
632	nb_aer = 10
633	ALLOCATE (aerosol_name(nb_aer))
634	aerosol_name(1) = id_ASSO4M
635	aerosol_name(2) = id_ASBCM
636	aerosol_name(3) = id_AIBCM
637	aerosol_name(4) = id_ASPOMM
638	aerosol_name(5) = id_AIPOMM
639	aerosol_name(6) = id_CSSSM
640	aerosol_name(7) = id_SSSSM
641	aerosol_name(8) = id_ASSSM
642	aerosol_name(9) = id_CIDUSTM
643	aerosol_name(10)= id_CSSO4M
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) THEN
681	soluble=.TRUE.
682	spsol=1
683	spss=0
684	ELSEIF (aerosol_name(m).EQ.id_ASPOMM) THEN
685	soluble=.TRUE.
686	spsol=2
687	spss=0
688	ELSEIF (aerosol_name(m).EQ.id_ASSO4M) 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) 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) THEN
699	soluble=.TRUE.
700	spsol=5
701	spss=3
702	ELSEIF (aerosol_name(m).EQ.id_CSSSM) THEN
703	soluble=.TRUE.
704	spsol=6
705	spss=2
706	ELSEIF (aerosol_name(m).EQ.id_ASSSM) THEN
707	soluble=.TRUE.
708	spsol=7
709	spss=1
710	ELSEIF (aerosol_name(m).EQ.id_CIDUSTM) THEN
711	soluble=.FALSE.
712	spinsol=1
713	spss=0
714	ELSEIF (aerosol_name(m).EQ.id_AIBCM) THEN
715	soluble=.FALSE.
716	spinsol=2
717	spss=0
718	ELSEIF (aerosol_name(m).EQ.id_AIPOMM) 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,inu)+tau_ae(i,k,id_CSSO4M,inu)+ &
947	tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ &
948	tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ &
949	tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+ &
950	tau_ae(i,k,id_SSSSM,inu)+ tau_ae(i,k,id_CIDUSTM,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,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
954	tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
955	tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
956	tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
957	tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
958	tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &
959	tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
960	tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
961	tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
962	tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,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,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu)+ &
967	tau_ae(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu)cg_ae(i,k,id_CSSO4M,inu)+ &
968	tau_ae(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu)cg_ae(i,k,id_ASBCM,inu)+ &
969	tau_ae(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu)+ &
970	tau_ae(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu)cg_ae(i,k,id_ASPOMM,inu)+ &
971	tau_ae(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu)+ &
972	tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu)+ &
973	tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu)+ &
974	tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu)+ &
975	tau_ae(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu)cg_ae(i,k,id_CIDUSTM,inu))/ &
976	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
977	ENDDO
978	ENDDO
979
980	ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT
981
982	DO k=1, KLEV
983	DO i=1, KLON
984	!RAF
985	tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M,inu)+ &
986	tau_ae_pi(i,k,id_CSSO4M,inu)+ &
987	tau_ae_pi(i,k,id_ASBCM,inu)+ &
988	tau_ae_pi(i,k,id_AIBCM,inu)+ &
989	tau_ae_pi(i,k,id_ASPOMM,inu)+ &
990	tau_ae_pi(i,k,id_AIPOMM,inu)+ &
991	tau_ae_pi(i,k,id_ASSSM,inu)+ &
992	tau_ae_pi(i,k,id_CSSSM,inu)+ &
993	tau_ae_pi(i,k,id_SSSSM,inu)+ &
994	tau_ae_pi(i,k,id_CIDUSTM,inu)
995	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
996
997	piz_allaer(i,k,mrfspecies,inu)=(tau_ae_pi(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
998	tau_ae_pi(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
999	tau_ae_pi(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
1000	tau_ae_pi(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
1001	tau_ae_pi(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
1002	tau_ae_pi(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &
1003	tau_ae_pi(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
1004	tau_ae_pi(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
1005	tau_ae_pi(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
1006	tau_ae_pi(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &
1007	/tau_allaer(i,k,mrfspecies,inu)
1008	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1009
1010	cg_allaer(i,k,mrfspecies,inu)=(&
1011	tau_ae_pi(i,k,id_ASSO4M,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu)+ &
1012	tau_ae_pi(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu)cg_ae(i,k,id_CSSO4M,inu)+ &
1013	tau_ae_pi(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu)cg_ae(i,k,id_ASBCM,inu)+ &
1014	tau_ae_pi(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu)+ &
1015	tau_ae_pi(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu)cg_ae(i,k,id_ASPOMM,inu)+ &
1016	tau_ae_pi(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu)+ &
1017	tau_ae_pi(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu)cg_ae(i,k,id_ASSSM,inu)+ &
1018	tau_ae_pi(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu)+ &
1019	tau_ae_pi(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu)+ &
1020	tau_ae_pi(i,k,id_CIDUSTM,inu)piz_ae(i,k,id_CIDUSTM,inu)&
1021	cg_ae(i,k,id_CIDUSTM,inu))/ &
1022	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1023	ENDDO
1024	ENDDO
1025
1026	ELSEIF (mrfspecies .EQ. 4) THEN ! = BC
1027	DO k=1, KLEV
1028	DO i=1, KLON
1029	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)
1030	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1031	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) &
1032	+tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ &
1033	tau_allaer(i,k,mrfspecies,inu)
1034	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1035	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)piz_ae(i,k,id_ASBCM,inu) cg_ae(i,k,id_ASBCM,inu)&
1036	+tau_ae(i,k,id_AIBCM,inu)piz_ae(i,k,id_AIBCM,inu)cg_ae(i,k,id_AIBCM,inu))/ &
1037	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1038	ENDDO
1039	ENDDO
1040
1041	ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4
1042
1043	DO k=1, KLEV
1044	DO i=1, KLON
1045	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)
1046	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1047	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) &
1048	+tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ &
1049	tau_allaer(i,k,mrfspecies,inu)
1050	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1051	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)piz_ae(i,k,id_CSSO4M,inu) cg_ae(i,k,id_CSSO4M,inu)&
1052	+tau_ae(i,k,id_ASSO4M,inu)piz_ae(i,k,id_ASSO4M,inu)cg_ae(i,k,id_ASSO4M,inu))/ &
1053	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1054	ENDDO
1055	ENDDO
1056
1057	ELSEIF (mrfspecies .EQ. 6) THEN ! = POM
1058
1059	DO k=1, KLEV
1060	DO i=1, KLON
1061	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)
1062	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1063	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) &
1064	+tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ &
1065	tau_allaer(i,k,mrfspecies,inu)
1066	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1067	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)piz_ae(i,k,id_ASPOMM,inu) cg_ae(i,k,id_ASPOMM,inu)&
1068	+tau_ae(i,k,id_AIPOMM,inu)piz_ae(i,k,id_AIPOMM,inu)cg_ae(i,k,id_AIPOMM,inu))/ &
1069	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1070	ENDDO
1071	ENDDO
1072
1073	ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST
1074
1075	DO k=1, KLEV
1076	DO i=1, KLON
1077	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu)
1078	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1079	piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu)
1080	cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu)
1081	ENDDO
1082	ENDDO
1083
1084	ELSEIF (mrfspecies .EQ. 8) THEN ! = SS
1085
1086	DO k=1, KLEV
1087	DO i=1, KLON
1088	tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)
1089	tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5)
1090	piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) &
1091	+tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) &
1092	+tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
1093	tau_allaer(i,k,mrfspecies,inu)
1094	piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1)
1095	cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)piz_ae(i,k,id_ASSSM,inu) cg_ae(i,k,id_ASSSM,inu)&
1096	+tau_ae(i,k,id_CSSSM,inu)piz_ae(i,k,id_CSSSM,inu)cg_ae(i,k,id_CSSSM,inu) &
1097	+tau_ae(i,k,id_SSSSM,inu)piz_ae(i,k,id_SSSSM,inu)cg_ae(i,k,id_SSSSM,inu))/ &
1098	(tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
1099	ENDDO
1100	ENDDO
1101
1102	ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3
1103
1104	DO k=1, KLEV
1105	DO i=1, KLON
1106	tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary
1107	piz_allaer(i,k,mrfspecies,inu)=0.
1108	cg_allaer(i,k,mrfspecies,inu)=0.
1109	ENDDO
1110	ENDDO
1111
1112	ELSE
1113
1114	DO k=1, KLEV
1115	DO i=1, KLON
1116	tau_allaer(i,k,mrfspecies,inu)=0.
1117	piz_allaer(i,k,mrfspecies,inu)=0.
1118	cg_allaer(i,k,mrfspecies,inu)=0.
1119	ENDDO
1120	ENDDO
1121
1122	ENDIF
1123
1124	ENDDO
1125	ENDDO
1126
1127
1128	inu=1
1129	DO i=1, KLON
1130	absvisaer(i)=SUM((1-piz_allaer(i,:,:,inu))*tau_allaer(i,:,:,inu))
1131	END DO
1132
1133	DEALLOCATE(aerosol_name)
1134
1135	END SUBROUTINE AEROPT_2BANDS

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

Download in other formats:

Original Format