Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

aeropt_6bands_rrtm.F90

Last change on this file was 2842, checked in by oboucher, 8 years ago
Removing absvisaer from calculations and outputs Adding abs550aer in rrtm/aeropt_5wv_rrtm and outputs Absorption aerosol optical depth at 550 nm
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: 48.2 KB

Rev	Line
[2004]	1	!
[2005]	2	! $Id: aeropt_6bands_rrtm.F90 2842 2017-04-01 16:09:13Z abarral $
[2004]	3	!
	4	SUBROUTINE AEROPT_6BANDS_RRTM ( &
[2737]	5	pdel, m_allaer, RHcl, &
	6	tau_allaer, piz_allaer, &
	7	cg_allaer, m_allaer_pi, &
	8	flag_aerosol, flag_bc_internal_mixture, zrho )
[2004]	9
	10	USE dimphy
	11	USE aero_mod
[2748]	12	USE YOMCST, ONLY: RG
[2004]	13
	14	! Yves Balkanski le 12 avril 2006
	15	! Celine Deandreis
	16	! Anne Cozic Avril 2009
	17	! a partir d'une sous-routine de Johannes Quaas pour les sulfates
	18	! Olivier Boucher février 2014 pour passage à RRTM
	19	! a partir des propriétés optiques fournies par Yves Balkanski
	20	!
	21	IMPLICIT NONE
[2748]	22	!
[2749]	23	INCLUDE "clesphys.h"
[2004]	24	!
	25	! Input arguments:
	26	!
[2748]	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
[2004]	34	!
	35	! Output arguments:
[2634]	36	! 2= total aerosols
	37	! 1= natural aerosols
[2004]	38	!
[2748]	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
[2004]	42	!
	43	! Local
	44	!
[2739]	45	LOGICAL :: soluble
[2004]	46	INTEGER :: i, k,n, inu, m
	47	INTEGER :: spsol, spinsol
	48	INTEGER :: RH_num(klon,klev)
	49
	50	INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
	51
[2738]	52	INTEGER, PARAMETER :: naero_soluble=9 ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
	53	! 5- seasalt super coarse 6- seasalt
	54	! 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
[2004]	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
[2634]	67	REAL :: zdh(klon,klev)
[2739]	68	INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name
[2004]	69	INTEGER :: nb_aer
	70
[2739]	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
[2644]	74
[2739]	75	REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae_pi
[2644]	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
[2004]	78	!
	79	! Proprietes optiques
	80	!
[2738]	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
[2644]	87	!
	88	! BC internal mixture
	89	!
[2739]	90	INTEGER, PARAMETER :: nbclassbc = 6 ! Added by Rong Wang/OB for the 5 fractions
[2738]	91	! of BC in the soluble mode:
[2739]	92	! bc_content/0.001, 0.01, 0.02, 0.05, 0.1, 0.2/
[2644]	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	!
[2738]	101	INTEGER :: aerindex, classbc, classbc_pi
[2644]	102	REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
[2748]	103	CHARACTER*20 :: modname
[2644]	104	!
[2739]	105	REAL, PARAMETER :: tau_min = 1.e-7
[2004]	106
	107	!***************************************************************************
	108	!--the order of the soluble species has to follow the spsol index below
	109	!--the order of the insoluble species has to follow the spinsol index below
	110
	111	DATA alpha_aers_6bands/ &
[2739]	112	! BC Accumulation Soluble (AS)
	113	5.276, 5.276, 5.276, 5.276, 5.276, 5.797, 6.353, 7.568, 8.779, 8.920, 9.062,10.403, &
	114	5.371, 5.371, 5.371, 5.371, 5.371, 5.870, 6.401, 7.559, 8.720, 8.857, 8.994,10.298, &
	115	5.126, 5.126, 5.126, 5.126, 5.126, 5.572, 6.043, 7.073, 8.113, 8.235, 8.359, 9.544, &
	116	4.225, 4.225, 4.225, 4.225, 4.225, 4.549, 4.887, 5.624, 6.370, 6.459, 6.548, 7.412, &
	117	2.667, 2.667, 2.667, 2.667, 2.667, 2.840, 3.015, 3.382, 3.745, 3.788, 3.832, 4.251, &
	118	1.113, 1.113, 1.113, 1.113, 1.113, 1.210, 1.307, 1.506, 1.699, 1.721, 1.744, 1.961, &
	119	! POM Accumulation Soluble (AS)
	120	5.962, 5.962, 5.962, 5.962, 5.962, 6.544, 7.153, 8.460, 9.732, 9.877,10.026,11.398, &
	121	5.609, 5.609, 5.609, 5.609, 5.609, 6.163, 6.752, 8.041, 9.324, 9.473, 9.624,11.044, &
	122	4.527, 4.527, 4.527, 4.527, 4.527, 4.982, 5.476, 6.586, 7.726, 7.861, 7.998, 9.306, &
	123	2.697, 2.697, 2.697, 2.697, 2.697, 2.980, 3.293, 4.017, 4.790, 4.884, 4.979, 5.902, &
	124	0.918, 0.918, 0.918, 0.918, 0.918, 1.025, 1.146, 1.435, 1.754, 1.794, 1.834, 2.232, &
	125	0.105, 0.105, 0.105, 0.105, 0.105, 0.154, 0.211, 0.344, 0.491, 0.509, 0.528, 0.710, &
	126	! Sulfate Coarse Soluble (CS)
	127	0.670, 0.670, 0.670, 0.670, 0.905, 0.981, 1.081, 1.217, 1.434, 1.615, 1.940, 2.773, &
	128	0.689, 0.689, 0.689, 0.689, 0.930, 1.009, 1.110, 1.247, 1.469, 1.655, 1.980, 2.825, &
	129	0.726, 0.726, 0.726, 0.726, 0.978, 1.060, 1.164, 1.306, 1.533, 1.723, 2.056, 2.915, &
	130	0.785, 0.785, 0.785, 0.785, 1.056, 1.144, 1.254, 1.406, 1.650, 1.852, 2.203, 3.103, &
	131	0.795, 0.795, 0.795, 0.795, 1.064, 1.160, 1.282, 1.449, 1.720, 1.945, 2.340, 3.348, &
	132	0.502, 0.502, 0.502, 0.502, 0.656, 0.728, 0.820, 0.950, 1.170, 1.360, 1.706, 2.666, &
	133	!-- Sulfate Accumulation (BC content=0)
	134	5.246, 5.246, 5.246, 5.246, 7.041, 7.658, 8.437, 9.503,11.204,12.608,15.030,21.080, &
	135	4.924, 4.924, 4.924, 4.924, 6.619, 7.258, 8.067, 9.186,10.998,12.507,15.142,21.822, &
	136	3.925, 3.925, 3.925, 3.925, 5.304, 5.888, 6.635, 7.686, 9.433,10.921,13.587,20.659, &
	137	2.248, 2.248, 2.248, 2.248, 3.079, 3.477, 3.993, 4.738, 6.022, 7.153, 9.261,15.301, &
	138	0.720, 0.720, 0.720, 0.720, 1.014, 1.169, 1.374, 1.678, 2.225, 2.726, 3.703, 6.779, &
	139	0.069, 0.069, 0.069, 0.069, 0.102, 0.121, 0.146, 0.185, 0.258, 0.328, 0.472, 0.974, &
	140	! Seasalt Super Coarse Soluble (SS)
	141	0.189, 0.232, 0.249, 0.270, 0.293, 0.322, 0.361, 0.425, 0.536, 0.634, 0.813, 1.250, &
	142	0.192, 0.235, 0.252, 0.274, 0.296, 0.325, 0.365, 0.429, 0.542, 0.639, 0.821, 1.259, &
	143	0.196, 0.240, 0.257, 0.279, 0.302, 0.331, 0.371, 0.437, 0.550, 0.649, 0.832, 1.274, &
	144	0.204, 0.249, 0.267, 0.289, 0.313, 0.342, 0.384, 0.450, 0.565, 0.666, 0.851, 1.300, &
	145	0.222, 0.270, 0.289, 0.313, 0.338, 0.369, 0.413, 0.482, 0.602, 0.705, 0.896, 1.355, &
	146	0.244, 0.293, 0.313, 0.338, 0.364, 0.397, 0.443, 0.517, 0.646, 0.757, 0.960, 1.443, &
	147	! Seasalt Coarse Soluble (CS)
	148	0.547, 0.657, 0.705, 0.754, 0.817, 0.896, 1.008, 1.169, 1.456, 1.724, 2.199, 3.358, &
	149	0.566, 0.679, 0.727, 0.776, 0.840, 0.920, 1.032, 1.196, 1.492, 1.760, 2.238, 3.416, &
	150	0.596, 0.714, 0.764, 0.816, 0.882, 0.965, 1.081, 1.250, 1.552, 1.828, 2.310, 3.509, &
	151	0.644, 0.771, 0.825, 0.880, 0.951, 1.040, 1.164, 1.345, 1.666, 1.957, 2.462, 3.700, &
	152	0.640, 0.772, 0.829, 0.887, 0.965, 1.061, 1.198, 1.398, 1.758, 2.085, 2.658, 4.031, &
	153	0.452, 0.562, 0.609, 0.659, 0.728, 0.813, 0.938, 1.125, 1.471, 1.797, 2.384, 3.855, &
	154	! Seasalt Accumulation Soluble (AS)
	155	4.270, 5.032, 5.395, 5.772, 6.158, 6.768, 7.510, 8.634,10.695,12.512,15.740,23.122, &
	156	3.952, 4.673, 5.030, 5.403, 5.781, 6.399, 7.156, 8.328,10.533,12.524,16.091,24.329, &
	157	3.090, 3.676, 3.979, 4.301, 4.622, 5.172, 5.853, 6.941, 9.077,11.083,14.830,23.951, &
	158	1.710, 2.069, 2.262, 2.471, 2.678, 3.048, 3.513, 4.285, 5.882, 7.467,10.618,19.102, &
	159	0.529, 0.663, 0.737, 0.819, 0.901, 1.050, 1.244, 1.575, 2.301, 3.067, 4.704, 9.712, &
	160	0.073, 0.142, 0.177, 0.215, 0.259, 0.328, 0.422, 0.579, 0.913, 1.263, 2.011, 4.358, &
	161	! Nitrate Coarse Soluble (CS)
	162	0.688, 0.688, 0.688, 0.755, 0.826, 0.900, 0.994, 1.190, 1.501, 1.804, 2.251, 2.374, &
	163	0.712, 0.712, 0.712, 0.781, 0.853, 0.929, 1.023, 1.225, 1.540, 1.846, 2.299, 2.421, &
	164	0.754, 0.754, 0.754, 0.826, 0.902, 0.980, 1.077, 1.286, 1.612, 1.925, 2.388, 2.509, &
	165	0.809, 0.809, 0.809, 0.887, 0.969, 1.053, 1.160, 1.384, 1.732, 2.069, 2.557, 2.685, &
	166	0.786, 0.786, 0.786, 0.867, 0.952, 1.043, 1.157, 1.403, 1.793, 2.171, 2.725, 2.871, &
	167	0.703, 0.703, 0.703, 0.767, 0.836, 0.909, 1.002, 1.210, 1.552, 1.901, 2.434, 2.579, &
	168	! Nitrate Accumulation Soluble (AS)
	169	5.251, 5.251, 5.251, 5.779, 6.336, 6.920, 7.655, 9.230,11.678,14.026,17.391,18.279, &
	170	4.643, 4.643, 4.643, 5.152, 5.696, 6.276, 7.016, 8.640,11.238,13.787,17.511,18.494, &
	171	3.373, 3.373, 3.373, 3.790, 4.245, 4.739, 5.383, 6.842, 9.289,11.804,15.637,16.675, &
	172	1.779, 1.779, 1.779, 2.029, 2.309, 2.619, 3.033, 4.008, 5.746, 7.645,10.730,11.601, &
	173	0.549, 0.549, 0.549, 0.639, 0.742, 0.859, 1.019, 1.412, 2.157, 3.029, 4.552, 5.002, &
	174	0.367, 0.367, 0.367, 0.416, 0.471, 0.532, 0.614, 0.808, 1.166, 1.575, 2.281, 2.489 /
