Discontinuous Galerkin method for coupling hydrostatic free surface flows to saturated subsurface systems. (1st May 2019)