Index: LMDZ6/trunk/libf/phylmd/cv_routines.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cv_routines.f90	(revision 5704)
+++ LMDZ6/trunk/libf/phylmd/cv_routines.f90	(revision 5705)
@@ -808,9 +808,10 @@
   ! Compute icbmax.
   ! -------------------------------------------------------------------
-
-  icbmax = 2
-  DO i = 1, ncum
-    icbmax = max(icbmax, icb(i))
-  END DO
+  
+!ym independance betwwen column
+!ym  icbmax = 2
+!ym  DO i = 1, ncum
+!ym    icbmax = max(icbmax, icb(i))
+!ym  END DO
 
   ! =====================================================================
@@ -836,10 +837,13 @@
 
   ! dtpbl = average of tvp-tv in the PBL (k=nk to icb-1).
-
-  DO k = minorig, icbmax
+!ym independance betwwen column
+!ym  DO k = minorig, icbmax
+  DO k = minorig, nd
     DO i = 1, ncum
-      IF ((k>=nk(i)) .AND. (k<=(icb(i)-1))) THEN
-        dtpbl(i) = dtpbl(i) + (tvp(i,k)-tv(i,k))*dph(i, k)
-      END IF
+      IF (k<=MAX(2,icb(i))) THEN
+        IF ((k>=nk(i)) .AND. (k<=(icb(i)-1))) THEN
+          dtpbl(i) = dtpbl(i) + (tvp(i,k)-tv(i,k))*dph(i, k)
+        END IF
+      ENDIF
     END DO
   END DO