A MATLAB code automatically performs partial curve matching of 1D apparent resistivity data recorded with the Schlumberger electrode array configuration. The two-layer master and auxiliary curves are used to systematically match through the branches of data extracting the corresponding model properties. Partial curve matching is a classical interpretation procedure of the sounding curve, which has been done manually. Results from the manual and automatic procedures are compared. The matched geoelectric models from the automatic process are retrieved more quickly, and these results are consistent because the process is digitalized and are not dependent on human numerical accuracy judgment. Magnitudes of random noise affect the final matched model parameters, yet these values are sufficient to be initial models for subsequent nonlinear inversion. It is hoped that for an inversion workflow, the code can be included to automatically find an initial resistivity model.