Statici: Heterogeneous Crust Deformation
Elastostatic Deformation in Heterogeneous Media
The deformation due to faulting in the presence of arbitrary spatial variations in elastic properties can in general be evaluated with fully numerical techniques like finite elements or in some simple cases with semi-analytic techniques, like the propagator matrix. Several investigators have proposed to tackle the general case that includes gradual variations as well as lateral and vertical changes in elastic moduli. Here we propose an iterative method that uses equivalent body forces and equivalent surface tractions to represent the effect of inclusions. The formulation for 2-D antiplane deformation can be found in (Barbot et al. 2008) and approach for the general 3-D case in (Barbot et al. 2009).
Heterogeneous Half Space Deformation Code: STATICI (last release: 08/2011)
The source code of STATICI, used to evaluate static deformation due to fault in a heterogeneous half space with spatially-varying elastic properties can be downloaded here:
Some programs for post processing and visualization with gmt can be downloaded here:
Some compilation options relate to the choice of FFT (fourt, FFTW, SGI FFT, MKL), output formats (text, GMT, Paraview) and parallelization method (MPI or OpenMP). These options are explained in the readme file contained in the tar ball. See also some related comments about the homogeneous version STATICU. Instructions and run examples can be found in the readme file.
The latest release includes a binary output format compatible with the 3D visualizer Paraview, provides commented examples and more installation instructions.