(вследствие образования встречных потоков или отрицательных
форм рельефа), условия, при которых получена система (6.126), не
выполняются. В результате исследований был предложен алгоритм,
позволяющий избежать возникающие в связи с этим допущения: 1)
отсутствуют замкнутые понижения и плато; 2) в локальных мини-
мумах разрезов рельефа вдоль осей координат соответствующая
составляющая скорости считается равной нулю.
Для численного интегрирования уравнений двумерной кинема-
тической волны может использоваться метод расщепления с приме-
нением неявной конечно-разностной схемы. Последовательно ре-
шается система однородных разностных уравнений:
й-Й ‚ а0 -4®) —
=—(Х -1),
Л: 
Ах 
2( 
)
‚1.!‹&‹.] _;'; ‚ 
(1) 
)
‚‚‚Ш ы АуЧ =%(Х_1)'‚ 
(6.128)
(=1,2,.3к, =Бо)=12, 
, -,
где Й, ‚ — значения глубин на первом полушаге; 94° ‹1)(‚0›‚ 4;” -
величины расходов, способ вычисления которых меняется в зави-
симости от способа задания уклона поверхности водосбора; & и /, —
число точек, заданных вдоль осей х и у.
Эта система решается способом итераций.
На границе водосбора и в точках локальных максимумов
рельефа задается условие й = 0. Считается, что вдоль границ водо-
сбора проходит канал (или река), и расчет движения в них осущест-
вляется по уравнениям одномерной кинематической волны.
Для расчета влагопереноса в почве к системе (6.126) добавля-
ется уравнение одномерной диффузии (учитывается только верти-
кальное движение влаги) и тогда решается система уравнений:
!
д 
.
_д..]./.!_‚{„ёз& Ч)' = Х 
- 1 _СИЗіЗ 
х 
›Ч_у = СЙЗП 
У 
‚
дг 
дх 
ду 
|втайт| 
|втаа т
® - 220% 
- @)+ -) 1=|- 20)8 + к@) 
› 6129
166