3D FD in spherical sections
FD3S is a finite-difference solver of the elastic wave equation in a spherical section. FD3S allows to model seismic wave attenuation as well as anisotropy with radial symmetry axis. The finite difference scheme is of fourth order in space and of second order in time. Arbitrarily complex Earth models can easily be included. This method has been developed by Tarje Nissen-Meyer, and extended by Andreas Fichtner, in collaboration with Heiner Igel.
Download
You may download this Fortran90/MPI package from the external pageSPICE software library pagescall_made.
References
Igel, H., Nissen-Meyer, T., Jahnke, G., 2002. Wave propagation in 3-D spherical sections: Effects of subduction zones. Phys. Earth Planet. Inter., 132, pp. 219-234.
Nissen-Meyer, T., 2001. Numerical Simulation of Seismic Wave Propagation through Subduction Zones, Diplom thesis, Ludwig-Maximilians-Universität Munich, Germany.