open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone.
I am a beginner in FORTRAN language.
I made mesh transformation type 3 in my 3D simulations, also I initialized my water temperature.
The code that I did was in the function of the plane, any suggestion to do the same but with the deep.
Thank you very much
! INITIALISES TRACERS
!
IF(NTRAC.NE.0) THEN
DO ITRAC=1,NTRAC
CALL OS('X=C ',X=TA%ADR(ITRAC)%P,C=TRAC0(ITRAC))
DO I=1,NPOIN2
DO IPLAN=1,NPLAN
J=I+((IPLAN-1)*NPOIN2)
IF(IPLAN.EQ.10.D0) THEN
TA%ADR(1)%P%R(J)=25.D0
ELSEIF(IPLAN.eq.9) THEN
TA%ADR(1)%P%R(J)=23.D0
ELSEIF(IPLAN.eq.8) THEN
TA%ADR(1)%P%R(J)=22.D0
ELSEIF(IPLAN.eq.7) THEN
TA%ADR(1)%P%R(J)=21.D0
ELSEIF(IPLAN.eq.6) THEN
TA%ADR(1)%P%R(J)=20.D0
ELSEIF(IPLAN.eq.5) THEN
TA%ADR(1)%P%R(J)=19.D0
ELSEIF(IPLAN.eq.4) THEN
TA%ADR(1)%P%R(J)=18.D0
ELSEIF(IPLAN.eq.3) THEN
TA%ADR(1)%P%R(J)=15.D0
ELSEIF(IPLAN.eq.2) THEN
TA%ADR(1)%P%R(J)=15.D0
ELSEIF(IPLAN.eq.1) THEN
TA%ADR(1)%P%R(J)=15.D0
ENDIF
ENDDO
ENDDO
ENDDO
ENDIF
!

My Fortran file

Hello,

Do you mean implementing an initial condition depending on the depth or elevation?

You can have a look at what is done for 3 TELEMAC-3D examples with FORTRAN FILE containing USER_CONDI3D_TRAC subroutine:
- V,
- bump_static,
- static_wind.

Hope this helps,

Chi-Tuan

Dear Chi-Tuan
I was forced to changed my TELEMAC version (from v7p2r3-v7p3r1) I am doing my real case, it is a 3D reservoir. I realized that in the new versions the subroutines were split.

1. I ran hydrodynamics with mesh type 2, with 10 planes, for now, the depth is approx 140 meters in the dam. It works!








2. I ran temperature in the real case, as classical adv-dif, not coupling with WQ yet, and how I made a test case before and I know for obtaining thermal stratification it is necessary to modify initial conditions.
So, I have two FORTRAN file, but the one for the tracer initialization condition is not working.





Thank you very much

Hello,

You should have a deeper look at the error messages. I think your FORTRAN file cannot compile:
IF(IPLAN.EQ.10.D0) THEN
IPLAN is declared as integer and you compare it to a float.

When getting an error while running a module of the TELEMAC-MASCARET system, please read the messages carefully, they may help. The error may not be at the last line, sometimes a few lines before, but read it from the beginning if you do not see the problem.

Hope this helps,

Chi-Tuan