open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I currently have a well-working model, allowing to model temperature diffusion within a fishing port (also using WAQTEL for air/water exchanges). For the moment, i consider a varying over time temperature/tide on liquid boundaries, from a LIQUID BOUNDARIES FILE, and i consider a constant temperature profil on the vertical.

I would like to know if it is possible to implement a vertical profile of temperatures, linked to a boundary condition file with different sea surface temperature and sea bottom temperature ? For exemple :

- Implement a sea surface temperature on my liquid boundaries, varying over time (ideally coming from a liquid boundaries file)
- Implement a sea bottom temperature on my liquid boundaries, also varying over time (also coming from a liquid boundaries file)
- Then make a linear stratification of the temperature depending on the depth

I tried to organize something using user_bord_3d and/or user_tra_prof_z, but I don't see how I can call two different temperatures, which are varying over time.

So my questions are :

- Is using user_bord_3d and/or user_tra_prof_z is the way to go?
- If yes, can I create a specific ascii/txt file for reading the temperature on the boundaries? How does it work?
- If no, is there any other solution to do it?

Hoping that I am making myself clear enough.

In advance thank you.

Regards,

Robin

Hi
The approach sounds good.
Possibilities exists to declare user formatted file in the steering file that could be used in your fortran development (look into the manuals)
You could also look at telemac2d source files read_fic_frliq.f to have an example of the management of such files.

Hope this helps

Hi,

Thank you for your answer.

I'm a newbie using TELEMAC (also using Fortran). Just to be sure of what to do :

- I must write a subroutine based on read_fic_frliq to read a FORMATTED DATA FILE which contains the temperature on the surface and on the bottom varying over time ;
- Then adapt the user_bord_3d, which will call the previous subroutine I will have written, to interpolate the temperatures at each plan of my 3D mesh ?

In advance thank you.

Hi
Yes this is the idea. Check if modification is also requested in user_tra_prof_z ...

Regards