--- trunk/Sources/phylmd/cltracrn.f	2015/07/16 17:39:10	156
+++ trunk/Sources/phylmd/cltracrn.f	2017/11/07 10:23:25	232
@@ -4,9 +4,9 @@
 
 contains 
 
-  SUBROUTINE cltracrn(itr, dtime, u1lay, v1lay, coef, t, ftsol, pctsrf, tr, &
-       trs, paprs, pplay, delp, masktr, fshtr, hsoltr, tautr, vdeptr, lat, &
-       d_tr, d_trs)
+  SUBROUTINE cltracrn(itr, dtime, u1lay, v1lay, coef, cdragh, t, ftsol, &
+       pctsrf, tr, trs, paprs, pplay, delp, masktr, fshtr, hsoltr, tautr, &
+       vdeptr, lat, d_tr, d_trs)
 
     ! From phylmd/cltracrn.F, version 1.2 2005/05/25 13:10:09
 
@@ -30,11 +30,11 @@
 
     REAL, intent(in):: dtime
     ! dtime----input-R- intervalle de temps (en second)
-    REAL u1lay(klon), v1lay(klon)
-    ! u1lay----input-R- vent u de la premiere couche (m/s)
-    ! v1lay----input-R- vent v de la premiere couche (m/s)
-    REAL coef(klon, klev)
+    REAL, intent(in):: u1lay(klon), v1lay(klon) ! vent de la premiere
+                                                ! couche (m/s)
+    REAL coef(:, 2:) ! (klon, 2:klev)
     ! coef-----input-R- le coefficient d'echange (m**2/s) l>1
+    real cdragh(:) ! klon
     REAL, intent(in):: t(klon, klev) ! temperature (K)
     real, intent(in):: ftsol(klon, nbsrf), pctsrf(klon, nbsrf) 
     ! ftsol----input-R- temperature du sol (en Kelvin)
@@ -116,7 +116,7 @@
     ENDDO
 
     DO i = 1, klon
-       zx_coef(i, 1) = coef(i, 1) &
+       zx_coef(i, 1) = cdragh(i) &
             * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) &
             * pplay(i, 1)/(RD*t(i, 1))
        zx_coef(i, 1) = zx_coef(i, 1) * dtime*RG