open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to implement input bedload as a function of flow discharge. I know that this needs to be done in CONLIT.F. Unfortunately the flow discharge variable, Q, is not declared as part of the subroutine.

I have done the following but haven't had any luck to get the subroutine to recognize Q.

1) In SISYPHE.f, I modified "CALL CONLIT(MESH%NBOR%I,AT0)" to "CALL CONLIT(MESH%NBOR%I,AT0,Q)"

2) In INTERFACE_SISYPHE.f, I modified "SUBROUTINE CONLIT(NBOR,AT)" to "SUBROUTINE CONLIT(NBOR,AT,Q)" and "DOUBLE PRECISION, INTENT(IN):: AT" to "DOUBLE PRECISION, INTENT(IN):: AT,Q"

3) In CONLIT.f, I modified "CONLIT(NBOR,AT)" to "CONLIT(NBOR,AT,Q)" and added "DOUBLE PRECISION, INTENT(IN):: Q"

The error message that I got is "Error: Symbol 'q' at (1) has no IMPLICIT type".

Any help will be greatly appreciated.

Hi,

conlit uses the declarations_sisyphe, which includes the declaration
TYPE(BIEF_OBJ), TARGET :: Q
for the flow rate, as well as QU and QV for the vector components.

So you should be able to use them without changing the subroutine calls.

Thanks for the quick reply. I tried to just use Q first but didn't work. I think the problem was that I used "Q(I)" instead of "Q%R(I)".

I made the following modification in CONLIT.f just to test it out.

I = NBOR(K)
QBOR%ADR(1)%P%R(K) = Q%R(I)

The results, however, seem to indicate 0 flow.

BOUNDARY 1 BEDLOAD FLUX = -0.000000 ( M3/S >0 = ENTERING )

did I not get the boundary node number assigned properly?

Hi,

I'm facing the same problem. In the listing, bedload fluxes are correct (I can see what I wanted to impose) but in my result file (.slf) I can't see my bedload discharge at the inlet.
Did you find a solution ?
Regards,