SPICE logo


SHaxi performs elastic and anelastic global SH wave propagation in the Earth's mantle for axi-symmetric geometries.

Authors: Gunnar Jahnke, Mike Thorne, Heiner Igel
General: Wave propagation
Code level: Research
Language: Fortran 90/95 or later
Associated groups: Numerical Methods
Supercomputing: Distributed Memory (e.g., MPI)
Grid dimensions: 2.5D
Coordinate system: Spherical
Method: Finite differences
Rheology: Acoustic scalar (SH)
Solution type: Numerical
Total hits: 8464

Project Description

The SHaxi program generates global high resolution SH wave seismograms and wavefield snapshots using high performance parallel computers or PC networks.

The algorithm is based on finite difference operators and staggered grids. For use on high performance parallel computers with distributed memory architecture or on PC clusters a domain decomposition of the model space was implemented, using Fortran 90 and MPI.

Possible models are restricted to axi-symmetric geometries which allows the computation of synthetic seismograms at dominant periods down to 2.5 seconds for global mantle models.

The code is verified against an analytical solution and the implemented ring-source is well understood. Appropriate routines for convenient definition of user-defined models are implemented and explained.

The SHaxi method has already proven to be a valuable method for computing global synthetic seismograms at high frequencies and for studying the seismic waveform effects from models where rotational symmetry may be assumed.

The code can be downloaded from: svn.geophysik.uni-muenchen.de/trac/shaxi

Detailed information abut the method can be found in Jahnke et al. (2008) as well as in the PhD thesis of Gunnar Jahnke.