Rquake (version 2.5-1)

XYSETUP: Set up matrix for hypocenter inversion

Description

Set up matrix for hypocenter inversion

Usage

XYSETUP(STAS, init, vel)

Value

matrix

Arguments

STAS

station information from pickfile

init

initial event location

vel

list, velocity

Author

Jonathan M. Lees<[e-mail]>

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.

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  )