open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I was debugging the fortran routine as pointed by the my post www.opentelemac.org/index.php/kunena/16-...me-dry?start=10#8941

I have noticed that when the 'INITIAL CONDITION' is set to 'CONSTANT DEPTH'
and 'INITIAL DEPTH' to 0, the water depth at the previous time step HN%R(K) at any K (node point) is always equal to zero. I checked even with H%R(K) which show the same values as HN%R(K).

This is the reason as why the program (www.opentelemac.org/media/kunena/attachm...2535/intake_fort.zip) cannot get a proper intake discharge, because at all time step HR%R(K) is equal to zero.
However on the completion of simulation the results in water depth, shows proper values which are not equal to zero.

Is the value of HN becoming zero and expected behavior when "INITIAL CONDITION" is set to "CONSTANT DEPTH".
If yes then which variable should I check for the water depths at previous time step.

Thanks
Gourish

Hello,

The answer to the question : Is the value of HN becoming zero and expected behavior when "INITIAL CONDITION" is set to "CONSTANT DEPTH" is no, the only use of "INITIAL CONDITION" is in condin.f for finding the value of depth. Then after this H and thus HN may evolve, but if your discharge depends on the depth, everything is blocked to 0.

If you give an initial depth equal to 0., of course there will be no water at the beginning, hence problems to enter a discharge.

Regards,

JMH

Hi Jean-Michel,

I have found the mistake in the program which I had done, and also verified that HN becoming zero is not an expected behavior. The problem was in the fortran code, which I have corrected and posted it at
www.opentelemac.org/media/kunena/attachm...35/intake_fortV2.zip

The fortran code did not take care of initial depth when it is equal to 0, so since the discharge was also dependent on the previous discharge value, there was a persisting value of 0 discharge, because anything multiplied by zero remains the same.

Thank and with Best regards,
Gourish