open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

does anybody know the syntax of the liquid boundaries file for multible suspended load diameters?

I tried

t CG(a,b1) CG(a,b2)

where a is the boundaries number and b the Diameter number.

It seems, that the twodimensional array is incorrect.

Thanks in advance

Matthias

Hi,

some more specification on my case:

Please see attached the Liquid boundaries file. It works wekk as long as I do not turn on "liquid boundaries file" in Sisyphe. Then CGL.f gives me an error and tells me about insufficient values:

"CG: MORE PRESCRIBED ELEVATIONS ARE REQUIRED, AT LEAST 3". Is seems to me, that this is a copy paste error, as the french translation talks about concentrations.

Kiind regards and thanks for help!

Matthias

Hi Matthias,

Unless I'm mistaken, this feature is not fully available in Sisyphe. Actually, with the liquid boundaries file, only one class of sediment concentration can be prescribed on several boundaries ( CGL(IFRLIQ,AT) ).
Some developments are missing at least in conlit.f, cgl.f and read_fic_conc.f

I don't know if Gaïa can support it.
By looking at your steering file, I saw that you're using Telemac-3D and Sisyphe for suspended load only (no bed load). So a solution (more accurate/suited) can be to consider these sediments directly in Telemac3d with its own sediment module SEDI-3D.
Moreover, I think you'll be able then to prescribed your 2 sediment class concentrations thanks to liquid boundary file (already implemented for tracers).

I hope it helps.

Regards,
Laurent

Hi,
here is a 'grep' of the JLM10's switch that I added for extended SISYPHE to a multi sediment discharge:

cgl.f: &(I,J,AT) !JLM10
cgl.f: FCT(5:5)=',' !JLM10
cgl.f: !FCT(5:5)=')' !JLM10
cgl.f: IF(J.LT.10)THEN !JLM10
cgl.f: WRITE(FCT(6:6),FMT='(I1)') J !JLM10
cgl.f: FCT(7:7)=')' !JLM10
cgl.f: ELSE !JLM10
cgl.f: WRITE(FCT(6:7),FMT='(I2)') J !JLM10
cgl.f: FCT(8:8)=')' !JLM10
cgl.f: ENDIF !JLM10
conlit.f: CBOR%ADR(I)%P%R(K)=CGL(IFRLIQ,I,AT)/XMVS !JLM10
interface_sisyphe.f: DOUBLE PRECISION FUNCTION CGL( I, J , AT) !JLM10 J
interface_sisyphe.f: INTEGER, INTENT(IN) :: I , J !JLM10 J
interface_sisyphe.f: LOGICAL, INTENT(INOUT) :: SEDCO(NSICLA) !JLM10

Easy ! Only three functions to change !

And then, you have to create a file like this:

#
# CONCENTRATION EN ENTREE
#
T CG(2,1) CG(2,2) CG(2,3) CG(2,4) CG(2,5) CG(2,6) CG(2,7) CG(2,8) CG(2,9) CG(2,10)
s g/l g/l g/l g/l g/l g/l g/l g/l g/l g/l
0 0 0 0 0 0 0 0 4.80E-03 4.80E-03 4.80E-03
900 0 0 0 0 0 0 0 4.80E-03 4.80E-03 4.80E-03
1800 0 0 0 0 0 0 0 4.80E-03 4.80E-03 4.80E-03
2700 0 0 0 0 0 0 0 4.80E-03 4.80E-03 4.80E-03



an example of the need to do it:
www.researchgate.net/publication/3287740..._boundary_conditions


I hope it helps !

Cheers !

Flanagan

Dear Flanagan,

thanks a lot for the great help. I inplemented the code. The only thing, I didn't know where to put the following line in interface_sisyphe.f

LOGICAL, INTENT(INOUT) :: SEDCO(NSICLA) !JLM10

I changed in cgl.f as well:

INTEGER, INTENT(IN) :: I , J !JLM10 J

I get the messages, that there are solid discharges missing. "QG: MORE PRESCRIBED SOLID DISCHARGES".

Do you have an idea?

Thanks a lot in advance

Matthias

You have to looking for the subroutine conlit or cgl in the SISYPHE interface.

Hi Flanagan,

thanks for the quick answer.

In the interface file I couldn´t find link to SEDCO(NSICLA) in the CGL or CONLIT function. When I turn it on, the code stops due to compilation error.

Regards

Hi,
it is probably not useful for your case: this variable is apparently a logical table referring to the nature cohesive or not of the sediment.

cheers

Hi,

thanks for the answer.

I fixed the bug! :)

The SEDCO(NSICLA) was not needed.

For the liquid Input file additional quantities for bed load is needed, because the flag does not know, if bed load or suspended load will enter.

The syntax is:

T QG(1) QG(2) CG(3) CG(1,1) CG(1,2) CG(2,1) CG(2,2) CG(3,1) CG(3,2)

Thanks once again

Matthias

Hi,
you probably need to change also the "qgl" function as well because you only can prescribe only one bedload flux for each boundary (same flux for every sediment class).

QG(1) => 1 is the number of the boundary

cheers