Computational Infrastructure
for Geodynamics (CIG)

Community-driven organization advancing Earth science by providing the infrastructure for the development and dissemination of software for geophysics and related fields.

2023 CIG Elections

Nominations are now being accepted for seats on the Executive and Science Steering Committees.  

Please contact the Nominations Committee chair with your suggestions by September 30th.

Elections open November 1.

Research Highlight

Evolution of the bound water content (top) and the free water (bottom) in a geodynamic one dimensional pipe flow model. Initially, all water is bound (left column). Since the middle layer has a zero solubility, water is almost instantaneously released from minerals and starts migrating upward with respect to the solid, leading to a higher bound water content when it reaches the top layer. Chemical reactions like these are typically much faster than the flow velocity in geodynamic models and are amongst the problems we want to solve by including dedicated ODE solver libraries.

Revisiting Open Source Libraries for Solving ODEs

Contributed by Wolfgang Bangerth, Colorado State University Fort Collins; Juliane Dannberg, University of Florida; Rene Gassmoeller, University of Florida

Many of CIG’s large flagship codes – for example Rayleigh, PyLith, SPECFEM, and ASPECT – solve partial differential equations. There is of course a long history of writing such codes and part of this history is that traditionally, PDE solvers have hand-rolled their own time integrators and nonlinear solvers: Exceptions – such as PyLith – notwithstanding, they generally use low-order time discretization methods (such as explicit or implicit Euler, Crank-Nicolson, or BDF-2) with time-step sizes based on CFL numbers rather than accuracy considerations; and they implement nonlinear solvers via fixed point iterations or relatively simple variations of Newton’s method, but without sophisticated line search, trust region, or acceleration methods.

This reliance on hand-written methods is perhaps surprising because we have all learned that at least for the discretization of PDEs, we should build on one of the widely used software libraries that provide everything one needs for the task: ... [full article]

Check out all Research Highlights

What Is happening?


CIG Business meeting





Posted 31 July 2023

2023-2024 Webinar Series

Fall 2023

Check back in late summer for our 2023-24 webinar schedule.

Posted 31 May 2023






Geochemistry, Geophysics, Geosystems Special Collection

Call for Papers for Frontiers in lithospheric dynamics: bridging scales through observations, experiments, and computations.   •  More Info

Deadline: 23 November 2023

Posted 17 January 2022

Get the Tools


A set of codes for MHD dynamo simulation in a rotating spherical shell using spherical harmonics expansion methods.
Current release: 2017-07-17
GNU GPL v2 or newer license


Finite-element code for dynamic and quasistatic simulations of crustal deformation, primarily earthquakes. and volcanoes
Current release: 2022-10-14
MIT license


A 3-D convection code designed for the study of dynamo behavior in spherical shell geometry. 
Current release: 2022-05-05
GNU GPL v3 or newer license



3-D seismic modeling, with a free surface condition on the top boundary, absorbing super-grid conditions on the far-field boundaries, and an arbitrary number of point force and/or point moment tensor source terms. 
Current release: 2023-08-30
GNU GPL v2 or newer license

Contribute Code

Find out more about contributing software to the community.
CIG on Zenodo
Join our forum to receive announcements and be part of community discussions. The forum has replaced our mailing lists.