Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90	(revision 1265)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90	(revision 1266)
@@ -144,12 +144,9 @@
   
   REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! ext. coeff. Soluble comp. units *** m2/g 
-   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
-  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! ext. coeff. Insoluble comp. 1- Dust: 2- BC; 3- POM
+  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! ext. coeff. Insoluble comp. 
   REAL :: cg_aers_5wv(nbre_RH,las,naero_soluble)      ! Asym. param. soluble comp. 
-   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
-  REAL :: cg_aeri_5wv(las,naero_insoluble)            ! Asym. param. insoluble comp. 1- Dust: 2- BC; 3- POM
+  REAL :: cg_aeri_5wv(las,naero_insoluble)            ! Asym. param. insoluble comp.
   REAL :: piz_aers_5wv(nbre_RH,las,naero_soluble)   
-   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
-  REAL :: piz_aeri_5wv(las,naero_insoluble)           ! Insoluble comp. 1- Dust: 2- BC; 3- POM
+  REAL :: piz_aeri_5wv(las,naero_insoluble)           ! Insoluble comp. 
 
   REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp
@@ -308,54 +305,54 @@
 !                                   le 12 AVRIL 2006 
 !  
- DATA alpha_aers_5wv/ & 
+  DATA alpha_aers_5wv/ & 
                                 ! bc soluble 
