Stanford Gibson запилил огненный видос по моделированию гидравлики водосливов. (см пост 1 и пост 2).

Главная мысль: если на водосливе у вас условие свободного падения воды (водопад) то конечно-разностная схема не может раскидать момент импульса между ячеек. Поэтому используйте аналитическое уравнение водослива в SA/2d. Но с некоторым трюком (см в видео)

https://youtu.be/Pe6v-C5bOQo?si=hycqHbjreimc826Y