--- trunk/phylmd/vdif_kcay.f	2014/09/04 10:40:24	105
+++ trunk/Sources/phylmd/vdif_kcay.f	2015/06/16 15:23:29	145
@@ -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, plev, 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
@@ -31,8 +33,7 @@
 
     REAL, intent(in):: dt
     real, intent(in):: g
-    real rconst
-    real plev(klon, klev+1), temp(klon, klev)
+    real plev(klon, klev+1)
     real ustar(klon), snstable
     REAL zlev(klon, klev+1)
     REAL zlay(klon, klev)
@@ -57,7 +58,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)
@@ -196,8 +197,7 @@
 
     ! Initialisation de q2
 
-    call yamada(ngrid, g, rconst, plev, temp, zlev, zlay, u, v, teta, &
-         q2diag, km, kn, ustar, l_mix)
+    call yamada(ngrid, g, zlev, zlay, u, v, teta, q2diag, km, kn)
     if (first.and.1.eq.1) then
        first=.false.
        q2=q2diag