Seismic anisotropy refers to the property of seismic waves to propagate at different wavespeeds depending on the direction of propagation. This property can be related to the coherent alignment of rock-forming minerals, which is thought to reflect the current dynamics or fossilized structure of Earth materials due to tectonic deformation. In the upper mantle, seismic anisotropy is most easily determined using the distortion of teleseismic body waves with a known initial polarity, typically core-refracted shear-waves (SKS, SKKS). SplitPy is a teleseismic shear-wave Splitting Toolbox based on the Matlab Tool SplitLab, but with modifications from Wustefeld et al (2008). Additional error surface implementation has been added, however these error surfaces have not been fully tested. The code produces output identical to those in Audet et al. (2016)