[2004]	175
	176	DATA piz_aers_6bands/&
[2739]	177	! BC Accumulation Soluble (AS)
	178	0.480, 0.480, 0.480, 0.480, 0.480, 0.487, 0.495, 0.513, 0.530, 0.532, 0.534, 0.553, &
	179	0.469, 0.469, 0.469, 0.469, 0.469, 0.474, 0.480, 0.497, 0.514, 0.517, 0.519, 0.538, &
	180	0.456, 0.456, 0.456, 0.456, 0.456, 0.464, 0.473, 0.495, 0.519, 0.522, 0.524, 0.550, &
	181	0.427, 0.427, 0.427, 0.427, 0.427, 0.437, 0.448, 0.476, 0.505, 0.509, 0.512, 0.543, &
	182	0.349, 0.349, 0.349, 0.349, 0.349, 0.358, 0.369, 0.396, 0.425, 0.429, 0.432, 0.465, &
	183	0.196, 0.196, 0.196, 0.196, 0.196, 0.198, 0.202, 0.213, 0.226, 0.227, 0.229, 0.243, &
	184	! POM Accumulation Soluble (AS)
	185	0.956, 0.956, 0.956, 0.956, 0.956, 0.961, 0.964, 0.970, 0.974, 0.975, 0.975, 0.978, &
	186	0.965, 0.965, 0.965, 0.965, 0.965, 0.969, 0.972, 0.977, 0.980, 0.980, 0.981, 0.983, &
	187	0.969, 0.969, 0.969, 0.969, 0.969, 0.972, 0.975, 0.979, 0.983, 0.983, 0.983, 0.986, &
	188	0.953, 0.953, 0.953, 0.953, 0.953, 0.958, 0.962, 0.969, 0.974, 0.975, 0.975, 0.979, &
	189	0.892, 0.892, 0.892, 0.892, 0.892, 0.902, 0.912, 0.928, 0.940, 0.942, 0.943, 0.952, &
	190	0.826, 0.826, 0.826, 0.826, 0.826, 0.658, 0.568, 0.479, 0.444, 0.442, 0.439, 0.426, &
	191	! Sulfate Coarse Soluble (CS)
[2004]	192	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	193	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	194	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	195	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	196	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	197	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
[2739]	198	!-- Sulfate Accumulation (BC content=0)
[2004]	199	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	200	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
[2739]	201	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	202	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	203	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	204	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	205	! Seasalt Super Coarse Soluble (SS)
	206	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	207	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	208	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	209	0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	210	0.983, 0.985, 0.986, 0.987, 0.987, 0.988, 0.988, 0.989, 0.989, 0.989, 0.989, 0.989, &
	211	0.957, 0.792, 0.765, 0.740, 0.721, 0.704, 0.687, 0.670, 0.652, 0.643, 0.631, 0.617, &
	212	! Seasalt Coarse Soluble (CS)
	213	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.001, 1.000, &
	214	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	215	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	216	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	217	0.994, 0.994, 0.995, 0.995, 0.995, 0.995, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
	218	0.976, 0.867, 0.837, 0.814, 0.796, 0.774, 0.754, 0.735, 0.713, 0.702, 0.690, 0.675, &
	219	! Seasalt Accumulation Soluble (AS)
	220	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	221	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	222	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	223	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	224	0.995, 0.996, 0.996, 0.996, 0.997, 0.997, 0.997, 0.998, 0.998, 0.998, 0.998, 0.999, &
	225	0.919, 0.607, 0.553, 0.515, 0.480, 0.455, 0.434, 0.421, 0.421, 0.429, 0.449, 0.492, &
	226	! Nitrate Coarse Soluble (CS)
	227	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	228	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	229	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	230	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	231	0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.998, 0.998, 0.998, 0.998, 0.998, 0.998, &
	232	0.707, 0.707, 0.707, 0.705, 0.702, 0.699, 0.697, 0.692, 0.686, 0.683, 0.679, 0.678, &
	233	! Nitrate Accumulation Soluble (AS)
	234	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	235	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	236	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	237	1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
	238	0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, &
	239	0.362, 0.362, 0.362, 0.362, 0.363, 0.365, 0.367, 0.375, 0.390, 0.406, 0.428, 0.434 /
