Index: LMDZ6/trunk/libf/phylmd/vdif_kcay.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/vdif_kcay.F90	(revision 4645)
+++ LMDZ6/trunk/libf/phylmd/vdif_kcay.F90	(revision 4649)
@@ -221,5 +221,5 @@
   END IF
 
-  DO ilev = 1, nlev
+  DO ilev = 2, nlay
     DO igrid = 1, ngrid
       q2(igrid, ilev) = amax1(q2(igrid,ilev), q2min)
@@ -591,4 +591,7 @@
   IF (1==1) THEN
 
+    sss=0.
+    sssq=0.
+    ! WARNING : travail sur le point ig=1 ????
     DO ilev = 2, klev - 1
       sss = sss + plev(1, ilev-1) - plev(1, ilev+1)
@@ -624,5 +627,5 @@
 
 
-    DO ilev = 1, nlev
+    DO ilev = 2, klev-1
       DO igrid = 1, ngrid
         q2(igrid, ilev) = max(q2(igrid,ilev), q2min)
Index: LMDZ6/trunk/libf/phylmd/yamada.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/yamada.F90	(revision 4645)
+++ LMDZ6/trunk/libf/phylmd/yamada.F90	(revision 4649)
@@ -128,5 +128,5 @@
     l0(ig) = 100.
   END DO
-  DO k = 2, klev - 1
+  DO k = 2, klev
     DO ig = 1, ngrid
       l(ig, k) = l0(ig)*kap*zlev(ig, k)/(kap*zlev(ig,k)+l0(ig))
@@ -139,5 +139,6 @@
       sqz(ig) = 1.E-10
     END DO
-    DO k = 2, klev - 1
+    !DO k = 2, klev - 1
+    DO k = 2, klev
       DO ig = 1, ngrid
         q2(ig, k) = l(ig, k)**2*zz(ig, k)