--- trunk/phylmd/clvent.f	2014/03/05 14:57:53	82
+++ trunk/Sources/phylmd/clvent.f	2017/10/16 12:35:41	225
@@ -15,12 +15,10 @@
     USE suphec_m, ONLY: rd, rg
 
     INTEGER knon
-    REAL, intent(in):: dtime
-    ! dtime----input-R- intervalle du temps (en second)
+    REAL, intent(in):: dtime ! intervalle de temps (en s)
 
-    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, intent(in):: u1lay(:), v1lay(:) ! (knon)
+    ! vent de la premiere couche (m/s)
 
     REAL, intent(in):: coef(:, :) ! (knon, klev)
     ! Coefficient d'echange (m**2/s) multiplié par le cisaillement du
@@ -36,8 +34,10 @@
     ! delp-----input-R- epaisseur de couche (Pa)
     REAL d_ven(klon, klev)
     ! d_ven----output-R- le changement de "ven"
-    REAL flux_v(klon, klev)
-    ! flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)
+
+    REAL, intent(out):: flux_v(:) ! (knon)
+    ! (diagnostic) flux du vent à la surface, en (kg m/s)/(m**2 s)
+    ! flux_v est le flux de moment angulaire (positif vers bas)
 
     ! Local:
     INTEGER i, k
@@ -99,19 +99,11 @@
        ENDDO
     ENDDO
 
-    ! flux_v est le flux de moment angulaire (positif vers bas) dont
-    ! l'unite est: (kg m/s)/(m**2 s)
     DO i = 1, knon
-       flux_v(i, 1) = zx_coef(i, 1)/(RG*dtime) &
+       flux_v(i) = zx_coef(i, 1)/(RG*dtime) &
             *(local_ven(i, 1)*zx_alf1(i) &
             +local_ven(i, 2)*zx_alf2(i))
     ENDDO
-    DO k = 2, klev
-       DO i = 1, knon
-          flux_v(i, k) = zx_coef(i, k)/(RG*dtime) &
-               * (local_ven(i, k)-local_ven(i, k-1))
-       ENDDO
-    ENDDO
 
     DO k = 1, klev
        DO i = 1, knon