Ray.time1D: Seismic 1D Travel Time and raypath

In RSEIS: Seismic Time Series Analysis Tools

Seismic 1D Travel Time and raypath

Description

Travel time and raypath from source to reciever in 1D local model.

Usage

Ray.time1D(indelta, inhpz, instaz, inlay, ztop, vel)

Arguments

indelta	distance in KM
inhpz	depth of hypocenter, km
instaz	elevation of station
inlay	number of layers
ztop	vector, tops of layers
vel	vector, velocities in layers

Details

Uses local 1D velocity model, not appropriate for spherical earth.

Value

list:

dtdr	derivative of t w.r.t. horizontal distance
dtdz	derivative of t w.r.t. z, depth
angle	incidence angle, degrees
tt	travel time, s
nnod	number of nodes
znod	node depths, km
rnod	node offset distances, km

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

See Also

travel.time1D, Get1Dvel

Examples

data(VELMOD1D)

v <-  VELMOD1D
indelta=23.;
 inhpz=7.;
 instaz=0.;
nz = length(v$zp)

tees  <-  travel.time1D(indelta, inhpz, instaz, nz , v$zp , v$vp)
rays  <-  Ray.time1D(indelta, inhpz, instaz, nz , v$zp , v$vp)



plot(rays$rnod[1:rays$nnod] , -rays$znod[1:rays$nnod],type="n",
 xlab="distance, km" , ylab="Depth, km")

abline(h=-v$zp, lty=2, col=grey(0.80) )
points(rays$rnod[1:rays$nnod] , -rays$znod[1:rays$nnod], pch=8, col='green')
lines(rays$rnod[1:rays$nnod] , -rays$znod[1:rays$nnod])
points(rays$rnod[rays$nnod] , -rays$znod[rays$nnod], pch=6, col='red', cex=2)
#####  to coordinate this in space, need to rotate about
#####    the line between source and receiver locations

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.