TELEMAC-MASCARET Forum: How to define non-erodable area within a channel (1/1)

Welcome, Guest

TOPIC: How to define non-erodable area within a channel

How to define non-erodable area within a channel 6 years 7 months ago #29688

Karki	Dear Friends, I am trying to run a coupling of Telemac2d and Sisyphe in a meandering channel. The inlet and outlet straight portion of the channel is a rigid bed (as shown in the attached jpeg file) and the curved portion is an erodible bed. In such case how/ where can I define those specific portions as a rigid bed or non-erodable bed? Thank you. Attachments: channel.jpg (64KB)
	The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29690

cyamin OFFLINE openTELEMAC Guru Posts: 997 Thank you received: 234	Hello, Your case is very simple to implement. Have a look at Sisyphe's example rigid-bed. User_fortran subroutine noerod.f does exactly what you want. The principle is to define the two extreme abscissae, X1,X2 where: DO I=1,NPOIN IF(X(I).LE.X1.OR.X(I).GE.X2) THEN ZR(I)=ZF(I) ENDIF ENDDO Regards, Costas
	The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29691

Karki	Dear Cyamin, Thank you very much for your help. I prepare the noerod.f file. Is this Ok? My Non-erodable domain are from X=-0.5 to X=0 and from X=4 to the end of the channel (X=5.5) Regards. Attachments: noerod.f (3KB)
	The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29692

cyamin OFFLINE openTELEMAC Guru Posts: 997 Thank you received: 234	It would be better if you grab a noerod.f file from the sources of your installation and work from there to adapt it like the example. The examples may define other parameters that may not apply to your case. But is your case it is just fine. Take note that in the example, the rest of the channel is set at -0.1m of sediment thickness. Costas
	The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29694

Karki	Dear Costas, I tried as your suggestion. But I encountered error while running the simulation. Can you please mention what is the problem here? The files are attached here. Thank you. Regards. Attachments: noerod_2018-04-19.f (3KB) error__.jpg (95KB) steer_case1.cas (3KB) steer_sisyphe.cas (4KB)
	The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29695

cyamin
OFFLINE
openTELEMAC Guru
Posts: 997
Thank you received: 234

First of all, when you define a sisyphe fortran file in a coupled computation, and there is no Telemac2D fortran file, then you need to define the fortran file in T2D cas file.

Regarding you noerod it self, XMAX is not declared (as in the example) and will give an error. Just replace it with its value:

      CALL OV('X=Y+C   ',ZR,ZF,ZF,-100.D0,NPOIN)         
        DO I=1,NPOIN
           IF(X(I).LE.0.D0.OR.X(I).GE.4.D0) THEN
             ZR(I)=ZF(I)
           ENDIF
        ENDDO

Costas

The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29729

Karki

Dear Costas,

I copied the file telemac2d.f from telemac directory and defined it as a fortran file in t2d.cas, but I got error.

I tried to see some examples for defining a fortran file in t2d.cas but couldn't find any.

Can you please mention how to define the fortran file and do we need to change any script in that fortran file?

The error msg is attached.

Thank you.

Attachments:

12456.jpg (67KB)

The administrator has disabled public write access.

How to define non-erodable area within a channel 6 years 7 months ago #29732

cyamin OFFLINE openTELEMAC Guru Posts: 997 Thank you received: 234	No, I meant to define noerod.f as FORTRAN FILE in the Telema2d cas file. This needs to be done when coupling Telemac2D with Sisyphe. All fortran files need to be defined from within the Telemac2d cas file. I hope I made my self clear. Costas
	The administrator has disabled public write access. The following user(s) said Thank You: Karki