open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone!

6 months ago I used ARTEMIS V5p6 to simulate the waves agitation with a gravity fondation in offshore. The results was good and i had been any error message.

Today I use ARTEMIS V6p0, the model is different with the type of fondation, "jacket". But I still using de same keys words and fortran file with BORH and CORFON subroutine. But when I run this model ARTEMIS send me a error message:
Here is the screen of artemis:

- FICHIER DES CONDITIONS AUX LIMITES : artemis.conlim

______________________________________________________________________________
*** EXECUTION ***






LISTING DE ARTEMIS

---

---

AAA RRRR TTTTT EEEEE M M IIIII SSSS
A A R R T E MM MM I S
AAAAA RRRR T EEEEE M M M I SSS
A A R R T E M M I S
A A R R T EEEEE M M IIIII SSSS

VERSION 6.0 FORTRAN 90








********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************


TITRE DE L'ETUDE : -
Cas: T=7.

OUVERTURE DES FICHIERS POUR ARTEMIS



*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************



READGEO1 : TITRE= V:\Telemac\Etudes\ARTEMIS_WPD_FECAMP\Maillage.MAT


NOMBRE D'ELEMENTS: 30230
NOMBRE REEL DE POINTS: 15317
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'UN POINT : 9
NOMBRE MAXIMUM DE POINTS VOISINS D'UNPOINT : 9
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES

MAILLAGE : MESH ALLOUE



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

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


IL Y A 20 ILE(S) DANS LE DOMAINE

IL Y A 2 FRONTIERE(S) LIQUIDE(S) :

FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 3
ET DE COORDONNEES : -2180.000 1006.750
FIN AU POINT DE BORD 87 , DE NUMERO GLOBAL 432
ET DE COORDONNEES : 2111.029 1000.000

FRONTIERE 2 :
DEBUT AU POINT DE BORD 92 , DE NUMERO GLOBAL 70
ET DE COORDONNEES : 2180.000 1180.902
FIN AU POINT DE BORD 199 , DE NUMERO GLOBAL 378
ET DE COORDONNEES : -2170.959 1204.872

IL Y A 22 FRONTIERE(S) SOLIDE(S) :

FRONTIERE 1 :
DEBUT AU POINT DE BORD 87 , DE NUMERO GLOBAL 432
ET DE COORDONNEES : 2111.029 1000.000
FIN AU POINT DE BORD 92 , DE NUMERO GLOBAL 70
ET DE COORDONNEES : 2180.000 1180.902

FRONTIERE 2 :
DEBUT AU POINT DE BORD 199 , DE NUMERO GLOBAL 378
ET DE COORDONNEES : -2170.959 1204.872
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 3
ET DE COORDONNEES : -2180.000 1006.750

FRONTIERE 3 :
DEBUT AU POINT DE BORD 311 , DE NUMERO GLOBAL 404
ET DE COORDONNEES : -9.447400 2140.185
FIN AU POINT DE BORD 311 , DE NUMERO GLOBAL 404
ET DE COORDONNEES : -9.447400 2140.185

FRONTIERE 4 :
DEBUT AU POINT DE BORD 359 , DE NUMERO GLOBAL 223
ET DE COORDONNEES : -3.665700 2140.367
FIN AU POINT DE BORD 359 , DE NUMERO GLOBAL 223
ET DE COORDONNEES : -3.665700 2140.367

FRONTIERE 5 :
DEBUT AU POINT DE BORD 358 , DE NUMERO GLOBAL 55
ET DE COORDONNEES : -8.948400 2146.151
FIN AU POINT DE BORD 358 , DE NUMERO GLOBAL 55
ET DE COORDONNEES : -8.948400 2146.151

FRONTIERE 6 :
DEBUT AU POINT DE BORD 323 , DE NUMERO GLOBAL 93
ET DE COORDONNEES : -1.650200 2143.816
FIN AU POINT DE BORD 323 , DE NUMERO GLOBAL 93
ET DE COORDONNEES : -1.650200 2143.816

