--- trunk/Sources/phylmd/phytrac.f	2017/11/02 15:47:03	227
+++ trunk/Sources/phylmd/phytrac.f	2017/11/07 10:23:25	232
@@ -77,7 +77,7 @@
 
     REAL pde_u(klon, llm) ! flux detraine dans le panache montant
     REAL pen_d(klon, llm) ! flux entraine dans le panache descendant
-    REAL coefh(klon, llm) ! coeff melange couche limite
+    REAL coefh(:, :) ! (klon, llm) coeff melange couche limite
     real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques
     REAL, intent(in):: yu1(:), yv1(:) ! (klon) vent au premier niveau
 
@@ -279,10 +279,10 @@
     DO it=1, nqmx - 2
        if (clsol(it)) then 
           ! couche limite avec quantite dans le sol calculee
-          CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, &
-               pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &
-               masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &
-               vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)
+          CALL cltracrn(it, pdtphys, yu1, yv1, coefh(:, 2:llm), coefh(:, 1), &
+               t_seri, ftsol, pctsrf, tr_seri(:, :, it), trs(:, it), paprs, &
+               pplay, delp, masktr(1, it), fshtr(1, it), hsoltr(it), &
+               tautr(it), vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)
           DO k = 1, llm
              DO i = 1, klon
                 tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)
@@ -297,7 +297,7 @@
              source(i) = 0. ! pas de source, pour l'instant
           ENDDO
 
-          CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &
+          CALL cltrac(pdtphys, coefh(:, 2:llm), t_seri, tr_seri(:, :, it), source, &
                paprs, pplay, delp, d_tr_cl(1, 1, it))
           DO k = 1, llm
              DO i = 1, klon