open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone,

I am doing Telemac 3D thermal stratification, I tried many things, every time I got a mixed temperature.
Getting steady-state conditions
I tried with more and fewer layers
with PPE
with Boussinesq
with all turbulence models
Coupled with waqtel
Changing COEFF. OF WIND INFLUENCE
Changing COEFFICIENT TO CALIBRATE THE ATMOSPHERE-WATER EXCHANGE MODEL
Changing diffusion coefficient
Finally, I am trying to modified free surface BC´s No Newman, but Dirichlet
How I can change in time because until now I have a constant value
Thank you very much

Hello Juli,

If you want to implement your own boundary conditions varying in time, as written in the TELEMAC-3D user manual, you have 2 possibilities:
- using the LIQUID BOUNDARIES FILE if your condition does not vary in space but only in time,
- implementing your boundary conditions in the USER_SL3 (or USER_Q3, or USER_VIT3 depending on the variable) if also varying in time but not in space,
- implementing your boundary conditions in the USER_BORD3D or the BORD3D subroutine.

Hope this helps,

Chi-Tuan

Dear Chi-Tuan
I was seeing USER_SL3, there is some examples to follow?
Thank you very much

Hello Juli,

In $HOMETEL/examples/telemac3d/vasque/user_fortran you have an example of implementing the USER_SL3 function.

Hope this helps,

Chi-Tuan

There are some examples for tracers?
Thank yo very much

Hello Juli,

To my knowledge, there is no example for tracers with USER_TR3 subroutine, but you can manage to write such a subroutine as a mix of USER_SL3 and USER_BORD3D. See in the TELEMAC-3D examples, e.g. stratification, tetra or Viollet examples.

Hope this helps,

Chi-Tuan

Hello,
I did what you suggest with example Stratification.

svn.opentelemac.org/svn/opentelemac/tags...ac3d/stratification/

I modified TR3 en v8p2r1 but is not changing with time as supposed to do.

When I used limi3d works but the problem is I don't have time in this subroutine



Any help, suggestion, I have to modify another subroutine?
Thanks in advance

Hello Juli,

In LIMI3D subroutine, as you have "USE DECLARATIONS_TELEMAC3D", you can use every variable declared in declarations_telemac3d module, such as AT for time.

To be sure that USER_TR3 is called, you have to have OKTR3(I,ITRAC) = .FALSE. in TR3 function or LIQUID BOUNDARIES FILE declared in your steering file.

Hope this helps,

Chi-Tuan

Dear Chi-Tuan,
Still, no success trying to modified surface boundary conditions for tracers.
With lateral works in BORD3d if I multiplied by AT, and with limid is taking only in the beginning.
Right now, I am sure that is reading the subroutines.
So this one is working:

I tried to do something similar but with lateral and does not work:

And for the last:

does not work either.
any idea what I am doing wrong
Thanks in advance

Dear Chi-Tuan,
I don't understand something: docs.opentelemac.org/doxydocs/v8p2r0/htm...rmics_8f_source.html
I was reading the code of thermics_3D and I was wondering why is used btabos

122 IF(atmosexch.EQ.1) THEN
123 a=(4.48d0+0.049d0*treel+2021.5d0*c_atmos*(1.d0+ww)*
124 & (1.12d0+0.018d0*treel+0.00158d0*treel**2))/lamb
125 atabos%ADR(ind_t)%P%R(ipoin2)=-a
126 btabos%ADR(ind_t)%P%R(ipoin2)= a*tair%R(ipoin2)
127 ELSEIF(atmosexch.EQ.2) THEN
128 !
129 ! SENSIBLE HEAT FLUXES
130 !
131 CALL evapo(treel,tair%R(ipoin2),ww2,patmos%R(ipoin2),
132 & hrel%R(ipoin2),ro,flux_evap,flux_sens,debevap,
133 & c_atmos,cfaer(1),cfaer(2))
134 !
135 ! LONGWAVE HEAT FLUXES
136 !
137 CALL shortrad(treel,tair%R(ipoin2),cldc%R(ipoin2),
138 & hrel%R(ipoin2),ray_atm,ray_eau)
139 !
140 ! BOUNDARY CONDITION FOR TEMPERATURE AT SURFACE
141 !
142 atabos%ADR(ind_t)%P%R(ipoin2) = 0.d0
143 btabos%ADR(ind_t)%P%R(ipoin2) = (ray_atm-ray_eau-flux_evap
144 & -flux_sens)/lamb
145 ENDIF

Thank you very much