NLSlocate {Rquake}

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)

Arguments

GH

List, RSEIS

vel

velocity model

init

initial guess for event location

PLOT

logical, TRUE=plot

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.

Value

vector, new location

Note

At this stage there are no weighting mechanisms or code to eliminate data that has residuals that are too large.

Author(s)

Jonathan M. Lees<[e-mail]>

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 )