Supplemental material: Untangling tectonic slip from the potentially misleading effects of landform geometry

Geosphere, August 2017, v. 13, p. 1310-1328, doi:10.1130/GES01386.1, Supplemental Item 3 - MATLAB script which performs two stages to estimate a 3D slip vector. First, parallel planes are fitted to a series of pairs of point sets, to measure orthogonal separations (diagnostic plots generated to assess quality of fit). Weighted or unweighted least squares minimization is then used to estimate a 3D slip vector, based on the orthogonal separations. Data is assumed to have the same units in all three dimensions (e.g., UTM meters).