open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi again!

For my previous problem, I was able to fix it by flatening the bathymetry near the output. I managed to simulate the duration I wanted to while using a certain size of sediments.

I want to study now two new classes of sediments, which sizes are 1E-3m and 2E-2m.
I have enabled the Bedload for these two classes, and everything seems to be going fine, except for one thing: the non-erodable domain.

When I try to set my dam as non-erodable with ZR = ZF-0.1 on it, I end up with the following error message:
Does anyone know how I could fix this?

I have attached to this post both my steering and my fortran file.
Thanks in advance!

Hi,

You can load the bathymetric points (txt) and the barrage_points_char.txt file.


Otherwise, try to limit the location of the file.

Example:
File_name = 'D: \ Travail \ MODSED_03_07 \ Telemac_char \ barrage_points_char.txt'

File_name = 'D: \ barrage_points_char.txt'

regards

Joaquin

Well your methods works, but when I set, for instance, ZR to ZF-1 by using a set of points, I end up with the error message I've written above.

That's only when I'm using the 2 classes and bedload. Otherwise, with the suspension class only, the evolution stops at the value written in the noerod file (for instance, if I wrote ZR = ZF-0.1, the evolution won't go beyond -0.1)

So it seems the NOEROD subroutine works well when you only have a single class of sediments....

I've tried with a class alone and had no problem, but when I try to simulate two different classes, it stops after a few iterations.

Here are the tests I did so far :
1) Simulation of 2 classes of the same diameter with an initial fraction of 50% each
2) Simulation of 2 different classes with an initial fraction of 100% for one and 0% for the other one
3) Removing the bathymetry of the dam and try the simulation in case 2)

The only thing was that I managed to make it work when I've set the initial fraction to 1 for the smallest class (1mm), but it ended up bugging anyway after a few hours.

Does anyone have any idea on this type of case? It seems that the fact the rigid bed is close to the initial bathymetry messes with the AVAIL value, if I understood this error message correctly in the layer.f subroutine :
Thanks in advance.

Hello,

I tried to do something like this several months ago, the problem could be that you are modeling non-uniform transport and in that case you should model BED LOAD or SUSPENDED LOAD, but not both simultaneously. (Please take a look at Page 40 of Sisyphe User Manual).

Try with BED LOAD = YES/true and SUSPENSION = NO/false, and vice versa, that could solve your problem.

I hope this helps.

Cheers,
José Andrés.

Oh, OK...

So that means that if you have more than a single class of sediment in your study, you can't have both Bedload and Suspension?

Thanks for your answer, I'll see how it goes on with only one transport mode.