open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I used a previous computation file in my Sisyphe simulation as a condition file (condinit.ser) to fix the values of my solid bed. I have a graded distribution of sediments (5 classes). When I ran my simulation setting up "Equilibrium sediment discharge" on the boundary, I realized that only the first class showed a solid discharge, the others 4 classes were at 0.

After un long search, I understood that it was the previous computation file that gave an error. It seems that the previous computation file erased the definition of the classes made on the steering file.

The solution had been ton replace the geo file by the condinit.ser file and delete the keywords linked ti the previous computation.

Hope it helps someone or change in the next version.

Hi again,

The solution exposed previously doesn't work correctly. I have now a right distribution between the classes in my inlet solid discharge but the rigid bed I had defined in my condinit.ser is ignored when I put it in my geo file.

Anyone have the solution or an idea?

I have two ideas maybe:
- Use a zonal variable in the selafin file (see 9.5 in the manual) to fix my erodable bed
- or change the subroutines init_compo, init_avail and/or init_sediment to use "INITIAL FRACTION FOR PARTICULAR SIZE CLASS" AVA0 with a computation continued...

Thank in advance for your help

I tried the first idea exposed in my previous post with a zonal variable (it could be set up also with a tickness variable, see 9.8 of the manual). It works properly

If we want to read the rigid bed only from the previous computation file and everything else from the steering file
you can do this by some small changes in init_avai.f and rescue_sisyphe.f.
In rescue_sisyphe you need to avoid a stop because the water depth is not in the file. And in init_avai the if-way for no previous computation file must be taken.
I add you a file containing init_avai and rescue_sisyphe in the trunk version (newer than v8p0r1) So you should check if both subroutines fit to your installation.
I hope this will work.

Best regards

Thank you for your response.
I had seen the note above in the ini-avail.f but I didn't understand how to change it.
"! NOTE: WHEN COMPUTATION CONTINUED INIT_COMPO MUST NOT
! CHANGE ES AND AVAIL

In my condinit.ser, I had the water depth, I was using the same condinit both for T2D and Sisyphe. Maybe it's not necessary to change the rescue subroutine in this cas?

I'll try your solution.

Best regards