--- trunk/Sources/phylmd/coefkz.f	2016/12/07 16:44:53	208
+++ trunk/Sources/phylmd/coefkz.f	2017/11/07 10:52:46	233
@@ -5,7 +5,7 @@
 contains
 
   SUBROUTINE coefkz(nsrf, paprs, pplay, ksta, ksta_ter, ts, rugos, u, v, t, &
-       q, qsurf, coefm, coefh)
+       q, qsurf, coefm, coefh, ycdragm, ycdragh)
 
     ! Authors: F. Hourdin, M. Forichon, Z. X. Li (LMD/CNRS)
     ! Date: September 22nd, 1993
@@ -15,7 +15,7 @@
     use clcdrag_m, only: clcdrag
     USE conf_phys_m, ONLY: iflag_pbl
     USE dimphy, ONLY: klev, klon
-    USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats
+    USE fcttre, ONLY: foede, foeew
     USE indicesol, ONLY: is_oce
     USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlstt, rlvtt, rtt
     USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2
@@ -29,23 +29,24 @@
     ! pression au milieu de chaque couche (en Pa)
 
     REAL, intent(in):: ksta, ksta_ter
-    REAL, intent(in):: ts(:) ! (klon) temperature du sol (en Kelvin)
+    REAL, intent(in):: ts(:) ! (knon) temperature du sol (en Kelvin)
     REAL, intent(in):: rugos(:) ! (klon) longeur de rugosite (en m)
     REAL, intent(in):: u(:, :), v(:, :) ! (klon, klev) wind
     REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K)
     real, intent(in):: q(:, :) ! (klon, klev) vapeur d'eau (kg/kg)
-    real, intent(in):: qsurf(:) ! (klon) 
-    REAL, intent(out):: coefm(:, :) ! (knon, klev) coefficient, vitesse
+    real, intent(in):: qsurf(:) ! (knon) 
+    REAL, intent(out):: coefm(:, 2:) ! (knon, 2:klev) coefficient, vitesse
 
-    real, intent(out):: coefh(:, :) ! (knon, klev) 
+    real, intent(out):: coefh(:, 2:) ! (knon, 2:klev) 
     ! coefficient, chaleur et humidité
 
+    real, intent(out):: ycdragm(:), ycdragh(:) ! (knon)
+
     ! Local:
 
     INTEGER knon ! nombre de points a traiter
-
-    INTEGER itop(size(coefm, 1))
-    ! (knon) numero de couche du sommet de la couche limite
+    INTEGER itop(size(coefm, 1)) ! (knon) numero de couche du sommet
+                                 ! de la couche limite
 
     ! Quelques constantes et options:
 
@@ -133,8 +134,7 @@
        z1(i) = zgeop(i, 1)
     ENDDO
 
-    CALL clcdrag(klon, knon, nsrf, .false., u1, v1, t1, q1, z1, ts, qsurf, &
-         rugos, coefm(:, 1), coefh(:, 1)) 
+    CALL clcdrag(nsrf, u1, v1, t1, q1, z1, ts, qsurf, rugos, ycdragm, ycdragh) 
 
     ! Calculer les coefficients turbulents dans l'atmosphere