open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

a question which I often asked myself,
how could I apply a stage hydrograph as internal boundary condition in Telemac-2D?
To clarify the question via an example:
The mesh consists of two reservoirs and the upstream reservoir releases water according to a water level draw down rule of 0.5 m/h. How can I model such a scenario? From a numerical point of view it should be possible?

Another maybe stupid question, how can I check if a subroutine is called every time step or only at the beginning of the simulation?

I would be very glad for any hints!
Clemens

Hi Clemens,

I am not sure to understand what you mean by internal boundary condition.

If you just want to apply a water level lowering in the upstream reservoir, you can use a time dependant stage condition on your liquid boundary using a liquid boundary file activated and called by your steering file LIQUID BOUNDARIES FILE.

If it is the link between the two reservoirs that is time varying, I dont really see how to do ?
How are the two reservoir linked ?
weir or a feature included in your mesh ?

To see the sequence of calls made by telemac you can add the keyword DEBUGGER = 1.

Hi Clemens

For your first question, I think it's possible with some manual programming. I think the easiest solution is to modify the weir subroutine where the discharge between two boundary is already implemented in order to compute the discharge not according the weir law but just with the emptying rule. if you want to follow a rule like 0.5m/h, you should probably fins a solution to compute the volume in your upstream basin...

For the second one, different options could be use:

• use the debugger mode (I hope all the call are well marked for the debugger mode)
• place the subroutine in your own fortran file and add a message in it (not really beautiful but efficient)
• look in the source code to find the different call of the subroutine and check if there is a conditionnal call with a test like (IF LT.LT.1) or something similar (LT is the iteration number)

Hope this helps
regards

Hello,

thank you very much for the answers!

Best regards,
Clemens