TauP.R (version 1.0)

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

data(model)

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

DistSummary(delta = 100, h = 40, model = model)