Gradual activation of SZ drainage

To prevent numerical oscillations the drainage constant may be adjusted between 0 and the actual drainage time constant defined in the input for SZ drainage. The option has been found to have a dampening effect when the groundwater table fluctuates around the drainage level between iterations (and does not entail reductions in the drain flow in the final solution).