[2004]	240
[2739]	241	DATA cg_aers_6bands/ &
	242	! BC Accumulation Soluble (AS)
	243	0.820, 0.820, 0.820, 0.820, 0.820, 0.831, 0.840, 0.852, 0.859, 0.860, 0.860, 0.864, &
	244	0.760, 0.760, 0.760, 0.760, 0.760, 0.777, 0.791, 0.812, 0.826, 0.827, 0.828, 0.837, &
	245	0.695, 0.695, 0.695, 0.695, 0.695, 0.716, 0.734, 0.762, 0.780, 0.782, 0.784, 0.797, &
	246	0.603, 0.603, 0.603, 0.603, 0.603, 0.628, 0.649, 0.685, 0.710, 0.712, 0.715, 0.734, &
	247	0.469, 0.469, 0.469, 0.469, 0.469, 0.494, 0.518, 0.559, 0.590, 0.593, 0.597, 0.622, &
	248	0.276, 0.276, 0.276, 0.276, 0.276, 0.292, 0.308, 0.337, 0.361, 0.364, 0.367, 0.389, &
	249	! POM Accumulation Soluble (AS)
	250	0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.711, 0.729, 0.740, 0.741, 0.742, 0.750, &
	251	0.685, 0.685, 0.685, 0.685, 0.685, 0.699, 0.712, 0.731, 0.744, 0.745, 0.747, 0.756, &
	252	0.673, 0.673, 0.673, 0.673, 0.673, 0.689, 0.702, 0.724, 0.739, 0.741, 0.742, 0.754, &
	253	0.637, 0.637, 0.637, 0.637, 0.637, 0.653, 0.667, 0.691, 0.709, 0.711, 0.712, 0.726, &
	254	0.553, 0.553, 0.553, 0.553, 0.553, 0.568, 0.582, 0.607, 0.628, 0.630, 0.632, 0.649, &
	255	0.340, 0.340, 0.340, 0.340, 0.340, 0.350, 0.361, 0.382, 0.400, 0.402, 0.404, 0.422, &
	256	! Sulfate Coarse Soluble (CS)
	257	0.746, 0.746, 0.746, 0.746, 0.779, 0.785, 0.790, 0.796, 0.804, 0.810, 0.816, 0.830, &
	258	0.730, 0.730, 0.730, 0.730, 0.764, 0.770, 0.775, 0.782, 0.791, 0.797, 0.806, 0.819, &
	259	0.707, 0.707, 0.707, 0.707, 0.746, 0.752, 0.758, 0.765, 0.775, 0.781, 0.790, 0.805, &
	260	0.691, 0.691, 0.691, 0.691, 0.736, 0.741, 0.747, 0.754, 0.763, 0.769, 0.778, 0.791, &
	261	0.696, 0.696, 0.696, 0.696, 0.748, 0.753, 0.759, 0.767, 0.776, 0.781, 0.789, 0.800, &
	262	0.711, 0.711, 0.711, 0.711, 0.766, 0.773, 0.782, 0.791, 0.803, 0.811, 0.822, 0.838, &
	263	!-- Sulfate Accumulation (BC content=0)
	264	0.674, 0.674, 0.674, 0.674, 0.724, 0.729, 0.735, 0.741, 0.749, 0.753, 0.760, 0.769, &
	265	0.677, 0.677, 0.677, 0.677, 0.729, 0.734, 0.740, 0.747, 0.756, 0.761, 0.768, 0.777, &
	266	0.671, 0.671, 0.671, 0.671, 0.723, 0.730, 0.737, 0.746, 0.757, 0.764, 0.773, 0.786, &
	267	0.638, 0.638, 0.638, 0.638, 0.691, 0.700, 0.710, 0.721, 0.736, 0.746, 0.760, 0.781, &
	268	0.553, 0.553, 0.553, 0.553, 0.604, 0.616, 0.629, 0.644, 0.665, 0.680, 0.700, 0.736, &
	269	0.344, 0.344, 0.344, 0.344, 0.388, 0.401, 0.416, 0.435, 0.462, 0.482, 0.512, 0.571, &
	270	! Seasalt Super Coarse Soluble (SS)
	271	0.795, 0.810, 0.814, 0.819, 0.824, 0.829, 0.832, 0.840, 0.846, 0.849, 0.854, 0.860, &
	272	0.788, 0.804, 0.808, 0.813, 0.818, 0.823, 0.828, 0.834, 0.842, 0.847, 0.851, 0.858, &
	273	0.775, 0.792, 0.797, 0.803, 0.808, 0.813, 0.820, 0.826, 0.835, 0.840, 0.846, 0.854, &
	274	0.754, 0.773, 0.779, 0.785, 0.791, 0.796, 0.803, 0.811, 0.821, 0.827, 0.834, 0.844, &
	275	0.735, 0.755, 0.761, 0.768, 0.774, 0.779, 0.786, 0.794, 0.804, 0.810, 0.817, 0.829, &
	276	0.712, 0.782, 0.795, 0.808, 0.819, 0.829, 0.840, 0.852, 0.865, 0.873, 0.882, 0.895, &
	277	! Seasalt Coarse Soluble (CS)
	278	0.754, 0.770, 0.776, 0.781, 0.784, 0.791, 0.797, 0.805, 0.815, 0.822, 0.828, 0.840, &
	279	0.736, 0.753, 0.759, 0.765, 0.771, 0.778, 0.785, 0.793, 0.804, 0.811, 0.820, 0.831, &
	280	0.716, 0.735, 0.742, 0.748, 0.754, 0.762, 0.769, 0.778, 0.789, 0.796, 0.807, 0.819, &
	281	0.704, 0.725, 0.733, 0.739, 0.745, 0.752, 0.759, 0.768, 0.778, 0.784, 0.792, 0.803, &
	282	0.716, 0.737, 0.744, 0.751, 0.756, 0.763, 0.770, 0.777, 0.786, 0.790, 0.795, 0.800, &
	283	0.688, 0.730, 0.741, 0.751, 0.761, 0.771, 0.782, 0.795, 0.810, 0.820, 0.833, 0.849, &
	284	! Seasalt Accumulation Soluble (AS)
	285	0.686, 0.711, 0.719, 0.726, 0.733, 0.740, 0.748, 0.755, 0.763, 0.767, 0.770, 0.774, &
	286	0.690, 0.716, 0.724, 0.731, 0.739, 0.747, 0.756, 0.764, 0.774, 0.778, 0.782, 0.784, &
	287	0.683, 0.710, 0.719, 0.727, 0.735, 0.744, 0.754, 0.765, 0.777, 0.784, 0.791, 0.797, &
	288	0.649, 0.676, 0.686, 0.695, 0.704, 0.714, 0.726, 0.739, 0.757, 0.768, 0.781, 0.797, &
	289	0.559, 0.586, 0.596, 0.606, 0.616, 0.629, 0.643, 0.661, 0.686, 0.703, 0.726, 0.758, &
	290	0.341, 0.358, 0.366, 0.375, 0.382, 0.395, 0.408, 0.428, 0.459, 0.484, 0.520, 0.581, &
	291	! Nitrate Coarse Soluble (CS)
	292	0.765, 0.765, 0.765, 0.773, 0.779, 0.785, 0.791, 0.803, 0.813, 0.819, 0.827, 0.828, &
	293	0.754, 0.754, 0.754, 0.762, 0.769, 0.774, 0.781, 0.790, 0.802, 0.810, 0.817, 0.818, &
	294	0.743, 0.743, 0.743, 0.751, 0.758, 0.763, 0.770, 0.779, 0.789, 0.797, 0.804, 0.806, &
	295	0.737, 0.737, 0.737, 0.745, 0.752, 0.758, 0.764, 0.773, 0.781, 0.786, 0.792, 0.793, &
	296	0.739, 0.739, 0.739, 0.749, 0.757, 0.763, 0.770, 0.779, 0.788, 0.793, 0.797, 0.798, &
	297	0.651, 0.651, 0.651, 0.672, 0.690, 0.707, 0.724, 0.751, 0.779, 0.796, 0.813, 0.817, &
	298	! Nitrate Accumulation Soluble (AS)
	299	0.713, 0.713, 0.713, 0.723, 0.732, 0.739, 0.746, 0.756, 0.765, 0.769, 0.772, 0.772, &
	300	0.724, 0.724, 0.724, 0.734, 0.742, 0.750, 0.757, 0.767, 0.777, 0.782, 0.785, 0.785, &
	301	0.719, 0.719, 0.719, 0.729, 0.738, 0.746, 0.754, 0.767, 0.780, 0.787, 0.793, 0.794, &
	302	0.673, 0.673, 0.673, 0.685, 0.696, 0.706, 0.717, 0.734, 0.753, 0.766, 0.778, 0.781, &
	303	0.565, 0.565, 0.565, 0.580, 0.593, 0.605, 0.619, 0.644, 0.673, 0.693, 0.716, 0.720, &
	304	0.328, 0.328, 0.328, 0.338, 0.348, 0.359, 0.371, 0.396, 0.431, 0.459, 0.494, 0.502 /
	305
	306	DATA alpha_aeri_6bands/ &
	307	! Dust insoluble
	308	0.735, 0.765, 0.817, 0.846, 0.718, 0.414, &
	309	! BC insoluble
	310	5.276, 5.371, 5.126, 4.225, 2.667, 1.113, &
	311	! POM insoluble
	312	5.962, 5.609, 4.527, 2.697, 0.918, 0.105, &
	313	! Nitrate insoluble
	314	0.688, 0.712, 0.754, 0.809, 0.786, 0.703 /
	315
