diff --git a/NEWS.txt b/NEWS.txt index 5a4f247d10d2738ac846259c04bc1818508eca20..88c525ed5e5ddcb57a13a9cf1802dc6b0d31e5f5 100644 --- a/NEWS.txt +++ b/NEWS.txt @@ -1,5 +1,7 @@ Latest changes ============== +TELEMAC-2D: Fix mass balance issue with finite volumes and tracer + MASCARET: Fix Bazin's friction Law KHIONE: Addition of a new formula to compute under ice cover friction diff --git a/sources/telemac2d/fv_balance.f b/sources/telemac2d/fv_balance.f index 90efc8819339dd51b471cb588a1b9004dbfbb0c1..dc836ad4bdc577cb434509a052e271280a6715c3 100644 --- a/sources/telemac2d/fv_balance.f +++ b/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. diff --git a/sources/telemac2d/telemac2d_init.F b/sources/telemac2d/telemac2d_init.F index f4c49cb669e826c49ef39501366b477887ca4658..7015132233031b4e059869b1759ce73c51b792de 100644 --- a/sources/telemac2d/telemac2d_init.F +++ b/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,