Sources and sinks

There is a source/sink term for each computational node. These sink terms are calculated from the root extraction due to transpiration in the upper part of the unsaturated zone. The integral of the root extraction over the entire root zone depth equals the total actual evapotranspiration. Direct evaporation from the soil is not allowed in the Gravity Flow method.