--- trunk/phylmd/vdif_kcay.f90	2013/11/15 18:45:49	76
+++ trunk/phylmd/vdif_kcay.f	2014/12/18 17:30:24	118
@@ -4,13 +4,15 @@
 
 contains
 
-  SUBROUTINE vdif_kcay(ngrid, dt, g, rconst, plev, temp, zlev, zlay, &
-       u, v, teta, cd, q2, q2diag, km, kn, ustar, l_mix)
+  SUBROUTINE vdif_kcay(ngrid, dt, g, rconst, plev, temp, zlev, zlay, u, v, &
+       teta, cd, q2, q2diag, km, kn, ustar, l_mix)
 
     ! From LMDZ4/libf/phylmd/vdif_kcay.F, version 1.1 2004/06/22 11:45:36
 
     USE dimphy, ONLY: klev, klon
+    use yamada_m, only: yamada
 
+    INTEGER ngrid
     ! dt : pas de temps
     ! g : g
     ! zlev : altitude a chaque niveau (interface inferieure de la couche
@@ -40,13 +42,13 @@
     REAL v(klon, klev)
     REAL teta(klon, klev)
     REAL, intent(in):: cd (:) ! (ngrid) cdrag, valeur au debut du pas de temps
-    REAL q2(klon, klev+1), q2s(klon, klev+1)
+    REAL q2(klon, klev+1)
     REAL q2diag(klon, klev+1)
     REAL km(klon, klev+1)
     REAL kn(klon, klev+1)
-    real sq(klon), sqz(klon), zz(klon, klev+1), zq, long0(klon)
+    real sq(klon), sqz(klon), zq, long0(klon)
 
-    integer l_mix, iii
+    integer l_mix
 
     ! nlay : nombre de couches 
     ! nlev : nombre de niveaux
@@ -57,7 +59,7 @@
     ! q : echelle de vitesse au bas de chaque couche
     ! (valeur a la fin du pas de temps)
 
-    INTEGER nlay, nlev, ngrid
+    INTEGER nlay, nlev
     REAL unsdz(klon, klev)
     REAL unsdzdec(klon, klev+1)
     REAL q(klon, klev+1)