We derive and compare several finite-difference frequency-domain (FD-FD) stencils for points on or near a planar dielectric interface. They are based on interface conditions or from modifying Helmholtz equation. We present a highly accurate formulation based on local plane wave expansion (LPWE). LPWE-based compact stencil is an extension of the analytically obtained LFE-9 stencil as used by the method of connected local fields [6]. We report that merely using five points per wavelength spatial sampling, LPWE coefficients achieve better than 0.01% local error near a planar interface. We numerically determine that we have fourth to eighth-order accuracy in the local errors for LPWE stencils.