-       7.930,7.930,7.930,7.930,7.930,7.930,     & 
-       7.930,7.930,10.893,12.618,14.550,16.613, & 
-       7.658,7.658,7.658,7.658,7.658,7.658,     & 
-       7.658,7.658,10.351,11.879,13.642,15.510, & 
-       7.195,7.195,7.195,7.195,7.195,7.195,     & 
-       7.195,7.195,9.551,10.847,12.381,13.994,  & 
-       6.736,6.736,6.736,6.736,6.736,6.736,     & 
-       6.736,6.736,8.818,9.938,11.283,12.687,   & 
-       6.277,6.277,6.277,6.277,6.277,6.277,     & 
-       6.277,6.277,8.123,9.094,10.275,11.501,   & 
-                                ! pom soluble 
-       6.676,6.676,6.676,6.676,6.710,6.934,   & 
-       7.141,7.569,8.034,8.529,9.456,10.511,  & 
-       5.109,5.109,5.109,5.109,5.189,5.535,   & 
-       5.960,6.852,8.008,9.712,12.897,19.676, & 
-       3.718,3.718,3.718,3.718,3.779,4.042,   & 
-       4.364,5.052,5.956,7.314,9.896,15.688,  & 
-       2.849,2.849,2.849,2.849,2.897,3.107,   & 
-       3.365,3.916,4.649,5.760,7.900,12.863,  & 
-       2.229,2.229,2.229,2.229,2.268,2.437,   & 
-       2.645,3.095,3.692,4.608,6.391,10.633,  & 
-                                ! Sulfate (Accumulation) 
-       5.751,6.215,6.690,7.024,7.599,8.195,      & 
-       9.156,10.355,12.660,14.823,18.908,24.508, & 
-       4.320,4.675,5.052,5.375,5.787,6.274,      & 
-       7.066,8.083,10.088,12.003,15.697,21.133,  & 
-       3.079,3.351,3.639,3.886,4.205,4.584,      & 
-       5.206,6.019,7.648,9.234,12.391,17.220,    & 
-       2.336,2.552,2.781,2.979,3.236,3.540,      & 
-       4.046,4.711,6.056,7.388,10.093,14.313,    & 
-       1.777,1.949,2.134,2.292,2.503,2.751,      & 
-       3.166,3.712,4.828,5.949,8.264,11.922,     & 
-                                ! Sulfate (Coarse) 
-       5.751,6.215,6.690,7.024,7.599,8.195,      & 
-       9.156,10.355,12.660,14.823,18.908,24.508, & 
-       4.320,4.675,5.052,5.375,5.787,6.274,      & 
-       7.066,8.083,10.088,12.003,15.697,21.133,  & 
-       3.079,3.351,3.639,3.886,4.205,4.584,      & 
-       5.206,6.019,7.648,9.234,12.391,17.220,    & 
-       2.336,2.552,2.781,2.979,3.236,3.540,      & 
-       4.046,4.711,6.056,7.388,10.093,14.313,    & 
-       1.777,1.949,2.134,2.292,2.503,2.751,      & 
-       3.166,3.712,4.828,5.949,8.264,11.922,     & 
-                                ! Seasalt soluble super_coarse (computed below for 550nm) 
-       0.50,0.90,1.05,1.21,1.40,2.41, &  
-       2.66,3.11,3.88,4.52,5.69,8.84, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
+     7.930,7.930,7.930,7.930,7.930,7.930,     & 
+     7.930,7.930,10.893,12.618,14.550,16.613, & 
+     7.658,7.658,7.658,7.658,7.658,7.658,     & 
+     7.658,7.658,10.351,11.879,13.642,15.510, & 
+     7.195,7.195,7.195,7.195,7.195,7.195,     & 
+     7.195,7.195,9.551,10.847,12.381,13.994,  & 
+     6.736,6.736,6.736,6.736,6.736,6.736,     & 
+     6.736,6.736,8.818,9.938,11.283,12.687,   & 
+     6.277,6.277,6.277,6.277,6.277,6.277,     & 
+     6.277,6.277,8.123,9.094,10.275,11.501,   & 
+     ! pom soluble 
+     6.676,6.676,6.676,6.676,6.710,6.934,   & 
+     7.141,7.569,8.034,8.529,9.456,10.511,  & 
+     5.109,5.109,5.109,5.109,5.189,5.535,   & 
+     5.960,6.852,8.008,9.712,12.897,19.676, & 
+     3.718,3.718,3.718,3.718,3.779,4.042,   & 
+     4.364,5.052,5.956,7.314,9.896,15.688,  & 
+     2.849,2.849,2.849,2.849,2.897,3.107,   & 
+     3.365,3.916,4.649,5.760,7.900,12.863,  & 
+     2.229,2.229,2.229,2.229,2.268,2.437,   & 
+     2.645,3.095,3.692,4.608,6.391,10.633,  & 
+     ! Sulfate (Accumulation) 
+     5.751,6.215,6.690,7.024,7.599,8.195,      & 
+     9.156,10.355,12.660,14.823,18.908,24.508, & 
+     4.320,4.675,5.052,5.375,5.787,6.274,      & 
+     7.066,8.083,10.088,12.003,15.697,21.133,  & 
+     3.079,3.351,3.639,3.886,4.205,4.584,      & 
+     5.206,6.019,7.648,9.234,12.391,17.220,    & 
+     2.336,2.552,2.781,2.979,3.236,3.540,      & 
+     4.046,4.711,6.056,7.388,10.093,14.313,    & 
+     1.777,1.949,2.134,2.292,2.503,2.751,      & 
+     3.166,3.712,4.828,5.949,8.264,11.922,     & 
+     ! Sulfate (Coarse) 
+     5.751,6.215,6.690,7.024,7.599,8.195,      & 
+     9.156,10.355,12.660,14.823,18.908,24.508, & 
+     4.320,4.675,5.052,5.375,5.787,6.274,      & 
+     7.066,8.083,10.088,12.003,15.697,21.133,  & 
+     3.079,3.351,3.639,3.886,4.205,4.584,      & 
+     5.206,6.019,7.648,9.234,12.391,17.220,    & 
+     2.336,2.552,2.781,2.979,3.236,3.540,      & 
+     4.046,4.711,6.056,7.388,10.093,14.313,    & 
+     1.777,1.949,2.134,2.292,2.503,2.751,      & 
+     3.166,3.712,4.828,5.949,8.264,11.922,     & 
+     ! Seasalt soluble super_coarse (computed below for 550nm) 
+     0.50,0.90,1.05,1.21,1.40,2.41, &  
+     2.66,3.11,3.88,4.52,5.69,8.84, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
      0.52,0.93,1.08,1.24,1.43,2.47, &  
      2.73,3.20,3.99,4.64,5.84,9.04, &  
@@ -364,9 +361,9 @@
      0.52,0.94,1.09,1.26,1.45,2.51, &  
      2.78,3.25,4.06,4.72,5.94,9.22, &  
-                                ! seasalt soluble coarse (computed below for 550nm) 
-       0.50,0.90,1.05,1.21,1.40,2.41, &  
-       2.66,3.11,3.88,4.52,5.69,8.84, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
+     ! seasalt soluble coarse (computed below for 550nm) 
+     0.50,0.90,1.05,1.21,1.40,2.41, &  
+     2.66,3.11,3.88,4.52,5.69,8.84, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
      0.52,0.93,1.08,1.24,1.43,2.47, &  
      2.73,3.20,3.99,4.64,5.84,9.04, &  
