Ray.time1D {RSEIS}

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