open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all,

I recently set up an Artemis model using Tomawac spectrum as boundary conditions. I did all the necessary configurations and Artemis seemed to abtain the wave conditions from Tomawac and run.

However, I noticed some strange behaviour. Some computation did manage to complete successfully, some did not and all of them exited with an error (exit code -107374094) relating to memory heap corruption. Xeon processors refused completely to run, giving the error straight away.

Debugging addressed the problem to twcale.f with the following error:Note that I used the default discretization of 5 direction and 5 frequencies (random waves mode).
I changed line 216 in twcale.f fromtoand now the computation seems to run smoothly on all cpu types.

Best Regards,
Costas

Thank you.

We will try to have that corrected in the next release for the system.

Sébastien.

Hello costas,

I want to do the same thing that you do, i have a result file of TOMAWAC and i want to make the specral wave from TOMAWAC to a boundary condition in ARTEMIS.

Can you help me to do this? i realy want to know the method that you use.

Hope that you can help me.

Thank you in advance.

Best Regards.
YASMINA.

Hello Yasmina,
You can study the ARTEMIS example chwac1, which demonstrates this functionality. The principle is simple:

1. Run a TOMAWAC model that outputs the spectrum on a single point on an appropriate location along the ARTEMIS boundary.
2. Then pass the file to ARTEMIS, which in turn applies the spectrum along the open boundary.

Just make sure that the spectrum discretisation in ARTEMIS lies within the one defined in the source TOMAWAC model.
Regards,
Costas

Hello Costas,

I am sorry for my late reply.
I had done what you suggest me, but in the first step of outputs the spectrum from TOMAWAC, think that i have make a mistake in my file with coordinates of spectra to write.

Because when i run the simulation it notice me the error below;



MASTER PROCESSOR NUMBER 0 OF THE GROUP OF 4
EXECUTABLE FILE: C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\A.EXE
BARRIER PASSED
LISTING OF TELEMAC2D

---

TTTTT EEEEE L EEEEE M M AAAAA CCCCC
T E L E MM MM A A C
T EEE L EEE M M M AAAAA C
T E L E M M A A C
T EEEEE LLLLL EEEEE M M A A CCCCC

2D VERSION V7P3 FORTRAN 2003

~^~^~^~^~^~^~^~^~^~^~^^~^~^~^~^~^~
~ ~
\ ' o '
/\ o \ o
>=)'> ' /\ '
\/ \ >=)'> ~
/ /\ \/
~ >=)'> / .
\/ )
/ (
~ ) )
} ~ ( ( (
{ ) ) )
} } . ( ( (
{ { /^^^^^^^^^^^^
^^^^^^^^^\ /
^^^^^^^^^

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:

OPENING FILES FOR TELEMAC2D
- FILE C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\T2DGEO00003-00000 OF TYPE BIG_ENDIAN
- FILE C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\T2DGEO00003-00000 IS SINGLE PRECISION

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

READ_MESH_INFO: TITLE= newSelafin
NUMBER OF ELEMENTS: 238530
NUMBER OF POINTS: 120731

FORMAT NOT INDICATED IN TITLE

MXPTEL (BIEF) : MAXIMUM NUMBER OF ELEMENTS AROUND A POINT: 9
MAXIMUM NUMBER OF POINTS AROUND A POINT: 10
(GLOBAL MESH)
SEGBOR (BIEF) : NUMBER OF BOUNDARY SEGMENTS = 2932
INCLUDING THOSE DUE TO DOMAIN DECOMPOSITION
CORRXY (BIEF):NO MODIFICATION OF COORDINATES

MESH: MESH ALLOCATED

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

LISTING OF TOMAWAC

---

TTTTT OOOOO M M AAAAA W W AAAAA CCCCC
T O O MM MM A A W W A A C
T O O M W M AAAAA W W W AAAAA C
T O O M M A A WW WW A A C
T OOOOO M M A A W W A A CCCCC

VERSION V7P3 FORTRAN 2003

| | |
)_) )_) )_) _
)___))___))___)\
)____)____)_____)\\
_____|____|____|____\\\__

---

\ /

---

^^^^^^^^^^^^^^^^^^^^^^^^^^^
^^^^ ^^^^ ^^^ ^^
^^^^ ^^^

COUPLED WITH TELEMAC2D INTERNALLY

********************************************
* SUBROUTINE LECDON_TOMAWAC *
* CALL OF DAMOCLES *
* VERIFICATION OF READ DATA *
* ON STEERING FILE *
********************************************

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

*************************************************
AT LINE 12 THE KEY-WORD : FREQUENCE DE PIC MAXIMALE INITIALE IS UNKNOWN...
*************************************************

