SiStER: Simple Stokes solver with Exotic Rheologies
SiStER is a Matlab-based, easy-to-use geodynamic modeling code for 2-D problems involving visco-elasto-plastic flow of the lithosphere and mantle. SiStER relies on a finite-difference / marker-in-cell discretization of the (incompressible) Stokes flow equations, with some additional components to handle spontaneous brittle strain localization (faulting) and elastic deformation. It can be used for a wide range of problems, from tectonic rifting to the dynamics of accretionary prisms. The implementation largely follows the methodology of Gerya , and relies on Matlab's capabilities for rapid vector operations and sparse matrix filling. It is partially described in Olive et al.  (please reference this study if you plan on using SiStER for research!), as well as in the attached documentation. SiStER is hosted on a GitHub repository and can be downloaded here. In addition to research applications, SiStER can be a great teaching tool to illustrate the various components of a modern geodynamic modeling code. Its fairly transparent and modular structure makes it easy to modify and improve... so new contributions are always welcome !
SiStER was developed by J.-A. Olive, E. Mittelstaedt, B. Klein, G. Ito, M. Behn, and S. Howell. The development was partially supported by NSF grants OCE-11-54238 (M. Behn); OCE-11-55098 (G. Ito); and an LDEO post-doctoral fellowship (J.-A. Olive).
Publications using SiStER (and topic) as of 2019:
Ice cliff collapse: Clerc, F., Minchew, B. M., and Behn, M. D., Marine ice cliff instability mitigated by slow removal of ice shelves, Geophys. Res. Lett., 46, 12108-12116, 2019. Oceanic transform faults: Morrow, T. A., Mittelstaedt, E., Kim, S.-S., Are segmented fracture zones weak? Analytical and numerical models constrain anomalous bathymetry at the Clarion and Murray fracture zones, Earth Planet. Sci Lett., 512, 214-226, 2019. Tectonics of icy worlds: Howell, S. M, and Pappalardo, R. T., Band formation and ocean-surface interaction on Europa and Ganymede, Geophys. Res. Lett., 45, 4701-4709, 2018. Accretionary prisms: Weiss, J. R., Ito, G., Brooks, B. A., Olive, J.-A., Moore, G. F., and Foster, J. H., Formation of the frontal thrust zone of accretionary wedges, Earth Planet. Sci Lett., 495, 87-100, 2018. Patient zero: Olive, J.-A., Behn, M. D., Mittelstaedt, E., Ito, G., and B. Z. Klein, The role of elasticity in simulating long-term tectonic extension, Geophys. J. Int, 205(2), 2016.
eiSPLIT - A suite of MATLAB tools for shear wave splitting
eiSPLIT is a suite of Matlab scripts aimed at conducting shear wave splitting measurements following the eigenvalue-minimization method of Silver and Chan . The compressed files can be downloaded here, and a user guide with some background on the method can be found here.
Below are individual elements of the eiSPLIT package that may prove helpful to become familiar with the method. The script split.m applies the eigenvalue minimization method to the east- and north- components of a split shear wave in order to estimate the best splitting parameters (fast azimuth and delay time) of the medium that generated it. This script can be tested by running synthetic.m in the same folder. This will generate a synthetic SKS waveform, split it according to prescribed parameters, then unsplit the resulting signals using split.m to estimate said parameters. Note that synthetic.m requires ttshift.m to carry out time-shifting operations.
2-D gravity and magnetic anomalies
The function gpoly.m calculates the vertical gravity anomaly generated by an infinite prism of polygonal section. It is a Matlab rewriting of the FORTRAN routine by Blakely . You can run buried_cylinder.m in the same folder to benchmark it against the analytical solution for a buried (infinite) cylinder.
The function quadmag.m calculates the magnetic anomaly generated by an infinite prism of quadrilateral section, following Talwani and Heirtzler . You can run buried_square.m in the same folder to test it with a right-square prism of prescribed magnetization.