Rquake (version 2.5-1)

NLSlocate: Nonlinear Least Squares Location

Description

Nonlinear Least Squares Location using Gieger's method

Usage

NLSlocate(GH, vel = list(), init = c(0, 0, 0, 0), PLOT = FALSE)

Value

vector, new location

Arguments

GH

List, RSEIS

vel

velocity model

init

initial guess for event location

PLOT

logical, TRUE=plot

Author

Jonathan M. Lees<[e-mail]>

Details

This is an adaptation of non-linear least squares inversion for earthquake location. A residual function is supplied, and iterations are performed until the location is determined.

References

Lee, W.H.K., and S.W. Stewart, Principles and Applications of Microearthquake Networks, Academic Press, New York, 1981.

See Also

swig

Examples

data(GH, package='RSEIS')
###  location is:
eqsol = NLSlocate(GH, vel=GH$velfile,  PLOT=TRUE )