open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I have a problem regarding multiple liquid boundary files, which are all open boundaries. How many liquid boundary files can actually be imposed within TELEMAC?

I want to give tidal prescribed elevations to 7 open boundaries (in total). In the cas file I have inserted the following information to account for the liquid boundaries (as described in the guidance):

LIQUID BOUNDARIES FILE ='Boundary_1.liq'
LIQUID BOUNDARIES FILE ='Boundary_2.liq'
LIQUID BOUNDARIES FILE ='Boundary_3.liq'
LIQUID BOUNDARIES FILE ='Boundary_4.liq'
LIQUID BOUNDARIES FILE ='Boundary_5.liq'
LIQUID BOUNDARIES FILE ='Boundary_6.liq'
LIQUID BOUNDARIES FILE ='Boundary_7.liq'

after running TELEMAC-2D, the message given is:

*** ALLOCATION OF USER FILES ***

- LIQUID BOUNDARIES FILE : Boundary_7.liq

THE KEY-WORD: LIQUID BOUNDARIES FILE
APPPEARS AT LEAST TWICE , THE LAST VALUE WILL BE KEPT...

The liquid boundary files are all imposed by elevation varying in space and time (tidal system), in which the boundary points are declared as prescribed elevation. For example in boundary_1.liq:

#
# boundary_1
#
T SL(1) SL(2) SL(3) SL(4) SL(5) SL(6)
s m m m m m m
100. -0.254 -0.288 -0.325 -0.371 -0.435 -0.523
200. -0.131 -0.161 -0.194 -0.235 -0.293 -0.376
300. -0.005 -0.032 -0.060 -0.096 -0.146 -0.224
400. 0.121 0.099 0.075 0.045 0.002 -0.070

SL function can be used for boundary condition which is constant in space but varies in time.

Have I prescribed the values/functions correctly?

I would really appreciate it if someone could offer me some assistance.

Many thanks,

Zheng

Hi
You could only reference one boundary file in which you describe the varaition in time af all the different boundaries (SL(1)... SL(7))
Each boundary contains a group of nodes which heve the same value imposed.

If you want to have a space varying condition you should adapt the program by manually programming in the subroutine bord.f

Hope this helps

Hi, Coulet,

Thank you very much for your replay.

First test, one liquid boundary file created and included in the cas file, which the variation in time of all the 7 boundaries (SL(1).. S(7)) and not varying in space. The boundaries.liq:

#
# Open boundaries file
# 7 boundaries managed
#
T SL(1) SL(2) SL(3) SL(4) SL(5) SL(6) SL(7)
s m m m m m m m
0. 0.000 0.000 0.000 0.000 0.000 0.000 0.000
200. -0.254 -0.288 -0.325 -0.371 -0.435 -0.523 -0.523
200. -0.131 -0.161 -0.194 -0.235 -0.293 -0.376 -0.376
300. -0.005 -0.032 -0.060 -0.096 -0.146 -0.224 -0.376

In my case, the number of the first open boundary from 1 to 9
the number of the second open boundary from 15 to 33
the number of the third open boundary from 1199 to 1207
the number of the four open boundary from 1208 to 1240

---

I insert print sentence for the IFRLIQ variable, which is number of boundary in subroutine bord.f , but after running TELEMAC_2D, the message given is:

ITERATION 0 TIME: 0.0000 S
k= 1 IFRLIQ= 1
THE LIQUID BOUNDARIES FILE CONTAINS
11 LINES WITH:
SL(1) SL(2) SL(3) SL(4) SL(5) SL(6) SL(7)
k= 2 IFRLIQ= 1
k= 3 IFRLIQ= 1
k= 4 IFRLIQ= 1
k= 5 IFRLIQ= 1
k= 6 IFRLIQ= 1
k= 7 IFRLIQ= 1
k= 8 IFRLIQ= 1
k= 9 IFRLIQ= 1
k= 15 IFRLIQ= 2
k= 16 IFRLIQ= 2
k= 17 IFRLIQ= 2
k= 18 IFRLIQ= 2
k= 19 IFRLIQ= 2
k= 20 IFRLIQ= 2
k= 21 IFRLIQ= 2
k= 22 IFRLIQ= 2
k= 23 IFRLIQ= 2
k= 24 IFRLIQ= 2
k= 25 IFRLIQ= 2
k= 26 IFRLIQ= 2
k= 27 IFRLIQ= 2
k= 28 IFRLIQ= 2
k= 29 IFRLIQ= 2
k= 30 IFRLIQ= 2
k= 31 IFRLIQ= 2
k= 32 IFRLIQ= 2
k= 33 IFRLIQ= 2
k= 1199 IFRLIQ= 1
k= 1200 IFRLIQ= 1
k= 1201 IFRLIQ= 1
k= 1202 IFRLIQ= 1

