We present a method for modelling seismic wave propagation in a whole-earth model by solving the elastodynamic equations in 2-D cylindrical coordinates (r, θ) using the Fourier pseudospectral method (PSM). In solving the ...