Seismic 1D Travel Time and raypath
Travel time and raypath from source to reciever in 1D local model.
Ray.time1D(indelta, inhpz, instaz, inlay, ztop, vel)
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 |
Uses local 1D velocity model, not appropriate for spherical earth.
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 |
Jonathan M. Lees<jonathan.lees.edu>
travel.time1D, Get1Dvel
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
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.