The first and second are right and after k=1199 IFRLIQ all equal to 1.

What is wrong?

Thanks

Zheng

Hi
Probably because the last boundary is connect to the first one

You cannot have boundary 3 from 1199 to 1207 and boundary 4 from 1208 to 1240
In the telemac rules boundariez 3 and 4 are only one boundary from 1199 to 1240.
I suppose this si similar for the next boundaries and the last node of the last boundary should be connected to the node 1...

regards

Hi, Coulet,

Thank you very much for your replay quickly.

I have total boundary nodes = 6575, in which has total 117 nodes of open boundaries and rest nodes are solid wall.

In my case, the first open boundary from 1 to 9
the second open boundary from 15 to 33
the third open boundary from 1199 to 1207
the four open boundary from 1208 to 1240
the five open boundary from 1241 to 1255
the six open boundary from 1256 to 1281
the seven open boundary from 1282 to 1287

You are right. The last boundary is connecting to the first one, the third boundary is connecting four, four connecting to five... six and last one. Only two is not connecting to one and third.

How I sort out the problem?

Zheng

Hi, Coulet,
May I reset up two boundaries for the whole 7 boundaries, which original two set up for the first and rest 6 boundaries to the second, which is right or not?

Hi, Coulet,

I am sorry to trouble you again.

I do the test with 2 open boundaries in one liquid boundary file as
#
# Open boundaries file
# 2 boundaries managed
#
T SL(1) SL(2)
s m m
0. 0.000 0.000
100. -0.254 -0.288
200. -0.131 -0.161
300. -0.005 -0.032
400. 0.121 0.099
500. 0.246 0.228
600. 0.367 0.353
700. 0.483 0.474
800. 0.591 0.587
900. 0.690 0.691
1000. 0.779 0.784

After the running TELEMAC-2D, the message given is:

=
VALIDATION PROCEDURE

---

-

1) READING THE REFERENCE FILE :

---

SUITE_SERAFIN : DIFFERENT MESH (OR WRONG STANDARD)



PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2
Duration of job : 264 seconds ( 0:4:24 ) (system=0.016 sec)

Two open boundaries I set for the test, in which node number from 1 to 9 and 1199 to 1287 (1287 connect to 1) for first one and node number from 15 to 33 for second one.

The case has multiple islands and set up solid wall.

What is wrong?

Thank you again,

Zheng

Hi
the last problem is different
Telemac try to made a validation of the result by a comparison between the result and a reference file. This is done for validation test case.
In your case, I suppose you don't have to validate your computation and it seems the reference file which is declare in you steering file is wrong.

Just change you steering file with Validation = NO or erase this keyword and everything goes well.

regards

Yes, that is right. Thanks Coulet!

Hello, Coulet,

I have read some code: bord.f,function fl and read_fic_frliq.f and ready to add some my fortran sentence to read multiple open boundaries files, which are all imposed by elevation varying in space and time as prescribed elevation declard. For example in boundary_1.txt:

100. -0.254 -0.288 -0.325 -0.371 -0.435 -0.523
200. -0.131 -0.161 -0.194 -0.235 -0.293 -0.376
300. -0.005 -0.032 -0.060 -0.096 -0.146 -0.224
400. 0.121 0.099 0.075 0.045 0.002 -0.070

I defined the format file in the steering file:

FORTRAN FILE ='BIEF_OPEN_FILES.f'

FORMATTED DATA FILE 1 ='boundary_1.txt'
FORMATTED DATA FILE 2 ='boundary_2.txt'
…
FORMATTED DATA FILE 7 ='boundary_7.txt'

add INTEGER T2DFO3,.. INTEGER T2dF07 to the MODULE DECLARATIONS_TELEMAC2D and then replaced to the source of TELEMAC 2d

add USE DECLARATIONS_TELEMAC2D to subroutine bief_open_files.f

After running by check out subroutine bief_open_files.f. It seems only upload to T2DF01

OPENING FILES FOR TELEMAC2D
…….
I= 26 LU= 6 T2DFO1 26

The error information is like this:

THE KEY-WORD: FORMATTED DATA FILE 1
APPPEARS AT LEAST TWICE , THE LAST VALUE WILL BE KEPT...

Could you please check it for me?

Cheers and thanks,

zheng