open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I'm using SISYPHE coupled with TELEMAC-2D. How can I use a sediment concentration for calculating additional terms in TELEMAC-2D? Is there an example for this?

Hello,

Basically, you can add an instruction :

USE DECLARATIONS_SISYPHE, ONLY : CS

in the subroutine where you want to use the concentration. Here CS is the concentration of suspended sediment and it is a block because there may be several classes, thus the double precision value of concentration of point I and of class J is : CS%ADR(J)%P%R(I)

This can be done for everything that is in DECLARATIONS_SISYPHE, unless there is a clash of names with Telemac-2D.

Then you must be careful that CS is already initialised when you use it.

With best regards,

Jean-Michel Hervouet

Thank you Jean-Michel.

I've read the declaration of variable CS in module DECLARATIONS_SISYPHE. It says that this is sediment concentration in time step N (or my interpretation due to lack of knowledge of French language drew me in wrong direction). When using it in TELEMAC, I'm curious if this is the value from previous time-step or value computed at the end of propagation for current time-step?

Best regards,
Budo

Hello Budo,

The basic principle is that both coupled programs start from values at time n and compute values at time n+1 (explicit coupling). So in your case it is indeed CS at time n that is needed on telemac-2d side.

Regards,

Jean-Michel Hervouet

Hello Jean-Michel.

I know that TELEMAC-2D can use fractional step method for its calculations and was uncertain whether it is the same thing with SISYPHE.

You mentioned previously that SISYPHE can use various classes for sediment transport. In the file DECLARATIONS_SISYPHE.F there are two values that define number of classes: NSICLA (number of classes) and NSICLM (maximum number of classes). Is the value of variable NSICLM default for NSICLA or is the context of these two different?

Hello Budo,

Actually NSICLA corresponds to the key word NUMBER OF CLASSES. NSICLM is an internal maximum for allocating some arrays, hence giving an NSICLA greater than NSICLM will (hopefully, not checked)trigger a message that NSICLM must be increased or to change the keyword NUMBER OF CLASSES.

Regards,

Jean-Michel