Theory
The two partial differential equations for water and heat are solved with an explicit forward difference method (Euler integration). This solution requires the soil profile to be approximated with a discrete number of internally homogeneous layers.
Slowly changing state variables are bypassed and changes of the integration time step are made during simulation to speed up execution times.
Difference approximation of soil heat and water flow
equations