open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear everyone:
I want to carry out a hydrodynamic simulation of a lake. The lake only have one inlet, which is assigned as a water level boundary. However, the inlet has a gate and I want to simulate gate operation. The operation rule is simple: If the boundary water level is bigger than lake level, the gate is open, Else the gate is close. Therefore, i change the bord.f as follows:
ELSEIF(NCOTE.GT.0.OR.NOMIMP(1:1).NE.' ') THEN
N=NBOR(K)
IF(NCSIZE.GT.1) N=MESH%KNOLG%I(N)
Z = SL(IFRLIQ,N)
IF(H%R(NBOR(K))<=MAX( 0.D0 , Z-ZF(NBOR(K)) )) THEN
HBOR(K) = MAX( 0.D0 , Z-ZF(NBOR(K)) )
H%R(NBOR(K))=HBOR(K)
ELSE

ENDIF
However, the simulation collapsed. Could anyone tell me how to fix this problem.
The case file and fortran file are attached

Best wieshes.

The boundary.liq and error information is attached

Dear everyone: the model is working now. I just change the time step to be 0.5s. Now the model is running. But the results may be not correct. Because the water still exit the domain. Could anyone give me some suggestions to change liquid boundary to wall boundary. Realy need some help, Best wishes.

Hello,

I would suggest you to read the TELEMAC-2D user manual available here:
wiki.opentelemac.org/doku.php?id=documentation_v8p4r0
or directly here:
gitlab.pam-retd.fr/otm/telemac-mascaret/...emac2d_user_v8p4.pdf
in particular section 4.2 "Prescribing boundary conditions"

Maybe you can also read the user manual of the mesh generator you use to build your geometry file (and boundary file?).

Hope this helps,

Chi-Tuan

Hi,
Sorry about that, but I'm complitly new in fortran .... I read your question and if I found tour answer I'll share it !

Until this happy end, could you tell me what's "NCSIZE" and where it's define in the fortran file ?

Thanks a lot,

Laurie

I think I found !

In the TELEMAC v7p0 user-manual :

FR version :

Noter que pour la version python de TELEMAC, le nombre de processeurs est donné comme un argument pour la commande de lancement et non pas comme un mot-clé codé en dur dans le fichier des paramètres (par exemple telemac2d.py --ncsize=4 cas.txt lancera TELEMAC-2D sur 4 processeurs).

EN version :
Note that for the python version of TELEMAC, the number of processors is given as an argument to the launch command and not as a hard-coded keyword in the parameter file (e.g. telemac2d.py --ncsize=4 cas.txt will launch TELEMAC-2D on 4 processors).

I understand I don't need to define it, because it should be a variable already defined somewhere in TELEMAC2D.