FRONTIERE 7 :
DEBUT AU POINT DE BORD 383 , DE NUMERO GLOBAL 285
ET DE COORDONNEES : -6.000600 2148.166
FIN AU POINT DE BORD 383 , DE NUMERO GLOBAL 285
ET DE COORDONNEES : -6.000600 2148.166

FRONTIERE 8 :
DEBUT AU POINT DE BORD 418 , DE NUMERO GLOBAL 415
ET DE COORDONNEES : 3.166400 2140.367
FIN AU POINT DE BORD 418 , DE NUMERO GLOBAL 415
ET DE COORDONNEES : 3.166400 2140.367

FRONTIERE 9 :
DEBUT AU POINT DE BORD 287 , DE NUMERO GLOBAL 324
ET DE COORDONNEES : -8.949300 2152.984
FIN AU POINT DE BORD 287 , DE NUMERO GLOBAL 324
ET DE COORDONNEES : -8.949300 2152.984

FRONTIERE 10 :
DEBUT AU POINT DE BORD 263 , DE NUMERO GLOBAL 262
ET DE COORDONNEES : 1.150400 2143.816
FIN AU POINT DE BORD 263 , DE NUMERO GLOBAL 262
ET DE COORDONNEES : 1.150400 2143.816

FRONTIERE 11 :
DEBUT AU POINT DE BORD 251 , DE NUMERO GLOBAL 25
ET DE COORDONNEES : -5.999900 2150.965
FIN AU POINT DE BORD 251 , DE NUMERO GLOBAL 25
ET DE COORDONNEES : -5.999900 2150.965

FRONTIERE 12 :
DEBUT AU POINT DE BORD 299 , DE NUMERO GLOBAL 380
ET DE COORDONNEES : 8.451400 2140.185
FIN AU POINT DE BORD 299 , DE NUMERO GLOBAL 380
ET DE COORDONNEES : 8.451400 2140.185

FRONTIERE 13 :
DEBUT AU POINT DE BORD 395 , DE NUMERO GLOBAL 286
ET DE COORDONNEES : -9.447400 2158.084
FIN AU POINT DE BORD 395 , DE NUMERO GLOBAL 286
ET DE COORDONNEES : -9.447400 2158.084

FRONTIERE 14 :
DEBUT AU POINT DE BORD 239 , DE NUMERO GLOBAL 76
ET DE COORDONNEES : -1.650800 2155.316
FIN AU POINT DE BORD 239 , DE NUMERO GLOBAL 76
ET DE COORDONNEES : -1.650800 2155.316

FRONTIERE 15 :
DEBUT AU POINT DE BORD 419 , DE NUMERO GLOBAL 147
ET DE COORDONNEES : 5.501100 2148.167
FIN AU POINT DE BORD 419 , DE NUMERO GLOBAL 147
ET DE COORDONNEES : 5.501100 2148.167

FRONTIERE 16 :
DEBUT AU POINT DE BORD 214 , DE NUMERO GLOBAL 154
ET DE COORDONNEES : 8.448800 2146.150
FIN AU POINT DE BORD 214 , DE NUMERO GLOBAL 154
ET DE COORDONNEES : 8.448800 2146.150

FRONTIERE 17 :
DEBUT AU POINT DE BORD 371 , DE NUMERO GLOBAL 414
ET DE COORDONNEES : -3.666800 2158.265
FIN AU POINT DE BORD 371 , DE NUMERO GLOBAL 414
ET DE COORDONNEES : -3.666800 2158.265

FRONTIERE 18 :
DEBUT AU POINT DE BORD 227 , DE NUMERO GLOBAL 235
ET DE COORDONNEES : 1.149400 2155.317
FIN AU POINT DE BORD 227 , DE NUMERO GLOBAL 235
ET DE COORDONNEES : 1.149400 2155.317

