open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello

I'm modelling a river with Telemac 2D (v6p3) coupled with Sisyphe (v6p3) for the bedload transport of a sediment of 6 class of diameters.
When I run only the hydrodinamic module (Telemac2d) everything goes well.
When I coupled Telemac2d with Sisyphe, and consider only 1 sediment class everything goes well.

But, when I include the 6 classes of sediments and modify the conlit.f to give the preescribed solid discharge for each class, the liquid flowrate discharge at the inlet start with the one I set but in the next time step the liquid flowrate is zero.
If I used the original conlit.f file, and set Q2BOR in the .cli file, the liquid discharge is ok, but the bedload flux at the inlet is zero.

I attached the result file so you can see what is happenning. (It is like a pulse of flux)
Also I attached the boundary condition .cli file for T2d and sisyphe, and the conlit.f and noerod.f file that where modify.


The modify lines of the conlit.f are:

DO K=1,NPTFR
!
I = NBOR(K)

.
.
.

if (BOUNDARY_COLOUR%I(K) .ge. 1957 .and. BOUNDARY_COLOUR%I(K) .le. 1988) then
LIQBOR%I(K)=KENT
QBOR%ADR(1)%P%R(K)=0.0000009D0
QBOR%ADR(2)%P%R(K)=0.0000009D0
QBOR%ADR(3)%P%R(K)=0.0000009D0
QBOR%ADR(4)%P%R(K)=0.0000009D0
QBOR%ADR(5)%P%R(K)=0.00000045D0
QBOR%ADR(6)%P%R(K)=0.00000045D0
end if

! EXAMPLE 2: IMPOSED BED EVOLUTON
!
! LIEBOR%I(K)=KENT
! (LIQBOR%I(K)=KSORT IS DONE IN SISYPHE.F)
! IF(LIEBOR%I(K).EQ.KENT) THEN
! EBOR%ADR(1)%P%R(K)=1.D-4
! EBOR%ADR(2)%P%R(K)=1.D-4.....
! ENDIF
!
ENDDO

.
.
.

The files attached are:

Geometry_Mf.slf = the geometry file containing the BOTTOM.
Bc_MF_Hydro.cli = Boundary Condition for Telemac2d.
BC_MF_THS.cli = Boundary Condition for Sisyphe.
Hydro_Mf.cas = T2d.cas
THS_MF.cas = Sisyphe.cas
liquid_MF.lqd = Liquid boundary file.

It is important to mantion that NPTFR = 1957 - 1987 Inlet boundary condition, and NPTFR=3959 - 3988 is the outlet BC.



Best regards,

Gonzalo

The missing files.

the missing files

geometry.slf file

geometry.slf missing file