open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I am working on a hydro-sedimentary model with one input and two outputs. I would like to impose at the imput a time dependent solid discharges. For that, I created a file of boundary conditions for sisyphe with :
4 5 5 1.000 0.000 0.000 0.000 4 0.000 0.000 0.000 9247 824 # Q

I also put the key word: PRESCRIBED SOLID DISCHARGES : 0;0.006;0
And I created a liquid boundaries file : (here an extract)
T QG(2)
s m3/s
0 0.005
20 0.1
40 0.3
60 0.3
80 0.6

But it seems to me that neither my liquid boundaries file neither the prescribed solid discharges is read.

I use the V7P2 so in the printout i have only the maximum and minimum evolution and no the bedload flux so I think it doesn't work because the evolution doesn't change.

Do you have an idea about what is not working?
Should I add or modify a FORTRAN file?

Thanks,
Claire

Dear Claire,

could you please show your listing file with the part where the boundaries are written (before the first time step) and then 2 or 3 time steps.

You don't need a fortran file for time dependent solid discharges.

Did you look at your results file? Is the velocity flow in the right direction? Is there a flow? Can you see an inflow at the inflow boundary and an outflow at the outflow boundary? What about the friction? Are the friction coefficients correct and fit to your chosen friction law?

Best regards,
Rebekka

Thank you for your reply !

I put the listing file and 2/3 time steps on attachments.

Yes i look my results. I first did a hydraulic study and all was ok, the flow, the direction, the boundaries conditions. I use the manning-stickler law (3) with a spatialized coefficient.
It's when i try to put a solid discharged. If i only use the keyword prescribed solid discharges it runs but when i change the value it doesn't change the evolution. And when I put a liquid boundaries file it doesn't run.

Dear Claire,

in the listing I miss something like "bedload flux". It seems that sisyphe doesn't see that there are 3 open boundaries. There is no open boundary at all and so your evolution will not change if you change the PRESCRIBED SOLID DISCHARGE...
In the sisyphe steering file the boundary file has to be given again. May be there is something wrong?

Try the following values in the sisyphe boundary file:
input sisyphe QS: 4 5 5 1. 0. 0. 0. 4
outlet sisyphe: 5 4 4 0. 0. 0. 0. 4

Best regards,
Rebekka

Dear Rebekka,

I tried the following value :
input sisyphe QS: 4 5 5 1. 0. 0. 0. 4
outlet sisyphe: 5 4 4 0. 0. 0. 0. 4

But it doesn't change...

In the printout I have a initialisation of QSL is it the beadload flux that you looking for ?

Best regards,
Claire

I found my mistake, it was in the boundary file, I forgot to change one value...

Now the constant solid discharges works, I'm trying to impose a time dependent solid discharges. It seems that i can't have two different liquid boundaries files (I have one for the hydrodynamic conditions) but when I put the values for the solid discharges in the liquid boudaries files it doesn't take them..

I added one column with :
QGL(2)
m3/s
0.01
0.03
...


Do you know why it doesn't work ?
Thanks,
Claire Couderc

Dear Claire,

good to hear that one problem is fixed

You need 2 files for liquid boundaries conditions. One for telemac which is set in the steering file for telemac and one for sisyphe which is in the sisyphe steering file. Both need the time for the first column.

Best regards,
Rebekka

Dear Rebekka,

Yes, thank you for your help !

I created my two liquid boundaries files: one for telemac with :
#
T SL(1) Q(2) SL(3)
s m m3/s m
0 91.39 36.21 95.5
....

And one for sisyphe with :
T QG(2)
s m3/s
0 0.002
...

But it doesn't run. I have this message :
ITERATION 0 TEMPS : 0.0000 S

BIEF_DESIMP
END BIEF_DESIMP
SISYPHE COUPLE AVEC : TELEMAC2D
END_INITIALIZATION
LE FICHIER DES FRONTIERES LIQUIDES CONTIENT
6 LIGNES AVEC :
SL(1) Q(2) SL(3)
USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS
TIME_LOOP
CONDIM_SISYPHE
END_CONDIM_SISYPHE
CONLIT
LA PREMIERE VARIABLE DOIT ETRE LE TEMPS T
DANS LE FICHIER DES FRONTIERES LIQUIDES

PLANTE : ARRET DU PROGRAMME APRES ERREUR
RETURNING EXIT CODE: 2

It seems that there is a problem with my liquid boundaries file in sisyphe but I put the time as the first variable.

Best regards,
Claire

Dear Claire,

you are right: somehow your liquid boundary file in sisyphe is not read correctly. In the file read_fic_conc.f the error message occurs if the "T" is not read.
Some versions are a little sensitive concerning ascii-files...
You can try the following things:
- avoid any comments in the file
- use the dos2unix command (or something similar in your editor) in order to be sure it is linux format (if you are working on linux)
- try to add a column for each liquid boundary (it will not be taken into account, but I doubt that this is the solution...)
- create a new file from the scratch with only 4 lines:
T QG(2)
s m3/s
0 0.002
1000000 0.002

If non of it works, I have no other idea than to go into the suboutine and print what the subroutine reads... but for this you need some fortran knowledge.

Best regards,
Rebekka

Dear Rebekka,

I have tried all the things but it doesn't work..
Is it difficult to print what the subroutine reads ?
Do you think that with an other version it will work ?

Best regards,
Claire