FRONTIERE 19 :
DEBUT AU POINT DE BORD 335 , DE NUMERO GLOBAL 294
ET DE COORDONNEES : 5.500600 2150.967
FIN AU POINT DE BORD 335 , DE NUMERO GLOBAL 294
ET DE COORDONNEES : 5.500600 2150.967

FRONTIERE 20 :
DEBUT AU POINT DE BORD 215 , DE NUMERO GLOBAL 225
ET DE COORDONNEES : 8.448700 2152.981
FIN AU POINT DE BORD 215 , DE NUMERO GLOBAL 225
ET DE COORDONNEES : 8.448700 2152.981

FRONTIERE 21 :
DEBUT AU POINT DE BORD 275 , DE NUMERO GLOBAL 164
ET DE COORDONNEES : 3.166200 2158.267
FIN AU POINT DE BORD 275 , DE NUMERO GLOBAL 164
ET DE COORDONNEES : 3.166200 2158.267

FRONTIERE 22 :
DEBUT AU POINT DE BORD 431 , DE NUMERO GLOBAL 427
ET DE COORDONNEES : 8.451800 2158.084
FIN AU POINT DE BORD 431 , DE NUMERO GLOBAL 427
ET DE COORDONNEES : 8.451800 2158.084
********************************************CORFON

================================================================================

PERIODE 1/ 2 : 6.5127 SECONDES

Tout s est bien passe

DIRECTION 1/ 2 : 263.1112 DEGRES

PHBOR : IL Y A : 1 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR
CNTPRE (ARTEMIS) : PRECONDITIONNEMENT DIAGONAL NON APPLIQUE (UN ELEMENT DIAGONA
L DE LA MATRICE EST NEGATIF OU NUL)

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 2 : 276.8813 DEGRES

PHBOR : IL Y A : 1 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

================================================================================

PERIODE 2/ 2 : 7.5078 SECONDES

Tout s est bien passe

DIRECTION 1/ 2 : 263.1112 DEGRES

PHBOR : IL Y A : 1 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 2 : 276.8813 DEGRES

PHBOR : IL Y A : 1 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE
OV (BIEF) : Y**C INTERDIT SI Y < 0



PLANTE : ARRET DU PROGRAMME APRES ERREUR
Returning exit code: 2
Duree du calcul : 3 secondes ( 0:0:3 ) (systeme=0.031 sec)
______________________________________________________________________________
*** RESTITUTION DES FICHIERS ***

- FICHIER DES RESULTATS : test.res


________________________________________________________
Execution terminee: artemis.bat
________________________________________________________
No compilation/linking/file errors detected.
No execution errors detected.
Returning exit status 0

=========================================================
Telemac System Freeware 6.0 - Perl scripts version V6.0
=========================================================
...stopping.

V:\Telemac\Etudes\artemis>


I check and recheck my file parameter, file fortran. But I don't understand why ARTEMIS send me "RESOLUTION DIRECTE OV (BIEF) : Y**C INTERDIT SI Y < 0".

I hope you can help me. Thanks in advance.

Romain_D

Hi,

After severals tests, my last error message has disappeared. But now I have a new one...

When I start a run for the first computation (period/direction) artemis send me:


PERIODE 1/ 2 : 6.5127 SECONDES

Tout s est bien passe

DIRECTION 1/ 2 : 263.1112 DEGRES

PHBOR : IL Y A : 1 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

SOUS-ITERATION NUMERO : 1


RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE
MATRICE AVEC PIVOT NUL A LA LIGNE
14



PLANTE : ARRET DU PROGRAMME APRES ERREUR
Returning exit code: 2
Duree du calcul : 2 secondes ( 0:0:2 ) (systeme=0.016 sec)

My model has 3 boundaries ,North 1 1 1 4, South 4 6 6 5, East/Weast/jacket 2 2 2 2. (img).

My mesh was built like this:
0.5m around the jacket and 20m at the boundaries.

I don't understand why the computation stop!

I can send you my files(cas conlim geo fortran), if you need it.

Thanks in advance.

Romain_D

Hello !

same problem for me If I find why I keep you in touch

