Index: trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90	(revision 3685)
+++ trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90	(revision 3686)
@@ -132,4 +132,5 @@
   dpr(:)        = 0.0
   lkcoef(:,:)   = 0.0
+  DTAUKV(:,:,:) = 0.0 
 
   do K=2,L_LEVELS
@@ -365,8 +366,9 @@
 
   ! Total extinction optical depths
-
-  DO NG=1,L_NGAUSS       ! full gauss loop
+  DO NG=1,L_NGAUSS ! full gauss loop
      DO NW=1,L_NSPECTV       
+        TAUV(1,NW,NG)=0.0D0
         TAUCUMV(1,NW,NG)=0.0D0
+        
         DO K=2,L_LEVELS
            TAUCUMV(K,NW,NG)=TAUCUMV(K-1,NW,NG)+DTAUKV(K,NW,NG)
@@ -378,8 +380,5 @@
         TAUV(L,NW,NG)=TAUCUMV(2*L_NLAYRAD+1,NW,NG)
      END DO            
-  END DO                 ! end full gauss loop
-
-
-
+  END DO ! end full gauss loop
 
 end subroutine optcv
Index: trunk/LMDZ.PLUTO/libf/phypluto/optcv_pluto_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/optcv_pluto_mod.F90	(revision 3685)
+++ trunk/LMDZ.PLUTO/libf/phypluto/optcv_pluto_mod.F90	(revision 3686)
@@ -260,36 +260,30 @@
 
         ! Last level
-
         L              = L_NLAYRAD
         K              = 2*L+1
-	DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
+	     DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
 
         WBARV(L,NW,NG) = ctemp(L,NW) / DTAUV(L,NW,NG)
-        !print*, 'TB22 : WBARV(L)=',WBARV(L,NW,NG),NW,NG
-        !print*, 'TB22 : ctemp(L)=',ctemp(L,NW),NW
-        !print*, 'TB22 : dtauv(L)=',DTAUV(L,NW,NG),NW,NG
-     END DO                 ! NW spectral loop
-  END DO                    ! NG Gauss loop
+     END DO ! NW spectral loop
+  END DO ! NG Gauss loop
 
   ! Total extinction optical depths
-
-  DO NG=1,L_NGAUSS       ! full gauss loop
+  DO NG=1,L_NGAUSS ! full gauss loop
      DO NW=1,L_NSPECTV
         TAUV(1,NW,NG)=0.0D0
+        TAUCUMV(1,NW,NG)=0.0D0
+        
         DO L=1,L_NLAYRAD
            TAUV(L+1,NW,NG)=TAUV(L,NW,NG)+DTAUV(L,NW,NG)
         END DO
 
-        TAUCUMV(1,NW,NG)=0.0D0
         DO K=2,L_LEVELS
            TAUCUMV(K,NW,NG)=TAUCUMV(K-1,NW,NG)+DTAUKV(K,NW,NG)
         END DO
      END DO
-  END DO                 ! end full gauss loop
-
+  END DO ! end full gauss loop
 
   return
 
-
 end subroutine optcv_pluto