@@ -375,9 +372,9 @@
      0.52,0.94,1.09,1.26,1.45,2.51, &  
      2.78,3.25,4.06,4.72,5.94,9.22, &  
-                                ! seasalt soluble accumulation (computed below for 550nm) 
+     ! seasalt soluble accumulation (computed below for 550nm) 
      4.28, 7.17, 8.44, 9.85,11.60,22.44,  &  
      25.34,30.54,39.38,46.52,59.33,91.77, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
-       0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
      2.48, 4.22, 5.02, 5.94, 7.11,15.29,  &  
      17.70,22.31,30.73,38.06,52.15,90.59, &  
@@ -386,93 +383,93 @@
      1.47, 2.59, 3.12, 3.74, 4.54, 10.42, &  
      12.24,15.82,22.66,28.91,41.54,79.33/ 
-
+  
   DATA alpha_aeri_5wv/ &
-                                 ! dust insoluble 
-        0.759, 0.770, 0.775, 0.775, 0.772, & 
-                                 !!jb bc insoluble 
-        11.536,10.033, 8.422, 7.234, 6.270, & 
-                                 ! pom insoluble 
-        5.042, 3.101, 1.890, 1.294, 0.934/ 
-   ! 
+     ! dust insoluble 
+     0.759, 0.770, 0.775, 0.775, 0.772, & 
+     !!jb bc insoluble 
+     11.536,10.033, 8.422, 7.234, 6.270, & 
+     ! pom insoluble 
+     5.042, 3.101, 1.890, 1.294, 0.934/ 
+  ! 
   DATA cg_aers_5wv/ &  
                                  ! bc soluble 
