Software

Short-Term Crustal Dynamics

PyLith

Developed, maintained, and supported.

Finite element code for the solution of dynamic and quasi-static tectonic deformation problems.

RELAX

Developed, maintained, and supported.

Relax implements a semi-analytic Fourier-domain solver and equivalent body forces to compute quasi-static relaxation of stress perturbation

SELEN

Accepting bug reports.

SELEN solves numerically the so-called "Sea Level Equation" (SLE) for a spherical, layered, non-rotating Earth with Maxwell viscoelastic rheology.

LithoMop

Accepting bug reports.

Finite element code for the solution of visco-elastic/plastic deformation that was designed for lithospheric modeling problems. A newer version of LithoMop, called PyLith, has been released by CIG.

Long-Term Tectonics

Gale

Accepting bug reports.

2D/3D code for the long-term tectonics community. The code solves problems related to orogenesis, rifting, and subduction with coupling to surface erosion models.

Plasti

Accepting bug reports.

2D ALE (Arbitrary Lagrangian Eulerian) code donated to CIG by Sean Willett and Chris Fuller of the University of Washington. The code originated at Dalhousie University in Canada.

SNAC

Accepting bug reports.

Updated Lagrangian explicit finite difference code for modeling a finitely deforming elasto-visco-plastic solid in 3D.

Mantle Convection

Aspect

Developed, maintained, and supported.

Finite element parallel code to simulate problems in thermal convection in both 2D and 3D models - currently in alpha testing.

CitcomCU

Developed, maintained, and supported.

Finite element parallel code capable of modeling thermochemical convection in a three dimensional domain appropriate for convection within the Earth's mantle.

CitcomS

Developed, maintained, and supported.

Finite element code designed to solve compressible thermochemical convection problems relevant to Earth's mantle.

ConMan

Maintained and supported.

Finite element program for the solution of the equations of incompressible, infinite-Prandtl number convection in two dimensions, originally written by Scott King, Arthur Raefsky, and Brad Hager.

Ellipsis3d

Accepting bug reports.

Three-dimensional version of the particle-in-cell finite element code Ellipsis, a solid modeling code for visco-elastoplastic materials. The particle-in-cell method combines the strengths of the Lagrangian and Eulerian formulations of mechanics while bypassing their limitations.

HC

Accepting bug reports.

Global mantle circulation solver following Hager & O'Connell (1981) which can compute velocities, tractions, and geoid for simple density distributions and plate velocities.

Seismology

Specfem3D

Developed, maintained, and supported.

Simulates seismic wave propagation in sedimentary basins. (formerly SPECFEM3D_BASIN)

Specfem3D Globe

Developed, maintained, and supported.

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

Specfem3D Geotech

Developed, maintained, and supported.

3D slope stability analysis and simulation of 3D multistage excavation based on the spectral-element method.

Specfem2D

Accepting bug reports.

Simulates forward and adjoint seismic wave propagation in two-dimensional acoustic, (an)elastic, poroelastic or coupled acoustic-(an)elastic-poroelastic media.

Specfem1D

Accepting bug reports.

Simulates seismic wave propagation in a one-dimensional heterogeneous medium. It is a small code that allows users to learn how a spectral-element program is written.

Mineos

Accepting bug reports.

1D code used to simulate synthetic seismograms in the spherical symmetric nonrotated Earth by normal mode summation.

Flexwin

Accepting bug reports.

Automates the time-window selection problem for seismologists. It operates on pairs of observed and synthetic single component seismograms, defining windows that cover as much of a given seismogram as possible, while avoiding portions of the waveform that are dominated by noise.

Seismic CPML

Accepting bug reports.

Eight open-source Fortran90 programs to solve the two-dimensional or three-dimensional isotropic or anisotropic elastic, viscoelastic or poroelastic wave equation using a finite-difference method with Convolutional Perfectly Matched Layer (C-PML) conditions, developed by Dimitri Komatitsch and Roland Martin from University of Pau, France.

Finite-Frequency Tomography Software

For the programs documented in Nolet, G., A Breviary of Seismic Tomography (CUP, 2008).

Geodynamo

MAG

Accepting bug reports.

Serial version of a rotating spherical convection/magnetoconvection/dynamo code, developed by Gary Glatzmaier and modified by Uli Christensen and Peter Olson.

Computational Science

Geodynamics AMR Suite (deal.II)

Developed.

C++ program library targeted at the computational solution of partial differential equations using adaptive finite elements. Its state-of-the-art programming techniques offer a modern interface to the complex data structures and algorithms required. CIG has sponsored Wolfgang Bangerth (Texas A&M) to produce tutorials that use deal.II in various implementations.

Cigma

Accepting bug reports.

Suite of tools that facilitates the comparison of numerical models, and performs error analysis, benchmarking, and code verification.

Exchanger

Accepting bug reports.

Package containing several C++ base classes. These classes, when customized for a solver, can provide communication channels between solvers. This packaged is used by CitcomS for solver coupling.

Pythia/Pyre

Accepting bug reports.

Pyre framework and a collection of packages that interact with it, such as an interface to the ACIS solid modeling package.

Donating Code

Interested in donating code to CIG? Take a look here.