Calculate Residuals

Description

given an earthquake hypocenter and a list of station information, retrieve the station residuals.

Usage

1
EQXYresid(XY, vel = list(), h1 = c(0, 0, 0, 0), PLOT = FALSE)

Arguments

XY

matrix of station location and arrival times.

vel

list, RSEIS velocity model

h1

hypocenter location, c(x,y,z,t)

PLOT

logical, TRUE=plot the residuals

Details

The XY mtrix is in cartesian coordinates, i.e. it has been projected into units of km. Only 1D velocity models are used at this time. Only residuals of P and S wave arrivals are estimated.

Value

vector, right hand side of the least squares problem.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

travel.time1D,UPdateEQLOC

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
## Not run: 
 XY = GEOmap::GLOB.XY(pstas$lat, pstas$lon, proj)
    elcor = rep(0, length(pstas$lat))
    
    DZ = pstas$z - mean(stas$z)
    
    elcor[pstas$phase=="P"] = DZ[pstas$phase=="P"]/v$vp[1]
    elcor[pstas$phase=="S"] = DZ[pstas$phase=="S"]/v$vs[1]
    
    XY$cor = elcor
    XY$phase = pstas$phase
    XY$sec = pstas$sec
    
    eqXY = GEOmap::GLOB.XY(sol[1], sol[2], proj)
    res =  EQXYresid(XY, vel=v , h1=c(eqXY$x, eqXY$y, sol[2],sol[4]) , PLOT=FALSE)
    
## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.