open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I'm still tying to use Artemis, but I've run into a new error... I can't understand where the problem is! In attachment my last command window.
Thank you for any help!

Luca


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

HOSTTYPE : win
PROJECT : C:\TELEMAC\V6P0
BASE DIRECTORY : C:\TELEMAC\VOLTRI\VOL1.MAT
LAUNCH DIRECTORY : C:\TELEMAC\VOLTRI\VOL1.MAT
WORK DIRECTORY : C:\TELEMAC\VOLTRI\VOL1.MAT\cas3236_tmp
PARAMETER FILE : cas


*** Fichier de configuration par defaut :
C:\TELEMAC\V6P0\config\systel.ini ***



*** ARTEMIS sur STATION ***


*** Lancement en interactif ***


*** VERSION v6p0 ***

________________________________________________________
Fichier des parametres : cas
________________________________________________________

________________________________________________________
Starting execution: artemis.bat
________________________________________________________
- FICHIER FORTRAN : borh.f
______________________________________________________________________________
*** EXECUTABLE LOCAL ***

borh_win_v6p0.exe

______________________________________________________________________________
*** ACQUISITION DES FICHIERS ***

- FICHIER DES PARAMETRES : cas
- DICTIONNAIRE : artemisv6p0.dico

- FICHIER DE GEOMETRIE : geo
- FICHIER DES CONDITIONS AUX LIMITES : cas.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






*************************************************
A LA LIGNE 12 LE MOT CLE SUIVANT : INITIAL CONDITIONS
EST INCONNU ...
*************************************************

************************
* ARRET DE DAMOCLES *
************************

---

- ERREUR DANS LE FICHIER DES PARAMETRES -

---

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


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


NAME OF THE STUDY : VOLTRINUM


OPENING FILES FOR ARTEMIS



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



READGEO1: TITLE= C:\TELEMAC\VOLTRI\VOL1.MAT


NUMBER OF ELEMENTS: 86390
NUMBER OF POINTS: 43732
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 1 LIQUID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 9
AND COORDINATES: 1.000000 2.700000
ENDS AT BOUNDARY POINT: 154 , WITH GLOBAL NUMBER: 4
AND COORDINATES: 1529.300 1.000000

THERE IS 1 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 154 , WITH GLOBAL NUMBER: 4
AND COORDINATES: 1529.300 1.000000
ENDS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 9
AND COORDINATES: 1.000000 2.700000

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

PERIOD 1/ 5 : 5.9707 SECONDS

PHBOR : THERE ARE : 1 BOUNDARIE(S) OF INCIDENT WAVE TYPE
END PHBOR

LINEAR SYSTEM SOLVING (SOLVE)

DIRECT SYSTEM SOLVER

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

PERIOD 2/ 5 : 8.1640 SECONDS

PHBOR : THERE ARE : 1 BOUNDARIE(S) OF INCIDENT WAVE TYPE
END PHBOR

LINEAR SYSTEM SOLVING (SOLVE)

DIRECT SYSTEM SOLVER

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

PERIOD 3/ 5 : 9.3032 SECONDS

PHBOR : THERE ARE : 1 BOUNDARIE(S) OF INCIDENT WAVE TYPE
END PHBOR

LINEAR SYSTEM SOLVING (SOLVE)

DIRECT SYSTEM SOLVER

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

PERIOD 4/ 5 : 10.2574 SECONDS

PHBOR : THERE ARE : 1 BOUNDARIE(S) OF INCIDENT WAVE TYPE
END PHBOR

LINEAR SYSTEM SOLVING (SOLVE)

DIRECT SYSTEM SOLVER

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

PERIOD 5/ 5 : 11.4296 SECONDS

PHBOR : THERE ARE : 1 BOUNDARIE(S) OF INCIDENT WAVE TYPE
END PHBOR

LINEAR SYSTEM SOLVING (SOLVE)

DIRECT SYSTEM SOLVER
OV (BIEF): Y**C FORBIDDEN IF Y < 0



PLANTE: PROGRAM STOPPED AFTER AN ERROR
Returning exit code: 2
Duree du calcul : 25 secondes ( 0:0:25 ) (systeme=0.015 sec)
______________________________________________________________________________
*** RESTITUTION DES FICHIERS ***

- FICHIER DES RESULTATS : resvol


________________________________________________________
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\VOLTRI\VOL1.MAT>

Dear Luca,

Well received. I have warned my colleagues in charge of Artémis and they will answer as soon as possible. Maybe you can anticipate that they will need your files to reproduce the problem. One thing you can try is to test an iterative solver like 3 (normal equation). It will require a lot of iterations but will possibly work.

With best regards,

Jean-Michel Hervouet

Hi,

I've solved the problem.
In the borh subroutine, I've first inizialized variables in this way:

TETABT( : ) = TETAH
TETAPT( : ) = 0.D0
ALFAPT( : ) = 0.D0
RPT( : ) = 0.D0
HBT( : ) = 0.D0

then, I've re-defined them in the following wrong way:

DO 100 I = 1,154
LIHBOR%I(I) = KINC
HB%I(I) = 5.50D0
TETAB%I(I) = 90.D0
TETAP%I(I) = 0.D0
100 CONTINUE

... that was a "copy and paste" error!

Note that, writing a value different to 0.D0 for HBT, the program works, even if the wave height is only defined in the inizialization.

LC

Yes Luca, I also prefer clear loops than e.g. TETABT%R(:)=0.D0 where the real dimension of the array is not obvious.

Users programming with BIEF structures like here must recall that the double precision values are stored in component R and integers in component I, hence in your case :

LIHBOR%I(I)=KINC

and HB%R(I)=5.5D0 HB%I(I) will go into a part of the memory which is not allocated (some compilers will not trigger an error and this is a problem).

Regards,

Jean-Michel