open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

You can look at subroutine trisou in library telemac3d to see how things are done (at the end of the subroutine). Though it is not programmed yet like this, I would say that when the source is in fact a sink, we should not add a modification of the right-hand sides of velocity equations CV1 and CV2. This is so for tracers, as when a tracer exits through a sink, its value is not a prescribed value, only the local value in the domain. So you can probably fix your problem by changing, in formulas updating CV1 and CV2 at the end of subroutine trisou :

QSCE(I) into MAX(QSCE(I),0.D0)

however it could be that some users will tell us that they do want this term.

With best regards,

Jean-Michel Hervouet