Learn R Programming
TauP.R (version 1.0)

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

  • extremalphaTakeoff angle for identified extreme epicentral distance (s/deg)
  • extrempRay parameter for extremalpha (s/deg)
  • extremdeltaIdentified 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.