open TELEMAC-MASCARET The mathematically superior suite of solvers

Could any one tell me how to prescribe the tracer value of all the surface node? I think this should be accomplished by modified the fortran code. I want to locate the paragraph of the code dealing it and the variable storing the value of the tracer concentration.

Greetings,

if I understand well, you want to prescribe the initial value of a specific tracer on the top layer of your TELEMAC-3D calculation.

To do so, I would suggest you to use the USER_CONDI3D_TRAC fortran file to setup your initial conditions.

You could look into the /examples/waqtel/waq3d_micropol/user_fortran_sed/ user_condi3d_trac.f file where an initial concentration of 1.0 was set to the top layer with the following paragraph:




Zied

Dear Zied,
Thank you for your reply. I think I did not express clearly. I meant that I need to set the value of tracer concentration on the surface every time step. could you give me a hand?

Unfortunately I am not sure about the way to do it at all time steps.

You could maybe look at the source_trac.f or user_bord3d.f files to impose source terms at the top layer, but I am not sure if this is the best way to go.

Sorry I cannot be of further help,

Zied

Zied,

Thank you for your help. I think Jean-Paul may have given a direction. Thank you for all of you.

Bert

Good afternoon,
I guess that it is possible to do it with Fortran subroutines as described by Zied.
Another way, if you find it too complicated with the subroutines might be to use TelApy python wrapper to interact with the code through Application Programming Interface. If you are not familiar with it you can found the user manual and some examples in the notebook.

You can directly change the state variables or retrieve them during the computation.

I would follow these steps:
1) Retrieve all the nodes number from the surface
2) Run the simulation with a loop one time-step at a time
3) At each time step set the tracer value for all nodes at the free-surface (look into set set_array functions or their parallel implementation mpi_set, mpi_set_array)

Good luck,

Jean-Paul

Jean-Paul,
Thank you for your reply. Using the API sounds good. I am struggling in this way.

Hello Bert,

You can read the following topic:
www.opentelemac.org/index.php/kunena/21-...g-in-time-telemac-3d

Anyway, you may encounter the same issue as Juli. Prescribing a Dirichlet condition for surface temperature may be too stiff.

Hope this helps,

Chi-Tuan

Chi-Tuan,

I think the post is very useful. Thank you very much.

Best wishes,

Bert