Travel time from source to reciever in 1D local model.
travel.time1D(indelta, inhpz, instaz, inlay, ztop, vel)
many.time1D(indelta, inhpz, instaz, inlay, ztop, vel)
list:
derivative of t w.r.t. horizontal distance
derivative of t w.r.t. z, depth
incidence angle, degrees
travel time, s
distance in KM
depth of hypocenter, km
elevation of station
number of layers
vector, tops of layers
vector, velocities in layers
Jonathan M. Lees<jonathan.lees.edu>
Uses local 1D velocity model, not appropriate for spherical earth. The many.time1D version will take a vector of distances (indelta) and either one station elevation or a vector.
The station elevation should be referenced to the top of the velocity model, not necessarily sea level. Usually this is set to zero and a station correction is used to take into account the topographic and other geologic effects.
Ray.time1D, Get1Dvel
data(VELMOD1D)
v <- VELMOD1D
tees <- travel.time1D(23, 7, 0, length(v$zs) , v$zp , v$vp)
print(tees)
Run the code above in your browser using DataLab