  We formulate a coupled surface/subsurface flow model that relies on
hydrostatic equations with free surface in the free flow domain and on the
Darcy model in the subsurface part. The model is discretized using the local
discontinuous Galerkin method, and a statement of discrete energy stability is
proved for the fully non-linear coupled system.