-      .651, .651, .651, .651, .651, .651, & 
-      .651, .651, .738, .764, .785, .800, & 
-      .597, .597, .597, .597, .597, .597, & 
-      .597, .597, .695, .725, .751, .770, & 
-      .543, .543, .543, .543, .543, .543, & 
-      .543, .543, .650, .684, .714, .736, &  
-      .504, .504, .504, .504, .504, .504, & 
-      .504, .504, .614, .651, .683, .708, &  
-      .469, .469, .469, .469, .469, .469, & 
-      .469, .469, .582, .620, .655, .681, & 
-                                 ! pom soluble 
-      .679, .679, .679, .679, .683, .691, & 
-      .703, .720, .736, .751, .766, .784, & 
-      .656, .656, .656, .656, .659, .669, & 
-      .681, .699, .717, .735, .750, .779, &  
-      .623, .623, .623, .623, .627, .637, & 
-      .649, .668, .688, .709, .734, .762, & 
-      .592, .592, .592, .592, .595, .605, & 
-      .618, .639, .660, .682, .711, .743, & 
-      .561, .561, .561, .561, .565, .575, & 
-      .588, .609, .632, .656, .688, .724, & 
-                                 ! Accumulation sulfate 
-      .671, .684, .697, .704, .714, .723, & 
-      .734, .746, .762, .771, .781, .789, & 
-      .653, .666, .678, .687, .697, .707, & 
-      .719, .732, .751, .762, .775, .789, & 
-      .622, .635, .648, .657, .667, .678, & 
-      .691, .705, .728, .741, .758, .777, & 
-      .591, .604, .617, .627, .638, .650, & 
-      .664, .679, .704, .719, .739, .761, & 
-      .560, .574, .587, .597, .609, .621, &  
-      .637, .653, .680, .697, .719, .745, & 
-                                 ! Coarse sulfate 
-      .671, .684, .697, .704, .714, .723, & 
-      .734, .746, .762, .771, .781, .789, & 
-      .653, .666, .678, .687, .697, .707, & 
-      .719, .732, .751, .762, .775, .789, & 
-      .622, .635, .648, .657, .667, .678, & 
-      .691, .705, .728, .741, .758, .777, & 
-      .591, .604, .617, .627, .638, .650, & 
-      .664, .679, .704, .719, .739, .761, & 
-      .560, .574, .587, .597, .609, .621, &  
-      .637, .653, .680, .697, .719, .745, & 
-                                 ! For super coarse seasalt (computed below for 550nm!) 
-      0.730,0.753,0.760,0.766,0.772,0.793, &  
-      0.797,0.802,0.809,0.813,0.820,0.830, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.721,0.744,0.750,0.756,0.762,0.784, &  
-      0.787,0.793,0.800,0.804,0.811,0.822, &  
-      0.717,0.741,0.747,0.753,0.759,0.780, &  
-      0.784,0.789,0.795,0.800,0.806,0.817, &  
-      0.715,0.739,0.745,0.751,0.757,0.777, &   
-      0.781,0.786,0.793,0.797,0.803,0.814, &  
-                                 ! For coarse-soluble seasalt (computed below for 550nm!) 
-      0.730,0.753,0.760,0.766,0.772,0.793, &  
-      0.797,0.802,0.809,0.813,0.820,0.830, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.721,0.744,0.750,0.756,0.762,0.784, &  
-      0.787,0.793,0.800,0.804,0.811,0.822, &  
-      0.717,0.741,0.747,0.753,0.759,0.780, &  
-      0.784,0.789,0.795,0.800,0.806,0.817, &  
-      0.715,0.739,0.745,0.751,0.757,0.777, &   
-      0.781,0.786,0.793,0.797,0.803,0.814, &  
-                                 ! accumulation-seasalt soluble (computed below for 550nm!)  
-      0.698,0.722,0.729,0.736,0.743,0.765, &  
-      0.768,0.773,0.777,0.779,0.781,0.779, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.000,0.000,0.000,0.000,0.000,0.000, &  
-      0.658,0.691,0.701,0.710,0.720,0.756, &  
-      0.763,0.771,0.782,0.788,0.795,0.801, &  
-      0.632,0.668,0.679,0.690,0.701,0.743, &  
-      0.750,0.762,0.775,0.783,0.792,0.804, &  
-      0.605,0.644,0.656,0.669,0.681,0.729, &  
-      0.737,0.750,0.765,0.775,0.787,0.803/
- !
-
+     .651, .651, .651, .651, .651, .651, & 
+     .651, .651, .738, .764, .785, .800, & 
+     .597, .597, .597, .597, .597, .597, & 
+     .597, .597, .695, .725, .751, .770, & 
+     .543, .543, .543, .543, .543, .543, & 
+     .543, .543, .650, .684, .714, .736, &  
+     .504, .504, .504, .504, .504, .504, & 
+     .504, .504, .614, .651, .683, .708, &  
+     .469, .469, .469, .469, .469, .469, & 
+     .469, .469, .582, .620, .655, .681, & 
+     ! pom soluble 
+     .679, .679, .679, .679, .683, .691, & 
+     .703, .720, .736, .751, .766, .784, & 
+     .656, .656, .656, .656, .659, .669, & 
+     .681, .699, .717, .735, .750, .779, &  
+     .623, .623, .623, .623, .627, .637, & 
+     .649, .668, .688, .709, .734, .762, & 
+     .592, .592, .592, .592, .595, .605, & 
+     .618, .639, .660, .682, .711, .743, & 
+     .561, .561, .561, .561, .565, .575, & 
+     .588, .609, .632, .656, .688, .724, & 
+     ! Accumulation sulfate 
+     .671, .684, .697, .704, .714, .723, & 
+     .734, .746, .762, .771, .781, .789, & 
+     .653, .666, .678, .687, .697, .707, & 
+     .719, .732, .751, .762, .775, .789, & 
+     .622, .635, .648, .657, .667, .678, & 
+     .691, .705, .728, .741, .758, .777, & 
+     .591, .604, .617, .627, .638, .650, & 
+     .664, .679, .704, .719, .739, .761, & 
+     .560, .574, .587, .597, .609, .621, &  
+     .637, .653, .680, .697, .719, .745, & 
+     ! Coarse sulfate 
+     .671, .684, .697, .704, .714, .723, & 
+     .734, .746, .762, .771, .781, .789, & 
+     .653, .666, .678, .687, .697, .707, & 
+     .719, .732, .751, .762, .775, .789, & 
+     .622, .635, .648, .657, .667, .678, & 
+     .691, .705, .728, .741, .758, .777, & 
+     .591, .604, .617, .627, .638, .650, & 
+     .664, .679, .704, .719, .739, .761, & 
+     .560, .574, .587, .597, .609, .621, &  
+     .637, .653, .680, .697, .719, .745, & 
+     ! For super coarse seasalt (computed below for 550nm!) 
+     0.730,0.753,0.760,0.766,0.772,0.793, &  
+     0.797,0.802,0.809,0.813,0.820,0.830, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.721,0.744,0.750,0.756,0.762,0.784, &  
+     0.787,0.793,0.800,0.804,0.811,0.822, &  
+     0.717,0.741,0.747,0.753,0.759,0.780, &  
+     0.784,0.789,0.795,0.800,0.806,0.817, &  
+     0.715,0.739,0.745,0.751,0.757,0.777, &   
+     0.781,0.786,0.793,0.797,0.803,0.814, &  
+     ! For coarse-soluble seasalt (computed below for 550nm!) 
+     0.730,0.753,0.760,0.766,0.772,0.793, &  
+     0.797,0.802,0.809,0.813,0.820,0.830, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.721,0.744,0.750,0.756,0.762,0.784, &  
+     0.787,0.793,0.800,0.804,0.811,0.822, &  
+     0.717,0.741,0.747,0.753,0.759,0.780, &  
+     0.784,0.789,0.795,0.800,0.806,0.817, &  
+     0.715,0.739,0.745,0.751,0.757,0.777, &   
+     0.781,0.786,0.793,0.797,0.803,0.814, &  
+     ! accumulation-seasalt soluble (computed below for 550nm!)  
+     0.698,0.722,0.729,0.736,0.743,0.765, &  
+     0.768,0.773,0.777,0.779,0.781,0.779, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.000,0.000,0.000,0.000,0.000,0.000, &  
+     0.658,0.691,0.701,0.710,0.720,0.756, &  
+     0.763,0.771,0.782,0.788,0.795,0.801, &  
+     0.632,0.668,0.679,0.690,0.701,0.743, &  
+     0.750,0.762,0.775,0.783,0.792,0.804, &  
+     0.605,0.644,0.656,0.669,0.681,0.729, &  
+     0.737,0.750,0.765,0.775,0.787,0.803/
+  !
+  
   DATA cg_aeri_5wv/&
      ! dust insoluble
@@ -485,82 +482,82 @@
   DATA piz_aers_5wv/&
                            ! bc soluble 
-  .445, .445, .445, .445, .445, .445, & 
-  .445, .445, .470, .487, .508, .531, & 
-  .442, .442, .442, .442, .442, .442, & 
-  .442, .442, .462, .481, .506, .533, & 
-  .427, .427, .427, .427, .427, .427, & 
-  .427, .427, .449, .470, .497, .526, & 
-  .413, .413, .413, .413, .413, .413, & 
-  .413, .413, .437, .458, .486, .516, & 
-  .399, .399, .399, .399, .399, .399, & 
-  .399, .399, .423, .445, .473, .506, & 
-                           ! pom soluble 
-  .975, .975, .975, .975, .975, .977, & 
-  .979, .982, .984, .987, .990, .994, & 
-  .972, .972, .972, .972, .973, .974, & 
-  .977, .980, .983, .986, .989, .993, & 
-  .963, .963, .963, .963, .964, .966, & 
-  .969, .974, .977, .982, .986, .991, & 
-  .955, .955, .955, .955, .955, .958, & 
-  .962, .967, .972, .977, .983, .989, & 
-  .944, .944, .944, .944, .944, .948, & 
-  .952, .959, .962, .972, .979, .987, & 
-                           ! sulfate soluble accumulation 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-                           ! sulfate soluble coarse 
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-                           ! seasalt super coarse (computed below for 550nm) 
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, & 
-                           ! seasalt coarse (computed below for 550nm) 
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-                           ! seasalt soluble accumulation (computed below for 550nm) 
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000, &  
-  1.000,1.000,1.000,1.000,1.000,1.000/ 
-
- !
+     .445, .445, .445, .445, .445, .445, & 
+     .445, .445, .470, .487, .508, .531, & 
+     .442, .442, .442, .442, .442, .442, & 
+     .442, .442, .462, .481, .506, .533, & 
+     .427, .427, .427, .427, .427, .427, & 
+     .427, .427, .449, .470, .497, .526, & 
+     .413, .413, .413, .413, .413, .413, & 
+     .413, .413, .437, .458, .486, .516, & 
+     .399, .399, .399, .399, .399, .399, & 
+     .399, .399, .423, .445, .473, .506, & 
+     ! pom soluble 
+     .975, .975, .975, .975, .975, .977, & 
+     .979, .982, .984, .987, .990, .994, & 
+     .972, .972, .972, .972, .973, .974, & 
+     .977, .980, .983, .986, .989, .993, & 
+     .963, .963, .963, .963, .964, .966, & 
+     .969, .974, .977, .982, .986, .991, & 
+     .955, .955, .955, .955, .955, .958, & 
+     .962, .967, .972, .977, .983, .989, & 
+     .944, .944, .944, .944, .944, .948, & 
+     .952, .959, .962, .972, .979, .987, & 
+     ! sulfate soluble accumulation 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     ! sulfate soluble coarse 
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     ! seasalt super coarse (computed below for 550nm) 
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, & 
+     ! seasalt coarse (computed below for 550nm) 
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     ! seasalt soluble accumulation (computed below for 550nm) 
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000, &  
+     1.000,1.000,1.000,1.000,1.000,1.000/ 
+  
+  !
   DATA piz_aeri_5wv/&
      ! dust insoluble
