This is a theoretical seismology paper which implements known algorithms for surface wave propagation across single, vertical discontinuities in 2-D media. Seismic surface waves undergo mode conversion at lateral discontinuities and this paper describes a Python software package for accurate computation of these conversion effects, a novel feature not available in the existing seismological literature. This work was done during my PhD research in order to validate a finite-difference wave solver.