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

list:
dtdrderivative of t w.r.t. horizontal distance
dtdzderivative of t w.r.t. z, depth
angleincidence angle, degrees
tttravel time, s
nnodnumber of nodes
znodnode depths, km
rnodnode 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