[2004]	316	DATA piz_aeri_6bands/ &
[2739]	317	! Dust insoluble
	318	0.894, 0.902, 0.937, 0.980, 0.993, 0.967, &
	319	! BC insoluble
	320	0.820, 0.760, 0.695, 0.603, 0.469, 0.276, &
	321	! POM insoluble
	322	0.687, 0.685, 0.673, 0.637, 0.553, 0.340, &
	323	! Nitrate insoluble
	324	1.000, 1.000, 1.000, 1.000, 0.999, 0.707 /
[2004]	325
[2739]	326	DATA cg_aeri_6bands/ &
	327	! Dust insoluble
	328	0.756, 0.737, 0.713, 0.701, 0.696, 0.672, &
	329	! BC insoluble
	330	0.480, 0.469, 0.456, 0.427, 0.349, 0.196, &
	331	! POM insoluble
	332	0.956, 0.965, 0.969, 0.953, 0.892, 0.826, &
	333	! Nitrate insoluble
	334	0.765, 0.754, 0.743, 0.737, 0.739, 0.651 /
	335
[2644]	336	! Added by R. Wang (July 31 2016)
	337	! properties for BC assuming Maxwell-Garnett rule and internal mixture
	338
[2739]	339	DATA alpha_MG_6bands/ &
	340	!--BC content=0.001
	341	3.546, 3.546, 3.546, 3.546, 3.640, 3.891, 3.883, 4.354, 4.967, 5.424, 6.412, 9.049, &
	342	3.995, 3.995, 3.995, 3.995, 4.073, 4.091, 4.111, 4.209, 4.502, 4.878, 5.591, 8.029, &
	343	4.650, 4.650, 4.650, 4.650, 4.661, 4.561, 4.472, 4.386, 4.373, 4.451, 4.726, 6.248, &
	344	4.445, 4.445, 4.445, 4.445, 4.393, 4.291, 4.171, 4.017, 3.812, 3.682, 3.514, 3.593, &
	345	2.837, 2.837, 2.837, 2.837, 2.824, 2.798, 2.760, 2.696, 2.579, 2.471, 2.265, 1.735, &
	346	1.087, 1.087, 1.087, 1.087, 1.095, 1.101, 1.107, 1.110, 1.109, 1.101, 1.074, 0.944, &
	347	!--BC content=0.010
	348	3.550, 3.550, 3.550, 3.550, 3.623, 3.852, 4.006, 4.342, 4.945, 5.381, 6.394, 8.972, &
	349	3.995, 3.995, 3.995, 3.995, 4.082, 4.094, 4.133, 4.235, 4.517, 4.871, 5.601, 8.001, &
	350	4.653, 4.653, 4.653, 4.653, 4.662, 4.564, 4.475, 4.393, 4.376, 4.447, 4.736, 6.221, &
	351	4.434, 4.434, 4.434, 4.434, 4.386, 4.286, 4.170, 4.021, 3.815, 3.680, 3.524, 3.583, &
	352	2.831, 2.831, 2.831, 2.831, 2.820, 2.795, 2.757, 2.695, 2.579, 2.472, 2.270, 1.734, &
	353	1.086, 1.086, 1.086, 1.086, 1.095, 1.101, 1.106, 1.110, 1.108, 1.101, 1.075, 0.944, &
	354	!--BC content=0.020
	355	3.549, 3.549, 3.549, 3.549, 3.630, 3.842, 4.016, 4.335, 4.928, 5.383, 6.369, 8.937, &
	356	4.004, 4.004, 4.004, 4.004, 4.096, 4.100, 4.145, 4.243, 4.529, 4.870, 5.595, 7.973, &
	357	4.650, 4.650, 4.650, 4.650, 4.662, 4.564, 4.478, 4.397, 4.379, 4.449, 4.734, 6.207, &
	358	4.422, 4.422, 4.422, 4.422, 4.378, 4.280, 4.167, 4.020, 3.817, 3.683, 3.527, 3.583, &
	359	2.824, 2.824, 2.824, 2.824, 2.815, 2.791, 2.754, 2.693, 2.579, 2.472, 2.272, 1.740, &
	360	1.085, 1.085, 1.085, 1.085, 1.094, 1.100, 1.105, 1.109, 1.108, 1.100, 1.075, 0.945, &
	361	!--BC content=0.050
	362	3.563, 3.563, 3.563, 3.563, 3.662, 3.834, 4.019, 4.318, 4.889, 5.360, 6.305, 8.842, &
	363	4.039, 4.039, 4.039, 4.039, 4.130, 4.126, 4.166, 4.261, 4.544, 4.863, 5.573, 7.895, &
	364	4.641, 4.641, 4.641, 4.641, 4.659, 4.565, 4.483, 4.407, 4.389, 4.455, 4.731, 6.161, &
	365	4.388, 4.388, 4.388, 4.388, 4.356, 4.264, 4.157, 4.017, 3.821, 3.691, 3.537, 3.582, &
	366	2.805, 2.805, 2.805, 2.805, 2.801, 2.779, 2.744, 2.687, 2.577, 2.474, 2.279, 1.757, &
	367	1.081, 1.081, 1.081, 1.081, 1.091, 1.098, 1.103, 1.107, 1.106, 1.099, 1.074, 0.949, &
	368	!--BC content=0.100
	369	3.607, 3.607, 3.607, 3.607, 3.701, 3.843, 4.021, 4.302, 4.841, 5.306, 6.214, 8.689, &
	370	4.098, 4.098, 4.098, 4.098, 4.175, 4.168, 4.201, 4.291, 4.554, 4.854, 5.529, 7.763, &
	371	4.627, 4.627, 4.627, 4.627, 4.651, 4.567, 4.490, 4.421, 4.403, 4.464, 4.723, 6.082, &
	372	4.333, 4.333, 4.333, 4.333, 4.319, 4.236, 4.139, 4.011, 3.828, 3.704, 3.554, 3.580, &
	373	2.773, 2.773, 2.773, 2.773, 2.778, 2.758, 2.727, 2.675, 2.573, 2.476, 2.291, 1.786, &
	374	1.076, 1.076, 1.076, 1.076, 1.087, 1.094, 1.099, 1.103, 1.103, 1.097, 1.073, 0.954, &
	375	!--BC content=0.200
	376	3.702, 3.702, 3.702, 3.702, 3.768, 3.878, 4.030, 4.278, 4.757, 5.190, 6.032, 8.370, &
	377	4.204, 4.204, 4.204, 4.204, 4.255, 4.246, 4.270, 4.345, 4.571, 4.837, 5.440, 7.494, &
	378	4.598, 4.598, 4.598, 4.598, 4.633, 4.565, 4.502, 4.445, 4.428, 4.481, 4.706, 5.923, &
	379	4.230, 4.230, 4.230, 4.230, 4.244, 4.179, 4.101, 3.996, 3.838, 3.728, 3.588, 3.581, &
	380	2.710, 2.710, 2.710, 2.710, 2.730, 2.716, 2.693, 2.651, 2.564, 2.480, 2.314, 1.846, &
	381	1.064, 1.064, 1.064, 1.064, 1.078, 1.085, 1.091, 1.096, 1.097, 1.091, 1.072, 0.965 /
	382
	383	DATA piz_MG_6bands/ &
	384	!--BC content=0.001
	385	-3.984, -3.984, -3.984, -3.984, -3.395, -3.075, -3.051, -2.567, -2.053, -1.807, -1.308, -0.592, &
	386	-3.014, -3.014, -3.014, -3.014, -2.520, -2.473, -2.411, -2.294, -2.027, -1.771, -1.368, -0.613, &
	387	-1.142, -1.142, -1.142, -1.142, -0.930, -0.961, -0.981, -1.001, -0.981, -0.934, -0.798, -0.337, &
	388	-0.202, -0.202, -0.202, -0.202, -0.129, -0.154, -0.183, -0.224, -0.284, -0.326, -0.383, -0.347, &
	389	0.133, 0.133, 0.133, 0.133, 0.158, 0.146, 0.131, 0.107, 0.061, 0.015, -0.081, -0.433, &
	390	0.137, 0.137, 0.137, 0.137, 0.148, 0.149, 0.149, 0.147, 0.139, 0.128, 0.098, -0.048, &
	391	!--BC content=0.010
	392	-3.640, -3.640, -3.640, -3.640, -3.219, -2.942, -2.761, -2.445, -1.983, -1.731, -1.262, -0.576, &
	393	-2.791, -2.791, -2.791, -2.791, -2.383, -2.352, -2.290, -2.180, -1.943, -1.707, -1.321, -0.595, &
	394	-1.069, -1.069, -1.069, -1.069, -0.889, -0.920, -0.945, -0.966, -0.954, -0.910, -0.775, -0.333, &
	395	-0.185, -0.185, -0.185, -0.185, -0.119, -0.144, -0.173, -0.213, -0.275, -0.318, -0.372, -0.346, &
	396	0.138, 0.138, 0.138, 0.138, 0.160, 0.149, 0.134, 0.110, 0.063, 0.018, -0.076, -0.430, &
	397	0.137, 0.137, 0.137, 0.137, 0.148, 0.149, 0.150, 0.148, 0.140, 0.129, 0.099, -0.046, &
	398	!--BC content=0.020
	399	-3.369, -3.369, -3.369, -3.369, -3.033, -2.799, -2.619, -2.339, -1.913, -1.657, -1.223, -0.558, &
	400	-2.592, -2.592, -2.592, -2.592, -2.254, -2.239, -2.185, -2.091, -1.869, -1.652, -1.284, -0.580, &
	401	-1.004, -1.004, -1.004, -1.004, -0.849, -0.882, -0.909, -0.932, -0.926, -0.885, -0.758, -0.326, &
	402	-0.168, -0.168, -0.168, -0.168, -0.108, -0.133, -0.163, -0.203, -0.265, -0.309, -0.364, -0.341, &
	403	0.142, 0.142, 0.142, 0.142, 0.163, 0.152, 0.137, 0.113, 0.066, 0.022, -0.072, -0.422, &
	404	0.138, 0.138, 0.138, 0.138, 0.149, 0.150, 0.150, 0.148, 0.140, 0.129, 0.100, -0.044, &
	405	!--BC content=0.050
	406	-2.743, -2.743, -2.743, -2.743, -2.582, -2.438, -2.292, -2.078, -1.729, -1.495, -1.123, -0.512, &
	407	-2.121, -2.121, -2.121, -2.121, -1.939, -1.951, -1.926, -1.865, -1.688, -1.511, -1.188, -0.541, &
	408	-0.839, -0.839, -0.839, -0.839, -0.742, -0.780, -0.812, -0.842, -0.848, -0.818, -0.707, -0.307, &
	409	-0.121, -0.121, -0.121, -0.121, -0.078, -0.103, -0.133, -0.174, -0.236, -0.281, -0.337, -0.325, &
	410	0.154, 0.154, 0.154, 0.154, 0.171, 0.160, 0.145, 0.121, 0.075, 0.032, -0.060, -0.399, &
	411	0.141, 0.141, 0.141, 0.141, 0.150, 0.151, 0.152, 0.150, 0.142, 0.131, 0.103, -0.038, &
	412	!--BC content=0.100
	413	-2.025, -2.025, -2.025, -2.025, -2.035, -1.965, -1.875, -1.731, -1.475, -1.285, -0.980, -0.447, &
	414	-1.569, -1.569, -1.569, -1.569, -1.540, -1.572, -1.580, -1.557, -1.444, -1.312, -1.052, -0.485, &
	415	-0.628, -0.628, -0.628, -0.628, -0.593, -0.634, -0.672, -0.710, -0.732, -0.716, -0.631, -0.279, &
	416	-0.054, -0.054, -0.054, -0.054, -0.032, -0.058, -0.088, -0.129, -0.191, -0.236, -0.295, -0.299, &
	417	0.173, 0.173, 0.173, 0.173, 0.184, 0.173, 0.158, 0.135, 0.090, 0.048, -0.039, -0.361, &
	418	0.145, 0.145, 0.145, 0.145, 0.153, 0.154, 0.154, 0.152, 0.144, 0.134, 0.107, -0.027, &
	419	!--BC content=0.200
	420	-1.155, -1.155, -1.155, -1.155, -1.300, -1.295, -1.272, -1.213, -1.077, -0.957, -0.747, -0.336, &
	421	-0.884, -0.884, -0.884, -0.884, -0.979, -1.028, -1.065, -1.086, -1.054, -0.986, -0.820, -0.388, &
	422	-0.332, -0.332, -0.332, -0.332, -0.360, -0.402, -0.444, -0.490, -0.531, -0.536, -0.492, -0.225, &
	423	0.050, 0.050, 0.050, 0.050, 0.048, 0.022, -0.007, -0.047, -0.107, -0.152, -0.213, -0.245, &
	424	0.204, 0.204, 0.204, 0.204, 0.208, 0.197, 0.182, 0.161, 0.120, 0.081, 0.002, -0.285, &
	425	0.153, 0.153, 0.153, 0.153, 0.159, 0.159, 0.159, 0.157, 0.150, 0.140, 0.115, -0.005 /
	426
	427	DATA cg_MG_6bands/ &
	428	!--BC content=0.001
	429	0.290, 0.290, 0.290, 0.290, 0.378, 0.381, 0.377, 0.378, 0.373, 0.344, 0.317, 0.118, &
	430	0.374, 0.374, 0.374, 0.374, 0.469, 0.474, 0.480, 0.480, 0.475, 0.462, 0.438, 0.242, &
	431	0.375, 0.375, 0.375, 0.375, 0.491, 0.502, 0.517, 0.527, 0.534, 0.529, 0.508, 0.255, &
	432	0.242, 0.242, 0.242, 0.242, 0.455, 0.516, 0.577, 0.631, 0.687, 0.715, 0.748, 0.772, &
	433	0.559, 0.559, 0.559, 0.559, 0.487, 0.461, 0.410, 0.303, -0.181, -3.594, 1.814, 1.164, &
	434	0.323, 0.323, 0.323, 0.323, 0.339, 0.339, 0.337, 0.328, 0.295, 0.243, 0.045, 2.778, &
	435	!--BC content=0.010
	436	0.300, 0.300, 0.300, 0.300, 0.385, 0.389, 0.390, 0.386, 0.377, 0.353, 0.319, 0.116, &
	437	0.384, 0.384, 0.384, 0.384, 0.478, 0.482, 0.488, 0.488, 0.483, 0.471, 0.442, 0.246, &
	438	0.380, 0.380, 0.380, 0.380, 0.496, 0.507, 0.520, 0.531, 0.538, 0.534, 0.512, 0.259, &
	439	0.233, 0.233, 0.233, 0.233, 0.454, 0.520, 0.581, 0.636, 0.691, 0.719, 0.752, 0.778, &
	440	0.553, 0.553, 0.553, 0.553, 0.487, 0.461, 0.414, 0.308, -0.144, -2.858, 1.889, 1.165, &
	441	0.323, 0.323, 0.323, 0.323, 0.338, 0.339, 0.337, 0.328, 0.295, 0.244, 0.051, 2.847, &
	442	!--BC content=0.020
	443	0.305, 0.305, 0.305, 0.305, 0.391, 0.394, 0.395, 0.390, 0.380, 0.358, 0.321, 0.112, &
	444	0.391, 0.391, 0.391, 0.391, 0.484, 0.488, 0.493, 0.493, 0.488, 0.476, 0.445, 0.246, &
	445	0.384, 0.384, 0.384, 0.384, 0.499, 0.511, 0.524, 0.534, 0.541, 0.537, 0.514, 0.260, &
	446	0.219, 0.219, 0.219, 0.219, 0.451, 0.523, 0.586, 0.642, 0.696, 0.724, 0.756, 0.783, &
	447	0.547, 0.547, 0.547, 0.547, 0.487, 0.461, 0.416, 0.315, -0.108, -2.308, 1.954, 1.170, &
	448	0.322, 0.322, 0.322, 0.322, 0.338, 0.339, 0.337, 0.328, 0.296, 0.245, 0.056, 2.949, &
	449	!--BC content=0.050
	450	0.312, 0.312, 0.312, 0.312, 0.400, 0.402, 0.403, 0.398, 0.385, 0.364, 0.321, 0.094, &
	451	0.402, 0.402, 0.402, 0.402, 0.495, 0.498, 0.503, 0.503, 0.496, 0.484, 0.451, 0.242, &
	452	0.388, 0.388, 0.388, 0.388, 0.506, 0.518, 0.530, 0.541, 0.548, 0.544, 0.521, 0.259, &
	453	0.152, 0.152, 0.152, 0.152, 0.434, 0.527, 0.599, 0.658, 0.711, 0.738, 0.768, 0.794, &
	454	0.534, 0.534, 0.534, 0.534, 0.486, 0.463, 0.423, 0.332, -0.018, -1.355, 2.203, 1.188, &
	455	0.321, 0.321, 0.321, 0.321, 0.337, 0.338, 0.336, 0.328, 0.297, 0.249, 0.071, 3.311, &
	456	!--BC content=0.100
	457	0.309, 0.309, 0.309, 0.309, 0.404, 0.405, 0.406, 0.401, 0.385, 0.364, 0.316, 0.057, &
	458	0.406, 0.406, 0.406, 0.406, 0.503, 0.507, 0.511, 0.512, 0.504, 0.491, 0.456, 0.232, &
	459	0.382, 0.382, 0.382, 0.382, 0.510, 0.523, 0.537, 0.548, 0.555, 0.551, 0.528, 0.255, &
	460	-0.178, -0.178, -0.178, -0.178, 0.312, 0.524, 0.626, 0.692, 0.741, 0.765, 0.791, 0.816, &
	461	0.519, 0.519, 0.519, 0.519, 0.486, 0.466, 0.432, 0.357, 0.091, -0.650, 2.968, 1.221, &
	462	0.319, 0.319, 0.319, 0.319, 0.336, 0.337, 0.335, 0.328, 0.300, 0.255, 0.094, 4.314, &
	463	!--BC content=0.200
	464	0.267, 0.267, 0.267, 0.267, 0.388, 0.391, 0.392, 0.387, 0.370, 0.345, 0.287, -0.053, &
	465	0.381, 0.381, 0.381, 0.381, 0.499, 0.506, 0.512, 0.514, 0.507, 0.494, 0.455, 0.198, &
	466	0.318, 0.318, 0.318, 0.318, 0.494, 0.515, 0.535, 0.551, 0.561, 0.559, 0.535, 0.237, &
	467	1.151, 1.151, 1.151, 1.151, 0.703, 0.643, 1.248, 0.874, 0.846, 0.846, 0.853, 0.868, &
	468	0.502, 0.502, 0.502, 0.502, 0.486, 0.472, 0.447, 0.393, 0.226, -0.108,-43.753, 1.313, &
	469	0.315, 0.315, 0.315, 0.315, 0.332, 0.334, 0.333, 0.327, 0.304, 0.266, 0.135, 21.079 /
