open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all,

I managed thanks to the forum to define uv initial conditions on the whole domain in 2D by modifying conlin.f

I would like now to do the same in 3D. I tried modifying condi3dh and condim and declaring it as my fortran file.
The code runs but my initial conditions do not appear in the result.

Any idea of the Reason (enclosed my fortran file).

ps: I am running v8p0r2

Thanks,

Thomas

Hello,

The initial condition you have implemented is only activated if using the keyword INITIAL CONDITIONS = SPECIAL, PARTICULAR or PARTICULIERES (see the condition IF...). Anyway, if you only do this modification, it is erased by the CONDI3DUVW subroutine called after in CONDIM subroutine.

What you should do is to edit the CONDI3DUVW subroutine and write your modified initial condition for U velocity component
CALL OS('X=C ',X=U,C=1.0D0)
rather than the
CALL OS('X=0 ' , X=U)
instruction.
Then, use INITIAL CONDITIONS = 'ZERO ELEVATION' to do the same as what you seem to want.

Hope this helps,

Chi-Tuan

Works perfectly,
great thanks!

Thomas