XYSETUP {Rquake}

Set up matrix for hypocenter inversion

Description

Set up matrix for hypocenter inversion

Usage

XYSETUP(STAS, init, vel)

Arguments

STAS

station information from pickfile

init

initial event location

vel

list, velocity

Details

This sets up the matrix used for nonlinear inversion. The code does not include information on the weighting. Station corrections are included.

The STAS are an internal component of the pickfile.

Value

matrix

Note

Need scheme for weighting according to errors in picks and distance weighting.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

setPROJ, GLOB.XY,NLSlocate

Examples

##  start with the location of the closest station
data(GH, package='RSEIS')

g1 = GH$pickfile
data(VELMOD1D, package='RSEIS')

vel= VELMOD1D

STAS = GH$pickfile$STAS
w1 = STAS$phase == 'P'
initz = 6
t0a = GH$pickfile$LOC$sec


XY = XYSETUP(STAS, c(STAS$lat[w1],STAS$lon[w1], initz,  STAS$sec[w1]-t0a  ) , vel  )