Before revision:  Our analysis of the stability properties of the proposed scheme reveals that a delicate scaling of the stabilization parameters is required for the equal higher-order case. This is also supported by numerical evidence from test experiments. Additionally, a geometrically robust estimate for the condition number of the stiffness matrix is provided. Numerical examples illustrate the implementation of the method and verify the theoretical findings .
After revision:  Numerical examples verify the theoretical findings .