You are here: Home / Resources / Software: Download / SPECFEM3D GLOBE / About


By Dimitri Komatitsch (primary-developer)1, ean-Pierre Vilotte (primary-developer)2, Jeroen Tromp (primary-developer)3, SPECFEM Development Team

1. CNRS, France 2. Institut de Physique du Globe 3. Princeton University, USA

Published on


SPECFEM3D_GLOBE simulates global and regional (continental-scale) seismic wave propagation.

Effects due to lateral variations in compressional-wave speed, shear-wave speed, density, a 3D crustal model, ellipticity, topography and bathymetry, the oceans, rotation, and self-gravitation are all included.

The version 7.0 release offers GPU graphics card support for both OpenCL and CUDA hardware accelerators, based on an automatic source-to-source transformation library (Videau et al. 2013). It offers additional support for ADIOS file I/O formats and contains important bug fixes related to 3D topography and geographic/geocentric transformations. Seismogram file names adapt a new naming convention, with better compatibility to the seismogram specifications by the Incorporated Research Institutions for Seismology (IRIS).

The version embeds non-blocking MPI communications and includes several performance improvements in mesher and solver. It provides a perfectly load-balanced mesh for 3D mantle models honoring shallow oceanic Moho (depths less than 15 km) and deep continental Moho (depths greater than 35 km). It also accommodates European crustal models EPcrust (Molinari & Morelli, 2011) and EuCrust07 (Tesauro et al., 2008), which may be combined with global crustal model Crust2.0. Sedimentary wavespeeds are superimposed on the mesh if sediment thickness exceeds 2 km.

Additional new model routines are provided for the Comprehensive Earth Model (CEM) project, generic point-profile models (PPM) and Gauss-Lobatto-Legendre based models (GLL), with complementary tools for postprocessing adjoint sensitivity kernels and gradient-based model updates. The structure of the software has been simplified to facilitate easier implementation of new 3D models. The code accommodates general moment tensor files, and provides complete information in the SAC headers, as explained in detail in the updated user manual. New matrix-matrix multiplication routines, adapted from the book of Deville et al. (2002), and loop-vectorization help reduce the total number of memory accesses performed in each spectral element and improve code vectorization, thus enhance numerical performance of the version.

Release Notes

For latest the version, download the software from the repository.

Sponsored by

This material is based upon work supported by the U.S. National Science Foundation under Grants No. EAR-0406751 and EAR-0711177, by the French CNRS, French INRIA Sud-Ouest MAGIQUE-3D, French ANR NUMASIS under Grant No. ANR-05-CIGC-002, and European FP6 Marie Curie International Reintegration Grant No. MIRG-CT-2005-017461. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the U.S. National Science Foundation, CNRS, INRIA, ANR or the European Marie Curie program.


The cover graphic of the manual was created by Santiago Lombeyda from Caltech’s Center for Advanced Computing Research (CACR), USA, with free satellite clipart pictures from and http: // added to it.