Fault contours are created by mapping the slip distribution along a fault surface. Fault contours have many useful applications in geophysics and geology including rock stress analysis, interpretation, and processing. Therefore, automatic calculation of fault contours from 3D data would be very valuable. As the first step in that direction, displacement is calculated from a simple 2D model of faulted seismic data. Making a stationary assumption, the displacement is calculated by fitting a smooth line to a cross-correlagram. The cross-correlagram is created by windowed cross-correlation across the fault. To fit this line, I use a non-linear optimization method that has similarities to Simulated Annealing. This removes most of the displacement which will later allow other methods to solve the non-stationary problem.