Rayfan: Ray Fan

In TauP.R: Earthquake Traveltime Calculations for 1-D Earth Models

Description

Calculate travel times and plot ray trajectories of phase(s) from focus to receiver(s).

Usage

Rayfan(phaselist, h, model, deltalist = 1:17 * 20, minp = 0.5, plot =
TRUE, add = TRUE, col = rep("black", length(phaselist)), verbose =
FALSE, mirror = FALSE)

Arguments

phaselist	Character vector of phases to plot (e.g. c('P','S','PP','SS') )
h	Focal depth (km)
model	Planet model
deltalist	Vector of epicentral distances (degrees)
minp	Smallest allowed ray parameter (s/deg) to prevent errors near the center of the planet.
plot	Logical: plot ray trajectories?
add	Add to existing Earthplot/Rayfan figure?
col	Color vector for 'image'
verbose	Print information as calculations are done?
mirror	Logical: should delta = x be considered equivalent to delta = 360 - x?

Details

It is useful to remember phases like PKKP that travel more than 180 degrees may physically arrive in the same place as a phase that travels less than 180 degrees like PKP, but this package does not recognize it unless 'mirror' is TRUE.

Value

Output from each Traveltime calculation is concatenated into the following list:

tt	vector of traveltimes (s)
p	vector of ray parameters (s/deg)
angles	vector of takeoff angles (degrees)
dists	vector of epicentral distances (degrees)

Author(s)

Jake Anderson

See Also

Earthplot, Traveltime, DistSummary

Examples

data(model)

Rayfan(c('S','ScS'),100,model)

TauP.R documentation built on May 2, 2019, 3:25 a.m.