--- trunk/Sources/phylmd/physiq.f	2017/10/16 13:04:05	226
+++ trunk/Sources/phylmd/physiq.f	2017/11/13 12:12:41	242
@@ -44,13 +44,14 @@
     USE indicesol, ONLY: clnsurf, epsfra, is_lic, is_oce, is_sic, is_ter, &
          nbsrf
     USE ini_histins_m, ONLY: ini_histins, nid_ins
+    use lift_noro_m, only: lift_noro
     use netcdf95, only: NF95_CLOSE
     use newmicro_m, only: newmicro
     use nr_util, only: assert
     use nuage_m, only: nuage
     USE orbite_m, ONLY: orbite
     USE ozonecm_m, ONLY: ozonecm
-    USE phyetat0_m, ONLY: phyetat0, rlat, rlon
+    USE phyetat0_m, ONLY: phyetat0
     USE phyredem_m, ONLY: phyredem
     USE phyredem0_m, ONLY: phyredem0
     USE phytrac_m, ONLY: phytrac
@@ -183,9 +184,7 @@
     REAL cdragh(klon) ! drag coefficient pour T and Q
     REAL cdragm(klon) ! drag coefficient pour vent
 
-    ! Pour phytrac :
-    REAL ycoefh(klon, llm) ! coef d'echange pour phytrac
-    REAL yu1(klon), yv1(klon) ! vent dans la premi\`ere couche
+    REAL ycoefh(klon, 2:llm) ! coef d'echange pour phytrac
 
     REAL, save:: ffonte(klon, nbsrf)
     ! flux thermique utilise pour fondre la neige
@@ -249,8 +248,9 @@
 
     REAL flux_q(klon, nbsrf) ! flux turbulent d'humidite à la surface
     REAL flux_t(klon, nbsrf) ! flux turbulent de chaleur à la surface
-    REAL flux_u(klon, nbsrf) ! flux turbulent de vitesse u à la surface
-    REAL flux_v(klon, nbsrf) ! flux turbulent de vitesse v à la surface
+
+    REAL flux_u(klon, nbsrf), flux_v(klon, nbsrf)
+    ! tension du vent (flux turbulent de vent) à la surface, en Pa
 
     ! Le rayonnement n'est pas calcul\'e tous les pas, il faut donc que
     ! les variables soient r\'emanentes.
@@ -370,7 +370,6 @@
 
     REAL zustrdr(klon), zvstrdr(klon)
     REAL zustrli(klon), zvstrli(klon)
-    REAL zustrph(klon), zvstrph(klon)
     REAL aam, torsfc
 
     REAL ve_lay(klon, llm) ! transport meri. de l'energie a chaque niveau vert.
@@ -582,9 +581,6 @@
     evap = - sum(flux_q * pctsrf, dim = 2)
     fder = dlw + dsens + devap
 
-    yu1 = u_seri(:, 1)
-    yv1 = v_seri(:, 1)
-
     DO k = 1, llm
        DO i = 1, klon
           t_seri(i, k) = t_seri(i, k) + d_t_vdf(i, k)
@@ -935,9 +931,9 @@
           ENDIF
        ENDDO
 
-       CALL lift_noro(klon, llm, dtphys, paprs, play, rlat, zmea, zstd, zpic, &
-            itest, t_seri, u_seri, v_seri, zulow, zvlow, zustrli, zvstrli, &
-            d_t_lif, d_u_lif, d_v_lif)
+       CALL lift_noro(dtphys, paprs, play, zmea, zstd, zpic, itest, t_seri, &
+            u_seri, v_seri, zulow, zvlow, zustrli, zvstrli, d_t_lif, &
+            d_u_lif, d_v_lif)
 
        ! Ajout des tendances :
        DO k = 1, llm
@@ -949,29 +945,16 @@
        ENDDO
     ENDIF
 
-    ! Stress n\'ecessaires : toute la physique
-
-    DO i = 1, klon
-       zustrph(i) = 0.
-       zvstrph(i) = 0.
-    ENDDO
-    DO k = 1, llm
-       DO i = 1, klon
-          zustrph(i) = zustrph(i) + (u_seri(i, k) - u(i, k)) / dtphys &
-               * zmasse(i, k)
-          zvstrph(i) = zvstrph(i) + (v_seri(i, k) - v(i, k)) / dtphys &
-               * zmasse(i, k)
-       ENDDO
-    ENDDO
-
-    CALL aaam_bud(rg, romega, rlat, rlon, pphis, zustrdr, zustrli, zustrph, &
-         zvstrdr, zvstrli, zvstrph, paprs, u, v, aam, torsfc)
+    CALL aaam_bud(rg, romega, pphis, zustrdr, zustrli, &
+         sum((u_seri - u) / dtphys * zmasse, dim = 2), zvstrdr, &
+         zvstrli, sum((v_seri - v) / dtphys * zmasse, dim = 2), paprs, u, v, &
+         aam, torsfc)
 
     ! 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, yu1, yv1, ftsol, &
-         pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, tr_seri, &
-         zmasse, ncid_startphy)
+         mfd, pde_u, pen_d, ycoefh, 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)