@@ -808,4 +805,9 @@
 !CDIR ON_ADB(alpha_aers_5wv)
             DO i=1, KLON
+              rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
+              RH_num = INT( rh(i,k)/10. + 1.)
+              IF (rh(i,k).GT.85.) RH_num=10
+              IF (rh(i,k).GT.90.) RH_num=11
+
               tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
                              (alpha_aers_5wv(RH_num+1,la,spsol) - & 
Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90	(revision 1265)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90	(revision 1266)
@@ -486,5 +486,7 @@
        ctetaSTD,dtime, ok_veget, &
        type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
-       ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz)
+       ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
+       new_aod, aerosol_couple)   
+
 
   USE iophy 
@@ -507,5 +509,5 @@
   logical                               :: ok_mensuel, ok_journe, ok_hf, ok_instan
   logical                               :: ok_LES,ok_ade,ok_aie
-
+  logical                               :: new_aod, aerosol_couple
   integer, intent(in)::  read_climoz ! read ozone climatology
   !     Allowed values are 0, 1 and 2
@@ -833,9 +835,11 @@
  CALL histdef2d(iff,o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
  CALL histdef2d(iff,o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day")
-     END DO
-
-     DO naero = 1, naero_spc
-CALL histdef2d(iff,o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
-     END DO
+END DO
+
+IF (new_aod .AND. (.NOT. aerosol_couple)) THEN
+  DO naero = 1, naero_spc
+  CALL histdef2d(iff,o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
+  END DO
+ENDIF
 
  IF (ok_ade) THEN
@@ -853,4 +857,5 @@
  CALL histdef2d(iff,o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2")
 
+ IF (.NOT. aerosol_couple) THEN 
  CALL histdef2d(iff,o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2")
  CALL histdef2d(iff,o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing  at SRF", "W/m2")
@@ -859,4 +864,5 @@
  CALL histdef2d(iff,o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2")
  CALL histdef2d(iff,o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2")
+ ENDIF
 
  ENDIF
Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h	(revision 1265)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h	(revision 1266)
@@ -834,5 +834,5 @@
 
 ! OD550 per species
-      IF (new_aod) THEN
+      IF (new_aod .and. (.not. aerosol_couple)) THEN
       DO naero = 1, naero_spc
           IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
@@ -898,4 +898,6 @@
 
 !cf
+
+        if (.not. aerosol_couple) then
         IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
         CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
@@ -927,4 +929,5 @@
      $      solswcf_aero(:,3))
         ENDIF
+        endif
 
 	endif ! new_aod
Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F	(revision 1265)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F	(revision 1266)
@@ -1457,5 +1457,6 @@
      &                        ctetaSTD,dtime,ok_veget,
      &                        type_ocean,iflag_pbl,ok_mensuel,ok_journe,
-     &              ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz)
+     &                        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, 
+     &                        read_climoz, new_aod, aerosol_couple)
 c$OMP END MASTER
 c$OMP BARRIER
@@ -1533,5 +1534,6 @@
      $                   pdtphys,
      $                   annee_ref,
-     $                   day_ini)
+     $                   day_ref, 
+     $                   itau_phy)
 
          CALL VTe(VTinca)
@@ -2826,5 +2828,5 @@
             CALL AEROSOL_METEO_CALC(
      $           calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs,
-     $           prfl,psfl,pctsrf,airephy,xjour,rlat,rlon,u10m,v10m)
+     $           prfl,psfl,pctsrf,airephy,rlat,rlon,u10m,v10m)
          END IF
 
@@ -3385,5 +3387,4 @@
      $                        day_ini,
      $                        airephy,
-     $                        xjour,
      $                        pphi,
      $                        pphis,
