TauP.R (version 1.5)

Rayfan: Ray Fan

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?

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)

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.

See Also

Earthplot, Traveltime, DistSummary

Examples

# NOT RUN {
data(model)

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