Skip to content
Snippets Groups Projects
Commit 4e06cdd8 authored by Florent Taccone's avatar Florent Taccone
Browse files

Merge branch 'fix_tracer_fv' into 'main'

fix for mass balance tracer in fv

See merge request !251
parents 8fa3dd1b 51f4b999
No related branches found
No related tags found
1 merge request!251fix for mass balance tracer in fv
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
......
......@@ -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.
......
......@@ -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,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment