Coseis is a toolkit for earthquake simulation featuring:
The Support Operator Rupture Dynamics (SORD) code for modeling spontaneous rupture and 3D wave propagation.
Utilities for mesh generation, coordinate projection, and visualization.
The primary interface is through a Python module which (for high-performance components) wraps Fortran parallelized with hybrid OpenMP and MPI.
Coseis is written by Geoffrey Ely with contributions from Steven Day, Bernard Minster, Feng Wang, Zheqiang Shi, and Jun Zhou. It is licensed under BSD terms.
WARNING: Coseis is a research code under active development. Changes are frequent and it has known bugs! The latest committed version may not be your best option or even working. Please contact me for guidance.
For MacOS only, install Xcode from the App Store followed by the Xcode the Command Line Tools with:
For MacOS only, install Homebrew and use it to install Fortran with:
brew install gfortran
Clone the source code from the Coseis GitHub repository:
git clone git://github.com/gely/coseis.git
Setup python to be able to find the
cd coseis python -m cst setup
To run the test suite interactively:
cd tests python test_runner.py --run=exec
Or, submit a job for batch processing:
python test_runner.py --run=submit
After completion, a report is printed to the screen (or saved in
PASSED: cst.tests.hello_mpi.test() PASSED: cst.tests.point_source.test() PASSED: cst.tests.pml_boundary.test() PASSED: cst.tests.kostrov.test()
Extract S-wave velocity at 500 meters depth. Plot using Matplotlib:
Build a fence diagram similar to Magistrale (2000) figure 10. Plot using Mayavi
Extract 3D mesh and search for the shallowest surface of Vs = 2.5 km/s. Plot over topography using Mayavi.