R: Seismic 1D Travel Time and raypath

Ray.time1D {RSEIS} | R Documentation

## 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

[Package RSEIS version 3.7-4 Index]