open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

For a project at work I have to make a model in which hydrodynamics, morphodynamics and waves are combined. The Telemac3d-Sisyphe coupling works fine, but I am having some troubles with the telemac3d-Tomawac coupling. Of course this coupling must work before I even try to combine all three modules.

I ran a single tomawac case with a constant current taken out of a selafin file from the hydrodynamic runs. This worked. However from the moment I try to run a coupled run it is suddenly stuck at the second time Tomawac is called.

When I ran the simulation with DEBUGGER = 1, I noticed the problem occurs when PREPRO is called. I added some debug information in that script and added it to my user-defined fortran file. In this way I found out the run stalls when CHARAC.F is called for the first frequency.

I can't believe there is something wrong with the code so I'm putting the question on the forum before I continue to dig in even more.

Could anyone check my .cas files and .sortie file for any discrepancies?

Thank you very much!

Kind regards

Tim

And here the end of the .sortie file

ITERATION 20 TIME 1 D 1 H 51 MN 0.0000 S ( 93060.0000 S)
================================================================================

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D_POS OPTION: 13 10 ITERATIONS
MURD3D_POS OPTION: 13 10 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GRACJG (BIEF) : 43 ITERATIONS, RELATIVE PRECISION: 0.9782214E-06
POSITIVE DEPTHS OBTAINED IN 10 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

MASS BALANCE

WATER
MASS AT THE PREVIOUS TIME STEP : 0.1532601E+10
MASS AT THE PRESENT TIME STEP : 0.1532573E+10
MASS LEAVING THE DOMAIN DURING THIS TIME STEP : 27877.85
ERROR ON THE MASS DURING THIS TIME STEP : -0.5324295E-05
FLUX BOUNDARY 1: -9292.618 M3/S ( >0 : ENTERING <0 : EXITING )

================================================================================
TEMPS : 93123.0000 SECONDES 1552IEME ITERATION

---

APPEL DE LIMWAC
RETOUR DE LIMWAC
APPEL DE ECRETE
RETOUR DE ECRETE
APPEL DE CORMAR
RETOUR DE CORMAR
APPEL DE INIPHY
RETOUR DE INIPHY

---

CALCUL DU CHAMP CONVECTEUR ET REMONTEE DES
CARACTERISTIQUES
APPEL DE PREPRO 2
PREPRO IS CALLED
CALL CONW4D 1
EXIT CONW4D 1
CALL CONW4D 2
EXIT CONW4D 2
CALL CONW4D 3
EXIT CONW4D 3
CALL CONW4D 4
EXIT CONW4D 4
CALL CONW4D 5
EXIT CONW4D 5
CALL CONW4D 6
EXIT CONW4D 6
CALL CONW4D 7
EXIT CONW4D 7
CALL CONW4D 8
EXIT CONW4D 8
CALL CONW4D 9
EXIT CONW4D 9
CALL CONW4D 10
EXIT CONW4D 10
CALL CHARAC 1

The simulations just stays there and doesn't continue