source: LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90 @ 2736

Last change on this file since 2736 was 2736, checked in by oboucher, 8 years ago

Complement to previous commit. Fixing an issue with mass array index.
This does not have any impact because ASSO4M and CSSO4M have same properties.
  • 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:keywords set to Author Date Id Revi
File size: 44.2 KB
Line 
1!
2! $Id: aeropt_6bands_rrtm.F90 2736 2016-12-17 19:48:54Z oboucher $
3!
4SUBROUTINE AEROPT_6BANDS_RRTM ( &
5     pdel, m_allaer, RHcl, tau_allaer, piz_allaer, &
6     cg_allaer, m_allaer_pi, flag_aerosol, flag_bc_internal_mixture, zrho )
7
8  USE dimphy
9  USE aero_mod
10  USE phys_local_var_mod, only: absvisaer
11
12  !    Yves Balkanski le 12 avril 2006
13  !    Celine Deandreis
14  !    Anne Cozic Avril 2009
15  !    a partir d'une sous-routine de Johannes Quaas pour les sulfates
16  !    Olivier Boucher février 2014 pour passage à RRTM
17  !    a partir des propriétés optiques fournies par Yves Balkanski
18  !
19  IMPLICIT NONE
20
21  INCLUDE "YOMCST.h"
22  INCLUDE "clesphys.h"
23
24  !
25  ! Input arguments:
26  !
27  REAL, DIMENSION(klon,klev),     INTENT(in)  :: pdel
28  REAL, DIMENSION(klon,klev,naero_tot),   INTENT(in)  :: m_allaer
29  REAL, DIMENSION(klon,klev,naero_tot),   INTENT(in)  :: m_allaer_pi
30  REAL, DIMENSION(klon,klev),     INTENT(in)  :: RHcl       ! humidite relative ciel clair
31  INTEGER,                        INTENT(in)  :: flag_aerosol
32  LOGICAL,                        INTENT(in)  :: flag_bc_internal_mixture
33  REAL, DIMENSION(klon,klev),     INTENT(in)  :: zrho
34  !
35  ! Output arguments:
36  ! 2= total aerosols
37  ! 1= natural aerosols
38  !
39  REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: tau_allaer ! epaisseur optique aerosol
40  REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: piz_allaer ! single scattering albedo aerosol
41  REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: cg_allaer  ! asymmetry parameter aerosol
42
43  !
44  ! Local
45  !
46  LOGICAL ::  soluble
47  INTEGER :: i, k,n, inu, m
48  INTEGER :: spsol, spinsol
49  INTEGER :: RH_num(klon,klev)
50
51  INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
52
53  INTEGER, PARAMETER :: naero_soluble=9    ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
54                                           ! 5- seasalt super coarse  6- seasalt coarse  7- seasalt acc.
55                                           ! 8- NO3 coarse 9- NO3 acc.
56  INTEGER, PARAMETER :: naero_insoluble=4  ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
57
58  INTEGER, PARAMETER :: nbre_RH=12
59  REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
60  REAL, PARAMETER :: RH_MAX=95.
61  REAL :: delta(klon,klev), rh(klon,klev)
62  REAL :: tau_ae2b_int   ! Intermediate computation of epaisseur optique aerosol
63  REAL :: piz_ae2b_int   ! Intermediate computation of Single scattering albedo
64  REAL :: cg_ae2b_int    ! Intermediate computation of Assymetry parameter
65  REAL :: Fact_RH(nbre_RH)
66  REAL :: fac
67  REAL :: zdh(klon,klev)
68  INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name
69  INTEGER :: nb_aer
70
71  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae
72  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: piz_ae
73  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: cg_ae
74
75  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae_pi
76  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: piz_ae_pi
77  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: cg_ae_pi
78  !
79  ! Proprietes optiques
80  !
81  REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit m2/g
82  REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)       !--unit m2/g
83  REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)    !--unitless
84  REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble)          !--unitless
85  REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !-- unitless
86  REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble)         !-- unitless
87  !
88  ! BC internal mixture
89  !
90  INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
91                                        ! of BC in the soluble mode:
92                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
93  ! for Maxwell-Garnet internal mixture
94  ! Detailed theory can be found in R. Wang Estimation of global black carbon direct
95  ! radiative forcing and its uncertainty constrained by observations. J.
96  ! Geophys. Res. Atmos. Added by R. Wang and OB
97  REAL :: alpha_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
98  REAL :: cg_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
99  REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
100  !
101  INTEGER :: aerindex, classbc, classbc_pi
102  REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
103
104  !
105  REAL, PARAMETER :: tau_min = 1.e-15
106!  REAL, PARAMETER :: tau_min = 1.e-7
107
108!***************************************************************************
109!--the order of the soluble   species has to follow the spsol   index below
110!--the order of the insoluble species has to follow the spinsol index below
111
112  DATA alpha_aers_6bands/  &
113       ! bc soluble AS
114  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
115  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
116  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
117  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
118  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
119  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, &
120       ! pom soluble AS
121  6.443, 6.443, 6.443, 6.443, 6.443, 7.100, 7.809, 9.277,10.721,10.882,11.056,12.628, &
122  6.443, 6.443, 6.443, 6.443, 6.443, 7.100, 7.809, 9.277,10.721,10.882,11.056,12.628, &
123  4.381, 4.381, 4.381, 4.381, 4.381, 4.828, 5.310, 6.309, 7.290, 7.400, 7.518, 8.587, &
124  1.846, 1.846, 1.846, 1.846, 1.846, 2.034, 2.237, 2.658, 3.072, 3.118, 3.168, 3.618, &
125  0.377, 0.377, 0.377, 0.377, 0.377, 0.415, 0.456, 0.542, 0.627, 0.636, 0.646, 0.738, &
126  0.052, 0.052, 0.052, 0.052, 0.052, 0.057, 0.063, 0.075, 0.087, 0.088, 0.089, 0.102, &
127       ! sulfate coarse CS
128  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
129  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
130  0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, &
131  0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, &
132  0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
133  0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
134       ! sulfate AS   
135  6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
136  6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
137  4.381, 4.381, 4.381, 4.828, 5.301, 5.792, 6.309, 7.176, 9.542,11.868,11.868,21.204, &
138  1.727, 1.727, 1.727, 1.903, 2.090, 2.283, 2.487, 2.829, 3.762, 4.679, 4.679, 8.359, &
139  0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, &
140  0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, &
141        ! Seasalt Super Coarse Soluble (SS)
142  0.214, 0.267, 0.287, 0.310, 0.337, 0.373, 0.421, 0.494, 0.625, 0.742, 0.956, 1.480, &
143  0.217, 0.270, 0.291, 0.314, 0.341, 0.377, 0.426, 0.499, 0.632, 0.751, 0.963, 1.490, &
144  0.221, 0.275, 0.297, 0.320, 0.348, 0.384, 0.434, 0.507, 0.642, 0.762, 0.976, 1.506, &
145  0.230, 0.285, 0.308, 0.331, 0.359, 0.396, 0.447, 0.522, 0.658, 0.780, 0.997, 1.536, &
146  0.250, 0.307, 0.330, 0.354, 0.384, 0.424, 0.477, 0.556, 0.696, 0.822, 1.044, 1.592, &
147  0.279, 0.347, 0.373, 0.401, 0.434, 0.478, 0.537, 0.625, 0.781, 0.918, 1.158, 1.744, &
148        ! Seasalt Coarse Soluble (CS)     
149  0.550, 0.676, 0.724, 0.779, 0.841, 0.925, 1.040, 1.213, 1.523, 1.803, 2.306, 3.535, &
150  0.568, 0.695, 0.744, 0.798, 0.863, 0.950, 1.066, 1.240, 1.556, 1.839, 2.342, 3.588, &
151  0.599, 0.727, 0.779, 0.834, 0.901, 0.993, 1.111, 1.290, 1.612, 1.899, 2.411, 3.669, &
152  0.647, 0.786, 0.841, 0.899, 0.969, 1.069, 1.193, 1.384, 1.717, 2.015, 2.543, 3.842, &
153  0.663, 0.825, 0.889, 0.957, 1.038, 1.126, 1.268, 1.483, 1.862, 2.196, 2.780, 4.183, &
154  0.479, 0.644, 0.711, 0.785, 0.873, 0.904, 1.048, 1.275, 1.692, 2.072, 2.753, 4.430, &
155        ! Seasalt Accumulation Soluble (AS)
156  6.128, 7.564, 8.127, 8.723, 9.421,10.014,11.250,13.113,16.320,19.061,23.675,34.158, &
157  5.080, 6.586, 7.197, 7.856, 8.645, 8.820,10.115,12.140,15.775,18.983,24.499,36.993, &
158  3.290, 4.541, 5.073, 5.663, 6.387, 6.227, 7.339, 9.161,12.645,15.918,21.908,36.673, &
159  1.389, 2.066, 2.371, 2.719, 3.161, 2.970, 3.623, 4.747, 7.064, 9.419,14.120,27.565, &
160  0.309, 0.497, 0.587, 0.693, 0.832, 0.777, 0.986, 1.364, 2.207, 3.136, 5.174,12.065, &
161  0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885, &
162        ! NO3 coarse (CS)
163  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
164  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
165  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
166  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
167  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
168  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
169        ! NO3 accumulation (AS)
170  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
171  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
172  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
173  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
174  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
175  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
176
177  DATA alpha_aeri_6bands/  &
178       ! dust insoluble CI
179  0.751, 0.751, 0.769, 0.772, 0.672, 0.437, &
180       ! BC insoluble AI
181  6.497, 6.497, 5.900, 4.284, 2.163, 0.966, &
182       ! POM insoluble AI
183  6.443, 6.443, 4.381, 1.846, 0.377, 0.052, &
184       ! NO3 insoluble CI
185  0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
186
187  DATA cg_aers_6bands/ &
188       ! bc soluble AS
189  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
190  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
191  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
192  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
193  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
194  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, &
195       ! pom soluble AS
196  0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.710, 0.726, 0.736, 0.737, 0.738, 0.745, &
197  0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.710, 0.726, 0.736, 0.737, 0.738, 0.745, &
198  0.658, 0.658, 0.658, 0.658, 0.658, 0.667, 0.674, 0.685, 0.692, 0.692, 0.693, 0.698, &
199  0.564, 0.564, 0.564, 0.564, 0.564, 0.566, 0.568, 0.571, 0.573, 0.573, 0.573, 0.574, &
200  0.363, 0.363, 0.363, 0.363, 0.363, 0.360, 0.357, 0.352, 0.350, 0.349, 0.349, 0.347, &
201  0.142, 0.142, 0.142, 0.142, 0.142, 0.139, 0.137, 0.133, 0.131, 0.131, 0.130, 0.129, &
202       ! sulfate coarse CS
203  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
204  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
205  0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, &
206  0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, &
207  0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
208  0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
209       ! sulfate AS
210  0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
211  0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
212  0.653, 0.653, 0.653, 0.662, 0.670, 0.676, 0.683, 0.692, 0.701, 0.706, 0.709, 0.716, &
213  0.563, 0.563, 0.563, 0.565, 0.567, 0.569, 0.570, 0.573, 0.575, 0.576, 0.577, 0.579, &
214  0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, &
215  0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, &
216        ! seasalt seasalt Super Coarse Soluble (SS)
217  0.797, 0.800, 0.801, 0.802, 0.804, 0.822, 0.825, 0.828, 0.832, 0.835, 0.838, 0.843, &
218  0.788, 0.792, 0.794, 0.795, 0.796, 0.815, 0.818, 0.822, 0.827, 0.829, 0.833, 0.838, &
219  0.773, 0.778, 0.780, 0.782, 0.783, 0.802, 0.806, 0.811, 0.817, 0.820, 0.825, 0.832, &
220  0.746, 0.753, 0.755, 0.759, 0.760, 0.781, 0.787, 0.792, 0.800, 0.805, 0.811, 0.820, &
221  0.706, 0.714, 0.716, 0.720, 0.722, 0.749, 0.753, 0.761, 0.769, 0.774, 0.783, 0.797, &
222  0.681, 0.682, 0.682, 0.683, 0.684, 0.723, 0.727, 0.732, 0.738, 0.741, 0.748, 0.757, &
223        ! seasalt seasalt Coarse Soluble (CS)     
224  0.756, 0.761, 0.764, 0.766, 0.769, 0.790, 0.793, 0.799, 0.805, 0.810, 0.815, 0.823, &
225  0.736, 0.743, 0.747, 0.749, 0.751, 0.773, 0.778, 0.784, 0.793, 0.797, 0.804, 0.815, &
226  0.712, 0.719, 0.721, 0.725, 0.726, 0.752, 0.758, 0.764, 0.773, 0.779, 0.786, 0.800, &
227  0.690, 0.694, 0.695, 0.698, 0.699, 0.731, 0.738, 0.742, 0.751, 0.756, 0.764, 0.776, &
228  0.682, 0.683, 0.683, 0.683, 0.684, 0.725, 0.729, 0.733, 0.737, 0.740, 0.744, 0.752, &
229  0.669, 0.673, 0.674, 0.675, 0.676, 0.718, 0.724, 0.730, 0.736, 0.739, 0.742, 0.746, &
230        ! seasalt seasalt Accumulation Soluble (AS)
231  0.694, 0.692, 0.692, 0.691, 0.689, 0.737, 0.740, 0.742, 0.742, 0.740, 0.737, 0.731, &
232  0.685, 0.690, 0.691, 0.692, 0.692, 0.735, 0.741, 0.746, 0.750, 0.751, 0.750, 0.744, &
233  0.650, 0.662, 0.666, 0.670, 0.673, 0.710, 0.719, 0.729, 0.741, 0.747, 0.754, 0.757, &
234  0.561, 0.585, 0.593, 0.601, 0.609, 0.637, 0.651, 0.669, 0.691, 0.705, 0.723, 0.745, &
235  0.392, 0.427, 0.439, 0.451, 0.464, 0.480, 0.500, 0.526, 0.563, 0.588, 0.621, 0.671, &
236  0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439, &
237        ! NO3 coarse (CS)
238  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
239  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
240  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
241  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
242  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
243  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
244        ! NO3 accumulation (AS)
245  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
246  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
247  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
248  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
249  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
250  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
251
252  DATA cg_aeri_6bands/ &
253       ! dust insoluble CI
254  0.718, 0.718, 0.699, 0.661, 0.676, 0.670, &
255       ! BC insoluble AI
256  0.721, 0.721, 0.643, 0.513, 0.321, 0.153, &
257       ! POM insoluble AI
258  0.687, 0.687, 0.658, 0.564, 0.363, 0.142, &
259       ! NO3 insoluble CI
260  0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
261
262  DATA piz_aers_6bands/&
263       ! BC soluble AS
264  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
265  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
266  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
267  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
268  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
269  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, &
270       ! pom soluble AS
271  0.973, 0.973, 0.973, 0.973, 0.973, 0.977, 0.980, 0.984, 0.987, 0.988, 0.988, 0.990, &
272  0.973, 0.973, 0.973, 0.973, 0.973, 0.977, 0.980, 0.984, 0.987, 0.988, 0.988, 0.990, &
273  0.972, 0.972, 0.972, 0.972, 0.972, 0.976, 0.979, 0.984, 0.987, 0.987, 0.988, 0.990, &
274  0.940, 0.940, 0.940, 0.940, 0.940, 0.948, 0.955, 0.965, 0.972, 0.973, 0.973, 0.978, &
275  0.816, 0.816, 0.816, 0.816, 0.816, 0.839, 0.859, 0.888, 0.908, 0.910, 0.911, 0.925, &
276  0.663, 0.663, 0.663, 0.663, 0.663, 0.607, 0.562, 0.492, 0.446, 0.441, 0.437, 0.404, &
277       ! sulfate coarse CS
278  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
279  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
280  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
281  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
282  0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
283  0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
284       ! sulfate AS
285  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
286  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
287  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
288  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
289  0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, &
290  0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, &
291        ! seasalt seasalt Super Coarse Soluble (SS)
292  0.983, 0.982, 0.982, 0.982, 0.982, 0.992, 0.993, 0.994, 0.995, 0.996, 0.997, 0.998, &
293  0.984, 0.983, 0.983, 0.983, 0.983, 0.993, 0.994, 0.995, 0.996, 0.996, 0.997, 0.998, &
294  0.986, 0.985, 0.985, 0.985, 0.984, 0.993, 0.994, 0.995, 0.996, 0.997, 0.997, 0.998, &
295  0.989, 0.988, 0.988, 0.988, 0.987, 0.995, 0.996, 0.996, 0.997, 0.997, 0.998, 0.999, &
296  0.994, 0.993, 0.992, 0.992, 0.992, 0.997, 0.997, 0.998, 0.998, 0.998, 0.999, 0.999, &
297  0.997, 0.997, 0.997, 0.996, 0.996, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, &
298        ! seasalt seasalt Coarse Soluble (CS)     
299  0.988, 0.987, 0.987, 0.987, 0.986, 0.995, 0.995, 0.996, 0.997, 0.997, 0.998, 0.999, &
300  0.990, 0.989, 0.989, 0.989, 0.989, 0.995, 0.996, 0.997, 0.997, 0.998, 0.998, 0.999, &
301  0.993, 0.992, 0.992, 0.992, 0.991, 0.997, 0.997, 0.997, 0.998, 0.998, 0.998, 0.999, &
302  0.996, 0.995, 0.995, 0.995, 0.994, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 0.999, &
303  0.998, 0.997, 0.997, 0.997, 0.997, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, &
304  0.999, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, &
305        ! seasalt seasalt Accumulation Soluble (AS)
306  0.999, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, &
307  0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
308  0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
309  0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
310  0.997, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, &
311  0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000,  &
312        ! NO3 coarse (CS)
313  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
314  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
315  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
316  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
317  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
318  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
319        ! NO3 accumulation (AS)
320  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
321  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
322  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
323  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
324  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
325  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /
326
327  DATA piz_aeri_6bands/ &
328       ! dust insoluble CI
329  0.954, 0.954, 0.965, 0.981, 0.996, 0.990, &
330       ! BC insoluble AI
331  0.460, 0.460, 0.445, 0.394, 0.267, 0.121, &
332       ! POM insoluble AI
333  0.973, 0.973, 0.972, 0.940, 0.816, 0.663, &
334       ! NO3 insoluble CI
335  1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /
336
337! Added by R. Wang (July 31 2016)
338! properties for BC assuming Maxwell-Garnett rule and internal mixture
339
340  DATA alpha_MG_6bands/ &
341     ! bc content = 0.001
342  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
343  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
344  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
345  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
346  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
347  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, &
348     ! bc content = 0.01
349  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
350  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
351  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
352  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
353  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
354  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, &
355     ! bc content = 0.02
356  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
357  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
358  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
359  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
360  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
361  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, &
362     ! bc content = 0.05
363  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
364  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
365  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
366  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
367  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
368  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, &
369     ! bc content = 0.10
370  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
371  6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, &
372  5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, &
373  4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, &
374  2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, &
375  0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894 /
376
377  DATA cg_MG_6bands/ &
378     ! bc content = 0.001
379  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
380  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
381  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
382  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
383  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
384  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, &
385     ! bc content = 0.01
386  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
387  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
388  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
389  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
390  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
391  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, &
392     ! bc content = 0.02
393  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
394  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
395  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
396  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
397  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
398  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, &
399     ! bc content = 0.05
400  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
401  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
402  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
403  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
404  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
405  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, &
406     ! bc content = 0.10
407  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
408  0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, &
409  0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, &
410  0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, &
411  0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, &
412  0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123 /
413
414  DATA piz_MG_6bands/ &
415     ! bc content = 0.001
416  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
417  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
418  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
419  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
420  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
421  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, &
422     ! bc content = 0.01
423  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
424  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
425  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
426  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
427  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
428  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, &
429     ! bc content = 0.02
430  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
431  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
432  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
433  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
434  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
435  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, &
436     ! bc content = 0.05
437  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
438  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
439  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
440  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
441  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
442  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, &
443     ! bc content = 0.10
444  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
445  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
446  0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, &
447  0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, &
448  0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, &
449  0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207 /
450
451!----BEGINNING OF CALCULATIONS
452
453  spsol = 0
454  spinsol = 0
455
456  IF (NSW.NE.nbands_sw_rrtm) THEN
457     print *,'Erreur NSW doit etre egal a 6 pour cette routine'
458     stop
459  ENDIF
460
461  zdh(:,:)=pdel(:,:)/(RG*zrho(:,:))      ! m
462
463  IF (flag_aerosol .EQ. 1) THEN
464     nb_aer = 2
465     ALLOCATE (aerosol_name(nb_aer))
466     aerosol_name(1) = id_ASSO4M_phy
467     aerosol_name(2) = id_CSSO4M_phy
468  ELSEIF (flag_aerosol .EQ. 2) THEN
469     nb_aer = 2
470     ALLOCATE (aerosol_name(nb_aer))
471     aerosol_name(1) = id_ASBCM_phy
472     aerosol_name(2) = id_AIBCM_phy
473  ELSEIF (flag_aerosol .EQ. 3) THEN
474     nb_aer = 2
475     ALLOCATE (aerosol_name(nb_aer))
476     aerosol_name(1) = id_ASPOMM_phy
477     aerosol_name(2) = id_AIPOMM_phy
478  ELSEIF (flag_aerosol .EQ. 4) THEN
479     nb_aer = 3
480     ALLOCATE (aerosol_name(nb_aer))
481     aerosol_name(1) = id_CSSSM_phy
482     aerosol_name(2) = id_SSSSM_phy
483     aerosol_name(3) = id_ASSSM_phy
484  ELSEIF (flag_aerosol .EQ. 5) THEN
485     nb_aer = 1
486     ALLOCATE (aerosol_name(nb_aer))
487     aerosol_name(1) = id_CIDUSTM_phy
488  ELSEIF (flag_aerosol .EQ. 6) THEN
489     nb_aer = 13
490     ALLOCATE (aerosol_name(nb_aer))
491     aerosol_name(1) = id_ASSO4M_phy     
492     aerosol_name(2) = id_ASBCM_phy
493     aerosol_name(3) = id_AIBCM_phy
494     aerosol_name(4) = id_ASPOMM_phy
495     aerosol_name(5) = id_AIPOMM_phy
496     aerosol_name(6) = id_CSSSM_phy
497     aerosol_name(7) = id_SSSSM_phy
498     aerosol_name(8) = id_ASSSM_phy
499     aerosol_name(9) = id_CIDUSTM_phy
500     aerosol_name(10)= id_CSSO4M_phy
501     aerosol_name(11)= id_CSNO3M_phy
502     aerosol_name(12)= id_ASNO3M_phy
503     aerosol_name(13)= id_CINO3M_phy
504  ENDIF
505
506  !
507  ! loop over modes, use of precalculated nmd and corresponding sigma
508  !    loop over wavelengths
509  !    for each mass species in mode
510  !      interpolate from Sext to retrieve Sext_at_gridpoint_per_species
511  !      compute optical_thickness_at_gridpoint_per_species
512
513  DO n=1,nbre_RH-1
514    fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n))
515  ENDDO
516   
517  DO k=1, klev
518    DO i=1, klon
519      rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
520      RH_num(i,k) = INT(rh(i,k)/10. + 1.)
521      IF (rh(i,k).GT.85.) RH_num(i,k)=10
522      IF (rh(i,k).GT.90.) RH_num(i,k)=11
523      delta(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k))
524    ENDDO
525  ENDDO
526
527  tau_ae(:,:,:,:)=0.
528  piz_ae(:,:,:,:)=0.
529  cg_ae(:,:,:,:)=0.
530
531  tau_ae_pi(:,:,:,:)=0.
532  piz_ae_pi(:,:,:,:)=0.
533  cg_ae_pi(:,:,:,:)=0.
534   
535  DO m=1,nb_aer   ! tau is only computed for each mass
536     fac=1.0
537     IF (aerosol_name(m).EQ.id_ASBCM_phy) THEN
538        soluble=.TRUE.
539        spsol=1
540     ELSEIF (aerosol_name(m).EQ.id_ASPOMM_phy) THEN
541        soluble=.TRUE.
542        spsol=2
543     ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
544        soluble=.TRUE.
545        spsol=3
546        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
547     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
548        soluble=.TRUE.
549        spsol=4
550        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
551     ELSEIF (aerosol_name(m).EQ.id_SSSSM_phy) THEN
552         soluble=.TRUE.
553         spsol=5
554     ELSEIF (aerosol_name(m).EQ.id_CSSSM_phy) THEN
555         soluble=.TRUE.
556         spsol=6
557     ELSEIF (aerosol_name(m).EQ.id_ASSSM_phy) THEN
558         soluble=.TRUE.
559         spsol=7
560     ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
561         soluble=.TRUE.
562         spsol=8
563         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
564     ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
565         soluble=.TRUE.
566         spsol=9
567         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
568     ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
569         soluble=.FALSE.
570         spinsol=1
571     ELSEIF  (aerosol_name(m).EQ.id_AIBCM_phy) THEN
572         soluble=.FALSE.
573         spinsol=2
574     ELSEIF (aerosol_name(m).EQ.id_AIPOMM_phy) THEN
575         soluble=.FALSE.
576         spinsol=3
577     ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
578         soluble=.FALSE.
579         spinsol=4
580         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
581     ELSE
582         CYCLE
583     ENDIF
584
585    !--shortname for aerosol index
586    aerindex=aerosol_name(m)
587
588    IF (soluble) THEN
589
590       !--here we treat the special case of soluble BC internal mixture with Maxwell-Garnett rule
591       IF (spsol.EQ.1 .AND. flag_bc_internal_mixture) THEN
592
593         DO k=1, klev
594           DO i=1, klon
595
596             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
597             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
598
599             ! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols
600             BC_massfra = m_allaer(i,k,id_ASBCM_phy)/(m_allaer(i,k,id_ASBCM_phy)+m_allaer(i,k,id_ASSO4M_phy))
601
602             IF (BC_massfra.GE.0.10) THEN
603               classbc = 5
604             ELSEIF  (BC_massfra.GE.0.05) THEN
605               classbc = 4
606             ELSEIF  (BC_massfra.GE.0.02) THEN
607               classbc = 3
608             ELSEIF  (BC_massfra.GE.0.01) THEN
609               classbc = 2
610             ELSE
611               classbc = 1
612             ENDIF
613
614             ! Calculate the dry BC/(BC+SUL) mass ratio for natural aerosols
615             BC_massfra_pi = m_allaer_pi(i,k,id_ASBCM_phy)/(m_allaer_pi(i,k,id_ASBCM_phy)+m_allaer_pi(i,k,id_ASSO4M_phy))
616
617             IF (BC_massfra_pi.GE.0.10) THEN
618               classbc_pi = 5
619             ELSEIF  (BC_massfra_pi.GE.0.05) THEN
620               classbc_pi = 4
621             ELSEIF  (BC_massfra_pi.GE.0.02) THEN
622               classbc_pi = 3
623             ELSEIF  (BC_massfra_pi.GE.0.01) THEN
624               classbc_pi = 2
625             ELSE
626               classbc_pi = 1
627             ENDIF
628
629             DO inu=1,NSW
630
631               !--all aerosols
632               tau_ae2b_int= alpha_MG_6bands(RH_num(i,k),inu,classbc)+                 &
633                             delta(i,k)* (alpha_MG_6bands(RH_num(i,k)+1,inu,classbc) - &
634                             alpha_MG_6bands(RH_num(i,k),inu,classbc))
635                   
636               piz_ae2b_int = piz_MG_6bands(RH_num(i,k),inu,classbc) +                 &
637                              delta(i,k)* (piz_MG_6bands(RH_num(i,k)+1,inu,classbc) -  &
638                              piz_MG_6bands(RH_num(i,k),inu,classbc))
639                   
640               cg_ae2b_int = cg_MG_6bands(RH_num(i,k),inu,classbc) +                   &
641                             delta(i,k)* (cg_MG_6bands(RH_num(i,k)+1,inu,classbc) -    &
642                             cg_MG_6bands(RH_num(i,k),inu,classbc))
643
644               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
645               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
646               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
647
648               !--natural aerosols
649               tau_ae2b_int= alpha_MG_6bands(RH_num(i,k),inu,classbc_pi)+                 &
650                             delta(i,k)* (alpha_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) - &
651                             alpha_MG_6bands(RH_num(i,k),inu,classbc_pi))
652                   
653               piz_ae2b_int = piz_MG_6bands(RH_num(i,k),inu,classbc_pi) +                 &
654                              delta(i,k)* (piz_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) -  &
655                              piz_MG_6bands(RH_num(i,k),inu,classbc_pi))
656                   
657               cg_ae2b_int = cg_MG_6bands(RH_num(i,k),inu,classbc_pi) +                   &
658                             delta(i,k)* (cg_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) -    &
659                             cg_MG_6bands(RH_num(i,k),inu,classbc_pi))
660
661               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
662               piz_ae_pi(i,k,aerindex,inu) = piz_ae2b_int
663               cg_ae_pi(i,k,aerindex,inu)  = cg_ae2b_int
664                       
665             ENDDO
666           ENDDO
667         ENDDO
668
669       !--else treat all other cases of soluble aerosols
670       ELSE
671
672         DO k=1, klev
673           DO i=1, klon
674             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
675             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
676
677             DO inu=1,NSW
678
679               tau_ae2b_int= alpha_aers_6bands(RH_num(i,k),inu,spsol)+ &
680                             delta(i,k)* (alpha_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
681                             alpha_aers_6bands(RH_num(i,k),inu,spsol))
682                   
683               piz_ae2b_int = piz_aers_6bands(RH_num(i,k),inu,spsol) + &
684                            delta(i,k)* (piz_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
685                              piz_aers_6bands(RH_num(i,k),inu,spsol))
686                   
687               cg_ae2b_int = cg_aers_6bands(RH_num(i,k),inu,spsol) + &
688                             delta(i,k)* (cg_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
689                             cg_aers_6bands(RH_num(i,k),inu,spsol))
690
691               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
692               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
693               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
694               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
695                       
696             ENDDO
697           ENDDO
698         ENDDO
699
700         !--external mixture case for soluble BC
701         IF (spsol.EQ.1) THEN
702           piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:)
703           cg_ae_pi(:,:,aerindex,:)  = cg_ae(:,:,aerindex,:)
704         ENDIF
705
706       ENDIF
707       
708     ELSE    ! For all aerosol insoluble components
709
710       DO k=1, klev
711         DO i=1, klon
712           tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
713           tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
714
715           DO inu=1,NSW
716             tau_ae2b_int = alpha_aeri_6bands(inu,spinsol)
717             piz_ae2b_int = piz_aeri_6bands(inu,spinsol)
718             cg_ae2b_int = cg_aeri_6bands(inu,spinsol)
719
720             tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
721             tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi*tau_ae2b_int
722             piz_ae(i,k,aerindex,inu) = piz_ae2b_int
723             cg_ae(i,k,aerindex,inu)= cg_ae2b_int
724           ENDDO
725         ENDDO
726       ENDDO
727
728     ENDIF ! soluble / insoluble
729
730  ENDDO  ! nb_aer 
731
732  DO inu=1, NSW
733     DO k=1, klev
734       DO i=1, klon
735!--all (natural + anthropogenic) aerosol
736         tau_allaer(i,k,2,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)+ &
737                               tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+   &
738                               tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ &
739                               tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+   &
740                               tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu)
741         tau_allaer(i,k,2,inu)=MAX(tau_allaer(i,k,2,inu),tau_min)
742
743         piz_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+   &
744                                tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+   &
745                                tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+     &
746                                tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+     &
747                                tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+   &
748                                tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+   &
749                                tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+     &
750                                tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+     &
751                                tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+     &
752                                tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
753                                /tau_allaer(i,k,2,inu)
754         piz_allaer(i,k,2,inu)=MIN(MAX(piz_allaer(i,k,2,inu),0.01),1.0)
755         IF (tau_allaer(i,k,2,inu).LE.tau_min) piz_allaer(i,k,2,inu)=1.0
756
757         cg_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ &
758                               tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)*cg_ae(i,k,id_CSSO4M_phy,inu)+ &
759                               tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)*cg_ae(i,k,id_ASBCM_phy,inu)+    &
760                               tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)*cg_ae(i,k,id_AIBCM_phy,inu)+    &
761                               tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)*cg_ae(i,k,id_ASPOMM_phy,inu)+ &
762                               tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)*cg_ae(i,k,id_AIPOMM_phy,inu)+ &
763                               tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)*cg_ae(i,k,id_ASSSM_phy,inu)+    &
764                               tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)*cg_ae(i,k,id_CSSSM_phy,inu)+    &
765                               tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)*cg_ae(i,k,id_SSSSM_phy,inu)+    &
766                               tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)*cg_ae(i,k,id_CIDUSTM_phy,inu))/ &
767                               (tau_allaer(i,k,2,inu)*piz_allaer(i,k,2,inu))
768         cg_allaer(i,k,2,inu)=MIN(MAX(cg_allaer(i,k,2,inu),0.0),1.0)
769
770!--natural aerosol
771!--ASBCM aerosols take _pi value because of internal mixture option
772         tau_allaer(i,k,1,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ &
773                               tau_ae_pi(i,k,id_ASBCM_phy,inu)+tau_ae_pi(i,k,id_AIBCM_phy,inu)+   &
774                               tau_ae_pi(i,k,id_ASPOMM_phy,inu)+tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ &
775                               tau_ae_pi(i,k,id_ASSSM_phy,inu)+tau_ae_pi(i,k,id_CSSSM_phy,inu)+   &
776                               tau_ae_pi(i,k,id_SSSSM_phy,inu)+ tau_ae_pi(i,k,id_CIDUSTM_phy,inu)
777         tau_allaer(i,k,1,inu)=MAX(tau_allaer(i,k,1,inu),tau_min)
778
779         piz_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+   &
780                                tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+   &
781                                tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae_pi(i,k,id_ASBCM_phy,inu)+  &
782                                tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+     &
783                                tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+   &
784                                tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+   &
785                                tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+     &
786                                tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+     &
787                                tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+     &
788                                tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
789                                /tau_allaer(i,k,1,inu)
790         piz_allaer(i,k,1,inu)=MIN(MAX(piz_allaer(i,k,1,inu),0.01),1.0)
791         IF (tau_allaer(i,k,1,inu).LE.tau_min) piz_allaer(i,k,1,inu)=1.0
792
793         cg_allaer(i,k,1,inu)=(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)+    &
794                               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)+    &
795                               tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae_pi(i,k,id_ASBCM_phy,inu)*cg_ae_pi(i,k,id_ASBCM_phy,inu)+ &
796                               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)+       &
797                               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)+    &
798                               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)+    &
799                               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)+       &
800                               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)+       &
801                               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)+       &
802                               tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)*cg_ae(i,k,id_CIDUSTM_phy,inu))/ &
803                               (tau_allaer(i,k,1,inu)*piz_allaer(i,k,1,inu))
804         cg_allaer(i,k,1,inu)=MIN(MAX(cg_allaer(i,k,1,inu),0.0),1.0)
805
806        ENDDO
807      ENDDO
808    ENDDO
809   
810!--waveband 2 and all aerosol (third index = 2)
811  inu=2
812  DO i=1, klon
813     absvisaer(i)=SUM((1-piz_allaer(i,:,2,inu))*tau_allaer(i,:,2,inu))
814  ENDDO
815
816  DEALLOCATE(aerosol_name)
817
818END SUBROUTINE AEROPT_6BANDS_RRTM
Note: See TracBrowser for help on using the repository browser.