TauP.R (version 1.5)

DistSummary: Arrival Summary

Description

Determine arrival times and information for all major phases arriving at a certain epicentral distance, and plot ray trajectories.

Usage

DistSummary(delta, h, model, phaselist = 'default', prop = "vp",
image.col = heat.colors(500), n = 200, ...)

Arguments

delta

Epicentral distance (degrees)

h

Focal depth (km)

model

Planet model

phaselist

Either 'default' for all available phases, or a character vector including desired phases

prop

Property by which to scale planet image: one of 'vp', 'vs', 'rho'.

image.col

Vector of colors for image

n

Resolution of image (pixels per side)

…

Other parameters for Rayfan

Value

Returns a list with the following elements:

p

Ray parameter for each arrival

t

Travel time for each arrival

dist

Epicentral distance (should be approximately the input dist)

phase

Phase of each arrival

Details

This function is really just a wrapper for Rayfan to calculate arrivals for many phases at just one epicentral distance. Since each phase must be calculated separately, the use of the default phaselist will result in a long calculation time (minutes), and the plot will probably be crowded. It is generally better to define phaselist as a smaller vector or use Rayfan instead.

See Also

Rayfan, Traveltime, Earthplot

Examples

# NOT RUN {
data(model)

# for an event occurring 100 degrees away at a depth of 40 km:

DistSummary(delta = 100, h = 40, model = model, phaselist = c('SKS', 'SKKS'))
# }