Ray.time1D: 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

Value

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

Details

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

Examples

Run this code



data(VELMOD1D)

v <-  VELMOD1D
tees  <-  travel.time1D(23, 7, 0, length(v$zp) , v$zp , v$vp)
rays  <-  Ray.time1D(23, 7, 0, length(v$zp) , 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