open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, everyone
I want to simulate one class of cohesive sediment. And I wonder if I must definite the LAW OF BOTTOM FRICTION and FRICTION COEFFICIENT both in the sisyphe steering file and the telema2d steering file?
1. If it is, why? In some examples, the LAW OF BOTTOM FRICTION and FRICTION COEFFICIENT in the sisyphe steering file and the telema2d steering file are different. I donot understand.

2. In my case, my bottom friction is read from external file. When I run it without coupling with sisyphe, it is work. But when coupling with sisyphe, itshows that the bottom friction is error.

Can anybody help me? Thanks!

Hi,

If you check the dictionary file for sisyphe, in the 'LAW OF BOTTOM FRICTION SECTION', you will find the note as mention below.


from which we can clearly see that while coupling telemac2d with sisyphe the friction coefficient is selected in the Telemac steering file, except when BED ROUGHNESS PREDICTION is set to YES

What is the error which you are getting when you run the simulation coupled with sisyphe.
Are you using the following statement to read the friction data file given as FORMATTED DATA FILE 1

With Best Regards,
Gourish

Yes, I use the SUBROUTINE STRCHE and add the fllowing statement in this file:

NFO1=T2D_FILES(T2DFO1)%LU
write(*,*) 'NPOIN=',NPOIN
DO I = 1,NPOIN
READ(NFO1,*) friction(I)
write(*,*) friction(I)
CHESTR%R (I) = friction(I)
enddo

I have not use BED ROUGHNESS PREDICTION=YES, and I know the friction coefficient is selected in the Telemac steering file while coupling telemac2d with sisyphe. The error and related files are attached, Thank you for help me.

Best Regards.

Hello,

A few remarks on your coding:

There is no need of an intermediate friction(I), you can write instead directly:

READ(NFO1,*) CHESTR%R(I)

this will avoid the problem of allocating the array friction.

Then the way you wrote the loop is meant for a file with one point given per line, is it what you have in the file ? If you have several values per line the program will read the first and will skip to the next line, and you must then change the coding.

I hope this brings you hints or solutions,

With best regards,

Jean-Michel Hervouet

Thank you,JMH
YEAH,I have only one value per line . And it is right without couplint with sisyphe. The problem occured when coupling with sisyphe, so I think problem is it failed in reading data file when compute the sediment.

Hello,

Could it be that you read the file twice ? In this case add a command:

REWIND(NFO1)

before the loop.

Regards,

Jean-Michel Hervouet

Dear, JMH
Yes, you are right. It is ok now. Thank you very much!