Index: /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3684)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3685)
@@ -104,11 +104,8 @@
       REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV)         ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1).
       REAL,INTENT(OUT) :: GSR_nu(ngrid,L_NSPECTV)         ! Surface SW radiation in each band (Normalized to the band width (W/m2/cm-1).
+      REAL,INTENT(OUT) :: int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
+      REAL,INTENT(OUT) :: int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
       REAL,INTENT(OUT) :: tau_col(ngrid)                  ! Diagnostic from aeropacity.
       REAL,INTENT(OUT) :: albedo_equivalent(ngrid)        ! Spectrally Integrated Albedo. For Diagnostic. By MT2015
-      REAL,INTENT(OUT) :: int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
-      REAL,INTENT(OUT) :: int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
-
-
-
 
 
Index: /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90	(revision 3684)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90	(revision 3685)
@@ -5,10 +5,11 @@
 CONTAINS
 
-    subroutine callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,   &
-          albedo,emis,mu0,pplev,pplay,pt,   &
-          zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,    &
-          dtlw,dtsw,fluxsurf_lw,    &
-          fluxsurf_sw,fluxtop_lw,fluxtop_sw,fluxtop_dn, &
-          reffrad,tau_col,ptime,pday,firstcall,lastcall)
+    subroutine callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,&
+          albedo,emis,mu0,pplev,pplay,pt,                      &
+          zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,          &
+          dtlw,dtsw,fluxsurf_lw,                               &
+          fluxsurf_sw,fluxtop_lw,fluxtop_sw,fluxtop_dn,        &
+          reffrad,int_dtaui,int_dtauv,tau_col,                 &
+          ptime,pday,firstcall,lastcall)
 
       use radinc_h
@@ -102,4 +103,6 @@
       REAL fluxtop_sw(ngrid)    ! outgoing LW flux to space (W/m2)
       REAL fluxtop_dn(ngrid)    ! incident top of atmosphere SW flux (W/m2)
+      REAL int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
+      REAL int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
 
 !-----------------------------------------------------------------------
@@ -745,4 +748,21 @@
          end do
 
+         ! Optical thickness diagnostics
+         ! Output exp(-tau) because gweight ponderates exp and not tau itself
+         int_dtauv(ig,:,:) = 0.0d0
+         int_dtaui(ig,:,:) = 0.0d0
+         do l=1,L_NLAYRAD
+            do nw=1,L_NSPECTV
+               do k=1,L_NGAUSS
+                  int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
+               enddo
+            enddo
+            do nw=1,L_NSPECTI
+               do k=1,L_NGAUSS
+                  int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
+               enddo
+            enddo
+         enddo
+
 ! **********************************************************
 !     NON NLTE correction in Pluto atmosphere
Index: /trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3684)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3685)
@@ -1022,6 +1022,7 @@
                                zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,       &
                                zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,   &
-                               fluxtop_sw,fluxtop_dn,reffrad,tau_col,ptime,pday, &
-                               firstcall,lastcall)
+                               fluxtop_sw,fluxtop_dn,reffrad,                    &
+                               int_dtaui,int_dtauv,tau_col,                      &
+                               ptime,pday,firstcall,lastcall)
                   albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
                   fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+       &