Index: LMDZ5/trunk/libf/phylmd/aero_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/aero_mod.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/aero_mod.F90	(revision 2550)
@@ -88,6 +88,7 @@
 
 ! Number of diagnostics wavelengths (5 SW + 1 LW @ 10 um)
-  INTEGER, PARAMETER :: nwave = 5
+  INTEGER, PARAMETER :: nwave_sw = 5
   INTEGER, PARAMETER :: nwave_lw = 1
+  INTEGER, PARAMETER :: nwave = nwave_sw + nwave_lw
 
 ! Number of modes spectral bands
Index: LMDZ5/trunk/libf/phylmd/aeropt_5wv.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/aeropt_5wv.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/aeropt_5wv.F90	(revision 2550)
@@ -74,5 +74,5 @@
   ! Local
   !
-  INTEGER, PARAMETER :: las = nwave
+  INTEGER, PARAMETER :: las = nwave_sw
   LOGICAL :: soluble
   
Index: LMDZ5/trunk/libf/phylmd/readaerosolstrato.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/readaerosolstrato.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/readaerosolstrato.F90	(revision 2550)
@@ -55,5 +55,5 @@
     data piz_strat  /0.9999998, 0.99762493/
     data cg_strat   /0.73107845,0.73229635/
-    real, dimension(nwave) :: alpha_strat_wave
+    real, dimension(nwave_sw) :: alpha_strat_wave
     data alpha_strat_wave/3.36780953,3.34667683,3.20444202,3.0293026,2.82108808/
 
@@ -153,5 +153,5 @@
 
 !--total vertical aod at the 6 wavelengths
-    DO wave=1, nwave
+    DO wave=1, nwave_sw
     DO k=1, klev
     tausum_aero(:,wave,id_STRAT_phy)=tausum_aero(:,wave,id_STRAT_phy)+tau_aer_strat(:,k)*alpha_strat_wave(wave)/alpha_strat_wave(2)
Index: LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 2550)
@@ -70,5 +70,5 @@
   ! Local
   !
-  INTEGER, PARAMETER :: las = nwave
+  INTEGER, PARAMETER :: las = nwave_sw
   LOGICAL :: soluble
   
Index: LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90	(revision 2550)
@@ -61,9 +61,9 @@
 !--diagnostics AOD in the SW
 ! alpha_sw_strat_wave is *not* normalised by the 550 nm extinction coefficient
-    real, dimension(nwave) :: alpha_sw_strat_wave
+    real, dimension(nwave_sw) :: alpha_sw_strat_wave
     data alpha_sw_strat_wave/3.708007,4.125824,4.136584,3.887478,3.507738/
 !
-!--diagnostics AOD in the LW at 10 um 
-    real :: alpha_lw_strat_wave
+!--diagnostics AOD in the LW at 10 um (not normalised by the 550 nm ext coefficient
+    real :: alpha_lw_strat_wave(nwave_lw)
     data alpha_lw_strat_wave/0.2746812/
 !
@@ -171,8 +171,8 @@
 
 !--total vertical aod at the 5 SW wavelengths
-    DO wave=1, nwave
+    DO wave=1, nwave_sw
     DO k=1, klev
-    tausum_aero(:,wave,id_STRAT_phy)=tausum_aero(:,wave,id_STRAT_phy)+ &
-       tau_aer_strat(:,k)*alpha_sw_strat_wave(wave)/alpha_sw_strat_wave(2)
+      tausum_aero(:,wave,id_STRAT_phy)=tausum_aero(:,wave,id_STRAT_phy)+ &
+          tau_aer_strat(:,k)*alpha_sw_strat_wave(wave)/alpha_sw_strat_wave(2)
     ENDDO
     ENDDO
@@ -210,4 +210,12 @@
     ENDIF 
 
+!--total vertical aod at the 1 LW wavelength
+    DO wave=1, nwave_lw
+    DO k=1, klev
+      tausum_aero(:,nwave_sw+wave,id_STRAT_phy)=tausum_aero(:,nwave_sw+wave,id_STRAT_phy)+ &
+         tau_aer_strat(:,k)*alpha_lw_strat_wave(wave)/alpha_sw_strat_wave(2)
+    ENDDO
+    ENDDO
+
     DO band=1, nbands_lw_rrtm
     tau_aero_lw_rrtm(:,:,2,band)  = tau_aero_lw_rrtm(:,:,2,band) + alpha_lw_abs_rrtm(band)*tau_aer_strat(:,:)
Index: LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90	(revision 2549)
+++ LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90	(revision 2550)
@@ -269,10 +269,11 @@
 !--total vertical aod at the 5 SW wavelengths
 !--for now use band 3 AOD into all 5 wavelengths
+!--it is only a reasonable approximation for 550 nm (wave=2)
     band=3
     DO i=1, klon
     DO k=1, klev
       IF (stratomask(i,k).GT.0.999999) THEN
-        DO wave=1, nwave
-          tausum_aero(:,wave,id_STRAT_phy)=tausum_aero(:,wave,id_STRAT_phy)+tau_aer_strat(:,k,band)
+        DO wave=1, nwave_sw
+          tausum_aero(i,wave,id_STRAT_phy)=tausum_aero(i,wave,id_STRAT_phy)+tau_aer_strat(i,k,band)
         ENDDO
       ENDIF
@@ -308,4 +309,17 @@
     ENDDO
 
+!--total vertical aod at 10 um
+!--this is approximated from band 7 of RRTM
+    band=7
+    DO i=1, klon
+    DO k=1, klev
+      IF (stratomask(i,k).GT.0.999999) THEN
+        DO wave=1, nwave_lw
+          tausum_aero(i,nwave_sw+wave,id_STRAT_phy)=tausum_aero(i,nwave_sw+wave,id_STRAT_phy)+taulw_aer_strat(i,k,band)
+        ENDDO
+      ENDIF
+    ENDDO
+    ENDDO
+
     DO band=1, NLW
     WHERE (stratomask.GT.0.999999)