Index: LMDZ6/trunk/libf/phylmd/cv3p2_closure.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cv3p2_closure.f90	(revision 5709)
+++ LMDZ6/trunk/libf/phylmd/cv3p2_closure.f90	(revision 5710)
@@ -23,6 +23,6 @@
   ! **************************************************************
 
-USE yomcst2_mod_h
-   USE lmdz_cv_ini, ONLY : alpha,alpha1,beta,flag_wb,minorig,nl,noconv_stop,pbcrit,rrd,wbmax,coef_peel
+  USE yomcst2_mod_h
+  USE lmdz_cv_ini, ONLY : alpha,alpha1,beta,flag_wb,minorig,nl,noconv_stop,pbcrit,rrd,wbmax,coef_peel
   USE conema3_mod_h
   USE cvflag_mod_h
@@ -167,19 +167,25 @@
   ! compute icbmax:
 
-  icbmax = 2
-  DO il = 1, ncum
-    icbmax = max(icbmax, icb(il))
-  END DO
+!ym break the column independance
+!ym  icbmax = 2
+!ym  DO il = 1, ncum
+!ym    icbmax = max(icbmax, icb(il))
+!ym  END DO
+
   ! if(prt.level.GE.20) print*,'cv3p2_closure apres 200'
 
   ! update sig and w0 below cloud base:
 
-  DO k = 1, icbmax
-    DO il = 1, ncum
-      IF (k<=icb(il)) THEN
-        sig(il, k) = beta*sig(il, k) - 2.*alpha*buoy(il, icb(il))*buoy(il,icb(il))
-        sig(il, k) = amax1(sig(il,k), 0.0)
-        w0(il, k) = beta*w0(il, k)
-      END IF
+!ym column independance
+!ym  DO k = 1, icbmax
+  DO k = 1, nd
+    DO il = 1, ncum
+      IF (k<=MAX(2,icb(il))) THEN
+        IF (k<=icb(il)) THEN
+          sig(il, k) = beta*sig(il, k) - 2.*alpha*buoy(il, icb(il))*buoy(il,icb(il))
+          sig(il, k) = amax1(sig(il,k), 0.0)
+          w0(il, k) = beta*w0(il, k)
+        END IF
+      ENDIF
     END DO
   END DO