Codes

I use a variety of codes, including a spectral boundary integral equation method (BIEM) code and a summation-by-parts finite difference (FD) code to model earthquake ruptures and volcanic eruptions. The BIEM code (version 4.1.7) includes input files that can be used to reproduce results in some of my papers on poroelastic bimaterial effects and ruptures with flash heating and thermal pressurization. The FD code has been used to study ruptures on nonplanar faults with strongly rate-weakening fault friction and off-fault plasticity. It currently only handles 2D problems.

MDSBI (MultiDimensional Spectral Boundary Integral code):
Language: Fortran 95 (with a few common extensions)
Use: Dynamic ruptures and interfacial elastodynamics
Features:
— Parallelized with MPI and OpenMP
— 2D (modes II and III and mixed mode) or 3D
— identical or dissimilar materials bounding the fault
— fully dynamic (displacement or velocity formulation) or quasidynamic elastic response
— automatic adaptive time stepping (using embedded Runge-Kutta error estimation method) for quasidynamic response
— adaptive substepping for arbitrary friction laws
— friction laws include slip weakening, rate and state, flash heating, and thermal pressurization (with Hiro Noda)
— poroelastic fault zone model

Versions:

Version 4.1.7, 20 September 2008 (MPI required, OpenMP optional)

Source code: mdsbi-v4.1.7.tgz (210 kB)
User Guide (with installation instructions): UserGuide-v4.1.7.pdf (75 kB)
Kernel tables: kernels.tgz (5.6 MB)


Version 3.9.10, 17 May 2007 (optional MPI/OpenMP support)

Source code: mdsbi-v3.9.10.tgz (213 kB)
User Guide (with installation instructions): UserGuide-v3.9.10.pdf (94 kB)
Kernel tables: kernels.tgz (5.6 MB)

FDMAP (Finite Difference code, uses coordinate transforms or MAPing to handle complex geometries):
Language: Fortran 95 (with a few common extensions)
Use: Dynamic ruptures and seismic wave propagation in complex geometries
Features:
— Parallelized with MPI
— summation-by-parts FD method with simultaneous-approximation-term method for provably stable enforcement of boundary conditions
— 2D only (antiplane shear or plane strain)
— viscoplastic or rate-independent plasticity
— friction laws include slip weakening, rate and state, flash heating

Versions:

Version 1.0, 10 January 2013 (MPI required)

Source code: fdmap-v1.0.tgz
User Guide: UserGuide-v1.0.pdf

homepage
Last updated: January 10, 2013