open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,
I have a problem with my model. I create a new variable to create erodable and no erodable zone in my domain.
When run simulation i have this error


VARIABLE NOER M FOUND IN THE GEOMETRY FILE
STRCHE (BIEF): NO MODIFICATION OF FRICTION

THERE IS 3 ISLAND(S) IN THE DOMAIN

THERE IS 2 LIQUID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1519 , WITH GLOBAL NUMBER: 15925
AND COORDINATES: 1700650. 4856366.
ENDS AT BOUNDARY POINT: 2 , WITH GLOBAL NUMBER: 22465
AND COORDINATES: 1700672. 4856327.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 989 , WITH GLOBAL NUMBER: 1
AND COORDINATES: 1700490. 4856750.
ENDS AT BOUNDARY POINT: 1032 , WITH GLOBAL NUMBER: 314
AND COORDINATES: 1700502. 4856708.

THERE IS 5 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 2 , WITH GLOBAL NUMBER: 22465
AND COORDINATES: 1700672. 4856327.
ENDS AT BOUNDARY POINT: 989 , WITH GLOBAL NUMBER: 1
AND COORDINATES: 1700490. 4856750.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 1032 , WITH GLOBAL NUMBER: 314
AND COORDINATES: 1700502. 4856708.
ENDS AT BOUNDARY POINT: 1519 , WITH GLOBAL NUMBER: 15925
AND COORDINATES: 1700650. 4856366.

BOUNDARY 3 :
BEGINS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.
ENDS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.

BOUNDARY 4 :
BEGINS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.
ENDS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.

BOUNDARY 5 :
BEGINS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
ENDS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
INITIALISATION ICURVE= 1
INITIALISATION ICURVE= 1

STAGE-DISCHARGE CURVE: 1

