OptimizeDist: Find Extrema in D(a)

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

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)

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.

Value

extremalpha	Takeoff angle for identified extreme epicentral distance (s/deg)
extremp	Ray parameter for extremalpha (s/deg)
extremdelta	Identified extreme epicentral distance

Author(s)

Jake Anderson

Examples

### not a user routine

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