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 calculated only for the first node in the soil column.