TauP.R (version 1.0)

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

data(model)

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