TauP.R (version 1.5)

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

# NOT RUN {
data(model)
v = model$vp
vsec = model$vs
z = model$z
rp = model$rp

AnalyzeLVZ(v, vsec, z, rp)
# }