open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I have a telemac 2d model and it's coupled with sisyphe in order to simulate the sediment movements in the domain. I also implement a roughness law in the model using the customized subroutine CORSTR_SISYPHE.

But recently I noticed that this subroutine doesn't change the roughness value CHESTR in Telemac 2D. And I tested two different parameters in CORSTR_SISYPHE but these two cases gave the same results and the values of CHESTR remian constant everywhere.

Then I checked the definition of CHESTR and I found it is defined in Telemac 2D and the same name is also defined in Sisyphe. So I'm wondering how does the CHESTR work when Telemac 2D is coupled with Sisyphe?

If I change it in subroutine CORSTR in Telemac 2D, it seems that Sisyphe uses this modified value too. But how about the opposite? If I want to change it in subroutine CORSTR_SISYPHE in Sisyphe, how can I pass the modified CHESTR from Sisyphe to Telemac 2D?

Thanks in advance!

Qilong

Hello,

I think I have found the solution for this problem. I add the line

USE DECLARATIONS_SISYPHE, ONLY: UNORM,HN,QSCL_C,CS,NPOIN
USE DECLARATIONS_TELEMAC2D, ONLY: CHESTR

in the subroutine CORSTR_SISYPHE then it only changes the value of CHESTR in Telemac 2D. And this modified value will be sent back to Sisyphe later.

Kind regards,
Qilong

Hello,

Yes, this is it.

JMH