TauP.R (version 1.5)

OptimizeDist: Find Extrema in D(a)

Description

Engine routine that identifies local extrema in the D(a) (epicentral distance/takeoff angle) function.

Usage

OptimizeDist(alphalimit, deltalimit, phase, h, imodel)

Arguments

alphalimit

Angle interval (2-element vector, deg)

deltalimit

Epicentral distances of alphalimit

phase

Arrival phase (e.g. 'P', 'PKIKP')

h

Focal depth (km)

imodel

Improved planet model (from ImproveModel)

Value

extremalpha

Takeoff angle for identified extreme epicentral distance (s/deg)

extremp

Ray parameter for extremalpha (s/deg)

extremdelta

Identified extreme epicentral distance

Details

OptimizeDist assumes that D(a) has only one extremum over the interval, and is finite and defined everywhere. It uses a Golden Section Search algorithm to find the extremum.

Examples

# NOT RUN {
### not a user routine
# }