open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi

I am trying to modelize the tide flow at Honikulu inlet in Wallis islands.

I process the bathymetric 3D model with BlueKenue and this seemed to be successful (?) this produced :

the geometry file
the boundary conditions file

I imported a steering file tide.cas from the Telemac2D example Tide

it was adaptated for tpxo database with the 2 binary tide files :

hf.PO_2009.out
uv.PO_2009.out

found in the tpxo lib in lib directories

the geographic system is : 3 = UTM south - zone number = 1

there are 3 open boundaries one entrance and two outputs

others options are standard.

whith most of the example I tryied the programm was running fine but not with my data

while launching the simulation it seems to start well but stops in BIEF module as you can see in the exec_file joined

ERREUR 5010 LORS DE L APPEL A GET_BND_CONNECTIVITY_SRF:READ
TEXTE DE L'ERROR : UNKONWN ERROR

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2


May be the geometry or some options are wrong ?

I'll be grateful if you can examine the outputs joinded and see what's wrong

Thank you for a reply

Best regards and thanks

T.H.

Hi
The boundary condition file should be a .cli file
You use the .bc2 file which is a Blue Kenue boundary file format which is not compatible with Telemac

Regards

Dear C. Coulet

thank you for your reply and help. I made a simplier model with better succes and improvments in the execution. Still some troubleshoutings

WARNING IN DICTIONARY:
FOR KEYWORD: VARIABLES FOR GRAPHIC PRINTOUTS
THE NUMBER OF DEFAULT VALUES 1 IS DIFFERENT FROM THE DECLARED SIZE 4

END OF FILE FOR DAMOCLES

********************************************
* LECDON: *
* AFTER CALLING DAMOCLES *
* CHECKING OF DATA READ *
* IN THE STEERING FILE *
********************************************

EXITING LECDON. NAME OF THE STUDY:
HONIKULU

OPENING FILES FOR TELEMAC2D
- FILE T2DGEO OF TYPE BIG_ENDIAN
- FILE T2DGEO IS SINGLE PRECISION

*****************************
* MEMORY ORGANIZATION *
*****************************

READ_MESH_INFO: TITLE= honik_geometry
NUMBER OF ELEMENTS: 9855
NUMBER OF POINTS: 5060

FORMAT NOT INDICATED IN TITLE

MXPTEL (BIEF) : MAXIMUM NUMBER OF ELEMENTS AROUND A POINT: 8
MAXIMUM NUMBER OF POINTS AROUND A POINT: 8
CORRXY (BIEF):NO MODIFICATION OF COORDINATES

MESH: MESH ALLOCATED

*************************************
* END OF MEMORY ORGANIZATION: *
*************************************

INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)
STRCHE (BIEF): NO MODIFICATION OF FRICTION

THERE IS 3 LIQUID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 13 , WITH GLOBAL NUMBER: 1638
AND COORDINATES: 584426.3 8519716.
ENDS AT BOUNDARY POINT: 38 , WITH GLOBAL NUMBER: 3673
AND COORDINATES: 585011.9 8519545.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 85 , WITH GLOBAL NUMBER: 5060
AND COORDINATES: 585660.0 8520070.
ENDS AT BOUNDARY POINT: 112 , WITH GLOBAL NUMBER: 4849
AND COORDINATES: 585518.8 8520662.

BOUNDARY 3 :
BEGINS AT BOUNDARY POINT: 192 , WITH GLOBAL NUMBER: 1854
AND COORDINATES: 584487.8 8521166.
ENDS AT BOUNDARY POINT: 221 , WITH GLOBAL NUMBER: 27
AND COORDINATES: 583842.4 8520961.

THERE IS 3 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 221 , WITH GLOBAL NUMBER: 27
AND COORDINATES: 583842.4 8520961.
ENDS AT BOUNDARY POINT: 13 , WITH GLOBAL NUMBER: 1638
AND COORDINATES: 584426.3 8519716.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 38 , WITH GLOBAL NUMBER: 3673
AND COORDINATES: 585011.9 8519545.
ENDS AT BOUNDARY POINT: 85 , WITH GLOBAL NUMBER: 5060
AND COORDINATES: 585660.0 8520070.

BOUNDARY 3 :
BEGINS AT BOUNDARY POINT: 112 , WITH GLOBAL NUMBER: 4849
AND COORDINATES: 585518.8 8520662.
ENDS AT BOUNDARY POINT: 192 , WITH GLOBAL NUMBER: 1854
AND COORDINATES: 584487.8 8521166.

OPTION FOR TIDAL BOUNDARY CONDITIONS:
GIVE THE SAME NUMBER OF VALUES AS THE NUMBER
OF LIQUID BOUNDARIES, I.E. 3

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2

maybe int he cas file ?

thank you for helping

Best regards

Thierry

Hello

The error message is explicit. For the keyword OPTION FOR TIDAL BOUNDARY CONDITIONS you need to give 3 values (the same numer of liquid boundaries) so write for instance
OPTION FOR TIDAL BOUNDARY CONDITIONS = 1;0;1
(1 for tidal and 0 no tide)
please read the user manual


with my best regards

Riadh

Dear Riadh Ata

thank you very much for your advices, the correction made great progress in the execution of the programm which can run until the end end. I made some corrections on the cas file but infortunately the results are still not valid. It seems that the variation of the sea level is not taken in account in the simulation even if the tide module TPXO seems to be integrated correctly. So there is no surface variation nor currents. there is probably some no coherent instruction in the cas or a mistake in the definitions of the boundary ? I send you the listing of the results so you could interpret the error code and find the error

thank you for your help

best regards

Thierry

Hello Thierry

Again, read the listing !
The solver is telling that
GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 20 RELATIVE PRECISION: NaN

when I looked at your steering file, I found that you are using a time step of 1200 s, which is the largest time-step I ever seen!
Conclusion: the solver is not converging because the time step is ... huge.
So please reduce the time-step, we often start with 1s and then we adapt it depending on the solver efficiency.

I hope that this helps.

with my kind regards

Riadh

Dear Riadh Ata

thank your very much for your reply. You were partially true concerning timestep. I changed to 60s but the problem still hanged the programm. I verify the params and changed option

// 1: ORIGINAL EQUATION - 2: WAVE EQ
TREATMENT OF THE LINEAR SYSTEM :

from 1 to 2 (WAVE EQUATION) and then all works good without planting, nice indeed

I'am a little deceipt by the results ; I thought the currents should be more importants in the entrance inlet (Honikulu) but it seems that currents are more importants at the other end of the inlet. (you can see the place on the picture joined). May be some more improvement to do in the choice of the params.

I go further with wider grids and let you informed.

thank you again for your help

Best regards and thanks

Thierry

Hello Thierry,

In addition, if you're simulating tidal flow in the inlet, maybe you should turn on TIDAL FLATS. Even if the tidal range is small in this area, some points can become dry during the simulation, especially if you're going to refine your mesh.

It could also be important to extend your computational domain, both in offshore direction and inside the lagoon to well reproduce tidal wave propagation in your domain.
Wave effects can also be important/dominant in this kind of place, in so far as they generally increase the net flow by breaking over the reef. An increase in the mean sea level is then observed inside the lagoon associated to a strengthening of the outgoing current by the pass ("ensachage" in french).

Best regards,
Laurent