Hello,
I have done what you adviced me, thank you for this.
So i had launched the simulation for 1000s just to verify that my subroutine works. The simulation works!!!! and the result seems correct juste i must launch it for more time.
I attached the output of the simulation:
MASTER PROCESSOR NUMBER 0 OF THE GROUP OF 4
EXECUTABLE FILE: C:\Stage_yasmina\Modélisation_hydrodynamique\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\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\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\T2DGEO00003-00000 OF TYPE BIG_ENDIAN
- FILE C:\Stage_yasmina\Modélisation_hydrodynamique\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\T2DGEO00003-00000 IS SINGLE PRECISION
*****************************
* MEMORY ORGANIZATION *
*****************************
READ_MESH_INFO: TITLE= newSelafin
NUMBER OF ELEMENTS: 102997
NUMBER OF POINTS: 53446
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 = 3911
INCLUDING THOSE DUE TO DOMAIN DECOMPOSITION
CORRXY (BIEF):NO MODIFICATION OF COORDINATES
MESH: MESH ALLOCATED
*************************************
* END OF MEMORY ORGANIZATION: *
*************************************
INITIALISING TELEMAC2D FOR
INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)
STRCHE (BIEF): NO MODIFICATION OF FRICTION
NUMBER OF LIQUID BOUNDARIES: 1
CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM
================================================================================
ITERATION 0 TIME: 0.0000 S
TELEMAC2D INITIALISED
USING STREAMLINE VERSION 7.3 FOR CHARACTERISTICS
PROSOU : IN CARTESIAN COORDINATES, THE CORIOLIS
PARAMETER IS READ IN THE STEERING FILE
IT IS THE KEY WORD 'CORIOLIS
COEFFICIENT', IT IS UNIFORM IN SPACE
GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 100 RELATIVE PRECISION: 0.1404821E-03
GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 100 RELATIVE PRECISION: 0.1539796E-03
GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 100 RELATIVE PRECISION: 0.1154087E-03
================================================================================
ITERATION 50 TIME: 8 MN 20.0000 S ( 500.0000 S)
--------------------------------------------------------------------------------
ADVECTION STEP
--------------------------------------------------------------------------------
DIFFUSION-PROPAGATION STEP
GRACJG (BIEF) : 34 ITERATIONS, RELATIVE PRECISION: 0.9383331E-04
--------------------------------------------------------------------------------
BALANCE OF WATER VOLUME
VOLUME IN THE DOMAIN : 0.9516157E+08 M3
FLUX BOUNDARY 1: -64721.91 M3/S ( >0 : ENTERING <0 : EXITING )
RELATIVE ERROR IN VOLUME AT T = 500.0 S : 0.2752341E-06
================================================================================
ITERATION 100 TIME: 16 MN 40.0000 S ( 1000.0000 S)
--------------------------------------------------------------------------------
ADVECTION STEP
--------------------------------------------------------------------------------
DIFFUSION-PROPAGATION STEP
GRACJG (BIEF) : 35 ITERATIONS, RELATIVE PRECISION: 0.9594290E-04
--------------------------------------------------------------------------------
BALANCE OF WATER VOLUME
VOLUME IN THE DOMAIN : 0.7236964E+08 M3
FLUX BOUNDARY 1: -31629.34 M3/S ( >0 : ENTERING <0 : EXITING )
RELATIVE ERROR IN VOLUME AT T = 1000. S : 0.1036767E-05
--------------------------------------------------------------------------------
FINAL BALANCE OF WATER VOLUME
RELATIVE ERROR CUMULATED ON VOLUME: -0.2310381E-04
INITIAL VOLUME : 0.1713690E+09 M3
FINAL VOLUME : 0.7236964E+08 M3
VOLUME THAT ENTERED THE DOMAIN: -0.9900332E+08 M3 ( IF <0 EXIT )
TOTAL VOLUME LOST : -3959.276 M3
END OF TIME LOOP
EXITING MPI
CORRECT END OF RUN
STOP 0
STOP 0
ELAPSE TIME :
STOP 0
30 SECONDS
STOP 0
The problem after the modification of the time step from 100 to 2000, the simulation don't work. In the begenning it's indicate me that the liquid boundary file has a tabulation replace it with space, wich i was done. But now it's indicate that there is some problems in the subroutine SL.
ITERATION 0 TIME: 0.0000 S
TELEMAC2D INITIALISED
THE LIQUID BOUNDARIES FILE CONTAINSapplication called MPI_Abort(MPI_COMM_WORLD,
2) - process 0
6 LINES WITH:
ETA(34331ETA(34288ETA(34240ETA(34197ETA(34147ETA(34094ETA(34043ETA(33989ETA(339
33ETA(33876ETA(33821ETA(33764ETA(33707ETA(33651ETA(33592ETA(33533ETA(33473ETA(33
414ETA(33353ETA(33291ETA(33230ETA(33168ETA(33104ETA(33041ETA(32974ETA(32907ETA(3
2839ETA(32773ETA(32707ETA(32639ETA(32573ETA(32506ETA(32439ETA(32372ETA(32306ETA(
32240ETA(32173ETA(32107ETA(32044ETA(31979ETA(31914ETA(31051ETA(31785ETA(31720ETA
(31054ETA(31588ETA(31521ETA(31459ETA(31398ETA(31338ETA(31278ETA(31220ETA(31162ET
A(31104ETA(31048ETA(30992ETA(30936ETA(30883ETA(30830ETA(30779ETA(30727ETA(30677E
TA(30629ETA(30582ETA(30535ETA(30490ETA(30447ETA(30403ETA(30361ETA(30321ETA(30281
ETA(30243ETA(30205ETA(30170ETA(30137ETA(30103ETA(30074ETA(30043ETA(30014ETA(2998
6ETA(29960ETA(29938ETA(29913ETA(20801ETA(29872ETA(29852ETA(29834ETA(29819ETA(298
44ETA(29826ETA(29809ETA(29802ETA(29806
SL: MORE PRESCRIBED ELEVATIONS ARE REQUIRED
IN THE PARAMETER FILE
AT LEAST 1 MUST BE GIVEN
PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2
job aborted:
rank: node: exit code[: error message]
Now i don't understand what it can be the solution.
I attached my liquid boundary file:
File Attachment:
File Name:
liq.txtFile Size: 12 KB
And this is my fortran file:
File Attachment:
File Name:
bord.fFile Size: 30 KB
Hope that you can help me in this error.
Best regards.
YASMINA.