open TELEMAC-MASCARET The mathematically superior suite of solvers

Good afternoon,

I have been trying to implement spatio-temporally varying boundary conditions in Telemac3D v7p1. In due course, I wish to have H,U,V all varying (these come from another model), but to start with I am just trying to get spatially varying H.

I'm new to both TELEMAC and Fortran, so I've followed the instructions given in a previous post:post here

The only instruction I haven't followed is:
- Do comment the line DEJA=.TRUE. (in the IF(.NOT.DEJA) THEN…) in SL( )
because that section of code has been removed.

The code compiles fine, however, when I try to run my model it seems to hang up. I get no error code but after 'telemac3d-v7p1 is displayed 'in dashes', nothing happens.

Would anyone have any suggestions as to where I've gone wrong or be kind enough to send me a version of the code thats been modified for variable boundaires that I can look at? I've attached my modified bord3.f, sl3.f and liquid boundary files (note that at present in the .liq file the water levels are not varying across the boundary in this simple test case).

Thanks in advance,
Iain

Hi
SL3 is wrong in the modified part.
I advise you to quickly start to read some fortran programming tutorial.
The right part of code should be :

---

K = N
IF (NCSIZE.GT.1) K = MESH2D%KNOLG%I(N)
FCT(1:4)='ETA('
IF(K.LT.10) THEN
WRITE(FCT(5:5),FMT='(I1)') K
FCT(6:11)=') '
ELSEIF(K.LT.100) THEN
WRITE(FCT(5:6),FMT='(I2)') K
FCT(7:11)=') '
ELSEIF(K.LT.1000) THEN
WRITE(FCT(5:7),FMT='(I3)') K
FCT(8:11)=') '
ELSEIF(K.LT.10000) THEN
WRITE(FCT(5:8),FMT='(I4)') K
FCT(10:11)=') '
ELSEIF(K.LT.100000) THEN
WRITE(FCT(5:9),FMT='(I5)') K
FCT(10:11)=') '

---

I didn't check bord3d for the moment

Regards

Having tried to understand the code a bit better and looking through some other examples, I thought I'd got the variable water level working, however, while it compiles and the model runs, my results are non-sensical. It seems like while the .liq file is being read (according to the sortie file), the boundaries are held at zero not read from the file.

Does anyone have any ideas? I've included my fortranfile, my .cas file, my .liq file, the 2-d results file, and my sortie file. I've been running a simple channel case with default parameters.
Do I need to modify any other files beyond bord3, sl3 (and vit3,q3,tr3 for the length of the WHAT variable), read_fic_frliq?

Thanks in advance.

Hello,

I looked but do not see any obvious mistake (which is not a proof...). You could add a print after the line :

HBOR%R(K) = SL3(ICOT,AT,N,INFOGR)-ZF%R(IPOIN2)

in bord3f.f to see if this line is reached and what value of HBOR is given.

With best regards,

Jean-Michel Hervouet

Only three files uploaded so the others attached again.

Also, I only modified bord3 in the section "PRESCRIBED ELEVATION GIVEN IN STEERING FILE (NCOTE<>0)," does it need to be modified elsewhere? it didnt look like it to me

Thanks, I realised in the end it was simply a typo in my .liq file