[2644]	470
[2004]	471	!----BEGINNING OF CALCULATIONS
	472
[2748]	473	modname='aeropt_6bands_rrtm'
	474
[2004]	475	spsol = 0
	476	spinsol = 0
[2634]	477
[2146]	478	IF (NSW.NE.nbands_sw_rrtm) THEN
[2750]	479	CALL abort_physic(modname,'Erreur NSW doit etre egal a 6 pour cette routine',1)
[2146]	480	ENDIF
[2004]	481
[2634]	482	zdh(:,:)=pdel(:,:)/(RG*zrho(:,:)) ! m
[2004]	483
	484	IF (flag_aerosol .EQ. 1) THEN
	485	nb_aer = 2
	486	ALLOCATE (aerosol_name(nb_aer))
[2146]	487	aerosol_name(1) = id_ASSO4M_phy
	488	aerosol_name(2) = id_CSSO4M_phy
[2004]	489	ELSEIF (flag_aerosol .EQ. 2) THEN
	490	nb_aer = 2
	491	ALLOCATE (aerosol_name(nb_aer))
[2146]	492	aerosol_name(1) = id_ASBCM_phy
	493	aerosol_name(2) = id_AIBCM_phy
[2004]	494	ELSEIF (flag_aerosol .EQ. 3) THEN
	495	nb_aer = 2
	496	ALLOCATE (aerosol_name(nb_aer))
