open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

I am running telemac2D simulation for a 30 km stretch of the river with a given free surface in the part of the river. The actual mesh encompasses the flood plain also. I have prescribed the initial spatially varying free surface with following lines of code in CONDIN subroutine.

DO I1 = 1, NPOIN
if(MESH%X%R(I1) <= 4453430) then
MESH%Z%R(I1) = 73.12
else
MESH%Z%R(I1) = 73.12 - (MESH%X%R(I1) - 4453430)*SLP
end if
END DO

My question is, do I need to do anything else if I want to incorporate spatially varying free-surface.

Best regards,
Sumit

Hello Sumit,

If I have understood, you only want to prescribe free surface for initial conditions, do you?
In that case, in CONDIN, as written in the header, you have to give H (see the example for 'CONSTANT ELEVATION' option). I would have written rather:

DO I1 = 1, NPOIN
H%R(I1) = 73.12D0 - ZF%R(I1) ! FREE SURFACE = ZF + H
IF (MESH%X%R(I1).GT.4453430D0) THEN
H%R(I1) = H%R(I1) - (MESH%X%R(I1) - 4453430D0)*SLP
ENDIF
ENDDO

I suppose that SLP is a slope defined by a double precision before?

Hope this helps,

Chi-Tuan

Hi Chi-Tuan,

Thanks for the help, I was thinking that if I prescribe the free-surface it will compute the depth by itself with the freesurface and bottom bathymetry information, silly of me :) .

Thanks again,
Sumit