Thank you very much!
In addition, I have encountered a problem. ve. I want to solve an advection–diffusion equation with a unit (negative) source term in the control region.
In a first step, I use the hydrodynamic model to generate the velocity and turbulence fields
as usual, and the intermediate results are stored. Then, in a second step,the equation I want to solve is integrated backward in time using the hydrodynamic fields computed in the first step.
So I read the velocity and turbulence fields in the telemac2d.f, and set the time step as negative,and add a unit (negative) source term via TIMP%ADR(1)%P%R(I)=-1.D0. The results is negative, and I wonder if I was wrong?