--- trunk/phylmd/vdif_kcay.f	2014/09/16 14:00:41	108
+++ trunk/Sources/phylmd/vdif_kcay.f	2016/03/11 18:47:26	178
@@ -4,14 +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
@@ -32,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)
@@ -58,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)
@@ -107,7 +107,6 @@
     ! gnmax : borne superieure de gn (0.0233)
     ! gninf : vrai si gn est en dessous de sa borne inferieure
     ! gnsup : vrai si gn est en dessus de sa borne superieure
-    ! gm : drole d'objet bien utile
     ! ri : nombre de Richardson
     ! sn : coefficient de stabilite pour n
     ! snq2 : premier terme du developement limite de sn en q2
@@ -119,8 +118,6 @@
     REAL gnmax
     LOGICAL gninf
     LOGICAL gnsup
-    REAL gm
-    ! REAL ri(klon, klev+1)
     REAL sn(klon, klev+1)
     REAL snq2(klon, klev+1)
     REAL sm(klon, klev+1)
@@ -197,8 +194,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
@@ -327,9 +323,6 @@
 
           gn=-long(igrid, ilev)**2 / q2(igrid, ilev) &
                * n2(igrid, ilev)
-          gm=long(igrid, ilev)**2 / q2(igrid, ilev) &
-               * m2(igrid, ilev)
-
           gninf=.false.
           gnsup=.false.
           long(igrid, ilev)=long(igrid, ilev)