Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
GAIA / SISYPHE
SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS
Flat Threaded

TOPIC: SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26099

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
Hi everyone,

I am trying to run a coupled telemac2d-sisyphe simulation from previous computation files. I am using 1 class of cohesive sediments, with only transport in suspension, and only one bottom layer.

In telemac2d, my PREVIOUS COMPUTATION FILE contains

- VELOCITY U
- VELOCITY V
- WATER DEPTH
- BOTTOM

In zizyphe, my PREVIOUS SEDIMENTOLOGICAL COMPUTATION FILE contains

- WATER DEPTH
- BOTTOM
- RIGID BED
- CONC VOL CL1
- LAYER1 THICKNE

I have the impression I have enough information, to start the next simulation. However, the latter stops with the following message:

TITLE OF PREVIOUS COMPUTATION: created with pputils

 NAME: VELOCITY U        UNIT: M/S
 NAME: VELOCITY V        UNIT: M/S
 NAME: WATER DEPTH       UNIT: M
 NAME: BOTTOM            UNIT: M
 NAME: TRANS. N-S        UNIT:
 NAME: TRANS. NW-SE      UNIT:
 NAME: TRANS. W-E        UNIT:
 NAME: TRANS. SW-NE      UNIT:
 NAME: COVERAGE SP. 1    UNIT:
 NAME: EXTRA DRAG COEF.  UNIT:

 READ_DATASET : READ OF RECORD     1

 TIME OF RECORD:     300.0000     S

================================================================================
 ITERATION        0    TIME:  5 MN   0.0000 S   (      300.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     1 LIQUID BOUNDARIES:

 BOUNDARY    1 :
  BEGINS AT BOUNDARY POINT:        1 , WITH GLOBAL NUMBER:         1
  AND COORDINATES:     0.000000          -280.0000
  ENDS AT BOUNDARY POINT:        6 , WITH GLOBAL NUMBER:         2
  AND COORDINATES:     100.0000          -280.0000

 THERE IS     1 SOLID BOUNDARIES:

 BOUNDARY    1 :
  BEGINS AT BOUNDARY POINT:        6 , WITH GLOBAL NUMBER:         2
  AND COORDINATES:     100.0000          -280.0000
  ENDS AT BOUNDARY POINT:        1 , WITH GLOBAL NUMBER:         1
  AND COORDINATES:     0.000000          -280.0000
 STRCHE (BIEF): NO MODIFICATION OF FRICTION

 TITLE OF PREVIOUS COMPUTATION: created with pputils

 NAME: WATER DEPTH       UNIT: M
 NAME: BOTTOM            UNIT: M
 NAME: RIGID BED         UNIT: M
 NAME: CONC VOL CL1      UNIT:
 NAME: LAYER1  THICKNE   UNIT: M

 READ_DATASET : READ OF RECORD     1

 TIME OF RECORD:     300.0000     S
--------------------------------------------------------------------------------
 RESCUE : HYDRODYNAMIC FILE
--------------------------------------------------------------------------------
 INITIALISATION EN CAS DE COUPLAGE : PART=           0
 ERROR IN INIT-MIXTE:
 THE SUM OF THICKNESS OF BED LAYERS                  IS DIFFERENT FROM ERODIBLE BED THICKNESS

 PLANTE: PROGRAM STOPPED AFTER AN ERROR
 RETURNING EXIT CODE:            2

Does somebody has an idea of what I could do to make it work?
Attachments:
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26103

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
I checked the problem a bit further this morning.

In the subroutine INIT_MIXTE, it seems that the thickness of layer 1 is still zero everywhere in the loop at line 180:

DO J=1,NOMBLAY
  EST=EST+ES(I,J)
  CONC(I,J) = CONC_VASE(J)
ENDDO

even if the variable 1ES is in the PREVIOUS SEDIMENTOLOGICAL COMPUTATION FILE.

I thought it could be because the corresponding item in the array ALIRE in sisyphe.F is 0. So I tried turning all items of ALIRE to 1, but it did not change the problem...

Does somebody has an idea?
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26104

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
As a work around, as I only have one class and one layer, I recompute the layer thickness from the bottom and the rigid bed positions. However, I would appreciate advices to solve that problem in a more elegant way.

Thank you all for your future inputs :-)
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26107

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi Olivier

Your analysis shows there is probably a problem in the reading of ES1...
I will have a look into your case to confirm this

Regards
Christophe
The administrator has disabled public write access.
The following user(s) said Thank You: o.gourgue

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 6 months ago #26311

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
Hi Christophe,

Any news on this? I'd like to model 2 layers now, and my work around does not work anymore...
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 6 months ago #26318

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi Olivier

Sorry for the delay, I started to look into your case but ...
I just look again and find the problem, I hope...

It's a little bit vicious because the problem I found is just a problem of variable name in your sisyphe geometry file...
You have LAYER THICKNE and sisyphe is looking for LAYER THICKNES

Kind regards
Christophe
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 6 months ago #26320

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
Now I feel a bit stupid :-s But thanks Christophe!!!
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26108

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
I have a final question about these simulations from previous computations.

As a test, I compared two cases:

case A: 1 only simulation of 30 M2 tidal cycles
case B: 30 simulations of 1 M2 tidal cycle each, where the initial condition of one simulation is the final result of the previous one.

In theory, I should have exactly the same results at the end.

However, if I compare them, they are slightly different. On the pictures, the differences can barely be seen (see attachments). However, if I compare point by point, the relative error on the final bathymetry ranges from 1.5e-7 to 0.13.

A difference of the order of 1e-8 is probably due to truncate errors when exporting to files, but a difference of the order of 1e-1 is more surprising.

Is it to be expected? Or is a problem in my setup?
The administrator has disabled public write access.

SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS 7 years 7 months ago #26109

  • o.gourgue
  • o.gourgue's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 155
  • Thank you received: 11
With attachments
Attachments:
The administrator has disabled public write access.
Flat Threaded
Moderators: Pablo, pavans
Forum
System Modules
GAIA / SISYPHE
SUM OF THICKNESS OF BED LAYERS DIFFERENT FROM ERODIBLE BED THICKNESS

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