Q= 0.0000000000000000 Z= 115.50000000000000
Q= 0.20000000000000001 Z= 115.55400000000000
Q= 0.40000000000000002 Z= 115.58300000000000
Q= 0.59999999999999998 Z= 115.60599999999999
Q= 0.80000000000000004 Z= 115.62600000000000
Q= 1.0000000000000000 Z= 115.64400000000001
Q= 1.2000000000000000 Z= 115.66000000000000
Q= 1.3999999999999999 Z= 115.67600000000000
Q= 1.6000000000000001 Z= 115.69100000000000
Q= 1.8000000000000000 Z= 115.70500000000000
Q= 2.0000000000000000 Z= 115.71800000000000
Q= 2.2000000000000002 Z= 115.73099999999999
Q= 2.3999999999999999 Z= 115.74400000000000
Q= 2.6000000000000001 Z= 115.75500000000000
Q= 2.7999999999999998 Z= 115.76700000000000
Q= 3.0000000000000000 Z= 115.77800000000001
Q= 3.2000000000000002 Z= 115.78900000000000
Q= 3.3999999999999999 Z= 115.80000000000000
Q= 3.6000000000000001 Z= 115.81100000000001
Q= 3.7999999999999998 Z= 115.82100000000000
Q= 4.0000000000000000 Z= 115.83100000000000
Q= 4.2000000000000002 Z= 115.84000000000000
Q= 4.4000000000000004 Z= 115.84999999999999
Q= 4.5999999999999996 Z= 115.86000000000000
Q= 4.7999999999999998 Z= 115.86900000000000
Q= 5.0000000000000000 Z= 115.87800000000000
Q= 5.2000000000000002 Z= 115.88700000000000
Q= 5.4000000000000004 Z= 115.89600000000000
Q= 5.5999999999999996 Z= 115.90500000000000
Q= 5.7999999999999998 Z= 115.91300000000000
Q= 6.0000000000000000 Z= 115.92200000000000
Q= 6.2000000000000002 Z= 115.93000000000001
Q= 6.4000000000000004 Z= 115.93800000000000
Q= 6.5999999999999996 Z= 115.94700000000000
Q= 6.7999999999999998 Z= 115.95500000000000
Q= 7.0000000000000000 Z= 115.96299999999999
Q= 7.2000000000000002 Z= 115.97100000000000
Q= 7.4000000000000004 Z= 115.97799999999999
Q= 7.5999999999999996 Z= 115.98699999999999
Q= 7.7999999999999998 Z= 115.99400000000000
Q= 8.0000000000000000 Z= 116.00100000000000
Q= 8.1999999999999993 Z= 116.00900000000000
Q= 8.4000000000000004 Z= 116.01600000000001
Q= 8.5999999999999996 Z= 116.02400000000000
Q= 8.8000000000000007 Z= 116.03100000000001
Q= 9.0000000000000000 Z= 116.03800000000000
Q= 9.1999999999999993 Z= 116.04500000000000
Q= 9.4000000000000004 Z= 116.05200000000001
Q= 9.5999999999999996 Z= 116.05900000000000
Q= 9.8000000000000007 Z= 116.06600000000000
Q= 10.000000000000000 Z= 116.07299999999999
Q= 10.199999999999999 Z= 116.08000000000000
Q= 10.400000000000000 Z= 116.08700000000000
Q= 10.600000000000000 Z= 116.09399999999999
Q= 10.800000000000001 Z= 116.10100000000000
Q= 11.000000000000000 Z= 116.10800000000000
Q= 11.199999999999999 Z= 116.11400000000000
Q= 11.400000000000000 Z= 116.12000000000000
Q= 11.600000000000000 Z= 116.12700000000000
Q= 11.800000000000001 Z= 116.13300000000000
Q= 12.000000000000000 Z= 116.14000000000000
Q= 12.199999999999999 Z= 116.14600000000000
Q= 12.400000000000000 Z= 116.15200000000000
Q= 12.600000000000000 Z= 116.15900000000001
Q= 12.800000000000001 Z= 116.16500000000001
Q= 13.000000000000000 Z= 116.17100000000001
Q= 13.199999999999999 Z= 116.17800000000000
Q= 13.400000000000000 Z= 116.18400000000000
Q= 13.600000000000000 Z= 116.19000000000000
Q= 13.800000000000001 Z= 116.19600000000000
Q= 14.000000000000000 Z= 116.20200000000000
Q= 14.199999999999999 Z= 116.20800000000000
Q= 14.400000000000000 Z= 116.21400000000000
Q= 14.600000000000000 Z= 116.22000000000000
Q= 14.800000000000001 Z= 116.22600000000000
Q= 15.000000000000000 Z= 116.23200000000000
Q= 15.199999999999999 Z= 116.23800000000000
Q= 15.400000000000000 Z= 116.24400000000000
Q= 15.600000000000000 Z= 116.24900000000000
Q= 15.800000000000001 Z= 116.25500000000000
Q= 16.000000000000000 Z= 116.26100000000000
Q= 16.199999999999999 Z= 116.26700000000000
Q= 16.399999999999999 Z= 116.27200000000001
Q= 16.600000000000001 Z= 116.27800000000001
Q= 16.800000000000001 Z= 116.28300000000000
Q= 17.000000000000000 Z= 116.28900000000000
Q= 17.199999999999999 Z= 116.29400000000000
Q= 17.399999999999999 Z= 116.30000000000000
Q= 17.600000000000001 Z= 116.30600000000000
Q= 17.800000000000001 Z= 116.31100000000001
Q= 18.000000000000000 Z= 116.31600000000000
Q= 18.199999999999999 Z= 116.32200000000000
Q= 18.399999999999999 Z= 116.32800000000000
Q= 18.600000000000001 Z= 116.33300000000000
Q= 18.800000000000001 Z= 116.33799999999999
Q= 19.000000000000000 Z= 116.34399999999999
Q= 19.199999999999999 Z= 116.34900000000000
Q= 19.399999999999999 Z= 116.35400000000000
Q= 19.600000000000001 Z= 116.35899999999999
Q= 19.800000000000001 Z= 116.36499999999999
Q= 20.000000000000000 Z= 116.37000000000000

CORFON (TELEMAC2D): 1 BOTTOM SMOOTHINGS