[2146]	497	aerosol_name(1) = id_ASPOMM_phy
	498	aerosol_name(2) = id_AIPOMM_phy
[2004]	499	ELSEIF (flag_aerosol .EQ. 4) THEN
	500	nb_aer = 3
	501	ALLOCATE (aerosol_name(nb_aer))
[2146]	502	aerosol_name(1) = id_CSSSM_phy
	503	aerosol_name(2) = id_SSSSM_phy
	504	aerosol_name(3) = id_ASSSM_phy
[2004]	505	ELSEIF (flag_aerosol .EQ. 5) THEN
	506	nb_aer = 1
	507	ALLOCATE (aerosol_name(nb_aer))
[2146]	508	aerosol_name(1) = id_CIDUSTM_phy
[2004]	509	ELSEIF (flag_aerosol .EQ. 6) THEN
[2738]	510	nb_aer = 13
[2004]	511	ALLOCATE (aerosol_name(nb_aer))
[2146]	512	aerosol_name(1) = id_ASSO4M_phy
	513	aerosol_name(2) = id_ASBCM_phy
	514	aerosol_name(3) = id_AIBCM_phy
	515	aerosol_name(4) = id_ASPOMM_phy
	516	aerosol_name(5) = id_AIPOMM_phy
	517	aerosol_name(6) = id_CSSSM_phy
	518	aerosol_name(7) = id_SSSSM_phy
	519	aerosol_name(8) = id_ASSSM_phy
	520	aerosol_name(9) = id_CIDUSTM_phy
	521	aerosol_name(10)= id_CSSO4M_phy
[2738]	522	aerosol_name(11)= id_CSNO3M_phy
	523	aerosol_name(12)= id_ASNO3M_phy
	524	aerosol_name(13)= id_CINO3M_phy
