Index: trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3694)
+++ trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3695)
@@ -927,4 +927,9 @@
 !-----------------------------------------------------------------------
 
+         ! Call everywhere for diagnostics.
+         call optcv(dtauv,tauv,taucumv,plevrad,                 &
+                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero,   &
+                 tmid,pmid,taugsurf,qvar,muvarrad)
+
          if(fract(ig) .ge. 1.0e-4) then ! Only during daylight.
             if((ngrid.eq.1).and.(global1d))then
@@ -938,8 +943,4 @@
             endif
 
-            call optcv(dtauv,tauv,taucumv,plevrad,                 &
-                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero,   &
-                 tmid,pmid,taugsurf,qvar,muvarrad)
-
             call sfluxv(dtauv,tauv,taucumv,albv,dwnv,wbarv,cosbv,   &
                  acosz,stel_fract,nfluxtopv,fluxtopvdn,nfluxoutv_nu,&
@@ -953,10 +954,8 @@
             nfluxoutv_nu(:) = 0.0d0
             nfluxgndv_nu(:) = 0.0d0
-            fmnetv_nu(:,:)=0.0
-            do l=1,L_NLAYRAD
-               fmnetv(l)=0.0d0
-               fluxupv(l)=0.0d0
-               fluxdnv(l)=0.0d0
-            end do
+            fmnetv_nu(:,:)  = 0.0d0
+            fmnetv(:)       = 0.0d0
+            fluxupv(:)      = 0.0d0
+            fluxdnv(:)      = 0.0d0
          end if
 
@@ -1068,10 +1067,10 @@
             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)
+                  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)
+                  int_dtaui(ig,l,nw) = int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
                enddo
             enddo
Index: trunk/LMDZ.PLUTO/libf/phypluto/optci.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/optci.F90	(revision 3694)
+++ trunk/LMDZ.PLUTO/libf/phypluto/optci.F90	(revision 3695)
@@ -133,4 +133,7 @@
   !     spectral interval, NW, and each Gauss point, NG.
 
+  dtaui(:,:,:) = 0.0
+  taucumi(:,:,:) = 0.0
+
   taugsurf(:,:) = 0.0
   dpr(:)        = 0.0
Index: trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90	(revision 3694)
+++ trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90	(revision 3695)
@@ -128,4 +128,8 @@
   !     Calculate the continuum opacities, i.e., those that do not depend on
   !     NG, the Gauss index.
+
+  dtauv(:,:,:)   = 0.0
+  tauv(:,:,:)    = 0.0
+  taucumv(:,:,:) = 0.0
 
   taugsurf(:,:) = 0.0
Index: trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3694)
+++ trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3695)
@@ -2201,5 +2201,9 @@
       !   call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw))
       !enddo
-
+      call write_output('dtauv_01','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,1))
+      call write_output('dtauv_23','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,23))
+      call write_output('dtaui_01','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,1))
+      call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17))
+      
       ! Temporary inclusions for heating diagnostics.
       if (.not.fast) then