open TELEMAC-MASCARET The mathematically superior suite of solvers

Thankyou. Unfortunately didn't seem to be the issue: head losses still aren't applied at the nodes...

Hi,

In addition of the mistake pointed out by Christophe, you just have to reverse the order of the DO and IF instructions in the 3rd DO loop.
Your case runs fine for me with the enclosed modifications.
I also did some optimisations on your code so as to remove unused variables and to keep only one DO loop instead of 3. I only look at the source.f subroutine.


Note also that you can stay in hydostatic version if you don't apply head loss to vertical component of momentum.

I hope it helps.

Laurent

Thankyou, Laurent!

That's great of you to do so. I see where I went wrong. It works well! Thanks for your help.