TauP.R (version 1.0)

AnalyzeLVZ: Analyze Low Velocity Zones

Description

Identifies low velocity zones and improves sampling to allow more accurate raypath calculation.

Usage

AnalyzeLVZ(v, vsec, z, rp)

Arguments

v

Velocity vector (km/s)

vsec

Other velocity vector (km/s)

z

Depth vector (km)

rp

Planet radius

Value

• List with following elements:
• newv
  Velocity (of the same type as input v) vector at new depths (km/s)
• newvsec
  Velocity (of the same type as input vsec) vector at new depths (km/s)
• newz
  New depths sampled (km)
• criticalz
  Critical depths requiring special treatment (km)

Details

Only v is checked for LVZs. However, since a velocity profile requires both P and S velocities, the other velocity vector is provided as vsec and interpolated within LVZs found in v.

Interpolated velocities might not match those returned by InterpModel because calculations are done after a flat earth transform here.

Examples

data(model)
v = model$vp
vsec = model$vs
z = model$z
rp = model$rp

AnalyzeLVZ(v, vsec, z, rp)