Ручной расчет ур-я неразрывности 1d Сен-Венана, замкнутого формулой Шези (без инерциальных членов).
Я использовал 3х точечную полуявную схему (некоторый апгрейд схемы Лакса). Которую вычитал у Грушевского. (в обеих книгах в формуле (5.11) есть очепятка. Предлагаю каждому вывести формулу самостоятельно и найти ошибку.

Выводы:
1. переход на следующий временной слой волновой плоскости по dA/dt, и здесь dQ/dx играет большую роль.
2. до того как расходы появятся на внутренних узлах, поток бежит по сухому руслу. Я тут проявил смекалку, как это обыграть в рамках схемы. :)
3. очень-очень качественная умственная тренировка на понимание Сен-Венана. Рекомендую всем.