Welcome, Guest
Username: Password: Remember me

TOPIC: Simulation stops after 1 iteration

Simulation stops after 1 iteration 11 years 6 months ago #8759

  • dsc1r12
  • dsc1r12's Avatar
Hi,

I am running a simple model in Telemac 2d with no subroutines where I have a channel flow with closed boundaries on each side and open boundaries at either end with prescribed elevations of 10 and 9.75 respectively, these are defined in the steering file. For some reason when I run the simulation it gets to iteration 0, then stops, and I am unable to type anything in the cmd window, as shown below:

C:\opentelemac\Practise\Practise5\Telemac files>runcode.py telemac2d -s Steering
.txt


Loading Options and Configurations
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

... parsing configuration file: c:\opentelemac\v6p2r2\config\systel.cfg


Running your CAS file for:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

+> configuration: wintel
+> root: C:\opentelemac\v6p2r2\
+> version v6p2


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



Running Steering.txt with telemac2d under C:\opentelemac\Practise\Practise5\Tele
mac files
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

... reading module dictionary
... running in English with Steering.txt
re-writing: C:\opentelemac\Practise\Practise5\Telemac files\Steering.txt_2013-
05-21-12h20min20s\T2DCAS
copying: Boundaryconditions.cli C:\opentelemac\Practise\Practise5\Telemac f
iles\Steering.txt_2013-05-21-12h20min20s\T2DCLI
copying: Geometry.slf C:\opentelemac\Practise\Practise5\Telemac files\Steer
ing.txt_2013-05-21-12h20min20s\T2DGEO
copying: telemac2dv6p2.dico C:\opentelemac\Practise\Practise5\Telemac files
\Steering.txt_2013-05-21-12h20min20s\T2DDICO


Running your simulation :
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


C:\opentelemac\Practise\Practise5\Telemac files\Steering.txt_2013-05-21-12h20min
20s\out_telemac2dv6p2.exe

LISTING OF TELEMAC-2D


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 6.2 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC

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

EXITING LECDON. NAME OF THE STUDY:
FinalHydroRuns

OPENING FILES FOR TELEMAC2D

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

READGEO1: TITLE= newSelafin

NUMBER OF ELEMENTS: 17359
NUMBER OF POINTS: 8907
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)

FONSTR (BIEF): NO BATHYMETRY IN THE GEOMETRY FILE
AND NO BATHYMETRY FILE. THE BOTTOM
LEVEL IS FIXED TO ZERO BUT STILL
CAN BE MODIFIED IN CORFON.

STRCHE (BIEF): NO MODIFICATION OF FRICTION

THERE IS 2 LIQUID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 207 , WITH GLOBAL NUMBER: 8887
AND COORDINATES: 200.0000 0.4000000
ENDS AT BOUNDARY POINT: 227 , WITH GLOBAL NUMBER: 8907
AND COORDINATES: 200.0000 19.60000

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 433 , WITH GLOBAL NUMBER: 21
AND COORDINATES: 0.000000 20.00000
ENDS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22
AND COORDINATES: 0.000000 0.000000

THERE IS 2 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22
AND COORDINATES: 0.000000 0.000000
ENDS AT BOUNDARY POINT: 207 , WITH GLOBAL NUMBER: 8887
AND COORDINATES: 200.0000 0.4000000

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 227 , WITH GLOBAL NUMBER: 8907
AND COORDINATES: 200.0000 19.60000
ENDS AT BOUNDARY POINT: 433 , WITH GLOBAL NUMBER: 21
AND COORDINATES: 0.000000 20.00000

CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM

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

ITERATION 0 TIME: 0.0000 S
USING STREAMLINE VERSION 6.2 FOR CHARACTERISTICS

Any help on why this is happening would be greatly appreciated.

Cheers

Danny
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8760

  • jmhervouet
  • jmhervouet's Avatar
Hello,

It could be that your case diverges, due to something wrong in your case, and an infinite velocity blocks the computation of characteristics. Put DEBUGGER = 1 in your steering file to have a better idea of the subroutine where it stops.

Regards,

JMH
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8761

  • dsc1r12
  • dsc1r12's Avatar
Hi JMH,

I do that and re run the simulation and it comes up with the following...

THERE IS 2 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22
AND COORDINATES: 0.000000 0.000000
ENDS AT BOUNDARY POINT: 207 , WITH GLOBAL NUMBER: 8887
AND COORDINATES: 200.0000 0.4000000

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 227 , WITH GLOBAL NUMBER: 8907
AND COORDINATES: 200.0000 19.60000
ENDS AT BOUNDARY POINT: 433 , WITH GLOBAL NUMBER: 21
AND COORDINATES: 0.000000 20.00000
BACK FROM FRONT2
CALLING CORNOR
BACK FROM CORNOR
CALLING MASBAS2D
BACK FROM MASBAS2D
CALLING CORFON

CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM

BACK FROM CORFON
CALLING CONDIN
BACK FROM CONDIN
CALLING PROPIN
BACK FROM PROPIN
CALLING FRICTION_CHOICE
BACK FROM FRICTION_CHOICE

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

ITERATION 0 TIME: 0.0000 S
CALLING HPROPA
BACK FROM HPROPA
CALLING DE CHPCON
BACK FROM CHPCON
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
USING STREAMLINE VERSION 6.2 FOR CHARACTERISTICS
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
BACK FROM CHARAC
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCON
BACK FROM CHPCON
CALLING UTIMP_TELEMAC2D
BACK FROM UTIMP_TELEMAC2D
CALLING CHPCON
BACK FROM CHPCON
CALLING HPROPA
BACK FROM HPROPA
CALLING PROPIN
BACK FROM PROPIN
CALLING CORSTR
BACK FROM CORSTR
CALLING BORD
BACK FROM BORD
CALLING PROPIN
BACK FROM PROPIN
CALLING WALL_FRICTION
BACK FROM WALL_FRICTION
BOUCLE 701 ISOUSI= 1
CALLING CHARAC
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8763

  • jmhervouet
  • jmhervouet's Avatar
Hello,

So it runs a few time steps (9 or 10 apparently), hence your initial or boundary conditions are wrong somewhere and the run stops with infinite velocities. Put a PERIOD FOR GRAPHIC PRINTOUTS = 1, ask only 9 time steps and you will get a result file that you can look at to uderstand what is wrong.

REgards,

JMH
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8764

  • dsc1r12
  • dsc1r12's Avatar
It seems I am unable to view the results when I do this, when I load the results file in Blue Kenue a WARNING message appears saying the following:


Reading SELAFIN file...
This appears to be a byte swapped file.
2D Grid Geometry:
8907 nodes, 17359 elements
Scanning Variable: VELOCITY UV
"ERROR: while loading first frame."
ERROR: while rewinding file:
C:\opentelemac\Practise\Practise5\Telemac files\Results.slf
WARNING
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8765

  • dsc1r12
  • dsc1r12's Avatar
I am able to open the ipobo in Blue Kenue, but not the results file.
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8769

  • dsc1r12
  • dsc1r12's Avatar
I have changed the boundary conditions to

PRESCRIBED ELEVATIONS = 0;0

INITIAL ELEVATION = 'CONSTANT ELEVATION'
INITIAL ELEVATION = 0

in the steering file, and now it seems to run, however at each iteration step there is a message saying SOLUTION X=0 BECAUSE L2-NORM OF B VERY SMALL: 0.000000. I have had a look on the forum for this and am wondering if it is anything to do with the 'Tracer'??
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8770

  • jmhervouet
  • jmhervouet's Avatar
Hello,

This is not a problem, in some configurations the linear system to solve leads to a right-hand side which is 0. and you get this warning.

Regards,

Jean-Michel Hervouet

P.S. Hello webmaster and dear colleagues ! This is my 1000th message, did I win something already or should I write an extra 1000 ?
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8785

  • dsc1r12
  • dsc1r12's Avatar
Ok thanks, I seem to have got it working by playing around with the boundary condition values.
The administrator has disabled public write access.

Simulation stops after 1 iteration 11 years 6 months ago #8791

  • Fauret
  • Fauret's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 44
  • Thank you received: 12
Jean Michel

I did noticed your 1000th post! Congratulation! and thank you for your faithfulness in monitoring the forum.

Thierry
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.