Index: LMDZ6/trunk/libf/phylmd/cv3_routines.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cv3_routines.f90	(revision 5705)
+++ LMDZ6/trunk/libf/phylmd/cv3_routines.f90	(revision 5706)
@@ -614,5 +614,6 @@
 
 ! Compute icbmax.
-
+  
+  !ym do not do that, independance between column !
   icbmax = 2
   DO i = 1, len
@@ -700,17 +701,21 @@
 
 ! Re-compute icbsmax (icbsmax2):                          !convect3
-!                                                         !convect3
-  icbsmax2 = 2                                            !convect3
-  DO i = 1, len                                           !convect3
-    icbsmax2 = max(icbsmax2, icbs(i))                     !convect3
-  END DO                                                  !convect3
+!
+!ym column independant, do not use reduction                                                         !convect3
+!ym  icbsmax2 = 2                                            !convect3
+!ym  DO i = 1, len                                           !convect3
+!ym    icbsmax2 = max(icbsmax2, icbs(i))                     !convect3
+!ym  END DO                                                  !convect3
 
 ! initialization outputs:
 
-  DO k = 1, icbsmax2                                      ! convect3
+!ym  DO k = 1, icbsmax2                                      ! convect3
+  DO k = 1, nd                                      ! convect3
     DO i = 1, len                                         ! convect3
-      tp(i, k) = 0.0                                      ! convect3
-      tvp(i, k) = 0.0                                     ! convect3
-      clw(i, k) = 0.0                                     ! convect3
+      IF (k<=MAX(2,icbs(i))) THEN
+        tp(i, k) = 0.0                                      ! convect3
+        tvp(i, k) = 0.0                                     ! convect3
+        clw(i, k) = 0.0                                     ! convect3
+      ENDIF
     END DO                                                ! convect3
   END DO                                                  ! convect3
@@ -718,8 +723,11 @@
 ! tp and tvp below cloud base:
 
-  DO k = minorig, icbsmax2 - 1
+!ym  DO k = minorig, icbsmax2 - 1
+  DO k = minorig, nd
     DO i = 1, len
-      tp(i, k) = tnk(i) - (gz(i,k)-gznk(i))*cpinv(i)
-      tvp(i, k) = tp(i, k)*(1.+qnk(i)/eps-qnk(i))        !whole thing (convect3)
+      IF (k<=MAX(2,icbs(i))-1) THEN
+        tp(i, k) = tnk(i) - (gz(i,k)-gznk(i))*cpinv(i)
+        tvp(i, k) = tp(i, k)*(1.+qnk(i)/eps-qnk(i))        !whole thing (convect3)
+      ENDIF
     END DO
   END DO