49
where
k
At K
i-1/2,1
A = ( ) (54)
i,l 2 k
Az C
1.1
k
At K
i+1/2,1
D = ( ) (55)
i ,1 2 k
A z C
U
k
At K
i,1+1/2
E = ( ) (56)
i,l 2 k
A r C
1,1
C = 1 A D E (57)
1,1 1.1 1,1 1,1
A t
k
F = h + _
1,1 1,1
A Z
k k
K k
i-1/2,1 i+1/2,1
)
k
C
1,1
At
k
S
k i ,1
C
1,1
(53)
Updating of the soil parameters, surface infiltration, mass balance
and time steps were implimented as described for the one-dimensional
soil water flow model.
The system of equations for the two dimensional model produces a
banded matrix. The matrix is a five banded matrix with the form of