open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone,

I am trying to simulate an overtopping experiment over an earthdam embankment. My simulation worked fine I was just having some issues at my boundary section due to some unwanted high oscillations that i was trying to resolve. By doing so i created i new problem which error is displayed in ( error.png). Ive changed some values in the condin.f ( X1/YI/X2/Y2) and thats when the problem occured. Changing them back to the initial value as displayed in the file condin.f ive loaded up the problem still remains. I hope somebody could help me out. Thank you
Kind regards
Christian

Hello,

in order to provide special initial conditions you should set the keyword:
INITIAL CONDITIONS ='SPECIAL'.
I don't know exactly why you need the corsui.f subroutine but in your case for a rectangular channel and if the initial water depth is depending only on the x-coordinate you could prescribe the water depth in condin.f as follows, where e.g. an x-ccordinate is 5.0 and the water surface level is 0.3m.

DO I=1,NPOIN
IF(X(I).LE.5.0D0) THEN
H%R(I) = 0.D0
ELSE
H%R(I) = 0.3D0-ZF%R(I)
ENDIF
END DO

By default, the velocity components U and V are zero anyway, see above in condin.f.

Hope this helps,
Clemens