================================================================================
ITERATION 0 TIME: 0.0000 S
TELEMAC2D COUPLED WITH: SISYPHE
TELEMAC-2D: INTERNAL COUPLING WITH SISYPHE
PART 0 : INITIALISING SISYPHE
INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)

THERE IS 3 ISLAND(S) IN THE DOMAIN

THERE IS 4 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22250
AND COORDINATES: 1700671. 4856328.
ENDS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22250
AND COORDINATES: 1700671. 4856328.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.
ENDS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.

BOUNDARY 3 :
BEGINS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.
ENDS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.

BOUNDARY 4 :
BEGINS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
ENDS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.

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.

VARIABLE NOER M FOUND IN THE GEOMETRY FILE
STRCHE (BIEF): NO MODIFICATION OF FRICTION

INITIALISATION EN CAS DE COUPLAGE : PART= 0
QSCL REINITIALISED IN INIT_TRANSPORT
FOR CLASS 1
QS REINITIALISED IN INIT_TRANSPORT

================================================================================
ITERATION 0 TIME: 0.0000 S

SISYPHE COUPLED WITH: TELEMAC2D
USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS

For create a file fortran to define erodable and non erodable zone i search an old discussion www.opentelemac.org/index.php/assistance...e-polygon-city-model

Hi Argirò,

I do not understand the problem you have because the computation seems to have started: "USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS").

Beware that the BOTTOM is set to 0.

Regards,
Luc

Thank you Lux,
the simulation not run because the message in the last row write:

SISYPHE COUPLED WITH: TELEMAC2D
USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS
_____________
runcode::main:
:
|runCode: Fail to run


I don't understand because not run. Maybe the formule in noerod.f not correct?
I set in bluekeneu a new variable in selafin file. This variable i call NOER and set 0 the portion of domain that I would not erodable and 100 the rest of domain. In noerod.f i set the formula X=Y-Z.
I attach the fortran file, maybe you see the error.
I'm sorry but I use telemac recently.




! *****************
SUBROUTINE NOEROD
! *****************
!
& (H , ZF , ZR , Z , X , Y , NPOIN , CHOIX , NLISS )
!
!***********************************************************************
! SISYPHE V7P1 21/07/2011
!***********************************************************************
!
!brief FIXES THE NON-ERODABLE BED ELEVATION ZR.
!
!note METHODS OF TREATMENT OF NON-ERODABLE BEDS CAN LEAD TO ZF.
!note CHOOSE TO SMOOTH THE SOLUTION WITH NLISS > 0.
!
!history C. LENORMANT
!+
!+ V5P1
!+
!
!history N.DURAND (HRW), S.E.BOURBAN (HRW)
!+ 13/07/2010
!+ V6P0
!+ Translation of French comments within the FORTRAN sources into
!+ English comments
!
!history N.DURAND (HRW), S.E.BOURBAN (HRW)
!+ 21/08/2010
!+ V6P0
!+ Creation of DOXYGEN tags for automated documentation and
!+ cross-referencing of the FORTRAN sources
!
!history J-M HERVOUET (EDF R&D, LNHE)
!+ 21/06/2013
!+ V6P3
!+ Now ZR=ZF-100.D0 by default
!+ previous versions was erronneously ZR=-100.D0
!history Costas (telemac forum) - Jose Díaz (Integral)
!+ 04/2016
!+ V7P1
!+ Using FONSTR as a base adapted to V7P1
!+ this version uses NOER as the DEPTH of the erodible bed
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CHOIX |-->| SELECTED METHOD FOR THE TREATMENT OF RIGID BEDS
!| H |-->| WATER DEPTH
!| NLISS |<->| NUMBER OF SMOOTHINGS
!| NPOIN |-->| NUMBER OF 2D POINTS
!| X,Y |-->| 2D COORDINATES
!| Z |-->| FREE SURFACE
!| ZF |-->| BED LEVEL
!| ZR |<--| RIGID BED LEVEL
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
USE BIEF
USE DECLARATIONS_TELEMAC
! USE DECLARATIONS_SISYPHE , ONLY : SIS_FILES, SISGEO, IELM, MESH
USE DECLARATIONS_TELEMAC2D , ONLY : T2D_FILES, T2DGEO, IELMT, MESH
! USE DECLARATIONS_SISYPHE, ONLY: PRIVE
USE DECLARATIONS_SPECIAL
USE INTERFACE_HERMES
!
IMPLICIT NONE
! INTEGER LNG,LU
! COMMON/INFO/LNG,LU
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
INTEGER, INTENT(IN):: NPOIN , CHOIX
INTEGER, INTENT(INOUT):: NLISS
!
DOUBLE PRECISION, INTENT(IN):: Z(NPOIN) , ZF(NPOIN)
DOUBLE PRECISION , INTENT(IN):: X(NPOIN) , Y(NPOIN), H(NPOIN)
DOUBLE PRECISION , INTENT(INOUT):: ZR(NPOIN)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
! USER VARIABLES
INTEGER ERR,IERR,RECORD,I
DOUBLE PRECISION BID
TYPE(BIEF_OBJ) :: NOER
CHARACTER(LEN=8) :: FFORMAT
! CHARACTER, INTENT(OUT) :: T2D_FILES(T2DGEO)%LU
! CHARACTER(LEN=16) :: NAME
! REAL, ALLOCATABLE :: W(:)
! LOGICAL OK
! INTEGER I,ERR
!
!

