Changeset 2787 for LMDZ5/branches/testing/libf/phylmd/rrtm
- Timestamp:
- Jan 30, 2017, 5:54:45 PM (8 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 8 edited
-
. (modified) (1 prop)
-
libf/phylmd/rrtm/aeropt_5wv_rrtm.F90 (modified) (12 diffs)
-
libf/phylmd/rrtm/aeropt_6bands_rrtm.F90 (modified) (24 diffs)
-
libf/phylmd/rrtm/aeropt_lw_rrtm.F90 (modified) (1 diff)
-
libf/phylmd/rrtm/read_rsun_rrtm.F90 (modified) (2 diffs)
-
libf/phylmd/rrtm/readaerosol_optic_rrtm.F90 (modified) (4 diffs)
-
libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90 (modified) (3 diffs)
-
libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2728,2730-2775,2777-2785
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90
r2669 r2787 12 12 USE DIMPHY 13 13 USE aero_mod 14 USE phys_local_var_mod, only: od550aer,od865aer,ec550aer,od550lt1aer 14 USE phys_local_var_mod, ONLY: od550aer,od865aer,ec550aer,od550lt1aer 15 USE YOMCST, ONLY: RD,RG 15 16 16 17 ! … … 20 21 ! a partir d'une sous-routine de Johannes Quaas pour les sulfates 21 22 ! Olivier Boucher mars 2014 pour adaptation RRTM 23 ! 22 24 ! 23 25 ! Refractive indices for seasalt come from Shettle and Fenn (1979) … … 50 52 ! 51 53 IMPLICIT NONE 52 INCLUDE "YOMCST.h"53 54 ! 54 55 ! Input arguments: 55 56 ! 56 REAL, DIMENSION(klon,klev), INTENT( in) :: pdel57 REAL, DIMENSION(klon,klev,naero_tot), INTENT( in) :: m_allaer58 REAL, DIMENSION(klon,klev), INTENT( in) :: RHcl ! humidite relative ciel clair59 INTEGER,INTENT( in) :: flag_aerosol60 LOGICAL,INTENT( in) :: flag_bc_internal_mixture61 REAL, DIMENSION(klon,klev), INTENT( in) :: pplay62 REAL, DIMENSION(klon,klev), INTENT( in) :: t_seri57 REAL, DIMENSION(klon,klev), INTENT(IN) :: pdel 58 REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer 59 REAL, DIMENSION(klon,klev), INTENT(IN) :: RHcl ! humidite relative ciel clair 60 INTEGER,INTENT(IN) :: flag_aerosol 61 LOGICAL,INTENT(IN) :: flag_bc_internal_mixture 62 REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay 63 REAL, DIMENSION(klon,klev), INTENT(IN) :: t_seri 63 64 ! 64 65 ! Output arguments: 65 66 ! 66 REAL, DIMENSION(klon), INTENT( out) :: ai ! POLDER aerosol index67 REAL, DIMENSION(klon,nwave,naero_tot), INTENT( out) :: tausum68 REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT( out) :: tau67 REAL, DIMENSION(klon), INTENT(OUT) :: ai ! POLDER aerosol index 68 REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT) :: tausum 69 REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT) :: tau 69 70 70 71 ! … … 83 84 INTEGER, PARAMETER :: la865 = 5 84 85 INTEGER, PARAMETER :: nbre_RH=12 85 INTEGER, PARAMETER :: naero_soluble=7 ! 1- BC soluble; 2- POM soluble; 3- SO4 acc. 86 ! 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc. 87 INTEGER, PARAMETER :: naero_insoluble=3 ! 1- Dust; 2- BC insoluble; 3- POM insoluble 86 INTEGER, PARAMETER :: naero_soluble=9 ! 1- BC soluble; 2- POM soluble; 3- SO4 coarse 87 ! 4- SO4 acc; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc. 88 ! 8- NO3 coarse; 9- NO3 accumulation 89 INTEGER, PARAMETER :: naero_insoluble=4 ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble 88 90 89 91 REAL :: zrho … … 101 103 REAL :: zdh(klon,klev) 102 104 103 ! Soluble components 1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-coarse; 6 seasalt coarse; 7 seasalt acc.104 105 REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble) ! Ext. coeff. ** m2/g 105 ! Insoluble components 1- Dust: 2- BC; 3- POM106 106 REAL :: alpha_aeri_5wv(las,naero_insoluble) ! Ext. coeff. ** m2/g 107 107 … … 109 109 ! BC internal mixture 110 110 ! 111 INTEGER, PARAMETER :: nbclassbc = 5! Added by Rong Wang/OB for the 5 fractions112 ! of BC in the soluble mode:113 ! bc_content/0.001, 0.01, 0.02, 0.05, !0.1/111 INTEGER, PARAMETER :: nbclassbc = 6 ! Added by Rong Wang/OB for the 5 fractions 112 ! of BC in the soluble mode: 113 ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/ 114 114 ! for Maxwell-Garnet internal mixture 115 115 ! Detailed theory can be found in R. Wang Estimation of global black carbon ! direct … … 124 124 INTEGER :: n, classbc 125 125 126 ! From here on we look at the optical parameters at 5 wavelengths: 127 ! 443nm, 550, 670, 765 and 865 nm 128 ! le 12 AVRIL 2006 129 ! 130 DATA alpha_aers_5wv/ & 131 ! bc soluble 132 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 133 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 134 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 135 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 136 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501, & 137 ! pom soluble 138 6.676,6.676,6.676,6.676,6.710,6.934,7.141,7.569,8.034,8.529,9.456,10.511, & 139 5.109,5.109,5.109,5.109,5.189,5.535,5.960,6.852,8.008,9.712,12.897,19.676, & 140 3.718,3.718,3.718,3.718,3.779,4.042,4.364,5.052,5.956,7.314,9.896,15.688, & 141 2.849,2.849,2.849,2.849,2.897,3.107,3.365,3.916,4.649,5.760,7.900,12.863, & 142 2.229,2.229,2.229,2.229,2.268,2.437,2.645,3.095,3.692,4.608,6.391,10.633, & 143 ! Sulfate (Accumulation) 144 5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, & 145 4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133, & 146 3.079,3.351,3.639,3.886,4.205,4.584,5.206,6.019,7.648,9.234,12.391,17.220, & 147 2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313, & 148 1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922, & 149 ! Sulfate (Coarse) 150 5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, & 151 4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133, & 152 3.079,3.351,3.639,3.886,4.205,4.584,5.206,6.019,7.648,9.234,12.391,17.220, & 153 2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313, & 154 1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922, & 155 ! seasalt seasalt Super Coarse Soluble (SS) 156 0.218, 0.272, 0.293, 0.316, 0.343, 0.380, 0.429, 0.501, 0.636, 0.755, 0.967, 1.495, & 157 0.221, 0.275, 0.297, 0.320, 0.348, 0.383, 0.432, 0.509, 0.640, 0.759, 0.972, 1.510, & 158 0.224, 0.279, 0.301, 0.324, 0.352, 0.388, 0.438, 0.514, 0.647, 0.768, 0.985, 1.514, & 159 0.227, 0.282, 0.303, 0.327, 0.356, 0.392, 0.441, 0.518, 0.652, 0.770, 0.987, 1.529, & 160 0.230, 0.285, 0.306, 0.330, 0.359, 0.396, 0.446, 0.522, 0.656, 0.777, 0.993, 1.539, & 161 ! seasalt seasalt Coarse Soluble (CS) 162 0.578, 0.706, 0.756, 0.809, 0.876, 0.964, 1.081, 1.256, 1.577, 1.858, 2.366, 3.613, & 163 0.598, 0.725, 0.779, 0.833, 0.898, 0.990, 1.109, 1.290, 1.609, 1.889, 2.398, 3.682, & 164 0.619, 0.750, 0.802, 0.857, 0.927, 1.022, 1.141, 1.328, 1.648, 1.939, 2.455, 3.729, & 165 0.633, 0.767, 0.820, 0.879, 0.948, 1.044, 1.167, 1.353, 1.683, 1.969, 2.491, 3.785, & 166 0.648, 0.785, 0.838, 0.896, 0.967, 1.066, 1.192, 1.381, 1.714, 2.006, 2.531, 3.836, & 167 ! seasalt seasalt Accumulation Soluble (AS) 168 4.432, 5.899, 6.505, 7.166, 7.964, 7.962, 9.232,11.257,14.979,18.337,24.223,37.811, & 169 3.298, 4.569, 5.110, 5.709, 6.446, 6.268, 7.396, 9.246,12.787,16.113,22.197,37.136, & 170 2.340, 3.358, 3.803, 4.303, 4.928, 4.696, 5.629, 7.198,10.308,13.342,19.120,34.296, & 171 1.789, 2.626, 2.999, 3.422, 3.955, 3.730, 4.519, 5.864, 8.593,11.319,16.653,31.331, & 172 1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043 / 173 174 DATA alpha_aeri_5wv/ & 175 ! dust insoluble 176 0.759, 0.770, 0.775, 0.775, 0.772, & 177 !!jb bc insoluble 178 11.536,10.033, 8.422, 7.234, 6.270, & 179 ! pom insoluble 180 5.042, 3.101, 1.890, 1.294, 0.934/ 181 182 ! internal mixture 183 DATA alpha_MG_5wv/ & 184 ! bc content = 0.001 185 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 186 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 187 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 188 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 189 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501, & 190 ! bc content = 0.01 191 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 192 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 193 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 194 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 195 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501, & 196 ! bc content = 0.02 197 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 198 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 199 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 200 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 201 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501, & 202 ! bc content = 0.05 203 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 204 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 205 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 206 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 207 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501, & 208 ! bc content = 0.10 209 7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, & 210 7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, & 211 7.195,7.195,7.195,7.195,7.195,7.195,7.195,7.195,9.551,10.847,12.381,13.994, & 212 6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687, & 213 6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501 / 214 126 ! From here on we look at the optical parameters at 5 wavelengths: 443, 550, 670, 765 and 865 nm 127 128 DATA alpha_aers_5wv/ & 129 ! BC Accumulation Soluble (AS) 130 5.342, 5.342, 5.342, 5.342, 5.342, 5.829, 6.344, 7.470, 8.603, 8.736, 8.870,10.149, & 131 5.159, 5.159, 5.159, 5.159, 5.159, 5.608, 6.083, 7.121, 8.169, 8.293, 8.418, 9.612, & 132 4.849, 4.849, 4.849, 4.849, 4.849, 5.251, 5.674, 6.598, 7.533, 7.644, 7.756, 8.829, & 133 4.573, 4.573, 4.573, 4.573, 4.573, 4.936, 5.318, 6.152, 6.996, 7.096, 7.198, 8.171, & 134 4.274, 4.274, 4.274, 4.274, 4.274, 4.600, 4.942, 5.686, 6.441, 6.530, 6.621, 7.495, & 135 ! POM Accumulation Soluble (AS) 136 5.300, 5.300, 5.300, 5.300, 5.300, 5.827, 6.392, 7.640, 8.898, 9.046, 9.195,10.606, & 137 4.569, 4.569, 4.569, 4.569, 4.569, 5.029, 5.528, 6.649, 7.802, 7.939, 8.077, 9.400, & 138 3.768, 3.768, 3.768, 3.768, 3.768, 4.152, 4.573, 5.533, 6.538, 6.658, 6.780, 7.955, & 139 3.210, 3.210, 3.210, 3.210, 3.210, 3.542, 3.909, 4.752, 5.644, 5.751, 5.860, 6.916, & 140 2.709, 2.709, 2.709, 2.709, 2.709, 2.994, 3.309, 4.041, 4.823, 4.917, 5.013, 5.949, & 141 ! Sulfate Coarse Soluble (CS) 142 0.702, 0.702, 0.702, 0.702, 0.947, 1.025, 1.127, 1.266, 1.490, 1.675, 2.003, 2.857, & 143 0.725, 0.725, 0.725, 0.725, 0.977, 1.057, 1.163, 1.304, 1.529, 1.718, 2.051, 2.914, & 144 0.751, 0.751, 0.751, 0.751, 1.011, 1.093, 1.200, 1.345, 1.576, 1.768, 2.110, 2.973, & 145 0.769, 0.769, 0.769, 0.769, 1.034, 1.120, 1.227, 1.375, 1.613, 1.811, 2.153, 3.032, & 146 0.786, 0.786, 0.786, 0.786, 1.056, 1.144, 1.254, 1.406, 1.646, 1.850, 2.202, 3.088, & 147 !-- Sulfate Accumulation (BC content=0) 148 4.639, 4.639, 4.639, 4.639, 6.244, 6.878, 7.684, 8.805,10.638,12.174,14.880,21.828, & 149 3.966, 3.966, 3.966, 3.966, 5.359, 5.950, 6.707, 7.771, 9.540,11.046,13.742,20.884, & 150 3.234, 3.234, 3.234, 3.234, 4.393, 4.914, 5.587, 6.543, 8.160, 9.556,12.101,19.072, & 151 2.721, 2.721, 2.721, 2.721, 3.712, 4.175, 4.774, 5.634, 7.101, 8.383,10.747,17.381, & 152 2.262, 2.262, 2.262, 2.262, 3.102, 3.505, 4.030, 4.789, 6.097, 7.251, 9.403,15.581, & 153 ! Seasalt Super Coarse Soluble (SS) 154 0.194, 0.237, 0.254, 0.275, 0.299, 0.327, 0.366, 0.432, 0.544, 0.642, 0.824, 1.265, & 155 0.196, 0.240, 0.257, 0.278, 0.303, 0.331, 0.371, 0.437, 0.550, 0.648, 0.831, 1.274, & 156 0.198, 0.243, 0.260, 0.283, 0.306, 0.335, 0.376, 0.442, 0.557, 0.654, 0.839, 1.285, & 157 0.201, 0.246, 0.263, 0.286, 0.308, 0.338, 0.380, 0.445, 0.559, 0.660, 0.846, 1.289, & 158 0.203, 0.249, 0.266, 0.289, 0.312, 0.341, 0.384, 0.449, 0.564, 0.665, 0.852, 1.297, & 159 ! Seasalt Coarse Soluble (CS) 160 0.576, 0.690, 0.738, 0.789, 0.855, 0.935, 1.046, 1.212, 1.512, 1.785, 2.258, 3.449, & 161 0.595, 0.713, 0.763, 0.814, 0.880, 0.963, 1.079, 1.248, 1.550, 1.826, 2.306, 3.507, & 162 0.617, 0.738, 0.789, 0.842, 0.911, 0.996, 1.113, 1.286, 1.592, 1.871, 2.369, 3.562, & 163 0.632, 0.755, 0.808, 0.862, 0.931, 1.018, 1.140, 1.316, 1.626, 1.909, 2.409, 3.622, & 164 0.645, 0.771, 0.825, 0.880, 0.951, 1.039, 1.164, 1.344, 1.661, 1.948, 2.455, 3.682, & 165 ! Seasalt Accumulation Soluble (AS) 166 3.684, 4.367, 4.711, 5.074, 5.438, 6.046, 6.793, 7.964,10.200,12.246,15.959,24.642, & 167 3.126, 3.717, 4.023, 4.349, 4.673, 5.229, 5.918, 7.018, 9.179,11.208,14.994,24.184, & 168 2.482, 2.973, 3.233, 3.511, 3.788, 4.272, 4.876, 5.858, 7.836, 9.739,13.393,22.658, & 169 2.086, 2.509, 2.735, 2.979, 3.220, 3.649, 4.186, 5.068, 6.874, 8.642,12.099,21.146, & 170 1.737, 2.097, 2.292, 2.503, 2.711, 3.086, 3.556, 4.337, 5.960, 7.571,10.779,19.427, & 171 ! Nitrate Coarse Soluble (CS) 172 0.726, 0.726, 0.726, 0.796, 0.868, 0.947, 1.041, 1.246, 1.563, 1.872, 2.328, 2.447, & 173 0.753, 0.753, 0.753, 0.825, 0.900, 0.979, 1.075, 1.285, 1.610, 1.922, 2.385, 2.503, & 174 0.780, 0.780, 0.780, 0.854, 0.932, 1.013, 1.113, 1.326, 1.656, 1.979, 2.447, 2.579, & 175 0.797, 0.797, 0.797, 0.874, 0.953, 1.035, 1.138, 1.356, 1.697, 2.020, 2.495, 2.621, & 176 0.811, 0.811, 0.811, 0.890, 0.971, 1.055, 1.160, 1.384, 1.733, 2.062, 2.547, 2.675, & 177 ! Nitrate Accumulation Soluble (AS) 178 4.208, 4.208, 4.208, 4.693, 5.217, 5.778, 6.502, 8.108,10.722,13.327,17.185,18.210, & 179 3.386, 3.386, 3.386, 3.808, 4.268, 4.768, 5.420, 6.897, 9.377,11.923,15.803,16.852, & 180 2.650, 2.650, 2.650, 2.997, 3.380, 3.801, 4.357, 5.638, 7.850,10.189,13.858,14.870, & 181 2.174, 2.174, 2.174, 2.471, 2.802, 3.167, 3.652, 4.784, 6.774, 8.917,12.345,13.302, & 182 1.776, 1.776, 1.776, 2.028, 2.309, 2.622, 3.040, 4.026, 5.787, 7.717,10.858,11.745 / 183 184 DATA alpha_aeri_5wv/ & 185 ! Dust insoluble 186 0.788, 0.818, 0.842, 0.851, 0.853, & 187 ! BC insoluble 188 5.342, 5.159, 4.849, 4.573, 4.274, & 189 ! POM insoluble 190 5.300, 4.569, 3.768, 3.210, 2.709, & 191 ! Nitrate insoluble 192 0.726, 0.753, 0.780, 0.797, 0.811 / 193 194 ! Added by R. Wang (July 31 2016) 195 ! properties for BC assuming Maxwell-Garnett rule and internal mixture 196 197 DATA alpha_MG_5wv/ & 198 !--BC content=0.001 199 4.293, 4.293, 4.293, 4.293, 4.320, 4.342, 4.271, 4.320, 4.476, 4.772, 5.310, 7.434, & 200 4.687, 4.687, 4.687, 4.687, 4.693, 4.602, 4.492, 4.413, 4.374, 4.462, 4.729, 6.274, & 201 4.802, 4.802, 4.802, 4.802, 4.776, 4.646, 4.516, 4.371, 4.231, 4.173, 4.217, 5.072, & 202 4.716, 4.716, 4.716, 4.716, 4.668, 4.548, 4.408, 4.249, 4.047, 3.951, 3.850, 4.259, & 203 4.520, 4.520, 4.520, 4.520, 4.461, 4.353, 4.230, 4.069, 3.850, 3.707, 3.524, 3.565, & 204 !--BC content=0.010 205 4.298, 4.298, 4.298, 4.298, 4.343, 4.333, 4.283, 4.325, 4.472, 4.751, 5.298, 7.402, & 206 4.692, 4.692, 4.692, 4.692, 4.695, 4.598, 4.499, 4.410, 4.383, 4.454, 4.739, 6.260, & 207 4.796, 4.796, 4.796, 4.796, 4.768, 4.644, 4.518, 4.376, 4.230, 4.172, 4.225, 5.048, & 208 4.708, 4.708, 4.708, 4.708, 4.659, 4.543, 4.411, 4.256, 4.053, 3.945, 3.855, 4.242, & 209 4.509, 4.509, 4.509, 4.509, 4.456, 4.351, 4.229, 4.072, 3.852, 3.707, 3.531, 3.560, & 210 !--BC content=0.020 211 4.301, 4.301, 4.301, 4.301, 4.353, 4.330, 4.291, 4.326, 4.478, 4.738, 5.288, 7.393, & 212 4.688, 4.688, 4.688, 4.688, 4.695, 4.596, 4.500, 4.412, 4.386, 4.454, 4.737, 6.248, & 213 4.787, 4.787, 4.787, 4.787, 4.761, 4.641, 4.516, 4.378, 4.231, 4.176, 4.226, 5.041, & 214 4.696, 4.696, 4.696, 4.696, 4.651, 4.538, 4.409, 4.256, 4.055, 3.948, 3.858, 4.240, & 215 4.497, 4.497, 4.497, 4.497, 4.448, 4.345, 4.225, 4.072, 3.854, 3.709, 3.535, 3.561, & 216 !--BC content=0.050 217 4.318, 4.318, 4.318, 4.318, 4.377, 4.337, 4.310, 4.334, 4.488, 4.724, 5.267, 7.342, & 218 4.678, 4.678, 4.678, 4.678, 4.693, 4.595, 4.506, 4.421, 4.396, 4.458, 4.734, 6.203, & 219 4.760, 4.760, 4.760, 4.760, 4.742, 4.631, 4.512, 4.381, 4.237, 4.185, 4.229, 5.015, & 220 4.662, 4.662, 4.662, 4.662, 4.629, 4.522, 4.401, 4.254, 4.062, 3.955, 3.867, 4.229, & 221 4.461, 4.461, 4.461, 4.461, 4.424, 4.328, 4.215, 4.068, 3.858, 3.718, 3.545, 3.562, & 222 !--BC content=0.100 223 4.348, 4.348, 4.348, 4.348, 4.404, 4.361, 4.337, 4.358, 4.503, 4.717, 5.240, 7.239, & 224 4.662, 4.662, 4.662, 4.662, 4.685, 4.596, 4.513, 4.437, 4.411, 4.468, 4.729, 6.123, & 225 4.716, 4.716, 4.716, 4.716, 4.713, 4.613, 4.505, 4.384, 4.249, 4.199, 4.235, 4.974, & 226 4.607, 4.607, 4.607, 4.607, 4.593, 4.497, 4.387, 4.252, 4.072, 3.969, 3.882, 4.212, & 227 4.403, 4.403, 4.403, 4.403, 4.385, 4.299, 4.196, 4.061, 3.865, 3.731, 3.564, 3.563, & 228 !--BC content=0.200 229 4.401, 4.401, 4.401, 4.401, 4.447, 4.409, 4.389, 4.405, 4.529, 4.715, 5.183, 7.007, & 230 4.631, 4.631, 4.631, 4.631, 4.666, 4.594, 4.526, 4.463, 4.439, 4.488, 4.714, 5.958, & 231 4.633, 4.633, 4.633, 4.633, 4.654, 4.575, 4.488, 4.387, 4.271, 4.224, 4.250, 4.894, & 232 4.505, 4.505, 4.505, 4.505, 4.520, 4.444, 4.356, 4.243, 4.089, 3.997, 3.912, 4.179, & 233 4.295, 4.295, 4.295, 4.295, 4.307, 4.239, 4.157, 4.045, 3.876, 3.757, 3.602, 3.569 / 215 234 ! 216 235 ! Initialisations … … 252 271 aerosol_name(1) = id_CIDUSTM_phy 253 272 ELSEIF (flag_aerosol .EQ. 6) THEN 254 nb_aer = 1 0273 nb_aer = 13 255 274 ALLOCATE (aerosol_name(nb_aer)) 256 275 aerosol_name(1) = id_ASSO4M_phy … … 264 283 aerosol_name(9) = id_CIDUSTM_phy 265 284 aerosol_name(10)= id_CSSO4M_phy 285 aerosol_name(11)= id_CSNO3M_phy 286 aerosol_name(12)= id_ASNO3M_phy 287 aerosol_name(13)= id_CINO3M_phy 266 288 ENDIF 267 289 … … 315 337 soluble=.TRUE. 316 338 spsol=7 339 ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN 340 soluble=.TRUE. 341 spsol=8 342 fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD 343 ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN 344 soluble=.TRUE. 345 spsol=9 346 fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD 317 347 ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN 318 348 soluble=.FALSE. … … 324 354 soluble=.FALSE. 325 355 spinsol=3 356 ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN 357 soluble=.FALSE. 358 spinsol=4 359 fac=1.2903 ! NO3NH4/NO3 / mass conversion factor for AOD 326 360 ELSE 327 361 CYCLE … … 345 379 BC_massfra = m_allaer(i,k,id_ASBCM_phy)/(m_allaer(i,k,id_ASBCM_phy)+m_allaer(i,k,id_ASSO4M_phy)) 346 380 347 IF (BC_massfra.GE.0.10) THEN 381 IF (BC_massfra.GE.0.20) THEN 382 classbc = 6 383 ELSEIF (BC_massfra.GE.0.10) THEN 348 384 classbc = 5 349 385 ELSEIF (BC_massfra.GE.0.05) THEN -
LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90
r2669 r2787 11 11 USE aero_mod 12 12 USE phys_local_var_mod, only: absvisaer 13 USE YOMCST, ONLY: RG 13 14 14 15 ! Yves Balkanski le 12 avril 2006 … … 20 21 ! 21 22 IMPLICIT NONE 22 23 INCLUDE "YOMCST.h" 23 ! 24 24 INCLUDE "clesphys.h" 25 26 25 ! 27 26 ! Input arguments: 28 27 ! 29 REAL, DIMENSION(klon,klev), INTENT( in) :: pdel30 REAL, DIMENSION(klon,klev,naero_tot), INTENT( in) :: m_allaer31 REAL, DIMENSION(klon,klev,naero_tot), INTENT( in) :: m_allaer_pi32 REAL, DIMENSION(klon,klev), INTENT( in) :: RHcl ! humidite relative ciel clair33 INTEGER, INTENT( in) :: flag_aerosol34 LOGICAL, INTENT( in) :: flag_bc_internal_mixture35 REAL, DIMENSION(klon,klev), INTENT( in) :: zrho28 REAL, DIMENSION(klon,klev), INTENT(IN) :: pdel 29 REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer 30 REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer_pi 31 REAL, DIMENSION(klon,klev), INTENT(IN) :: RHcl ! humidite relative ciel clair 32 INTEGER, INTENT(IN) :: flag_aerosol 33 LOGICAL, INTENT(IN) :: flag_bc_internal_mixture 34 REAL, DIMENSION(klon,klev), INTENT(IN) :: zrho 36 35 ! 37 36 ! Output arguments: … … 39 38 ! 1= natural aerosols 40 39 ! 41 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: tau_allaer ! epaisseur optique aerosol 42 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: piz_allaer ! single scattering albedo aerosol 43 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: cg_allaer ! asymmetry parameter aerosol 44 40 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: tau_allaer ! epaisseur optique aerosol 41 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: piz_allaer ! single scattering albedo aerosol 42 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: cg_allaer ! asymmetry parameter aerosol 45 43 ! 46 44 ! Local 47 45 ! 48 LOGICAL :: 46 LOGICAL :: soluble 49 47 INTEGER :: i, k,n, inu, m 50 48 INTEGER :: spsol, spinsol … … 53 51 INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA 54 52 55 INTEGER, PARAMETER :: naero_soluble=7 ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse 56 ! 5- seasalt super coarse 6- seasalt coarse 7- seasalt acc. 57 INTEGER, PARAMETER :: naero_insoluble=3 ! 1- Dust; 2- BC insoluble; 3- POM insoluble 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 55 ! coarse 7- seasalt acc. 56 ! 8- NO3 coarse 9- NO3 acc. 57 INTEGER, PARAMETER :: naero_insoluble=4 ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble 58 58 59 59 INTEGER, PARAMETER :: nbre_RH=12 … … 67 67 REAL :: fac 68 68 REAL :: zdh(klon,klev) 69 INTEGER, ALLOCATABLE, DIMENSION(:) 69 INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name 70 70 INTEGER :: nb_aer 71 71 72 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: 73 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: 74 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: 75 76 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: 72 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae 73 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: piz_ae 74 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: cg_ae 75 76 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae_pi 77 77 REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: piz_ae_pi 78 78 REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: cg_ae_pi … … 80 80 ! Proprietes optiques 81 81 ! 82 REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit m2/g SO483 REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble) 84 REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit 85 REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble) 86 REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !-- unit87 REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble) !-- unit82 REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit m2/g 83 REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble) !--unit m2/g 84 REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unitless 85 REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble) !--unitless 86 REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unitless 87 REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble) !--unitless 88 88 ! 89 89 ! BC internal mixture 90 90 ! 91 INTEGER, PARAMETER :: nbclassbc = 5! Added by Rong Wang/OB for the 5 fractions92 ! of BC in the soluble mode:93 ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/91 INTEGER, PARAMETER :: nbclassbc = 6 ! Added by Rong Wang/OB for the 5 fractions 92 ! of BC in the soluble mode: 93 ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1, 0.2/ 94 94 ! for Maxwell-Garnet internal mixture 95 95 ! Detailed theory can be found in R. Wang Estimation of global black carbon direct … … 100 100 REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc) 101 101 ! 102 INTEGER :: id, classbc, classbc_pi102 INTEGER :: aerindex, classbc, classbc_pi 103 103 REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi 104 105 ! 106 REAL, PARAMETER :: tau_min = 1.e-15 107 ! REAL, PARAMETER :: tau_min = 1.e-7 104 CHARACTER*20 :: modname 105 ! 106 REAL, PARAMETER :: tau_min = 1.e-7 108 107 109 108 !*************************************************************************** … … 112 111 113 112 DATA alpha_aers_6bands/ & 114 ! bc soluble AS 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 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 117 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 118 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 119 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 120 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, & 121 ! pom soluble AS 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 6.443, 6.443, 6.443, 6.443, 6.443, 7.100, 7.809, 9.277,10.721,10.882,11.056,12.628, & 124 4.381, 4.381, 4.381, 4.381, 4.381, 4.828, 5.310, 6.309, 7.290, 7.400, 7.518, 8.587, & 125 1.846, 1.846, 1.846, 1.846, 1.846, 2.034, 2.237, 2.658, 3.072, 3.118, 3.168, 3.618, & 126 0.377, 0.377, 0.377, 0.377, 0.377, 0.415, 0.456, 0.542, 0.627, 0.636, 0.646, 0.738, & 127 0.052, 0.052, 0.052, 0.052, 0.052, 0.057, 0.063, 0.075, 0.087, 0.088, 0.089, 0.102, & 128 ! sulfate AS 129 6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, & 130 6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, & 131 4.381, 4.381, 4.381, 4.828, 5.301, 5.792, 6.309, 7.176, 9.542,11.868,11.868,21.204, & 132 1.727, 1.727, 1.727, 1.903, 2.090, 2.283, 2.487, 2.829, 3.762, 4.679, 4.679, 8.359, & 133 0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, & 134 0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, & 135 ! sulfate coarse CS 136 0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, & 137 0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, & 138 0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, & 139 0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, & 140 0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, & 141 0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, & 142 ! seasalt seasalt Super Coarse Soluble (SS) 143 0.214, 0.267, 0.287, 0.310, 0.337, 0.373, 0.421, 0.494, 0.625, 0.742, 0.956, 1.480, & 144 0.217, 0.270, 0.291, 0.314, 0.341, 0.377, 0.426, 0.499, 0.632, 0.751, 0.963, 1.490, & 145 0.221, 0.275, 0.297, 0.320, 0.348, 0.384, 0.434, 0.507, 0.642, 0.762, 0.976, 1.506, & 146 0.230, 0.285, 0.308, 0.331, 0.359, 0.396, 0.447, 0.522, 0.658, 0.780, 0.997, 1.536, & 147 0.250, 0.307, 0.330, 0.354, 0.384, 0.424, 0.477, 0.556, 0.696, 0.822, 1.044, 1.592, & 148 0.279, 0.347, 0.373, 0.401, 0.434, 0.478, 0.537, 0.625, 0.781, 0.918, 1.158, 1.744, & 149 ! seasalt seasalt Coarse Soluble (CS) 150 0.550, 0.676, 0.724, 0.779, 0.841, 0.925, 1.040, 1.213, 1.523, 1.803, 2.306, 3.535, & 151 0.568, 0.695, 0.744, 0.798, 0.863, 0.950, 1.066, 1.240, 1.556, 1.839, 2.342, 3.588, & 152 0.599, 0.727, 0.779, 0.834, 0.901, 0.993, 1.111, 1.290, 1.612, 1.899, 2.411, 3.669, & 153 0.647, 0.786, 0.841, 0.899, 0.969, 1.069, 1.193, 1.384, 1.717, 2.015, 2.543, 3.842, & 154 0.663, 0.825, 0.889, 0.957, 1.038, 1.126, 1.268, 1.483, 1.862, 2.196, 2.780, 4.183, & 155 0.479, 0.644, 0.711, 0.785, 0.873, 0.904, 1.048, 1.275, 1.692, 2.072, 2.753, 4.430, & 156 ! seasalt seasalt Accumulation Soluble (AS) 157 6.128, 7.564, 8.127, 8.723, 9.421,10.014,11.250,13.113,16.320,19.061,23.675,34.158, & 158 5.080, 6.586, 7.197, 7.856, 8.645, 8.820,10.115,12.140,15.775,18.983,24.499,36.993, & 159 3.290, 4.541, 5.073, 5.663, 6.387, 6.227, 7.339, 9.161,12.645,15.918,21.908,36.673, & 160 1.389, 2.066, 2.371, 2.719, 3.161, 2.970, 3.623, 4.747, 7.064, 9.419,14.120,27.565, & 161 0.309, 0.497, 0.587, 0.693, 0.832, 0.777, 0.986, 1.364, 2.207, 3.136, 5.174,12.065, & 162 0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885 / 113 ! BC Accumulation Soluble (AS) 114 5.276, 5.276, 5.276, 5.276, 5.276, 5.797, 6.353, 7.568, 8.779, 8.920, 9.062,10.403, & 115 5.371, 5.371, 5.371, 5.371, 5.371, 5.870, 6.401, 7.559, 8.720, 8.857, 8.994,10.298, & 116 5.126, 5.126, 5.126, 5.126, 5.126, 5.572, 6.043, 7.073, 8.113, 8.235, 8.359, 9.544, & 117 4.225, 4.225, 4.225, 4.225, 4.225, 4.549, 4.887, 5.624, 6.370, 6.459, 6.548, 7.412, & 118 2.667, 2.667, 2.667, 2.667, 2.667, 2.840, 3.015, 3.382, 3.745, 3.788, 3.832, 4.251, & 119 1.113, 1.113, 1.113, 1.113, 1.113, 1.210, 1.307, 1.506, 1.699, 1.721, 1.744, 1.961, & 120 ! POM Accumulation Soluble (AS) 121 5.962, 5.962, 5.962, 5.962, 5.962, 6.544, 7.153, 8.460, 9.732, 9.877,10.026,11.398, & 122 5.609, 5.609, 5.609, 5.609, 5.609, 6.163, 6.752, 8.041, 9.324, 9.473, 9.624,11.044, & 123 4.527, 4.527, 4.527, 4.527, 4.527, 4.982, 5.476, 6.586, 7.726, 7.861, 7.998, 9.306, & 124 2.697, 2.697, 2.697, 2.697, 2.697, 2.980, 3.293, 4.017, 4.790, 4.884, 4.979, 5.902, & 125 0.918, 0.918, 0.918, 0.918, 0.918, 1.025, 1.146, 1.435, 1.754, 1.794, 1.834, 2.232, & 126 0.105, 0.105, 0.105, 0.105, 0.105, 0.154, 0.211, 0.344, 0.491, 0.509, 0.528, 0.710, & 127 ! Sulfate Coarse Soluble (CS) 128 0.670, 0.670, 0.670, 0.670, 0.905, 0.981, 1.081, 1.217, 1.434, 1.615, 1.940, 2.773, & 129 0.689, 0.689, 0.689, 0.689, 0.930, 1.009, 1.110, 1.247, 1.469, 1.655, 1.980, 2.825, & 130 0.726, 0.726, 0.726, 0.726, 0.978, 1.060, 1.164, 1.306, 1.533, 1.723, 2.056, 2.915, & 131 0.785, 0.785, 0.785, 0.785, 1.056, 1.144, 1.254, 1.406, 1.650, 1.852, 2.203, 3.103, & 132 0.795, 0.795, 0.795, 0.795, 1.064, 1.160, 1.282, 1.449, 1.720, 1.945, 2.340, 3.348, & 133 0.502, 0.502, 0.502, 0.502, 0.656, 0.728, 0.820, 0.950, 1.170, 1.360, 1.706, 2.666, & 134 !-- Sulfate Accumulation (BC content=0) 135 5.246, 5.246, 5.246, 5.246, 7.041, 7.658, 8.437, 9.503,11.204,12.608,15.030,21.080, & 136 4.924, 4.924, 4.924, 4.924, 6.619, 7.258, 8.067, 9.186,10.998,12.507,15.142,21.822, & 137 3.925, 3.925, 3.925, 3.925, 5.304, 5.888, 6.635, 7.686, 9.433,10.921,13.587,20.659, & 138 2.248, 2.248, 2.248, 2.248, 3.079, 3.477, 3.993, 4.738, 6.022, 7.153, 9.261,15.301, & 139 0.720, 0.720, 0.720, 0.720, 1.014, 1.169, 1.374, 1.678, 2.225, 2.726, 3.703, 6.779, & 140 0.069, 0.069, 0.069, 0.069, 0.102, 0.121, 0.146, 0.185, 0.258, 0.328, 0.472, 0.974, & 141 ! Seasalt Super Coarse Soluble (SS) 142 0.189, 0.232, 0.249, 0.270, 0.293, 0.322, 0.361, 0.425, 0.536, 0.634, 0.813, 1.250, & 143 0.192, 0.235, 0.252, 0.274, 0.296, 0.325, 0.365, 0.429, 0.542, 0.639, 0.821, 1.259, & 144 0.196, 0.240, 0.257, 0.279, 0.302, 0.331, 0.371, 0.437, 0.550, 0.649, 0.832, 1.274, & 145 0.204, 0.249, 0.267, 0.289, 0.313, 0.342, 0.384, 0.450, 0.565, 0.666, 0.851, 1.300, & 146 0.222, 0.270, 0.289, 0.313, 0.338, 0.369, 0.413, 0.482, 0.602, 0.705, 0.896, 1.355, & 147 0.244, 0.293, 0.313, 0.338, 0.364, 0.397, 0.443, 0.517, 0.646, 0.757, 0.960, 1.443, & 148 ! Seasalt Coarse Soluble (CS) 149 0.547, 0.657, 0.705, 0.754, 0.817, 0.896, 1.008, 1.169, 1.456, 1.724, 2.199, 3.358, & 150 0.566, 0.679, 0.727, 0.776, 0.840, 0.920, 1.032, 1.196, 1.492, 1.760, 2.238, 3.416, & 151 0.596, 0.714, 0.764, 0.816, 0.882, 0.965, 1.081, 1.250, 1.552, 1.828, 2.310, 3.509, & 152 0.644, 0.771, 0.825, 0.880, 0.951, 1.040, 1.164, 1.345, 1.666, 1.957, 2.462, 3.700, & 153 0.640, 0.772, 0.829, 0.887, 0.965, 1.061, 1.198, 1.398, 1.758, 2.085, 2.658, 4.031, & 154 0.452, 0.562, 0.609, 0.659, 0.728, 0.813, 0.938, 1.125, 1.471, 1.797, 2.384, 3.855, & 155 ! Seasalt Accumulation Soluble (AS) 156 4.270, 5.032, 5.395, 5.772, 6.158, 6.768, 7.510, 8.634,10.695,12.512,15.740,23.122, & 157 3.952, 4.673, 5.030, 5.403, 5.781, 6.399, 7.156, 8.328,10.533,12.524,16.091,24.329, & 158 3.090, 3.676, 3.979, 4.301, 4.622, 5.172, 5.853, 6.941, 9.077,11.083,14.830,23.951, & 159 1.710, 2.069, 2.262, 2.471, 2.678, 3.048, 3.513, 4.285, 5.882, 7.467,10.618,19.102, & 160 0.529, 0.663, 0.737, 0.819, 0.901, 1.050, 1.244, 1.575, 2.301, 3.067, 4.704, 9.712, & 161 0.073, 0.142, 0.177, 0.215, 0.259, 0.328, 0.422, 0.579, 0.913, 1.263, 2.011, 4.358, & 162 ! Nitrate Coarse Soluble (CS) 163 0.688, 0.688, 0.688, 0.755, 0.826, 0.900, 0.994, 1.190, 1.501, 1.804, 2.251, 2.374, & 164 0.712, 0.712, 0.712, 0.781, 0.853, 0.929, 1.023, 1.225, 1.540, 1.846, 2.299, 2.421, & 165 0.754, 0.754, 0.754, 0.826, 0.902, 0.980, 1.077, 1.286, 1.612, 1.925, 2.388, 2.509, & 166 0.809, 0.809, 0.809, 0.887, 0.969, 1.053, 1.160, 1.384, 1.732, 2.069, 2.557, 2.685, & 167 0.786, 0.786, 0.786, 0.867, 0.952, 1.043, 1.157, 1.403, 1.793, 2.171, 2.725, 2.871, & 168 0.703, 0.703, 0.703, 0.767, 0.836, 0.909, 1.002, 1.210, 1.552, 1.901, 2.434, 2.579, & 169 ! Nitrate Accumulation Soluble (AS) 170 5.251, 5.251, 5.251, 5.779, 6.336, 6.920, 7.655, 9.230,11.678,14.026,17.391,18.279, & 171 4.643, 4.643, 4.643, 5.152, 5.696, 6.276, 7.016, 8.640,11.238,13.787,17.511,18.494, & 172 3.373, 3.373, 3.373, 3.790, 4.245, 4.739, 5.383, 6.842, 9.289,11.804,15.637,16.675, & 173 1.779, 1.779, 1.779, 2.029, 2.309, 2.619, 3.033, 4.008, 5.746, 7.645,10.730,11.601, & 174 0.549, 0.549, 0.549, 0.639, 0.742, 0.859, 1.019, 1.412, 2.157, 3.029, 4.552, 5.002, & 175 0.367, 0.367, 0.367, 0.416, 0.471, 0.532, 0.614, 0.808, 1.166, 1.575, 2.281, 2.489 / 176 177 DATA piz_aers_6bands/& 178 ! BC Accumulation Soluble (AS) 179 0.480, 0.480, 0.480, 0.480, 0.480, 0.487, 0.495, 0.513, 0.530, 0.532, 0.534, 0.553, & 180 0.469, 0.469, 0.469, 0.469, 0.469, 0.474, 0.480, 0.497, 0.514, 0.517, 0.519, 0.538, & 181 0.456, 0.456, 0.456, 0.456, 0.456, 0.464, 0.473, 0.495, 0.519, 0.522, 0.524, 0.550, & 182 0.427, 0.427, 0.427, 0.427, 0.427, 0.437, 0.448, 0.476, 0.505, 0.509, 0.512, 0.543, & 183 0.349, 0.349, 0.349, 0.349, 0.349, 0.358, 0.369, 0.396, 0.425, 0.429, 0.432, 0.465, & 184 0.196, 0.196, 0.196, 0.196, 0.196, 0.198, 0.202, 0.213, 0.226, 0.227, 0.229, 0.243, & 185 ! POM Accumulation Soluble (AS) 186 0.956, 0.956, 0.956, 0.956, 0.956, 0.961, 0.964, 0.970, 0.974, 0.975, 0.975, 0.978, & 187 0.965, 0.965, 0.965, 0.965, 0.965, 0.969, 0.972, 0.977, 0.980, 0.980, 0.981, 0.983, & 188 0.969, 0.969, 0.969, 0.969, 0.969, 0.972, 0.975, 0.979, 0.983, 0.983, 0.983, 0.986, & 189 0.953, 0.953, 0.953, 0.953, 0.953, 0.958, 0.962, 0.969, 0.974, 0.975, 0.975, 0.979, & 190 0.892, 0.892, 0.892, 0.892, 0.892, 0.902, 0.912, 0.928, 0.940, 0.942, 0.943, 0.952, & 191 0.826, 0.826, 0.826, 0.826, 0.826, 0.658, 0.568, 0.479, 0.444, 0.442, 0.439, 0.426, & 192 ! Sulfate Coarse Soluble (CS) 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, & 198 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 199 !-- Sulfate Accumulation (BC content=0) 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, & 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 206 ! Seasalt Super Coarse Soluble (SS) 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 210 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 211 0.983, 0.985, 0.986, 0.987, 0.987, 0.988, 0.988, 0.989, 0.989, 0.989, 0.989, 0.989, & 212 0.957, 0.792, 0.765, 0.740, 0.721, 0.704, 0.687, 0.670, 0.652, 0.643, 0.631, 0.617, & 213 ! Seasalt Coarse Soluble (CS) 214 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.001, 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 218 0.994, 0.994, 0.995, 0.995, 0.995, 0.995, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, & 219 0.976, 0.867, 0.837, 0.814, 0.796, 0.774, 0.754, 0.735, 0.713, 0.702, 0.690, 0.675, & 220 ! Seasalt Accumulation Soluble (AS) 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 225 0.995, 0.996, 0.996, 0.996, 0.997, 0.997, 0.997, 0.998, 0.998, 0.998, 0.998, 0.999, & 226 0.919, 0.607, 0.553, 0.515, 0.480, 0.455, 0.434, 0.421, 0.421, 0.429, 0.449, 0.492, & 227 ! Nitrate Coarse Soluble (CS) 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 232 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.998, 0.998, 0.998, 0.998, 0.998, 0.998, & 233 0.707, 0.707, 0.707, 0.705, 0.702, 0.699, 0.697, 0.692, 0.686, 0.683, 0.679, 0.678, & 234 ! Nitrate Accumulation Soluble (AS) 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 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 239 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, & 240 0.362, 0.362, 0.362, 0.362, 0.363, 0.365, 0.367, 0.375, 0.390, 0.406, 0.428, 0.434 / 241 242 DATA cg_aers_6bands/ & 243 ! BC Accumulation Soluble (AS) 244 0.820, 0.820, 0.820, 0.820, 0.820, 0.831, 0.840, 0.852, 0.859, 0.860, 0.860, 0.864, & 245 0.760, 0.760, 0.760, 0.760, 0.760, 0.777, 0.791, 0.812, 0.826, 0.827, 0.828, 0.837, & 246 0.695, 0.695, 0.695, 0.695, 0.695, 0.716, 0.734, 0.762, 0.780, 0.782, 0.784, 0.797, & 247 0.603, 0.603, 0.603, 0.603, 0.603, 0.628, 0.649, 0.685, 0.710, 0.712, 0.715, 0.734, & 248 0.469, 0.469, 0.469, 0.469, 0.469, 0.494, 0.518, 0.559, 0.590, 0.593, 0.597, 0.622, & 249 0.276, 0.276, 0.276, 0.276, 0.276, 0.292, 0.308, 0.337, 0.361, 0.364, 0.367, 0.389, & 250 ! POM Accumulation Soluble (AS) 251 0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.711, 0.729, 0.740, 0.741, 0.742, 0.750, & 252 0.685, 0.685, 0.685, 0.685, 0.685, 0.699, 0.712, 0.731, 0.744, 0.745, 0.747, 0.756, & 253 0.673, 0.673, 0.673, 0.673, 0.673, 0.689, 0.702, 0.724, 0.739, 0.741, 0.742, 0.754, & 254 0.637, 0.637, 0.637, 0.637, 0.637, 0.653, 0.667, 0.691, 0.709, 0.711, 0.712, 0.726, & 255 0.553, 0.553, 0.553, 0.553, 0.553, 0.568, 0.582, 0.607, 0.628, 0.630, 0.632, 0.649, & 256 0.340, 0.340, 0.340, 0.340, 0.340, 0.350, 0.361, 0.382, 0.400, 0.402, 0.404, 0.422, & 257 ! Sulfate Coarse Soluble (CS) 258 0.746, 0.746, 0.746, 0.746, 0.779, 0.785, 0.790, 0.796, 0.804, 0.810, 0.816, 0.830, & 259 0.730, 0.730, 0.730, 0.730, 0.764, 0.770, 0.775, 0.782, 0.791, 0.797, 0.806, 0.819, & 260 0.707, 0.707, 0.707, 0.707, 0.746, 0.752, 0.758, 0.765, 0.775, 0.781, 0.790, 0.805, & 261 0.691, 0.691, 0.691, 0.691, 0.736, 0.741, 0.747, 0.754, 0.763, 0.769, 0.778, 0.791, & 262 0.696, 0.696, 0.696, 0.696, 0.748, 0.753, 0.759, 0.767, 0.776, 0.781, 0.789, 0.800, & 263 0.711, 0.711, 0.711, 0.711, 0.766, 0.773, 0.782, 0.791, 0.803, 0.811, 0.822, 0.838, & 264 !-- Sulfate Accumulation (BC content=0) 265 0.674, 0.674, 0.674, 0.674, 0.724, 0.729, 0.735, 0.741, 0.749, 0.753, 0.760, 0.769, & 266 0.677, 0.677, 0.677, 0.677, 0.729, 0.734, 0.740, 0.747, 0.756, 0.761, 0.768, 0.777, & 267 0.671, 0.671, 0.671, 0.671, 0.723, 0.730, 0.737, 0.746, 0.757, 0.764, 0.773, 0.786, & 268 0.638, 0.638, 0.638, 0.638, 0.691, 0.700, 0.710, 0.721, 0.736, 0.746, 0.760, 0.781, & 269 0.553, 0.553, 0.553, 0.553, 0.604, 0.616, 0.629, 0.644, 0.665, 0.680, 0.700, 0.736, & 270 0.344, 0.344, 0.344, 0.344, 0.388, 0.401, 0.416, 0.435, 0.462, 0.482, 0.512, 0.571, & 271 ! Seasalt Super Coarse Soluble (SS) 272 0.795, 0.810, 0.814, 0.819, 0.824, 0.829, 0.832, 0.840, 0.846, 0.849, 0.854, 0.860, & 273 0.788, 0.804, 0.808, 0.813, 0.818, 0.823, 0.828, 0.834, 0.842, 0.847, 0.851, 0.858, & 274 0.775, 0.792, 0.797, 0.803, 0.808, 0.813, 0.820, 0.826, 0.835, 0.840, 0.846, 0.854, & 275 0.754, 0.773, 0.779, 0.785, 0.791, 0.796, 0.803, 0.811, 0.821, 0.827, 0.834, 0.844, & 276 0.735, 0.755, 0.761, 0.768, 0.774, 0.779, 0.786, 0.794, 0.804, 0.810, 0.817, 0.829, & 277 0.712, 0.782, 0.795, 0.808, 0.819, 0.829, 0.840, 0.852, 0.865, 0.873, 0.882, 0.895, & 278 ! Seasalt Coarse Soluble (CS) 279 0.754, 0.770, 0.776, 0.781, 0.784, 0.791, 0.797, 0.805, 0.815, 0.822, 0.828, 0.840, & 280 0.736, 0.753, 0.759, 0.765, 0.771, 0.778, 0.785, 0.793, 0.804, 0.811, 0.820, 0.831, & 281 0.716, 0.735, 0.742, 0.748, 0.754, 0.762, 0.769, 0.778, 0.789, 0.796, 0.807, 0.819, & 282 0.704, 0.725, 0.733, 0.739, 0.745, 0.752, 0.759, 0.768, 0.778, 0.784, 0.792, 0.803, & 283 0.716, 0.737, 0.744, 0.751, 0.756, 0.763, 0.770, 0.777, 0.786, 0.790, 0.795, 0.800, & 284 0.688, 0.730, 0.741, 0.751, 0.761, 0.771, 0.782, 0.795, 0.810, 0.820, 0.833, 0.849, & 285 ! Seasalt Accumulation Soluble (AS) 286 0.686, 0.711, 0.719, 0.726, 0.733, 0.740, 0.748, 0.755, 0.763, 0.767, 0.770, 0.774, & 287 0.690, 0.716, 0.724, 0.731, 0.739, 0.747, 0.756, 0.764, 0.774, 0.778, 0.782, 0.784, & 288 0.683, 0.710, 0.719, 0.727, 0.735, 0.744, 0.754, 0.765, 0.777, 0.784, 0.791, 0.797, & 289 0.649, 0.676, 0.686, 0.695, 0.704, 0.714, 0.726, 0.739, 0.757, 0.768, 0.781, 0.797, & 290 0.559, 0.586, 0.596, 0.606, 0.616, 0.629, 0.643, 0.661, 0.686, 0.703, 0.726, 0.758, & 291 0.341, 0.358, 0.366, 0.375, 0.382, 0.395, 0.408, 0.428, 0.459, 0.484, 0.520, 0.581, & 292 ! Nitrate Coarse Soluble (CS) 293 0.765, 0.765, 0.765, 0.773, 0.779, 0.785, 0.791, 0.803, 0.813, 0.819, 0.827, 0.828, & 294 0.754, 0.754, 0.754, 0.762, 0.769, 0.774, 0.781, 0.790, 0.802, 0.810, 0.817, 0.818, & 295 0.743, 0.743, 0.743, 0.751, 0.758, 0.763, 0.770, 0.779, 0.789, 0.797, 0.804, 0.806, & 296 0.737, 0.737, 0.737, 0.745, 0.752, 0.758, 0.764, 0.773, 0.781, 0.786, 0.792, 0.793, & 297 0.739, 0.739, 0.739, 0.749, 0.757, 0.763, 0.770, 0.779, 0.788, 0.793, 0.797, 0.798, & 298 0.651, 0.651, 0.651, 0.672, 0.690, 0.707, 0.724, 0.751, 0.779, 0.796, 0.813, 0.817, & 299 ! Nitrate Accumulation Soluble (AS) 300 0.713, 0.713, 0.713, 0.723, 0.732, 0.739, 0.746, 0.756, 0.765, 0.769, 0.772, 0.772, & 301 0.724, 0.724, 0.724, 0.734, 0.742, 0.750, 0.757, 0.767, 0.777, 0.782, 0.785, 0.785, & 302 0.719, 0.719, 0.719, 0.729, 0.738, 0.746, 0.754, 0.767, 0.780, 0.787, 0.793, 0.794, & 303 0.673, 0.673, 0.673, 0.685, 0.696, 0.706, 0.717, 0.734, 0.753, 0.766, 0.778, 0.781, & 304 0.565, 0.565, 0.565, 0.580, 0.593, 0.605, 0.619, 0.644, 0.673, 0.693, 0.716, 0.720, & 305 0.328, 0.328, 0.328, 0.338, 0.348, 0.359, 0.371, 0.396, 0.431, 0.459, 0.494, 0.502 / 163 306 164 307 DATA alpha_aeri_6bands/ & 165 ! dust insoluble CI 166 0.751, 0.751, 0.769, 0.772, 0.672, 0.437, & 167 ! bc insoluble AI 168 6.497, 6.497, 5.900, 4.284, 2.163, 0.966, & 169 ! pom insoluble AI 170 6.443, 6.443, 4.381, 1.846, 0.377, 0.052 / 171 172 DATA cg_aers_6bands/ & 173 ! bc soluble AS 174 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 175 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 176 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 177 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 178 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 179 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, & 180 ! pom soluble AS 181 0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.710, 0.726, 0.736, 0.737, 0.738, 0.745, & 182 0.687, 0.687, 0.687, 0.687, 0.687, 0.700, 0.710, 0.726, 0.736, 0.737, 0.738, 0.745, & 183 0.658, 0.658, 0.658, 0.658, 0.658, 0.667, 0.674, 0.685, 0.692, 0.692, 0.693, 0.698, & 184 0.564, 0.564, 0.564, 0.564, 0.564, 0.566, 0.568, 0.571, 0.573, 0.573, 0.573, 0.574, & 185 0.363, 0.363, 0.363, 0.363, 0.363, 0.360, 0.357, 0.352, 0.350, 0.349, 0.349, 0.347, & 186 0.142, 0.142, 0.142, 0.142, 0.142, 0.139, 0.137, 0.133, 0.131, 0.131, 0.130, 0.129, & 187 ! sulfate AS 188 0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, & 189 0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, & 190 0.653, 0.653, 0.653, 0.662, 0.670, 0.676, 0.683, 0.692, 0.701, 0.706, 0.709, 0.716, & 191 0.563, 0.563, 0.563, 0.565, 0.567, 0.569, 0.570, 0.573, 0.575, 0.576, 0.577, 0.579, & 192 0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, & 193 0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, & 194 ! sulfate coarse CS 195 0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, & 196 0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, & 197 0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, & 198 0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, & 199 0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, & 200 0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, & 201 ! seasalt seasalt Super Coarse Soluble (SS) 202 0.797, 0.800, 0.801, 0.802, 0.804, 0.822, 0.825, 0.828, 0.832, 0.835, 0.838, 0.843, & 203 0.788, 0.792, 0.794, 0.795, 0.796, 0.815, 0.818, 0.822, 0.827, 0.829, 0.833, 0.838, & 204 0.773, 0.778, 0.780, 0.782, 0.783, 0.802, 0.806, 0.811, 0.817, 0.820, 0.825, 0.832, & 205 0.746, 0.753, 0.755, 0.759, 0.760, 0.781, 0.787, 0.792, 0.800, 0.805, 0.811, 0.820, & 206 0.706, 0.714, 0.716, 0.720, 0.722, 0.749, 0.753, 0.761, 0.769, 0.774, 0.783, 0.797, & 207 0.681, 0.682, 0.682, 0.683, 0.684, 0.723, 0.727, 0.732, 0.738, 0.741, 0.748, 0.757, & 208 ! seasalt seasalt Coarse Soluble (CS) 209 0.756, 0.761, 0.764, 0.766, 0.769, 0.790, 0.793, 0.799, 0.805, 0.810, 0.815, 0.823, & 210 0.736, 0.743, 0.747, 0.749, 0.751, 0.773, 0.778, 0.784, 0.793, 0.797, 0.804, 0.815, & 211 0.712, 0.719, 0.721, 0.725, 0.726, 0.752, 0.758, 0.764, 0.773, 0.779, 0.786, 0.800, & 212 0.690, 0.694, 0.695, 0.698, 0.699, 0.731, 0.738, 0.742, 0.751, 0.756, 0.764, 0.776, & 213 0.682, 0.683, 0.683, 0.683, 0.684, 0.725, 0.729, 0.733, 0.737, 0.740, 0.744, 0.752, & 214 0.669, 0.673, 0.674, 0.675, 0.676, 0.718, 0.724, 0.730, 0.736, 0.739, 0.742, 0.746, & 215 ! seasalt seasalt Accumulation Soluble (AS) 216 0.694, 0.692, 0.692, 0.691, 0.689, 0.737, 0.740, 0.742, 0.742, 0.740, 0.737, 0.731, & 217 0.685, 0.690, 0.691, 0.692, 0.692, 0.735, 0.741, 0.746, 0.750, 0.751, 0.750, 0.744, & 218 0.650, 0.662, 0.666, 0.670, 0.673, 0.710, 0.719, 0.729, 0.741, 0.747, 0.754, 0.757, & 219 0.561, 0.585, 0.593, 0.601, 0.609, 0.637, 0.651, 0.669, 0.691, 0.705, 0.723, 0.745, & 220 0.392, 0.427, 0.439, 0.451, 0.464, 0.480, 0.500, 0.526, 0.563, 0.588, 0.621, 0.671, & 221 0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439 / 308 ! Dust insoluble 309 0.735, 0.765, 0.817, 0.846, 0.718, 0.414, & 310 ! BC insoluble 311 5.276, 5.371, 5.126, 4.225, 2.667, 1.113, & 312 ! POM insoluble 313 5.962, 5.609, 4.527, 2.697, 0.918, 0.105, & 314 ! Nitrate insoluble 315 0.688, 0.712, 0.754, 0.809, 0.786, 0.703 / 316 317 DATA piz_aeri_6bands/ & 318 ! Dust insoluble 319 0.894, 0.902, 0.937, 0.980, 0.993, 0.967, & 320 ! BC insoluble 321 0.820, 0.760, 0.695, 0.603, 0.469, 0.276, & 322 ! POM insoluble 323 0.687, 0.685, 0.673, 0.637, 0.553, 0.340, & 324 ! Nitrate insoluble 325 1.000, 1.000, 1.000, 1.000, 0.999, 0.707 / 222 326 223 327 DATA cg_aeri_6bands/ & 224 ! dust insoluble CI 225 0.718, 0.718, 0.699, 0.661, 0.676, 0.670, & 226 ! bc insoluble AI 227 0.721, 0.721, 0.643, 0.513, 0.321, 0.153, & 228 ! pom insoluble AI 229 0.687, 0.687, 0.658, 0.564, 0.363, 0.142 / 230 231 DATA piz_aers_6bands/& 232 ! bc soluble AS 233 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 234 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 235 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 236 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 237 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 238 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, & 239 ! pom soluble AS 240 0.973, 0.973, 0.973, 0.973, 0.973, 0.977, 0.980, 0.984, 0.987, 0.988, 0.988, 0.990, & 241 0.973, 0.973, 0.973, 0.973, 0.973, 0.977, 0.980, 0.984, 0.987, 0.988, 0.988, 0.990, & 242 0.972, 0.972, 0.972, 0.972, 0.972, 0.976, 0.979, 0.984, 0.987, 0.987, 0.988, 0.990, & 243 0.940, 0.940, 0.940, 0.940, 0.940, 0.948, 0.955, 0.965, 0.972, 0.973, 0.973, 0.978, & 244 0.816, 0.816, 0.816, 0.816, 0.816, 0.839, 0.859, 0.888, 0.908, 0.910, 0.911, 0.925, & 245 0.663, 0.663, 0.663, 0.663, 0.663, 0.607, 0.562, 0.492, 0.446, 0.441, 0.437, 0.404, & 246 ! sulfate AS 247 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 248 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 249 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 250 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 251 0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, & 252 0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, & 253 ! sulfate coarse CS 254 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 255 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 256 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 257 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 258 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, & 259 0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, & 260 ! seasalt seasalt Super Coarse Soluble (SS) 261 0.983, 0.982, 0.982, 0.982, 0.982, 0.992, 0.993, 0.994, 0.995, 0.996, 0.997, 0.998, & 262 0.984, 0.983, 0.983, 0.983, 0.983, 0.993, 0.994, 0.995, 0.996, 0.996, 0.997, 0.998, & 263 0.986, 0.985, 0.985, 0.985, 0.984, 0.993, 0.994, 0.995, 0.996, 0.997, 0.997, 0.998, & 264 0.989, 0.988, 0.988, 0.988, 0.987, 0.995, 0.996, 0.996, 0.997, 0.997, 0.998, 0.999, & 265 0.994, 0.993, 0.992, 0.992, 0.992, 0.997, 0.997, 0.998, 0.998, 0.998, 0.999, 0.999, & 266 0.997, 0.997, 0.997, 0.996, 0.996, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, & 267 ! seasalt seasalt Coarse Soluble (CS) 268 0.988, 0.987, 0.987, 0.987, 0.986, 0.995, 0.995, 0.996, 0.997, 0.997, 0.998, 0.999, & 269 0.990, 0.989, 0.989, 0.989, 0.989, 0.995, 0.996, 0.997, 0.997, 0.998, 0.998, 0.999, & 270 0.993, 0.992, 0.992, 0.992, 0.991, 0.997, 0.997, 0.997, 0.998, 0.998, 0.998, 0.999, & 271 0.996, 0.995, 0.995, 0.995, 0.994, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 0.999, & 272 0.998, 0.997, 0.997, 0.997, 0.997, 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, & 273 0.999, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, & 274 ! seasalt seasalt Accumulation Soluble (AS) 275 0.999, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, & 276 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 277 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 278 0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, & 279 0.997, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, & 280 0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000 / 281 282 DATA piz_aeri_6bands/ & 283 ! dust insoluble CI 284 0.954, 0.954, 0.965, 0.981, 0.996, 0.990, & 285 ! bc insoluble AI 286 0.460, 0.460, 0.445, 0.394, 0.267, 0.121, & 287 ! pom insoluble AI 288 0.973, 0.973, 0.972, 0.940, 0.816, 0.663 / 328 ! Dust insoluble 329 0.756, 0.737, 0.713, 0.701, 0.696, 0.672, & 330 ! BC insoluble 331 0.480, 0.469, 0.456, 0.427, 0.349, 0.196, & 332 ! POM insoluble 333 0.956, 0.965, 0.969, 0.953, 0.892, 0.826, & 334 ! Nitrate insoluble 335 0.765, 0.754, 0.743, 0.737, 0.739, 0.651 / 289 336 290 337 ! Added by R. Wang (July 31 2016) 291 338 ! properties for BC assuming Maxwell-Garnett rule and internal mixture 292 339 293 DATA alpha_MG_6bands/ & 294 ! bc content = 0.001 295 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 296 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 297 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 298 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 299 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 300 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, & 301 ! bc content = 0.01 302 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 303 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 304 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 305 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 306 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 307 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, & 308 ! bc content = 0.02 309 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 310 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 311 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 312 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 313 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 314 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, & 315 ! bc content = 0.05 316 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 317 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 318 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 319 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 320 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 321 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894, & 322 ! bc content = 0.10 323 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 324 6.497, 6.497, 6.497, 6.497, 6.497, 7.160, 7.875, 9.356,10.811,10.974,11.149,12.734, & 325 5.900, 5.900, 5.900, 5.900, 5.900, 6.502, 7.151, 8.496, 9.818, 9.965,10.124,11.564, & 326 4.284, 4.284, 4.284, 4.284, 4.284, 4.721, 5.193, 6.169, 7.129, 7.236, 7.352, 8.397, & 327 2.163, 2.163, 2.163, 2.163, 2.163, 2.384, 2.622, 3.115, 3.600, 3.654, 3.712, 4.240, & 328 0.966, 0.966, 0.966, 0.966, 0.966, 1.065, 1.171, 1.392, 1.608, 1.632, 1.658, 1.894 / 329 330 DATA cg_MG_6bands/ & 331 ! bc content = 0.001 332 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 333 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 334 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 335 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 336 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 337 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, & 338 ! bc content = 0.01 339 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 340 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 341 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 342 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 343 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 344 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, & 345 ! bc content = 0.02 346 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 347 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 348 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 349 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 350 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 351 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, & 352 ! bc content = 0.05 353 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 354 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 355 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 356 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 357 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 358 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123, & 359 ! bc content = 0.10 360 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 361 0.721, 0.721, 0.721, 0.729, 0.735, 0.741, 0.746, 0.754, 0.762, 0.766, 0.769, 0.775, & 362 0.643, 0.643, 0.643, 0.654, 0.662, 0.670, 0.677, 0.688, 0.698, 0.704, 0.707, 0.715, & 363 0.513, 0.513, 0.513, 0.522, 0.530, 0.536, 0.542, 0.552, 0.560, 0.565, 0.568, 0.575, & 364 0.321, 0.321, 0.321, 0.323, 0.325, 0.327, 0.328, 0.331, 0.333, 0.334, 0.335, 0.337, & 365 0.153, 0.153, 0.153, 0.149, 0.145, 0.142, 0.139, 0.135, 0.130, 0.128, 0.127, 0.123 / 366 367 DATA piz_MG_6bands/ & 368 ! bc content = 0.001 369 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 370 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 371 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 372 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 373 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 374 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, & 375 ! bc content = 0.01 376 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 377 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 378 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 379 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 380 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 381 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, & 382 ! bc content = 0.02 383 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 384 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 385 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 386 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 387 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 388 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, & 389 ! bc content = 0.05 390 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 391 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 392 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 393 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 394 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 395 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207, & 396 ! bc content = 0.10 397 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 398 0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, & 399 0.445, 0.445, 0.445, 0.445, 0.445, 0.521, 0.583, 0.679, 0.741, 0.747, 0.753, 0.798, & 400 0.394, 0.394, 0.394, 0.394, 0.394, 0.477, 0.545, 0.649, 0.718, 0.724, 0.730, 0.779, & 401 0.267, 0.267, 0.267, 0.267, 0.267, 0.365, 0.446, 0.571, 0.652, 0.660, 0.667, 0.725, & 402 0.121, 0.121, 0.121, 0.121, 0.121, 0.139, 0.155, 0.178, 0.193, 0.195, 0.196, 0.207 / 340 DATA alpha_MG_6bands/ & 341 !--BC content=0.001 342 3.546, 3.546, 3.546, 3.546, 3.640, 3.891, 3.883, 4.354, 4.967, 5.424, 6.412, 9.049, & 343 3.995, 3.995, 3.995, 3.995, 4.073, 4.091, 4.111, 4.209, 4.502, 4.878, 5.591, 8.029, & 344 4.650, 4.650, 4.650, 4.650, 4.661, 4.561, 4.472, 4.386, 4.373, 4.451, 4.726, 6.248, & 345 4.445, 4.445, 4.445, 4.445, 4.393, 4.291, 4.171, 4.017, 3.812, 3.682, 3.514, 3.593, & 346 2.837, 2.837, 2.837, 2.837, 2.824, 2.798, 2.760, 2.696, 2.579, 2.471, 2.265, 1.735, & 347 1.087, 1.087, 1.087, 1.087, 1.095, 1.101, 1.107, 1.110, 1.109, 1.101, 1.074, 0.944, & 348 !--BC content=0.010 349 3.550, 3.550, 3.550, 3.550, 3.623, 3.852, 4.006, 4.342, 4.945, 5.381, 6.394, 8.972, & 350 3.995, 3.995, 3.995, 3.995, 4.082, 4.094, 4.133, 4.235, 4.517, 4.871, 5.601, 8.001, & 351 4.653, 4.653, 4.653, 4.653, 4.662, 4.564, 4.475, 4.393, 4.376, 4.447, 4.736, 6.221, & 352 4.434, 4.434, 4.434, 4.434, 4.386, 4.286, 4.170, 4.021, 3.815, 3.680, 3.524, 3.583, & 353 2.831, 2.831, 2.831, 2.831, 2.820, 2.795, 2.757, 2.695, 2.579, 2.472, 2.270, 1.734, & 354 1.086, 1.086, 1.086, 1.086, 1.095, 1.101, 1.106, 1.110, 1.108, 1.101, 1.075, 0.944, & 355 !--BC content=0.020 356 3.549, 3.549, 3.549, 3.549, 3.630, 3.842, 4.016, 4.335, 4.928, 5.383, 6.369, 8.937, & 357 4.004, 4.004, 4.004, 4.004, 4.096, 4.100, 4.145, 4.243, 4.529, 4.870, 5.595, 7.973, & 358 4.650, 4.650, 4.650, 4.650, 4.662, 4.564, 4.478, 4.397, 4.379, 4.449, 4.734, 6.207, & 359 4.422, 4.422, 4.422, 4.422, 4.378, 4.280, 4.167, 4.020, 3.817, 3.683, 3.527, 3.583, & 360 2.824, 2.824, 2.824, 2.824, 2.815, 2.791, 2.754, 2.693, 2.579, 2.472, 2.272, 1.740, & 361 1.085, 1.085, 1.085, 1.085, 1.094, 1.100, 1.105, 1.109, 1.108, 1.100, 1.075, 0.945, & 362 !--BC content=0.050 363 3.563, 3.563, 3.563, 3.563, 3.662, 3.834, 4.019, 4.318, 4.889, 5.360, 6.305, 8.842, & 364 4.039, 4.039, 4.039, 4.039, 4.130, 4.126, 4.166, 4.261, 4.544, 4.863, 5.573, 7.895, & 365 4.641, 4.641, 4.641, 4.641, 4.659, 4.565, 4.483, 4.407, 4.389, 4.455, 4.731, 6.161, & 366 4.388, 4.388, 4.388, 4.388, 4.356, 4.264, 4.157, 4.017, 3.821, 3.691, 3.537, 3.582, & 367 2.805, 2.805, 2.805, 2.805, 2.801, 2.779, 2.744, 2.687, 2.577, 2.474, 2.279, 1.757, & 368 1.081, 1.081, 1.081, 1.081, 1.091, 1.098, 1.103, 1.107, 1.106, 1.099, 1.074, 0.949, & 369 !--BC content=0.100 370 3.607, 3.607, 3.607, 3.607, 3.701, 3.843, 4.021, 4.302, 4.841, 5.306, 6.214, 8.689, & 371 4.098, 4.098, 4.098, 4.098, 4.175, 4.168, 4.201, 4.291, 4.554, 4.854, 5.529, 7.763, & 372 4.627, 4.627, 4.627, 4.627, 4.651, 4.567, 4.490, 4.421, 4.403, 4.464, 4.723, 6.082, & 373 4.333, 4.333, 4.333, 4.333, 4.319, 4.236, 4.139, 4.011, 3.828, 3.704, 3.554, 3.580, & 374 2.773, 2.773, 2.773, 2.773, 2.778, 2.758, 2.727, 2.675, 2.573, 2.476, 2.291, 1.786, & 375 1.076, 1.076, 1.076, 1.076, 1.087, 1.094, 1.099, 1.103, 1.103, 1.097, 1.073, 0.954, & 376 !--BC content=0.200 377 3.702, 3.702, 3.702, 3.702, 3.768, 3.878, 4.030, 4.278, 4.757, 5.190, 6.032, 8.370, & 378 4.204, 4.204, 4.204, 4.204, 4.255, 4.246, 4.270, 4.345, 4.571, 4.837, 5.440, 7.494, & 379 4.598, 4.598, 4.598, 4.598, 4.633, 4.565, 4.502, 4.445, 4.428, 4.481, 4.706, 5.923, & 380 4.230, 4.230, 4.230, 4.230, 4.244, 4.179, 4.101, 3.996, 3.838, 3.728, 3.588, 3.581, & 381 2.710, 2.710, 2.710, 2.710, 2.730, 2.716, 2.693, 2.651, 2.564, 2.480, 2.314, 1.846, & 382 1.064, 1.064, 1.064, 1.064, 1.078, 1.085, 1.091, 1.096, 1.097, 1.091, 1.072, 0.965 / 383 384 DATA piz_MG_6bands/ & 385 !--BC content=0.001 386 -3.984, -3.984, -3.984, -3.984, -3.395, -3.075, -3.051, -2.567, -2.053, -1.807, -1.308, -0.592, & 387 -3.014, -3.014, -3.014, -3.014, -2.520, -2.473, -2.411, -2.294, -2.027, -1.771, -1.368, -0.613, & 388 -1.142, -1.142, -1.142, -1.142, -0.930, -0.961, -0.981, -1.001, -0.981, -0.934, -0.798, -0.337, & 389 -0.202, -0.202, -0.202, -0.202, -0.129, -0.154, -0.183, -0.224, -0.284, -0.326, -0.383, -0.347, & 390 0.133, 0.133, 0.133, 0.133, 0.158, 0.146, 0.131, 0.107, 0.061, 0.015, -0.081, -0.433, & 391 0.137, 0.137, 0.137, 0.137, 0.148, 0.149, 0.149, 0.147, 0.139, 0.128, 0.098, -0.048, & 392 !--BC content=0.010 393 -3.640, -3.640, -3.640, -3.640, -3.219, -2.942, -2.761, -2.445, -1.983, -1.731, -1.262, -0.576, & 394 -2.791, -2.791, -2.791, -2.791, -2.383, -2.352, -2.290, -2.180, -1.943, -1.707, -1.321, -0.595, & 395 -1.069, -1.069, -1.069, -1.069, -0.889, -0.920, -0.945, -0.966, -0.954, -0.910, -0.775, -0.333, & 396 -0.185, -0.185, -0.185, -0.185, -0.119, -0.144, -0.173, -0.213, -0.275, -0.318, -0.372, -0.346, & 397 0.138, 0.138, 0.138, 0.138, 0.160, 0.149, 0.134, 0.110, 0.063, 0.018, -0.076, -0.430, & 398 0.137, 0.137, 0.137, 0.137, 0.148, 0.149, 0.150, 0.148, 0.140, 0.129, 0.099, -0.046, & 399 !--BC content=0.020 400 -3.369, -3.369, -3.369, -3.369, -3.033, -2.799, -2.619, -2.339, -1.913, -1.657, -1.223, -0.558, & 401 -2.592, -2.592, -2.592, -2.592, -2.254, -2.239, -2.185, -2.091, -1.869, -1.652, -1.284, -0.580, & 402 -1.004, -1.004, -1.004, -1.004, -0.849, -0.882, -0.909, -0.932, -0.926, -0.885, -0.758, -0.326, & 403 -0.168, -0.168, -0.168, -0.168, -0.108, -0.133, -0.163, -0.203, -0.265, -0.309, -0.364, -0.341, & 404 0.142, 0.142, 0.142, 0.142, 0.163, 0.152, 0.137, 0.113, 0.066, 0.022, -0.072, -0.422, & 405 0.138, 0.138, 0.138, 0.138, 0.149, 0.150, 0.150, 0.148, 0.140, 0.129, 0.100, -0.044, & 406 !--BC content=0.050 407 -2.743, -2.743, -2.743, -2.743, -2.582, -2.438, -2.292, -2.078, -1.729, -1.495, -1.123, -0.512, & 408 -2.121, -2.121, -2.121, -2.121, -1.939, -1.951, -1.926, -1.865, -1.688, -1.511, -1.188, -0.541, & 409 -0.839, -0.839, -0.839, -0.839, -0.742, -0.780, -0.812, -0.842, -0.848, -0.818, -0.707, -0.307, & 410 -0.121, -0.121, -0.121, -0.121, -0.078, -0.103, -0.133, -0.174, -0.236, -0.281, -0.337, -0.325, & 411 0.154, 0.154, 0.154, 0.154, 0.171, 0.160, 0.145, 0.121, 0.075, 0.032, -0.060, -0.399, & 412 0.141, 0.141, 0.141, 0.141, 0.150, 0.151, 0.152, 0.150, 0.142, 0.131, 0.103, -0.038, & 413 !--BC content=0.100 414 -2.025, -2.025, -2.025, -2.025, -2.035, -1.965, -1.875, -1.731, -1.475, -1.285, -0.980, -0.447, & 415 -1.569, -1.569, -1.569, -1.569, -1.540, -1.572, -1.580, -1.557, -1.444, -1.312, -1.052, -0.485, & 416 -0.628, -0.628, -0.628, -0.628, -0.593, -0.634, -0.672, -0.710, -0.732, -0.716, -0.631, -0.279, & 417 -0.054, -0.054, -0.054, -0.054, -0.032, -0.058, -0.088, -0.129, -0.191, -0.236, -0.295, -0.299, & 418 0.173, 0.173, 0.173, 0.173, 0.184, 0.173, 0.158, 0.135, 0.090, 0.048, -0.039, -0.361, & 419 0.145, 0.145, 0.145, 0.145, 0.153, 0.154, 0.154, 0.152, 0.144, 0.134, 0.107, -0.027, & 420 !--BC content=0.200 421 -1.155, -1.155, -1.155, -1.155, -1.300, -1.295, -1.272, -1.213, -1.077, -0.957, -0.747, -0.336, & 422 -0.884, -0.884, -0.884, -0.884, -0.979, -1.028, -1.065, -1.086, -1.054, -0.986, -0.820, -0.388, & 423 -0.332, -0.332, -0.332, -0.332, -0.360, -0.402, -0.444, -0.490, -0.531, -0.536, -0.492, -0.225, & 424 0.050, 0.050, 0.050, 0.050, 0.048, 0.022, -0.007, -0.047, -0.107, -0.152, -0.213, -0.245, & 425 0.204, 0.204, 0.204, 0.204, 0.208, 0.197, 0.182, 0.161, 0.120, 0.081, 0.002, -0.285, & 426 0.153, 0.153, 0.153, 0.153, 0.159, 0.159, 0.159, 0.157, 0.150, 0.140, 0.115, -0.005 / 427 428 DATA cg_MG_6bands/ & 429 !--BC content=0.001 430 0.290, 0.290, 0.290, 0.290, 0.378, 0.381, 0.377, 0.378, 0.373, 0.344, 0.317, 0.118, & 431 0.374, 0.374, 0.374, 0.374, 0.469, 0.474, 0.480, 0.480, 0.475, 0.462, 0.438, 0.242, & 432 0.375, 0.375, 0.375, 0.375, 0.491, 0.502, 0.517, 0.527, 0.534, 0.529, 0.508, 0.255, & 433 0.242, 0.242, 0.242, 0.242, 0.455, 0.516, 0.577, 0.631, 0.687, 0.715, 0.748, 0.772, & 434 0.559, 0.559, 0.559, 0.559, 0.487, 0.461, 0.410, 0.303, -0.181, -3.594, 1.814, 1.164, & 435 0.323, 0.323, 0.323, 0.323, 0.339, 0.339, 0.337, 0.328, 0.295, 0.243, 0.045, 2.778, & 436 !--BC content=0.010 437 0.300, 0.300, 0.300, 0.300, 0.385, 0.389, 0.390, 0.386, 0.377, 0.353, 0.319, 0.116, & 438 0.384, 0.384, 0.384, 0.384, 0.478, 0.482, 0.488, 0.488, 0.483, 0.471, 0.442, 0.246, & 439 0.380, 0.380, 0.380, 0.380, 0.496, 0.507, 0.520, 0.531, 0.538, 0.534, 0.512, 0.259, & 440 0.233, 0.233, 0.233, 0.233, 0.454, 0.520, 0.581, 0.636, 0.691, 0.719, 0.752, 0.778, & 441 0.553, 0.553, 0.553, 0.553, 0.487, 0.461, 0.414, 0.308, -0.144, -2.858, 1.889, 1.165, & 442 0.323, 0.323, 0.323, 0.323, 0.338, 0.339, 0.337, 0.328, 0.295, 0.244, 0.051, 2.847, & 443 !--BC content=0.020 444 0.305, 0.305, 0.305, 0.305, 0.391, 0.394, 0.395, 0.390, 0.380, 0.358, 0.321, 0.112, & 445 0.391, 0.391, 0.391, 0.391, 0.484, 0.488, 0.493, 0.493, 0.488, 0.476, 0.445, 0.246, & 446 0.384, 0.384, 0.384, 0.384, 0.499, 0.511, 0.524, 0.534, 0.541, 0.537, 0.514, 0.260, & 447 0.219, 0.219, 0.219, 0.219, 0.451, 0.523, 0.586, 0.642, 0.696, 0.724, 0.756, 0.783, & 448 0.547, 0.547, 0.547, 0.547, 0.487, 0.461, 0.416, 0.315, -0.108, -2.308, 1.954, 1.170, & 449 0.322, 0.322, 0.322, 0.322, 0.338, 0.339, 0.337, 0.328, 0.296, 0.245, 0.056, 2.949, & 450 !--BC content=0.050 451 0.312, 0.312, 0.312, 0.312, 0.400, 0.402, 0.403, 0.398, 0.385, 0.364, 0.321, 0.094, & 452 0.402, 0.402, 0.402, 0.402, 0.495, 0.498, 0.503, 0.503, 0.496, 0.484, 0.451, 0.242, & 453 0.388, 0.388, 0.388, 0.388, 0.506, 0.518, 0.530, 0.541, 0.548, 0.544, 0.521, 0.259, & 454 0.152, 0.152, 0.152, 0.152, 0.434, 0.527, 0.599, 0.658, 0.711, 0.738, 0.768, 0.794, & 455 0.534, 0.534, 0.534, 0.534, 0.486, 0.463, 0.423, 0.332, -0.018, -1.355, 2.203, 1.188, & 456 0.321, 0.321, 0.321, 0.321, 0.337, 0.338, 0.336, 0.328, 0.297, 0.249, 0.071, 3.311, & 457 !--BC content=0.100 458 0.309, 0.309, 0.309, 0.309, 0.404, 0.405, 0.406, 0.401, 0.385, 0.364, 0.316, 0.057, & 459 0.406, 0.406, 0.406, 0.406, 0.503, 0.507, 0.511, 0.512, 0.504, 0.491, 0.456, 0.232, & 460 0.382, 0.382, 0.382, 0.382, 0.510, 0.523, 0.537, 0.548, 0.555, 0.551, 0.528, 0.255, & 461 -0.178, -0.178, -0.178, -0.178, 0.312, 0.524, 0.626, 0.692, 0.741, 0.765, 0.791, 0.816, & 462 0.519, 0.519, 0.519, 0.519, 0.486, 0.466, 0.432, 0.357, 0.091, -0.650, 2.968, 1.221, & 463 0.319, 0.319, 0.319, 0.319, 0.336, 0.337, 0.335, 0.328, 0.300, 0.255, 0.094, 4.314, & 464 !--BC content=0.200 465 0.267, 0.267, 0.267, 0.267, 0.388, 0.391, 0.392, 0.387, 0.370, 0.345, 0.287, -0.053, & 466 0.381, 0.381, 0.381, 0.381, 0.499, 0.506, 0.512, 0.514, 0.507, 0.494, 0.455, 0.198, & 467 0.318, 0.318, 0.318, 0.318, 0.494, 0.515, 0.535, 0.551, 0.561, 0.559, 0.535, 0.237, & 468 1.151, 1.151, 1.151, 1.151, 0.703, 0.643, 1.248, 0.874, 0.846, 0.846, 0.853, 0.868, & 469 0.502, 0.502, 0.502, 0.502, 0.486, 0.472, 0.447, 0.393, 0.226, -0.108,-43.753, 1.313, & 470 0.315, 0.315, 0.315, 0.315, 0.332, 0.334, 0.333, 0.327, 0.304, 0.266, 0.135, 21.079 / 403 471 404 472 !----BEGINNING OF CALCULATIONS 473 474 modname='aeropt_6bands_rrtm' 405 475 406 476 spsol = 0 … … 408 478 409 479 IF (NSW.NE.nbands_sw_rrtm) THEN 410 print *,'Erreur NSW doit etre egal a 6 pour cette routine' 411 stop 480 CALL abort_physic(modname,'Erreur NSW doit etre egal a 6 pour cette routine',1) 412 481 ENDIF 413 482 … … 440 509 aerosol_name(1) = id_CIDUSTM_phy 441 510 ELSEIF (flag_aerosol .EQ. 6) THEN 442 nb_aer = 1 0511 nb_aer = 13 443 512 ALLOCATE (aerosol_name(nb_aer)) 444 513 aerosol_name(1) = id_ASSO4M_phy … … 452 521 aerosol_name(9) = id_CIDUSTM_phy 453 522 aerosol_name(10)= id_CSSO4M_phy 523 aerosol_name(11)= id_CSNO3M_phy 524 aerosol_name(12)= id_ASNO3M_phy 525 aerosol_name(13)= id_CINO3M_phy 454 526 ENDIF 455 527 … … 491 563 soluble=.TRUE. 492 564 spsol=2 493 ELSEIF (aerosol_name(m).EQ.id_ ASSO4M_phy) THEN565 ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN 494 566 soluble=.TRUE. 495 567 spsol=3 496 568 fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD 497 ELSEIF (aerosol_name(m).EQ.id_ CSSO4M_phy) THEN569 ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN 498 570 soluble=.TRUE. 499 571 spsol=4 … … 508 580 soluble=.TRUE. 509 581 spsol=7 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 510 590 ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN 511 591 soluble=.FALSE. … … 517 597 soluble=.FALSE. 518 598 spinsol=3 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 519 603 ELSE 520 604 CYCLE … … 522 606 523 607 !--shortname for aerosol index 524 id=aerosol_name(m)608 aerindex=aerosol_name(m) 525 609 526 610 IF (soluble) THEN … … 532 616 DO i=1, klon 533 617 534 tmp_var=m_allaer(i,k, spsol)/1.e6*zdh(i,k)*fac535 tmp_var_pi=m_allaer_pi(i,k, spsol)/1.e6*zdh(i,k)*fac618 tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac 619 tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac 536 620 537 621 ! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols 538 622 BC_massfra = m_allaer(i,k,id_ASBCM_phy)/(m_allaer(i,k,id_ASBCM_phy)+m_allaer(i,k,id_ASSO4M_phy)) 539 623 540 IF (BC_massfra.GE.0.10) THEN 624 IF (BC_massfra.GE.0.20) THEN 625 classbc = 6 626 ELSEIF (BC_massfra.GE.0.10) THEN 541 627 classbc = 5 542 628 ELSEIF (BC_massfra.GE.0.05) THEN … … 553 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)) 554 640 555 IF (BC_massfra_pi.GE.0.10) THEN 641 IF (BC_massfra_pi.GE.0.20) THEN 642 classbc_pi = 6 643 ELSEIF (BC_massfra_pi.GE.0.10) THEN 556 644 classbc_pi = 5 557 645 ELSEIF (BC_massfra_pi.GE.0.05) THEN … … 565 653 ENDIF 566 654 567 DO inu=1, NSW655 DO inu=1, NSW 568 656 569 657 !--all aerosols … … 580 668 cg_MG_6bands(RH_num(i,k),inu,classbc)) 581 669 582 tau_ae(i,k, id,inu) = tmp_var*tau_ae2b_int583 piz_ae(i,k, id,inu) = piz_ae2b_int584 cg_ae(i,k, id,inu) = cg_ae2b_int670 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 585 673 586 674 !--natural aerosols … … 597 685 cg_MG_6bands(RH_num(i,k),inu,classbc_pi)) 598 686 599 tau_ae_pi(i,k, id,inu) = tmp_var_pi* tau_ae2b_int600 piz_ae_pi(i,k, id,inu) = piz_ae2b_int601 cg_ae_pi(i,k, id,inu) = cg_ae2b_int687 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 602 690 603 691 ENDDO … … 610 698 DO k=1, klev 611 699 DO i=1, klon 612 tmp_var=m_allaer(i,k, spsol)/1.e6*zdh(i,k)*fac613 tmp_var_pi=m_allaer_pi(i,k, spsol)/1.e6*zdh(i,k)*fac614 615 DO inu=1, NSW700 tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac 701 tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac 702 703 DO inu=1, NSW 616 704 617 705 tau_ae2b_int= alpha_aers_6bands(RH_num(i,k),inu,spsol)+ & … … 627 715 cg_aers_6bands(RH_num(i,k),inu,spsol)) 628 716 629 tau_ae(i,k, id,inu) = tmp_var*tau_ae2b_int630 tau_ae_pi(i,k, id,inu) = tmp_var_pi* tau_ae2b_int631 piz_ae(i,k, id,inu) = piz_ae2b_int632 cg_ae(i,k, id,inu) = cg_ae2b_int717 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 633 721 634 722 ENDDO … … 638 726 !--external mixture case for soluble BC 639 727 IF (spsol.EQ.1) THEN 640 piz_ae_pi(:,:, id,:) = piz_ae(:,:,id,:)641 cg_ae_pi(:,:, id,:) = cg_ae(:,:,id,:)728 piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:) 729 cg_ae_pi(:,:,aerindex,:) = cg_ae(:,:,aerindex,:) 642 730 ENDIF 643 731 … … 648 736 DO k=1, klev 649 737 DO i=1, klon 650 tmp_var=m_allaer(i,k, naero_soluble+spinsol)/1.e6*zdh(i,k)*fac651 tmp_var_pi=m_allaer_pi(i,k, naero_soluble+spinsol)/1.e6*zdh(i,k)*fac652 653 DO inu=1, NSW738 tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac 739 tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac 740 741 DO inu=1, NSW 654 742 tau_ae2b_int = alpha_aeri_6bands(inu,spinsol) 655 743 piz_ae2b_int = piz_aeri_6bands(inu,spinsol) 656 744 cg_ae2b_int = cg_aeri_6bands(inu,spinsol) 657 745 658 tau_ae(i,k, id,inu) = tmp_var*tau_ae2b_int659 tau_ae_pi(i,k, id,inu) = tmp_var_pi*tau_ae2b_int660 piz_ae(i,k, id,inu) = piz_ae2b_int661 cg_ae(i,k, id,inu)= cg_ae2b_int746 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 662 750 ENDDO 663 751 ENDDO -
LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_lw_rrtm.F90
r2542 r2787 2 2 ! aeropt_lw_rrtm.F90 2014-05-13 C. Kleinschmitt 3 3 ! 2016-05-03 O. Boucher 4 ! 2016-12-17 O. Boucher 4 5 ! 5 6 ! This routine feeds aerosol LW properties to RRTM 6 7 ! we only consider absorption (not scattering) 8 ! we only consider dust for now 7 9 8 SUBROUTINE AEROPT_LW_RRTM( aerosol_couple,paprs,tr_seri)10 SUBROUTINE AEROPT_LW_RRTM(ok_alw, pdel, zrho, flag_aerosol, m_allaer, m_allaer_pi) 9 11 10 12 USE dimphy 11 13 USE aero_mod 12 USE infotrac_phy13 USE phys_state_var_mod, ONLY : tau_aero_lw_rrtm14 USE YO ERAD, ONLY : NLW14 USE phys_state_var_mod, ONLY: tau_aero_lw_rrtm 15 USE YOERAD, ONLY: NLW 16 USE YOMCST, ONLY: RG 15 17 16 18 IMPLICIT NONE 17 19 18 INCLUDE "YOMCST.h"19 20 INCLUDE "clesphys.h" 20 21 ! 21 22 ! Input arguments: 22 23 ! 23 LOGICAL, INTENT(IN) :: aerosol_couple 24 REAL, DIMENSION(klon,klev+1), INTENT(in) :: paprs 25 REAL, DIMENSION(klon,klev,nbtr), INTENT(in) :: tr_seri 26 ! 27 REAL, DIMENSION(klon,klev) :: zdp, mass_temp 28 ! 24 LOGICAL, INTENT(IN) :: ok_alw 25 INTEGER, INTENT(IN) :: flag_aerosol 26 REAL, DIMENSION(klon,klev), INTENT(IN) :: pdel, zrho 27 REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer, m_allaer_pi 29 28 ! 30 29 INTEGER inu, i, k 31 INTEGER :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM32 INTEGER :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M33 30 REAL :: zdh(klon,klev) 31 REAL :: tmp_var, tmp_var_pi 32 CHARACTER*20 modname 34 33 ! 35 34 !--absorption coefficient for CIDUST 36 35 REAL:: alpha_abs_CIDUST_16bands(nbands_lw_rrtm) !--unit m2/g 37 DATA alpha_abs_CIDUST_16bands / & 38 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & 39 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / 36 DATA alpha_abs_CIDUST_16bands / & 37 0.001, 0.003, 0.005, 0.006, 0.012, 0.030, 0.148, 0.098, & 38 0.017, 0.053, 0.031, 0.008, 0.010, 0.011, 0.013, 0.015 / 39 ! 40 modname='aeropt_lw_rrtm' 40 41 ! 41 42 IF (NLW.NE.nbands_lw_rrtm) THEN 42 print *,'Erreur NLW doit etre egal a 16 pour cette routine' 43 stop 43 CALL abort_physic(modname,'Erreur NLW doit etre egal a 16 pour cette routine',1) 44 44 ENDIF 45 45 ! 46 IF (aerosol_couple) THEN 47 ! 48 do i=1,nbtr 49 select case(trim(solsym(i))) 50 case ("ASBCM") 51 id_ASBCM = i 52 case ("ASPOMM") 53 id_ASPOMM = i 54 case ("ASSO4M") 55 id_ASSO4M = i 56 case ("ASMSAM") 57 id_ASMSAM = i 58 case ("CSSO4M") 59 id_CSSO4M = i 60 case ("CSMSAM") 61 id_CSMSAM = i 62 case ("SSSSM") 63 id_SSSSM = i 64 case ("CSSSM") 65 id_CSSSM = i 66 case ("ASSSM") 67 id_ASSSM = i 68 case ("CIDUSTM") 69 id_CIDUSTM = i 70 case ("AIBCM") 71 id_AIBCM = i 72 case ("AIPOMM") 73 id_AIPOMM = i 74 case ("ASNO3M") 75 id_ASNO3M = i 76 case ("CSNO3M") 77 id_CSNO3M = i 78 case ("CINO3M") 79 id_CINO3M = i 80 end select 81 enddo 46 IF (ok_alw) THEN !--aerosol LW effects 47 ! 48 IF (flag_aerosol.EQ.5.OR.flag_aerosol.EQ.6) THEN !-Dust 49 ! 50 zdh(:,:)=pdel(:,:)/(RG*zrho(:,:)) ! m 82 51 ! 83 52 DO k=1, klev 84 zdp(:,k) = (paprs(:,k)-paprs(:,k+1))/RG !--kg/m2 85 ENDDO 86 ! 87 !--for now only consider dust 88 mass_temp(:,:)=tr_seri(:,:,id_CIDUSTM) !--kg/kg 89 ! 90 DO inu=1,NLW 91 ! 92 !--total aerosol 93 tau_aero_lw_rrtm(:,:,2,inu) = mass_temp(:,:)*zdp(:,:)*1000.*alpha_abs_CIDUST_16bands(inu) 94 !--no aerosol at all 95 tau_aero_lw_rrtm(:,:,1,inu) = 0.0 53 DO i=1, klon 54 ! 55 tmp_var =m_allaer(i,k,id_CIDUSTM_phy) /1.e6*zdh(i,k) !--g/m2 56 tmp_var_pi=m_allaer_pi(i,k,id_CIDUSTM_phy)/1.e6*zdh(i,k) !--g/m2 57 ! 58 DO inu=1, NLW 59 ! 60 !--total aerosol 61 tau_aero_lw_rrtm(i,k,2,inu) = MAX(1.e-15,tmp_var*alpha_abs_CIDUST_16bands(inu)) 62 !--natural aerosol 63 ! tau_aero_lw_rrtm(:,:,1,inu) = MAX(1.e-15,tmp_var_pi*alpha_abs_CIDUST_16bands(inu)) 64 tau_aero_lw_rrtm(i,k,1,inu) = 1.e-15 !--test 65 ! 66 ENDDO 67 ENDDO 96 68 ! 97 69 ENDDO 70 ! 71 ENDIF 72 ! 73 ELSE !--no aerosol LW effects 98 74 ! 99 !--avoid very small values100 tau_aero_lw_rrtm = MAX(tau_aero_lw_rrtm,1.e-15)101 !102 ELSE !--not aerosol_couple103 !104 !--no LW effects if not coupled to INCA105 75 tau_aero_lw_rrtm = 1.e-15 106 76 ENDIF -
LMDZ5/branches/testing/libf/phylmd/rrtm/read_rsun_rrtm.F90
r2729 r2787 16 16 USE mod_phys_lmdz_para 17 17 18 USE YOESW 18 USE YOESW, ONLY : RSUN 19 19 20 20 IMPLICIT NONE … … 33 33 !$OMP THREADPRIVATE(TSI) 34 34 35 INTEGER, SAVE :: day_pre 35 INTEGER, SAVE :: day_pre=1 36 !$OMP THREADPRIVATE(day_pre) 36 37 37 38 !--only one processor reads -
LMDZ5/branches/testing/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90
r2720 r2787 1 1 ! $Id$ 2 2 ! 3 SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple, &3 SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple, ok_alw, & 4 4 new_aod, flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, & 5 5 pdtphys, pplay, paprs, t_seri, rhcl, presnivs, & … … 30 30 LOGICAL, INTENT(IN) :: debut 31 31 LOGICAL, INTENT(IN) :: aerosol_couple 32 LOGICAL, INTENT(IN) :: ok_alw 32 33 LOGICAL, INTENT(IN) :: new_aod 33 34 INTEGER, INTENT(IN) :: flag_aerosol … … 281 282 ! Calculate the total mass of all soluble aersosols 282 283 ! to be revisited for AR6 283 mass_solu_aero(:,:) = sulfacc(:,:) + bcsol(:,:) + pomsol(:,:) + nitracc(:,:)+ ssacu(:,:)284 mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:) 284 mass_solu_aero(:,:) = sulfacc(:,:) + bcsol(:,:) + pomsol(:,:) + nitracc(:,:) + ssacu(:,:) 285 mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:) 285 286 286 287 !**************************************************************************************** … … 310 311 pplay, t_seri, & 311 312 tausum_aero, tau3d_aero ) 313 314 !--call LW optical properties for tropospheric aerosols 315 CALL aeropt_lw_rrtm(ok_alw, pdel, zrho, flag_aerosol, m_allaer, m_allaer_pi) 312 316 313 317 ! Diagnostics calculation for CMIP5 protocol -
LMDZ5/branches/testing/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90
r2720 r2787 35 35 REAL, POINTER:: lev(:) 36 36 INTEGER k, band, wave, i 37 INTEGER, SAVE :: mth_pre 37 INTEGER, SAVE :: mth_pre=1 38 !$OMP THREADPRIVATE(mth_pre) 38 39 39 40 REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tau_aer_strat … … 44 45 REAL, ALLOCATABLE:: tauaerstrat_mois(:, :, :) 45 46 REAL, ALLOCATABLE:: tauaerstrat_mois_glo(:, :) 46 47 REAL, ALLOCATABLE:: sum_tau_aer_strat(:)48 47 49 48 ! For NetCDF: … … 77 76 78 77 IF (.not.ALLOCATED(tau_aer_strat)) ALLOCATE(tau_aer_strat(klon,klev)) 79 IF (.not.ALLOCATED(sum_tau_aer_strat)) ALLOCATE(sum_tau_aer_strat(klon))80 78 81 79 !--we only read monthly strat aerosol data -
LMDZ5/branches/testing/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90
r2729 r2787 41 41 REAL, POINTER:: wav(:) 42 42 INTEGER i,k,wave,band 43 INTEGER, SAVE :: mth_pre 43 INTEGER, SAVE :: mth_pre=1 44 !$OMP THREADPRIVATE(mth_pre) 44 45 45 46 REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: tau_aer_strat
Note: See TracChangeset
for help on using the changeset viewer.