[2004]	525	ENDIF
	526
	527	!
	528	! loop over modes, use of precalculated nmd and corresponding sigma
	529	! loop over wavelengths
	530	! for each mass species in mode
	531	! interpolate from Sext to retrieve Sext_at_gridpoint_per_species
	532	! compute optical_thickness_at_gridpoint_per_species
	533
	534	DO n=1,nbre_RH-1
	535	fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n))
	536	ENDDO
	537
[2634]	538	DO k=1, klev
	539	DO i=1, klon
[2004]	540	rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
[2146]	541	RH_num(i,k) = INT(rh(i,k)/10. + 1.)
[2004]	542	IF (rh(i,k).GT.85.) RH_num(i,k)=10
	543	IF (rh(i,k).GT.90.) RH_num(i,k)=11
	544	delta(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k))
	545	ENDDO
	546	ENDDO
	547
[2634]	548	tau_ae(:,:,:,:)=0.
	549	piz_ae(:,:,:,:)=0.
	550	cg_ae(:,:,:,:)=0.
[2644]	551
	552	tau_ae_pi(:,:,:,:)=0.
	553	piz_ae_pi(:,:,:,:)=0.
	554	cg_ae_pi(:,:,:,:)=0.
[2004]	555
	556	DO m=1,nb_aer ! tau is only computed for each mass
	557	fac=1.0
[2146]	558	IF (aerosol_name(m).EQ.id_ASBCM_phy) THEN
[2004]	559	soluble=.TRUE.
	560	spsol=1
[2146]	561	ELSEIF (aerosol_name(m).EQ.id_ASPOMM_phy) THEN
[2004]	562	soluble=.TRUE.
	563	spsol=2
[2738]	564	ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
[2004]	565	soluble=.TRUE.
	566	spsol=3
[2810]	567	!fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
	568	fac=0.0 !--6 March 2017 - OB as Didier H said CSSO4 should not be used
[2738]	569	ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
[2004]	570	soluble=.TRUE.
	571	spsol=4
[2644]	572	fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
[2146]	573	ELSEIF (aerosol_name(m).EQ.id_SSSSM_phy) THEN
[2004]	574	soluble=.TRUE.
	575	spsol=5
[2146]	576	ELSEIF (aerosol_name(m).EQ.id_CSSSM_phy) THEN
[2004]	577	soluble=.TRUE.
	578	spsol=6
[2146]	579	ELSEIF (aerosol_name(m).EQ.id_ASSSM_phy) THEN
[2004]	580	soluble=.TRUE.
	581	spsol=7
[2738]	582	ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
	583	soluble=.TRUE.
	584	spsol=8
	585	fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD
	586	ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
	587	soluble=.TRUE.
	588	spsol=9
	589	fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD
[2146]	590	ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
[2004]	591	soluble=.FALSE.
	592	spinsol=1
[2146]	593	ELSEIF (aerosol_name(m).EQ.id_AIBCM_phy) THEN
[2004]	594	soluble=.FALSE.
	595	spinsol=2
[2146]	596	ELSEIF (aerosol_name(m).EQ.id_AIPOMM_phy) THEN
[2004]	597	soluble=.FALSE.
	598	spinsol=3
[2738]	599	ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
	600	soluble=.FALSE.
	601	spinsol=4
	602	fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD
[2004]	603	ELSE
	604	CYCLE
	605	ENDIF
	606
[2644]	607	!--shortname for aerosol index
[2738]	608	aerindex=aerosol_name(m)
[2004]	609
	610	IF (soluble) THEN
	611
[2644]	612	!--here we treat the special case of soluble BC internal mixture with Maxwell-Garnett rule
	613	IF (spsol.EQ.1 .AND. flag_bc_internal_mixture) THEN
[2004]	614
[2644]	615	DO k=1, klev
	616	DO i=1, klon
[2004]	617
[2738]	618	tmp_var=m_allaer(i,k,aerindex)/1.e6zdh(i,k)fac
	619	tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6zdh(i,k)fac
[2644]	620
	621	! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols
	622	BC_massfra = m_allaer(i,k,id_ASBCM_phy)/(m_allaer(i,k,id_ASBCM_phy)+m_allaer(i,k,id_ASSO4M_phy))
	623
[2739]	624	IF (BC_massfra.GE.0.20) THEN
	625	classbc = 6
	626	ELSEIF (BC_massfra.GE.0.10) THEN
[2644]	627	classbc = 5
	628	ELSEIF (BC_massfra.GE.0.05) THEN
	629	classbc = 4
	630	ELSEIF (BC_massfra.GE.0.02) THEN
	631	classbc = 3
	632	ELSEIF (BC_massfra.GE.0.01) THEN
	633	classbc = 2
	634	ELSE
	635	classbc = 1
	636	ENDIF
	637
	638	! Calculate the dry BC/(BC+SUL) mass ratio for natural aerosols
	639	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))
	640
[2739]	641	IF (BC_massfra_pi.GE.0.20) THEN
	642	classbc_pi = 6
	643	ELSEIF (BC_massfra_pi.GE.0.10) THEN
[2644]	644	classbc_pi = 5
	645	ELSEIF (BC_massfra_pi.GE.0.05) THEN
	646	classbc_pi = 4
	647	ELSEIF (BC_massfra_pi.GE.0.02) THEN
	648	classbc_pi = 3
	649	ELSEIF (BC_massfra_pi.GE.0.01) THEN
	650	classbc_pi = 2
	651	ELSE
	652	classbc_pi = 1
	653	ENDIF
	654
[2739]	655	DO inu=1, NSW
[2644]	656
	657	!--all aerosols
	658	tau_ae2b_int= alpha_MG_6bands(RH_num(i,k),inu,classbc)+ &
	659	delta(i,k)* (alpha_MG_6bands(RH_num(i,k)+1,inu,classbc) - &
	660	alpha_MG_6bands(RH_num(i,k),inu,classbc))
[2004]	661
[2644]	662	piz_ae2b_int = piz_MG_6bands(RH_num(i,k),inu,classbc) + &
	663	delta(i,k)* (piz_MG_6bands(RH_num(i,k)+1,inu,classbc) - &
	664	piz_MG_6bands(RH_num(i,k),inu,classbc))
	665
	666	cg_ae2b_int = cg_MG_6bands(RH_num(i,k),inu,classbc) + &
	667	delta(i,k)* (cg_MG_6bands(RH_num(i,k)+1,inu,classbc) - &
	668	cg_MG_6bands(RH_num(i,k),inu,classbc))
	669
[2738]	670	tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
	671	piz_ae(i,k,aerindex,inu) = piz_ae2b_int
	672	cg_ae(i,k,aerindex,inu) = cg_ae2b_int
[2644]	673
	674	!--natural aerosols
	675	tau_ae2b_int= alpha_MG_6bands(RH_num(i,k),inu,classbc_pi)+ &
	676	delta(i,k)* (alpha_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) - &
	677	alpha_MG_6bands(RH_num(i,k),inu,classbc_pi))
	678
	679	piz_ae2b_int = piz_MG_6bands(RH_num(i,k),inu,classbc_pi) + &
	680	delta(i,k)* (piz_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) - &
	681	piz_MG_6bands(RH_num(i,k),inu,classbc_pi))
	682
	683	cg_ae2b_int = cg_MG_6bands(RH_num(i,k),inu,classbc_pi) + &
	684	delta(i,k)* (cg_MG_6bands(RH_num(i,k)+1,inu,classbc_pi) - &
	685	cg_MG_6bands(RH_num(i,k),inu,classbc_pi))
	686
[2738]	687	tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
	688	piz_ae_pi(i,k,aerindex,inu) = piz_ae2b_int
	689	cg_ae_pi(i,k,aerindex,inu) = cg_ae2b_int
[2644]	690
	691	ENDDO
	692	ENDDO
	693	ENDDO
	694
	695	!--else treat all other cases of soluble aerosols
	696	ELSE
	697
	698	DO k=1, klev
	699	DO i=1, klon