bye

hello

here is my problem wich is the same than ROMAIN. I attached my files, perhaps you can see where I made a mistake




TITRE DE L'ETUDE : PAS DE TITRE DANS LE FICHIER CAS


OUVERTURE DES FICHIERS POUR ARTEMIS



*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************



READGEO1 : TITRE= C:\TELEMAC\jedi\agitation\test_2\test2.MAT


NOMBRE D'ELEMENTS: 42181
NOMBRE REEL DE POINTS: 21394
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UNPOINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES

MAILLAGE : MESH ALLOUE



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

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


IL Y A 1 FRONTIERE(S) SOLIDE(S) :

FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 332
ET DE COORDONNEES : 2827.700 2042.600
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 332
ET DE COORDONNEES : 2827.700 2042.600

================================================================================

PERIODE 1/ 4 : 6.7546 SECONDES


DIRECTION 1/ 5 : 76.9860 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 5 : 90.5580 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 3/ 5 : 99.9900 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 4/ 5 : 109.4220 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 5/ 5 : 122.9940 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

================================================================================

PERIODE 2/ 4 : 7.5964 SECONDES


DIRECTION 1/ 5 : 76.9860 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 5 : 90.5580 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 3/ 5 : 99.9900 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 4/ 5 : 109.4220 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 5/ 5 : 122.9940 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

================================================================================

PERIODE 3/ 4 : 8.1139 SECONDES


DIRECTION 1/ 5 : 76.9860 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 5 : 90.5580 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 3/ 5 : 99.9900 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 4/ 5 : 109.4220 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 5/ 5 : 122.9940 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

================================================================================

PERIODE 4/ 4 : 8.9858 SECONDES


DIRECTION 1/ 5 : 76.9860 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 2/ 5 : 90.5580 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 3/ 5 : 99.9900 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 4/ 5 : 109.4220 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE

DIRECTION 5/ 5 : 122.9940 DEGRES

PHBOR : IL Y A : 0 FRONTIERE(S) DE TYPE ONDE INCIDENTE
END PHBOR

RESOLUTION DU SYSTEME LINEAIRE (SOLVE)

RESOLUTION DIRECTE
OV (BIEF) : Y**C INTERDIT SI Y < 0



PLANTE : ARRET DU PROGRAMME APRES ERREUR
Returning exit code: 2
Duree du calcul : 22 secondes ( 0:0:22 ) (systeme=0.031 sec)
______________________________________________________________________________
*** RESTITUTION DES FICHIERS ***

- FICHIER DES RESULTATS : test2_res


________________________________________________________
Execution terminee: artemis.bat
________________________________________________________
No compilation/linking/file errors detected.
No execution errors detected.
Returning exit status 0

=========================================================
Telemac System Freeware 6.0 - Perl scripts version V6.0
=========================================================
...stopping.

C:\TELEMAC\jedi\agitation\test_2>

Hello,

Yesterday I called the assistance and the problem comes from Artemis V6P0. The problem doesn't fixed yet.
I will try with a old version of TELEMAC like V5P9.

Regards.

Romain_D

Hello,

To overcome the problem, I propose that you use the enclosed version of ov.f in your Fortran file (look for operation X=Y**C)

Basically the bug is just that someone used the formula X**2.D0 instead of X**2, X is positive but can be slightly negative due to truncation errors, the second formula will work but the first one crashes. In ov.f I have just put MAX(Y,1.D-15) instead of Y.

With best regards,

Jean-Michel Hervouet

hello

Thank you very much for your help ! It works in several case but I have a new problem with a new project. It is working well when I launch simulations without wave breaking. When I try to launch it with breakin it crashes with "pivot nul à la ligne xxx". Is it a way to fix it ?

thanks a lot

hello,

Please can you send your artemis files?

Regards,

Christophe

hello
I finally solve the problem. I made a mistake on the direction of the swell.
Hope this help !

hello

Romain, can you tell me what was your solution to remove the first error ? Thanks a lot !