Learn R Programming
RSEIS (version 2.6-0)

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

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.

See Also

travel.time1D, Get1Dvel

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

Run the code above in your browser using DataLab