NEWS.txt

 Latest changes
 ==============
+
+TELEMAC-2D: Fix mass balance issue with finite volumes and tracer
+Initialization of MASSOU for finite volumes kernel
+Correction of mass balance computation for parellel runs
+
+PYTHON: TelemacFile.print_info now correctly displays the number of interface
+points if any, as well as the number of boundary points and elements regardless
+of the presence or not of a boundary conditions file.
+
 MASCARET: Fix Bazin's friction Law
 
 KHIONE: Addition of a new formula to compute under ice cover friction
@@ -18,10 +27,6 @@ TELEMAC-3D: New keywords SIGMA DISTRIBUTION OF PLANES and
 PRESCRIBED ELEVATION OF PLANES to set vertical planes directly in the steering
 file
 
-PYTHON: TelemacFile.print_info now correctly displays the number of interface
-points if any, as well as the number of boundary points and elements regardless
-of the presence or not of a boundary conditions file.
-
 TOMAWAC: fix time series for wave directions
 
 MASCARET: Reference Documentation has been added

sources/telemac2d/fv_balance.f

@@ -117,7 +117,12 @@
             ENDIF
 !
 !           FINAL BALANCE OF TRACER
-            MASSOU(ITRAC) = MASSOU(ITRAC) + DT*SMTR%ADR(ITRAC)%P%R(I)
+            IF(NCSIZE.GT.1) THEN
+              MASSOU(ITRAC) = MASSOU(ITRAC) + DT*SMTR%ADR(ITRAC)%P%R(I)
+     &                      * MESH%IFAC%I(I)
+            ELSE
+              MASSOU(ITRAC) = MASSOU(ITRAC) + DT*SMTR%ADR(ITRAC)%P%R(I)
+            ENDIF
           ENDDO
 !
           IF(INCLUS(COUPLING,'GAIA').AND.

sources/telemac2d/telemac2d_init.F

@@ -1368,6 +1368,7 @@
             FLUTENT = 0.D0
             IF(DEBUG.GT.0) WRITE(LU,*) 'CALLING BILANT1'
             DO ITRAC=1,NTRAC
+            MASSOU(ITRAC) = 0.D0
             CALL BILANT1(H,UCONV,VCONV,HPROP,T3,T4,T5,T6,
      &                   LT,TOTAL_ITER,ENTET,MASKTR,
      &                   T%ADR(ITRAC)%P,TN%ADR(ITRAC)%P,TETAT,