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

fix for mass balance tracer in fv

parent 8fa3dd1b
No related branches found
No related tags found
1 merge request!251fix for mass balance tracer in fv
Latest changes Latest changes
============== ==============
TELEMAC-2D: Fix mass balance issue with finite volumes and tracer
MASCARET: Fix Bazin's friction Law MASCARET: Fix Bazin's friction Law
KHIONE: Addition of a new formula to compute under ice cover friction KHIONE: Addition of a new formula to compute under ice cover friction
......
...@@ -117,7 +117,12 @@ ...@@ -117,7 +117,12 @@
ENDIF ENDIF
! !
! FINAL BALANCE OF TRACER ! 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 ENDDO
! !
IF(INCLUS(COUPLING,'GAIA').AND. IF(INCLUS(COUPLING,'GAIA').AND.
......
...@@ -1368,6 +1368,7 @@ ...@@ -1368,6 +1368,7 @@
FLUTENT = 0.D0 FLUTENT = 0.D0
IF(DEBUG.GT.0) WRITE(LU,*) 'CALLING BILANT1' IF(DEBUG.GT.0) WRITE(LU,*) 'CALLING BILANT1'
DO ITRAC=1,NTRAC DO ITRAC=1,NTRAC
MASSOU(ITRAC) = 0.D0
CALL BILANT1(H,UCONV,VCONV,HPROP,T3,T4,T5,T6, CALL BILANT1(H,UCONV,VCONV,HPROP,T3,T4,T5,T6,
& LT,TOTAL_ITER,ENTET,MASKTR, & LT,TOTAL_ITER,ENTET,MASKTR,
& T%ADR(ITRAC)%P,TN%ADR(ITRAC)%P,TETAT, & 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