---

!
! SET VARIABLES
!
RECORD = 0
FFORMAT = 'SERAFIN '

!
!

---

!

---

! RIGID BEDS POSITION
!

---

!
! DEFAULT VALUE: ZR=ZF-100.D0
!
! Read rigid bed level from geometry file (variable 'RIGID BED' which is
! erodible bed)
! Inspired by subroutine FONSTR
!
! LOOKS FOR THE RIGID BOTTOM ELEVATION IN THE GEOMETRY FILE
!
CALL BIEF_ALLVEC( 1 , NOER , 'NOER ' , IELMT , 1 , 2 , MESH )
! ALLOCATE(W(NPOIN),STAT=ERR)


IF(LNG.EQ.1) CALL FIND_VARIABLE(FFORMAT, T2D_FILES(T2DGEO)%LU,
& 'NOER ',NOER%R, NPOIN,
& IERR,RECORD=RECORD,TIME_RECORD=BID)
IF(LNG.EQ.2) CALL FIND_VARIABLE(FFORMAT, T2D_FILES(T2DGEO)%LU,
& 'NOER ',NOER%R, NPOIN,
& IERR,RECORD=RECORD,TIME_RECORD=BID)
IF((IERR.EQ.HERMES_VAR_UNKNOWN_ERR).AND.LNG.EQ.1) THEN
CALL FIND_VARIABLE(FFORMAT, T2D_FILES(T2DGEO)%LU,
& 'NOER ',NOER%R,NPOIN,
& IERR,RECORD=RECORD,TIME_RECORD=BID)
ENDIF
IF((IERR.EQ.HERMES_VAR_UNKNOWN_ERR).AND.LNG.EQ.2) THEN
CALL FIND_VARIABLE(FFORMAT, T2D_FILES(T2DGEO)%LU,
& 'NOER ',NOER%R,NPOIN,
& IERR,RECORD=RECORD,TIME_RECORD=BID)
ENDIF
!
! INITIALISES THE RIGID BOTTOM ELEVATION
!
! CALL OV('X=Y ',ZR,NOER%R,ZF,-100.D0,NPOIN)
! CALL OV('X=Y-C ',ZR,ZF,ZF,0.D0,NPOIN)

CALL OV('X=Y-Z ',ZR,ZF,NOER%R,0.D0,NPOIN)

! DO I=1,NPOIN
! IF(ZR(I).GT.ZF(I)) THEN
! WRITE(LU,*) I, ZF(I)-ZR(I)
! ENDIF
! ENDDO
! CALL OV('X=Y+C ',ZR,ZF,ZF,-100.D0,NPOIN)
!
!

---

! SMOOTHING OPTION
!

