GETpsTT {Rquake} | R Documentation |
Get Pand S travel times and derivatives
Description
Get Pand S travel times and derivatives
Usage
GETpsTT(phase, eqz = 6, staz = 0, delx = 1, dely = 1, deltadis = 6, vel)
Arguments
phase |
character vector, phase |
eqz |
event depth |
staz |
station elevation |
delx |
km, delta X |
dely |
km, delta Y |
deltadis |
km, distance |
vel |
velocity models (P and S) |
Details
Creates a vector of travel times, and a matrix and derivatives used for inversion.
Value
list:
TT |
travel time vector |
Derivs |
matrix of derivatives, dtdx, dtdy, dtdz |
Author(s)
Jonathan M. Lees<jonathan.lees@unc.edu>
See Also
many.time1D
Examples
library(RSEIS)
library(GEOmap)
data(GH, package='RSEIS')
data(VELMOD1D, package='RSEIS')
vel = VELMOD1D
p1 = GH$pickfile$STAS
loc = GH$pickfile$LOC
proj = GEOmap::setPROJ(type = 2, LAT0 =loc$lat, LON0 = loc$lon)
XYsta = GEOmap::GLOB.XY(p1$lat, p1$lon, proj)
XYq = GEOmap::GLOB.XY(loc$lat, loc$lon, proj)
delx = XYq$x-XYsta$x
dely = XYq$y-XYsta$y
dists = sqrt(delx^2+dely^2)
G1 = GETpsTT(p1$phase, eqz=loc$z, staz=0, delx=delx, dely=dely, deltadis=dists , vel)
[Package Rquake version 2.5-1 Index]