open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to model sieche that occur in a reservoir during an earthquake. For this i need to add peak horizontal velocity to the actual velocity components. If I change condin, only the initial condition I can alter. I would like to know in which subroutine I have to update this additional horizontal velocity.
Pls help in this regards


Thank you

Is there a way to add velocity to the computed velocity inside time loop?

Hi
You could add a source point which allow you to add velocity (and mass) locally.
Otherwise you should add it by programming it (prosou.f ?)

Regards

Hi,
Programming Prosou.f, will it work with Telemac3d? As you said, adding velocity as a source, we can add locally, but I want to add velocity at all nodes at an instant.

Maybe as an initial condition?

Initial condition means, I tried programming condiuvw.f and it worked well. But I want to add intermittently inside the time loop.

is it right to program velres.f?

Hello,

To add a source term for momentum equations in TELEMAC-3D, you can have a look at TRISOU subroutine.

Hope this helps,

Chi-Tuan