Okada: Surface deformation due to a finite rectangular source (Scripts) Publisher's description
from FranГ§ois Beauducel
The Okada  model allows to calculate analytic solution for surface deformation due to shear and tensile faults in an elastic half-space.
The Okada  model allows to calculate analytic solution for surface deformation due to shear and tensile faults in an elastic half-space. This model is widely used to simulate ground deformation produced by local perturbation like tectonic faults (earthquakes) or volcanic dykes. Given rectangular fault geometry (length, width, depth, strike, dip) and 3-component dislocation amplitude, it computes the displacements, tilt and horizontal strain on the free-surface.
The proposed Matlab script is a literal transcription of the Okada's equations, except that it proposes a strike angle of the fault, dislocation parameters are given by: rake, slip and opening (instead of U1, U2, U3), following Aki & Richards  definition, and coordinates are relative to fault centroid. LamГ©'s constants О» and Ој are also replaced by Poisson's ratio ОЅ, since the equations are independent of other elastic parameters. The equations are also vectorized for (x,y) coordinates. See help for syntax, and script comments for details.
System Requirements:MATLAB 6.5 (R13)
Program Release Status: New Release
Program Install Support: Install and Uninstall