---

!
! NLISS : NUMBER OF SMOOTHING IF (ZF - ZR ) NEGATIVE
! DEFAULT VALUE : NLISS = 0 (NO SMOOTHING)
!
NLISS = 0
!
!

---

! CONTROL (CAN BE ACTIVATED IF ZR USER DEFINED...)
!

---

!
! DO I=1,NPOIN
! IF(ZR(I).GT.ZF(I)) THEN
! WRITE(LU,*) 'POINT ',I,' NON ERODABLE BED HIGHER THAN BED'
! CALL PLANTE(1)
! STOP
! ENDIF
! ENDDO
!
!

---

!
RETURN
END

You say that bottom is 0, but in selafin file in bluekeneu I read value

Ok, there is indeed a crash and the error is probably not written in the master processor listing.
Try to launch the simulation in a single processor or have a look to the listing of the guilty processor.

The bottom is set to 0 because Telemac does not recognize it in GEOMETRY FILE (it should have the correct name (BOTTOM) and unit(M)), see the master processor listing:

Thank you for reply,
I do a new geometry, i define variables and now the error about bathymetry there isn't.

But the simulation is not start. This is the last part of cmd processor.

THERE IS 4 SOLID BOUNDARIES:

BOUNDARY 1 :
BEGINS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22250
AND COORDINATES: 1700671. 4856328.
ENDS AT BOUNDARY POINT: 1 , WITH GLOBAL NUMBER: 22250
AND COORDINATES: 1700671. 4856328.

BOUNDARY 2 :
BEGINS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.
ENDS AT BOUNDARY POINT: 1554 , WITH GLOBAL NUMBER: 6422
AND COORDINATES: 1700603. 4856778.

BOUNDARY 3 :
BEGINS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.
ENDS AT BOUNDARY POINT: 1585 , WITH GLOBAL NUMBER: 21167
AND COORDINATES: 1700668. 4856794.

BOUNDARY 4 :
BEGINS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
ENDS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
VARIABLE NOER M FOUND IN THE GEOMETRY FILE
STRCHE (BIEF): NO MODIFICATION OF FRICTION

INITIALISATION EN CAS DE COUPLAGE : PART= 0
QSCL REINITIALISED IN INIT_TRANSPORT
FOR CLASS 1
QS REINITIALISED IN INIT_TRANSPORT

================================================================================
ITERATION 0 TIME: 0.0000 S

SISYPHE COUPLED WITH: TELEMAC2D
USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS
_____________
runcode::main:
:
|runCode: Fail to run

Have you tried to launch the simulation on a single processor?

To have more hints about the crash, you could try to set DEBUGGER=1 in the CAS file.

Hope it helps,
Luc

Yes, i lanch as single processor, now try with debbuger=1.
Thank you

With debbuger=1 the error is:

BOUNDARY 4 :
BEGINS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
ENDS AT BOUNDARY POINT: 1586 , WITH GLOBAL NUMBER: 48413
AND COORDINATES: 1700731. 4856812.
STRCHE (BIEF): NO MODIFICATION OF FRICTION

INITIALISATION EN CAS DE COUPLAGE : PART= 0
QSCL REINITIALISED IN INIT_TRANSPORT
FOR CLASS 1
QS REINITIALISED IN INIT_TRANSPORT

================================================================================
ITERATION 0 TIME: 0.0000 S

SISYPHE COUPLED WITH: TELEMAC2D
BACK FROM SISYPHE
CALLING CHPCONV
BACK FROM CHPCONV
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 OPTION STRONG
USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS
BACK FROM CHARAC OPTION STRONG
CALLING SIPHON
BACK FROM SIPHON
CALLING PROSOU
BACK FROM PROSOU
CALLING PROPAG
BACK FROM PROPAG
TRAITEMENT BANCS DECOUVRANTS
FIN DU TRAITEMENT BANCS DECOUVRANTS
CALLING CHPCONV
BACK FROM CHPCONV
CALLING SISYPHE, CHARRIAGE
_____________
runcode::main:
:
|runCode: Fail to run