open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all,
I am modeling a dam break in a flume with the reservoir. upstream boundary condition is closed boundary condition(wall) . right after the reservoir is a gate.The gate located at 21.5 from the upstream end of the flume.upstream of the gate is non erodible. So the model consists of two sections. when X is less than 21.5 meter it is non erodible and the other area (when X is greater than 21.5) is erodible.I just add a few lines in the subroutine. I would appreciate if you could correct me if I am wrong.Thanks

! RIGID BEDS POSITION
!

---

!
! DEFAULT VALUE: ZR=ZF-100.D0
!
CALL OV('X=Y+C ',ZR,ZF,ZF,-100.D0,NPOIN)
XMAX=21.5D0
DO I=1,NPOIN
IF(X(I).LE.XMAX) THEN
ZR(I)=ZF(I)
ENDIF
ENDDO

Hello,

Yes it is correct, X is here declared as a double precision array.

With best regards,

Jean-Michel Hervouet