************************
* DAMOCLES STOPPED *
************************

---

- ERROR IN THE STEERING FILE -

---

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

DAMOCLE: TRYING ANOTHER LANGUAGE

ATTENTION DANS LE DICTIONNAIRE :
POUR LE MOT-CLE : VARIABLES POUR LES SORTIES GRAPHIQUES 2D
LE NOMBRE DE VALEURS PAR DEFAUT 1 EST DIFFERENT DE LA TAILLE ANNONCEE 2

OUVERTURE DES FICHIERS POUR TOMAWAC
- FICHIER C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\WACGEO00003-00000 DE TYPE BIG_ENDIAN
- FICHIER C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\WACGEO00003-00000 EST SIMPLE PRECISION

*******************************
* CONSTRUCTION DES POINTEURS: *
*******************************

READ_MESH_INFO : TITRE= newSelafin
NOMBRE D'ELEMENTS: 238530
NOMBRE REEL DE POINTS: 120731

FORMAT NON PRECISE DANS LE TITRE

MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D' UN POINT : 9
NOMBRE MAXIMUM DE POINTS VOISINS D'UN POINT : 10
(MAILLAGE GLOBAL)
SEGBOR (BIEF) : NOMBRE DE SEGMENTS DE BORD = 2932
EN COMPTANT CEUX DUS A LA DECOMPOSITION DE DOMAINE
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES

MAILLAGE : MESH ALLOUE

READ_MESH_INFO : TITRE= newSelafin
NOMBRE D'ELEMENTS: 238530
NOMBRE REEL DE POINTS: 120731

FORMAT NON PRECISE DANS LE TITRE

MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D' UN POINT : 9
NOMBRE MAXIMUM DE POINTS VOISINS D'UN POINT : 10
(MAILLAGE GLOBAL)
SEGBOR (BIEF) : NOMBRE DE SEGMENTS DE BORD = 2932
EN COMPTANT CEUX DUS A LA DECOMPOSITION DE DOMAINE
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES

MAILLAGE : MESH3D ALLOUE

****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************

INITIALISATION DE TELEMAC2D POUR
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT

NOMBRE DE FRONTIERES LIQUIDES : 1

CORFON (TELEMAC2D) : PAS DE MODIFICATION DU FOND

INITIALISATION BASEE SUR TPXO :
C_ID( 1 ) = m2
C_ID( 2 ) = s2
C_ID( 3 ) = n2
C_ID( 4 ) = k2
C_ID( 5 ) = k1
C_ID( 6 ) = o1
C_ID( 7 ) = p1
C_ID( 8 ) = q1
C_ID( 9 ) = mf
C_ID( 10 ) = mm
C_ID( 11 ) = m4
C_ID( 12 ) = ms4
C_ID( 13 ) = mn4
- OBTENTION DES NIVEAUX
- INTERPOLATION DES NIVEAUX
- OBTENTION DES VITESSES
- INTERPOLATION DES VITESSES
FIN DE L'INITIALISATION TPXO

================================================================================
ITERATION 0 TEMPS : 0.0000 S
TELEMAC2D COUPLE AVEC : TOMAWAC
TELEMAC-2D : COUPLAGE INTERNE AVEC TOMAWAC
At line 68 of file C:\opentelemac-mascaret\v7p3r1\sources\tomawac\read_spectra_coords.f (unit = 15, file = 'C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\LEOWXY00003-00003')
Fortran runtime error: Bad integer for item 1 in list input
At line 68 of file C:\opentelemac-mascaret\v7p3r1\sources\tomawac\read_spectra_coords.f (unit = 15, file = 'C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\LEOWXY00003-00001')
Fortran runtime error: Bad integer for item 1 in list input
At line 68 of file C:\opentelemac-mascaret\v7p3r1\sources\tomawac\read_spectra_coords.f (unit = 15, file = 'C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\LEOWXY00003-00002')
Fortran runtime error: Bad integer for item 1 in list input
At line 68 of file C:\opentelemac-mascaret\v7p3r1\sources\tomawac\read_spectra_coords.f (unit = 15, file = 'C:\Stage_yasmina\Modélisation_hydrodynamique\large_box\Toma_tele\init2_Golf.cas_2018-06-13-11h00min53s\LEOWXY00003-00000')
Fortran runtime error: Bad integer for item 1 in list input.


I atached my file with coordinates of spectra to write.



Hope that you can clear me more.

Thank you.
Best Regards.
YASMINA.

Hello yasmina,
First of all, you only need to define one point, otherwise Artemis won't accept the spectrum file. Secondly, in your spectrum coordinates file, you should add another (unused) integer right after the number of points (1948 in your case). 0 will do fine.
Regards,
Costas