Add support for discharge scaling depth in drainage package

[Huite]

MODFLOW 6 supports drainage discharge scaling: when heads near the drainage elevation, they may be reduced by a scaling factor. See also the supplementary info chapter 3.

Cubic smoothing is used to ensure continuous discharge derivatives when the Newton-Rhapson formulation is used, which helps convergence.