open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I'm modelling a river with bedload transport of 6 clasess of sediments diamters using Sisyphe coupled with telemac2d.
I'm trying to set zero evolution at the inlet boundary condition, but i haven't had success.

I tried with:
- Setting in the boundary condition file, LIEBOR=5, and EBOR=0 at the boundary condition.

It didn't works so I modify de conlite.f file:

LIEBOR%I(K)=KENT
LIQBOR%I(K)=KENT
if (BOUNDARY_COLOUR%I(K) .ge. 5510 .and. BOUNDARY_COLOUR%I(K) .le. 5556) then

QBOR%ADR(1)%P%R(K)=0.013D0
QBOR%ADR(2)%P%R(K)=0.013D0
QBOR%ADR(3)%P%R(K)=0.013D0
QBOR%ADR(4)%P%R(K)=0.013D0
QBOR%ADR(5)%P%R(K)=0.0065106D0
QBOR%ADR(6)%P%R(K)=0.0065106D0
EBOR%ADR(1)%P%R(K)=0.D0
EBOR%ADR(2)%P%R(K)=0.D0
EBOR%ADR(3)%P%R(K)=0.D0
EBOR%ADR(4)%P%R(K)=0.D0
EBOR%ADR(5)%P%R(K)=0.D0
EBOR%ADR(6)%P%R(K)=0.D0
.
.
.

It didn't work. For that reason I modify the noerod.f file:

CALL OV('X=Y+C ',ZR,ZF,ZF,-100.D0,NPOIN)
do K=1,NPOIN
if (BOUNDARY_COLOUR%I(K) .ge. 5510 .and. BOUNDARY_COLOUR%I(K) .le. 5556) then
ZR(K)=ZF(K)
end if
end do

It didn't work too.

what could be happening?

I attached the boundary condition file, noerod.f, conlit.f file and the .cas file.

Best regards,

Gonzalo

Hello,

I think you should do the loop along the boundary nodes,
so instead of DO I=1,NPOIN -> DO I=1,NPTFR and adding NPTFR in
USE DECLARATIONS_SISYPHE, ONLY : BOUNDARY_COLOUR, NPTFR

Have also a look in the forum, there are some nice examples around for the noerod.f subroutine e.g. for directly editing in BlueKenue and reading.

Hope this helps.
Clemens

Hello,

Yes, Clemens is right, the loop should be:

USE DECLARATIONS_SISYPHE, ONLY : BOUNDARY_COLOUR,NPTFR,MESH

do K=1,NPTFR
if (BOUNDARY_COLOUR%I(K) .ge. 5510 .and. BOUNDARY_COLOUR%I(K) .le. 5556) then
I=MESH%NBOR%I(K)
ZR(I)=ZF(I)
end if
end do

With best regards,

Jean-Michel Hervouet

Dear Jean-Michel Hervouet,

Did you mean that add those line in the Noerod.f file than it is OK? your mean use 'DECLARATIONS_SISYPHE.f' it is mean in steering file. add 'Fortran File = 'DECLARATIONS_SISYPHE.f' '

do K=1,NPTFR
if (BOUNDARY_COLOUR%I(K) .ge. 5510 .and. BOUNDARY_COLOUR%I(K) .le. 5556) then
I=MESH%NBOR%I(K)
ZR(I)=ZF(I)
end if
end do



I have another question related to conlit.f

For the time series solid discharge, how should we set in conlit.f

For me, I use

IF(NIT.GE.500.0D0 .AND.NIT.LT. 1000.0D0) THEN

QBOR%ADR(1)%P%R(K)=1.D-3
...
END IF

but this is not suitable for the long term time series.

Do you have suggestion on this point-

Best Regards

Yao

Thanks!!!!
That solve the problem.