--- trunk/Sources/dyn3d/fluxstokenc.f	2015/04/29 15:47:56	134
+++ trunk/Sources/dyn3d/fluxstokenc.f	2016/09/01 10:30:53	207
@@ -8,13 +8,15 @@
 
     ! Author: F. Hourdin
 
-    USE histwrite_m, ONLY: histwrite
-    USE dimens_m, ONLY: jjm, llm, nqmx
-    USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1, jjp1
     USE comgeom, ONLY: aire
+    USE dimens_m, ONLY: jjm, llm
+    use groupe_m, only: groupe
+    USE histwrite_m, ONLY: histwrite
+    use initfluxsto_m, only: initfluxsto
+    USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1
     USE tracstoke, ONLY: istdyn, istphy
 
-    REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
+    REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
     REAL, intent(in):: masse(ip1jmp1, llm)
     real, intent(in):: phi(ip1jmp1, llm)
     real, intent(in):: teta(ip1jmp1, llm)
@@ -22,7 +24,7 @@
     REAL, intent(in):: time_step
     INTEGER, INTENT (IN):: itau
 
-    ! Variables local to the procedure:
+    ! Local:
 
     REAL, SAVE:: pbaruc(ip1jmp1, llm), pbarvc(ip1jm, llm)
     REAL, SAVE:: massem(ip1jmp1, llm)
@@ -37,7 +39,7 @@
     !-------------------------------------------------------------
 
     IF (itau == 0) THEN
-       CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, &
+       CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, &
             fluxid, fluxvid, fluxdid)
        CALL histwrite(fluxid, 'phis', 1, phis)
        CALL histwrite(fluxid, 'aire', 1, aire)
@@ -54,7 +56,7 @@
        CALL initial0(ijmllm, pbarvc)
     END IF
 
-    !   accumulation des flux de masse horizontaux
+    ! Accumulation des flux de masse horizontaux
     DO l = 1, llm
        DO ij = 1, ip1jmp1
           pbaruc(ij, l) = pbaruc(ij, l) + pbaru(ij, l)
@@ -66,12 +68,12 @@
        END DO
     END DO
 
-    !   selection de la masse instantannee des mailles avant le transport.
+    ! S\'election de la masse instantan\'ee des mailles avant le transport.
     IF (itau == 0) massem = masse
 
     IF (mod(itau + 1, istdyn) == 0) THEN
-       ! on advecte a ce pas de temps
-       !    normalisation
+       ! On advecte \`a ce pas de temps
+       ! normalisation
        DO l = 1, llm
           DO ij = 1, ip1jmp1
              pbaruc(ij, l) = pbaruc(ij, l)/float(istdyn)
@@ -83,11 +85,10 @@
           END DO
        END DO
 
-       !   traitement des flux de masse avant advection.
-       !     1. calcul de w
-       !     2. groupement des mailles pres du pole.
-
-       CALL groupe(massem, pbaruc, pbarvc, pbarug, pbarvg, wg)
+       ! Traitement des flux de masse avant advection.
+       ! 1. Calcul de w
+       ! 2. Groupement des mailles pr\`es du p\^ole.
+       CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg)
 
        CALL histwrite(fluxid, 'masse', itau, massem)
        CALL histwrite(fluxid, 'pbaru', itau, pbarug)