Microscopic Modeling of Metasurfaces by the Mixed Finite Element Numerical Mode-Matching Method