open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

I am trying to put a water quality model under T3D framework, I have two non-conservative tracers T1 and T2. My questions are as follows

a) I don't want T2 to be advected or diffused it just decays in time but T1 will get advected and diffused with the flow, how can I do some thing like this ?

b) Also if some one could remind me which is the user-subroutine that is called at every time step that can be used to put in some user calculations.

Thanks for all the help in advance.

Kind regards,
Sumit

Hello,

a) This is not implemented for a single tracer among others. I would suggest that you change subroutine diff3d.f for your case: on line 259 or 261, depending on the case, you put FN instead of FC (thus cancelling the effect of advection). To do this only for one tracer use FN%NAME which will be TAN1, TAN2, etc., depending on the tracer rank (in version 6.3 it will be TRN1, TRN2, we changed to avoid a conflict with the Tangent function).

b) the name is utimp.f

With best regards,

Jean-Michel Hervouet

Thanks a ton JMH

kind regards,
Sumit