--- trunk/Sources/phylmd/physiq.f	2017/11/06 17:20:45	229
+++ trunk/Sources/phylmd/physiq.f	2018/01/05 14:45:45	247
@@ -170,7 +170,7 @@
     REAL, save:: zval(klon) ! Minimum de l'OESM
     REAL, save:: rugoro(klon) ! longueur de rugosite de l'OESM
     REAL zulow(klon), zvlow(klon)
-    INTEGER igwd, itest(klon)
+    INTEGER ktest(klon)
 
     REAL, save:: agesno(klon, nbsrf) ! age de la neige
     REAL, save:: run_off_lic_0(klon)
@@ -184,7 +184,7 @@
     REAL cdragh(klon) ! drag coefficient pour T and Q
     REAL cdragm(klon) ! drag coefficient pour vent
 
-    REAL ycoefh(klon, llm) ! coef d'echange pour phytrac
+    REAL coefh(klon, 2:llm) ! coef d'echange pour phytrac
 
     REAL, save:: ffonte(klon, nbsrf)
     ! flux thermique utilise pour fondre la neige
@@ -571,7 +571,7 @@
          paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, &
          snow_fall, fsolsw, fsollw, frugs, agesno, rugoro, d_t_vdf, d_q_vdf, &
          d_u_vdf, d_v_vdf, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, &
-         cdragm, q2, dsens, devap, ycoefh, t2m, q2m, u10m_srf, v10m_srf, &
+         cdragm, q2, dsens, devap, coefh, t2m, q2m, u10m_srf, v10m_srf, &
          pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, trmb3, plcl, &
          fqcalving, ffonte, run_off_lic_0)
 
@@ -897,18 +897,16 @@
 
     IF (ok_orodr) THEN
        ! S\'election des points pour lesquels le sch\'ema est actif :
-       igwd = 0
        DO i = 1, klon
-          itest(i) = 0
+          ktest(i) = 0
           IF (zpic(i) - zmea(i) > 100. .AND. zstd(i) > 10.) THEN
-             itest(i) = 1
-             igwd = igwd + 1
+             ktest(i) = 1
           ENDIF
        ENDDO
 
-       CALL drag_noro(klon, llm, dtphys, paprs, play, zmea, zstd, zsig, zgam, &
-            zthe, zpic, zval, itest, t_seri, u_seri, v_seri, zulow, zvlow, &
-            zustrdr, zvstrdr, d_t_oro, d_u_oro, d_v_oro)
+       CALL drag_noro(dtphys, paprs, play, zmea, zstd, zsig, zgam, zthe, &
+            zpic, zval, ktest, t_seri, u_seri, v_seri, zulow, zvlow, zustrdr, &
+            zvstrdr, d_t_oro, d_u_oro, d_v_oro)
 
        ! ajout des tendances
        DO k = 1, llm
@@ -922,16 +920,14 @@
 
     IF (ok_orolf) THEN
        ! S\'election des points pour lesquels le sch\'ema est actif :
-       igwd = 0
        DO i = 1, klon
-          itest(i) = 0
+          ktest(i) = 0
           IF (zpic(i) - zmea(i) > 100.) THEN
-             itest(i) = 1
-             igwd = igwd + 1
+             ktest(i) = 1
           ENDIF
        ENDDO
 
-       CALL lift_noro(dtphys, paprs, play, zmea, zstd, zpic, itest, t_seri, &
+       CALL lift_noro(dtphys, paprs, play, zmea, zstd, zpic, ktest, t_seri, &
             u_seri, v_seri, zulow, zvlow, zustrli, zvstrli, d_t_lif, &
             d_u_lif, d_v_lif)
 
@@ -952,9 +948,9 @@
 
     ! Calcul des tendances traceurs
     call phytrac(julien, time, firstcal, lafin, dtphys, t, paprs, play, mfu, &
-         mfd, pde_u, pen_d, ycoefh, fm_therm, entr_therm, u(:, 1), v(:, 1), &
-         ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, &
-         tr_seri, zmasse, ncid_startphy)
+         mfd, pde_u, pen_d, coefh, cdragh, fm_therm, entr_therm, u(:, 1), &
+         v(:, 1), ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, &
+         dnwd, tr_seri, zmasse, ncid_startphy)
 
     ! Calculer le transport de l'eau et de l'energie (diagnostique)
     CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq)