open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone!

Still working on my case, I managed to find ways to solve some of my previous problems. Here's a new one I got lately.


To help you have a better look at my problem, here's the context: I have a dam located at the end of a small portion of a river. I've set, using the NOEROD subroutine, the dam as non-erodable (ZR = ZF-0.1 at this place, whereas ZR = ZF-100 otherwise), and I'm using a single class of sediment with the following characteristics: 30E-6 m of diameter, transported by bedload and suspension.

After more than 150 000 time steps, I end up with an error on a certain node:(If I understood correctly, the English version would be "BEWARE ALL LAYERS EMPTY: NODE I = 2526")

Sometimes, with a slightly different bathymetry on the input, I manage to go a bit beyond those 150k time steps, but this message shows up at some point anyway.
I am a bit puzzled by this error as I can't seem to find how to correct it. You can find, attached below, the steering files ("sis.cas" and "t2d_raf.cas"), the fortran file I use("conlit_inverse.f") as well as the boundary files ("t2d_conlim.txt" and "sis_conlim.txt").

I have also attached a picture of the bathymetry at the glitch's localization.

If anyone has any idea on this problem, I'd love to hear it! :)
Thanks in advance!

I'm re-upping this topic, because I still have that problem when I use a single class transported by suspension only.

The simulation runs fine until I end up with the following error:


It seems this issue occures on points surrounded by two highly different bathymetry, as shown in the example of my first post and the new one I'm joining here today (129m to 120.5m).

Does anyone knows what is the cause of this error and how can I fix it?
Thanks in advance!

Hi lpilavoine, sorry for the re-upping, but did you manage to find which error occurred ? And did you fix it ?
Thanks you