Elastic wave calculations by the Fourier method

We introduce a two-dimensional forward modeling algorithim based on a Fourier method. In order to be able to handle the free surface boundary condition with the Fourier method, a new set of wave equations are derived which contain the stresses as unknowns instead of the displacements. The solution algorithm includes a discretization in both space and time. Spatial derivatives are approximated with the use of the Fast Fourier Transform, whereas temporal derivatives are calculated with second order differencing. The numerical method is tested against the analytic solution for Lamb's problem in two dimensions.