гроу — |СКВЫ В, 
› при АМ <,
О, при АЛ > Л, 
(5.10)
й = й)- АИ,.
где Иро — начальная глубина, а В, Го — ширина и уклон реки; С —
коэффициент Шези.
Внутренний итерационный цикл завершается при согласовании
элементов водного баланса и морфометрических характеристик
водоема, когда разница расчетных глубин между итерациями ме-
нее 0.01 м.
3. 
Расчет уровня водоема Н, = НГ0 —- ЛА/ и переход к п. 1 - для рас-
четов на слеДуЮЩИй момент времени 1 = 1 + Аг.
В результате, в ходе моделирования осуществляется пошаговая
корректировка морфометрических параметров акватории озера по
мере его осушения на основе оценки изменений водного баланса во-
доема.
ТбСТОВЬ1Й анализ такого алгоритма показал, что на сходимость ите-
рационного процесса подбора глубины и площади водоема сущест-
венно влияет расчетный шаг по времени А/. При этом, в случае бес-
сточного водоема, когда Р` = (, его гидрологический режим характе-
ризуются гораздо большей динамичностью и требует использования
меньшего расчетного шага по времени.
Рассмотренный алгоритм моделирования взаимодействия подзем-
ных вод с водоемами реализован в программе МССО (Гриневский,
2002).
Принципиально похожий алгоритм используется в расчетном па-
кете программы МойНоу БАКЗ [СБепе, Апйегзоп, 1993; Мегии, Копт-
Ко\у, 2000]. Здесь водоем может занимать несколько вертикальных
блоков модели в зависимости от его глубины, однако осушение аква-
тории водоема учитывается только с блочной дискретностью и осу-
ществляется «отключением» блоков модели с водоемом, когда его
уровень опускается ниже подошвы блока. Авторы также отмечают
возможность возникновения счетных осцилляций итерационного ре-
шения в зависимости от расчетного шага по времени и для их мини-
мизации используют средневзвешенную между временными интерва-
лами глубину водоема, рассчитываемую с использованием задаваемо-
го пользователем весового коэффициента.
Особого учета в расчетном алгоритме моделирования гидрологи-
ческого режима водоема требует случай «внутреннего» осушения ак-
ватории и деления единого водоема на части [Мегп&, КомКом, 2000]
— рис. 5.3. В этом случае итерационный расчет интегральной глубины
водоема, аналогичный рассмотренному в п. 2, должен производиться
по каждому вновь образованному водоему в отдельности, что может
приводить к различию уровней в них. В случае их повторного объе-
106