XYSETUP: Set up matrix for hypocenter inversion

In Rquake: Seismic Hypocenter Determination

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  )

Rquake documentation built on June 22, 2024, 11:54 a.m.