The lower boundary is calculated from the analytical
solution of the sine variation at the soil surface and a mean value of
damping depth for the whole soil profile.
Constant heat flow
A constant heat flow is given by the value of the
parameter GeothermalFlow.