This article describes the programs included in the COMPLOC computer program package that are designed to apply the source-specific station term (SSST) method to solve for local earthquake locations using P- and S- wave phase data. These programs can greatly improve the relative location accuracy of nearby events by applying empirical corrections for the biasing effects of three-dimensional velocity structure. They have been tested on data from both the Southern California Seismic Network (SCSN) and Northern California Seismic Network (NCSN) The SSST method (Richards-Dinger and Shearer 2000; Lin and Shearer 2005) works by assigning each station a travel-time correction that varies as a function of source position. This approach provides relative location accuracy comparable to master event or hypocentroidal decomposition (Jordan and Sverdrup 1981) methods within compact event clusters, but is applicable to distributed seismicity. It has some similarities to the double-difference algorithm (Waldhauser and Ellsworth 2000, 2002) and can be shown to give comparable results in tests on synthetic data (Lin and Shearer 2005). However, the SSST approach has computational advantages for big data sets because the location and station term parts of the computation are separate, so that large matrix inversions are not necessary. In addition, our implementation of SSST provides the option to use L1-norm misfit measures, which are more robust than least squares in the case of occasional timing errors or bad phase picks. We implement the SSST approach by selecting nearby events located within a sphere of specified radius r max around the target event. The station term for the target event is then computed as the median (or mean) residual of these events. Different results will be obtained depending upon the size of the cutoff distance r max. If r max is set to a large enough distance, then the SSST …