136
k+1/2 k+1/2
BR h CR h = DR
i ,1 i ,2 i ,1
where
k+1/2
1 At
K
i+1/2,1
k+1/2
CR =
(
)
h
2
1,2
2 Ar
C
1.1
BR = 1 +
CR
and
k+1
k+1
k+1
- AV h
+ BV h
- CV h
= DV
1-1,1
1,1
1+1,1
i
where
k+1/2
k+1/2
1 At
K K
1 At
k
1-1/2,1
i+1/2,1
DV = h + (
)
1.1
)
4 AZ
C
2 C
1,1
1,1
k+1/2
1 At K
i,1+1/2 k+1/2 k+1/2
+ ( ) ( h h )
2 1.2 1.1
2 Ar C
1,1
Updating of the soil parameters, surface infiltration, mass
(99)
(100)
(101)
(102)
(103)
balance
and time steps were implimented as described for the one-dimensional
soil water flow model.