  In this paper, we propose a new numerical scheme for the coupled Stokes-Darcy
model with Beavers-Joseph-Saffman interface condition. We use the weak Galerkin
method to discretize the Stokes equation and the mixed finite element method to
the Darcy equation. A discrete inf-sup condition is proved and optimal error
estimates are also derived. Numerical experiments validate the theoretical
analysis.