[2738]	700	tmp_var=m_allaer(i,k,aerindex)/1.e6zdh(i,k)fac
	701	tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6zdh(i,k)fac
[2644]	702
[2739]	703	DO inu=1, NSW
[2644]	704
	705	tau_ae2b_int= alpha_aers_6bands(RH_num(i,k),inu,spsol)+ &
	706	delta(i,k)* (alpha_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
	707	alpha_aers_6bands(RH_num(i,k),inu,spsol))
	708
	709	piz_ae2b_int = piz_aers_6bands(RH_num(i,k),inu,spsol) + &
[2004]	710	delta(i,k)* (piz_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
[2644]	711	piz_aers_6bands(RH_num(i,k),inu,spsol))
[2004]	712
[2644]	713	cg_ae2b_int = cg_aers_6bands(RH_num(i,k),inu,spsol) + &
	714	delta(i,k)* (cg_aers_6bands(RH_num(i,k)+1,inu,spsol) - &
	715	cg_aers_6bands(RH_num(i,k),inu,spsol))
[2004]	716
[2738]	717	tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
	718	tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
	719	piz_ae(i,k,aerindex,inu) = piz_ae2b_int
	720	cg_ae(i,k,aerindex,inu) = cg_ae2b_int
[2644]	721
	722	ENDDO
[2004]	723	ENDDO
	724	ENDDO
[2644]	725
	726	!--external mixture case for soluble BC
	727	IF (spsol.EQ.1) THEN
[2738]	728	piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:)
	729	cg_ae_pi(:,:,aerindex,:) = cg_ae(:,:,aerindex,:)
[2644]	730	ENDIF
	731
	732	ENDIF
[2004]	733
	734	ELSE ! For all aerosol insoluble components
	735
[2634]	736	DO k=1, klev
	737	DO i=1, klon
[2738]	738	tmp_var=m_allaer(i,k,aerindex)/1.e6zdh(i,k)fac
	739	tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6zdh(i,k)fac
[2004]	740
[2739]	741	DO inu=1, NSW
[2004]	742	tau_ae2b_int = alpha_aeri_6bands(inu,spinsol)
	743	piz_ae2b_int = piz_aeri_6bands(inu,spinsol)
	744	cg_ae2b_int = cg_aeri_6bands(inu,spinsol)
	745
[2738]	746	tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
	747	tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi*tau_ae2b_int
	748	piz_ae(i,k,aerindex,inu) = piz_ae2b_int
	749	cg_ae(i,k,aerindex,inu)= cg_ae2b_int
[2004]	750	ENDDO
	751	ENDDO
	752	ENDDO
	753
	754	ENDIF ! soluble / insoluble
	755
	756	ENDDO ! nb_aer
	757
	758	DO inu=1, NSW
[2634]	759	DO k=1, klev
	760	DO i=1, klon
[2644]	761	!--all (natural + anthropogenic) aerosol
[2146]	762	tau_allaer(i,k,2,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)+ &
	763	tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+ &
	764	tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ &
	765	tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+ &
	766	tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu)
[2644]	767	tau_allaer(i,k,2,inu)=MAX(tau_allaer(i,k,2,inu),tau_min)
[2004]	768
[2644]	769	piz_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
	770	tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ &
	771	tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ &
	772	tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ &
	773	tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ &
	774	tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ &
	775	tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ &
	776	tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ &
	777	tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ &
[2146]	778	tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
[2004]	779	/tau_allaer(i,k,2,inu)
[2644]	780	piz_allaer(i,k,2,inu)=MIN(MAX(piz_allaer(i,k,2,inu),0.01),1.0)
	781	IF (tau_allaer(i,k,2,inu).LE.tau_min) piz_allaer(i,k,2,inu)=1.0
[2004]	782
[2146]	783	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)+ &
	784	tau_ae(i,k,id_CSSO4M_phy,inu)piz_ae(i,k,id_CSSO4M_phy,inu)cg_ae(i,k,id_CSSO4M_phy,inu)+ &
[2644]	785	tau_ae(i,k,id_ASBCM_phy,inu)piz_ae(i,k,id_ASBCM_phy,inu)cg_ae(i,k,id_ASBCM_phy,inu)+ &
	786	tau_ae(i,k,id_AIBCM_phy,inu)piz_ae(i,k,id_AIBCM_phy,inu)cg_ae(i,k,id_AIBCM_phy,inu)+ &
[2146]	787	tau_ae(i,k,id_ASPOMM_phy,inu)piz_ae(i,k,id_ASPOMM_phy,inu)cg_ae(i,k,id_ASPOMM_phy,inu)+ &
	788	tau_ae(i,k,id_AIPOMM_phy,inu)piz_ae(i,k,id_AIPOMM_phy,inu)cg_ae(i,k,id_AIPOMM_phy,inu)+ &
[2644]	789	tau_ae(i,k,id_ASSSM_phy,inu)piz_ae(i,k,id_ASSSM_phy,inu)cg_ae(i,k,id_ASSSM_phy,inu)+ &
	790	tau_ae(i,k,id_CSSSM_phy,inu)piz_ae(i,k,id_CSSSM_phy,inu)cg_ae(i,k,id_CSSSM_phy,inu)+ &
	791	tau_ae(i,k,id_SSSSM_phy,inu)piz_ae(i,k,id_SSSSM_phy,inu)cg_ae(i,k,id_SSSSM_phy,inu)+ &
[2146]	792	tau_ae(i,k,id_CIDUSTM_phy,inu)piz_ae(i,k,id_CIDUSTM_phy,inu)cg_ae(i,k,id_CIDUSTM_phy,inu))/ &
[2004]	793	(tau_allaer(i,k,2,inu)*piz_allaer(i,k,2,inu))
[2644]	794	cg_allaer(i,k,2,inu)=MIN(MAX(cg_allaer(i,k,2,inu),0.0),1.0)
[2004]	795
	796	!--natural aerosol
[2644]	797	!--ASBCM aerosols take _pi value because of internal mixture option
[2146]	798	tau_allaer(i,k,1,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ &
	799	tau_ae_pi(i,k,id_ASBCM_phy,inu)+tau_ae_pi(i,k,id_AIBCM_phy,inu)+ &
	800	tau_ae_pi(i,k,id_ASPOMM_phy,inu)+tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ &
	801	tau_ae_pi(i,k,id_ASSSM_phy,inu)+tau_ae_pi(i,k,id_CSSSM_phy,inu)+ &
	802	tau_ae_pi(i,k,id_SSSSM_phy,inu)+ tau_ae_pi(i,k,id_CIDUSTM_phy,inu)
[2644]	803	tau_allaer(i,k,1,inu)=MAX(tau_allaer(i,k,1,inu),tau_min)
[2004]	804
[2644]	805	piz_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
	806	tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ &
	807	tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae_pi(i,k,id_ASBCM_phy,inu)+ &
	808	tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ &
	809	tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ &
	810	tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ &
	811	tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ &
	812	tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ &
	813	tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ &
[2146]	814	tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
[2004]	815	/tau_allaer(i,k,1,inu)
[2644]	816	piz_allaer(i,k,1,inu)=MIN(MAX(piz_allaer(i,k,1,inu),0.01),1.0)
	817	IF (tau_allaer(i,k,1,inu).LE.tau_min) piz_allaer(i,k,1,inu)=1.0
[2004]	818
[2644]	819	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)+ &
	820	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)+ &
	821	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)+ &
	822	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)+ &
	823	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)+ &
	824	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)+ &
	825	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)+ &
	826	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)+ &
	827	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)+ &
[2146]	828	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))/ &
[2004]	829	(tau_allaer(i,k,1,inu)*piz_allaer(i,k,1,inu))
[2644]	830	cg_allaer(i,k,1,inu)=MIN(MAX(cg_allaer(i,k,1,inu),0.0),1.0)
[2004]	831
	832	ENDDO
	833	ENDDO
	834	ENDDO
	835
	836	DEALLOCATE(aerosol_name)
	837
	838	END SUBROUTINE AEROPT_6BANDS_RRTM

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

Download in other formats:

Original Format