open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I am preparing a Sedimentation model (For Both bedload and suspended load and Two Sediment grain size) for the coastal region. In my model, There are Four boundaries.
1. Oceanside Boundary
Input - Tidal height input and no condition for Sediment
for that, I have used following boundary condition-
5 4 4 0.000 0.000 0.000 0.000 4 0.000 0.000 0.000 1368 1

2. Two boundaries for two rivers
Input - Height of water level and Time varying Sediment Input
Que 1 - Can I use a liquid boundary file for Sediment input?
Que 2 - If Ans 1 is YES, What is the format of liquid boundary file?
for that, I have used following boundary condition-
5 4 4 0.000 0.000 0.000 0.000 5 0.000 0.000 0.000 1368 1

3. Fixed wall
for that, I have used following boundary condition-
2 2 2 0.000 0.000 0.000 0.000 2 0.000 0.000 0.000 1405 38

Que 3 - Should I need to make two separate boundary condition file for Hydrodynamic model and Sisyphe model?

---

I am attaching both control files for constant Sediment input at river boundary.
Boundary condition for river boundary (Constant sediment, value 0.03)-
5 4 4 0.000 0.000 0.000 0.000 5 0.030 0.000 0.000 179 4932

Que 4 - For that attached model, Simulation is terminated after some time. Please let me know if I am doing something wrong.

Thanks In advance. Please help me.

Hi,have you solved your problem? I face the same problem. I still don't know the fomat of a time-series suspended load. I am using v8p1r1 coupling with gaia.
Can anyone give some advises? thanks.

Hi,

I am not very familiar with suspended load, but can advice you concerning bed load:

Que 1: yes you can use the LIQUID BOUNDARIES FILE for sisyphe. You need to have the keyword in your sisyphe steering file and you need the keyword PRESCRIBED SOLID DISCHARGES even if it is not used.
The liquid boundaries file should looks like the following:
# comment line
#
T QG(1) QG(2)
t m3/s m3/s
0 5910 0 # 0h
900 5870 0 # 1h
1800 5840 0 # 2h
2700 5810 0 # 3h

QG(1)=bed load in m3/s for the first open boundary,
CG(1)=concentration in g/l for the first open boundary. The number in brackets is the number of the open boundary. You have to be sure that LIQBOR (for bedload) and LICBOR (for suspended load) is set to kent (=5) for the wanted boundary nodes.

In case of multigrain the initial fraction for each class (given in the steering file) is used. So If you have 2 classes with initially 50 % each and your QG from the liquid boundaries file is 100m3/s you will have an input of 50 m3/s each.

The type ob boundaries
2nd column is LIUBOR/LIQBOR 8th column is LITBOR/LIEBOR/LICBOR

Inlet Sisyphe QS : 4 5 5 1.0 0.0 0.0 0.0 4
(beware the 1.0 in the 4th column)
Outlet Sisyphe : 5 4 4 0.0 0.0 0.0 0.0 4
LIEBOR need to be 4 (no setting of the evolution at the inlet). But LICBOR needs to be 5 for the inlet (setting of a concentration of at the inlet).
You need to program that in the conlit.f.

Que 3: you can use different cli-files for hydrodynamic and sediment. But I think in your case it is not needed.

Que 4: If a simulation run is terminated, it is best to look in the listing. Ususally there is a message why the simulation broke down. If you run in parallel it can happen that the message is only in one of the PE-files.

Best regards,
Rebekka