CPS¶
surf96¶
- disba.surf96(t, d, a, b, rho, mode=0, itype=0, ifunc=2, dc=0.005, dt=0.025)[source]¶
Get phase or group velocity dispersion curve.
- Parameters
t (array_like) – Periods (in s).
d (array_like) – Layer thickness (in km).
a (array_like) – Layer P-wave velocity (in km/s).
b (array_like) – Layer S-wave velocity (in km/s).
rho (array_like) – Layer density (in g/cm3).
mode (int, optional, default 0) – Mode number (0 if fundamental).
itype (int, optional, default 0) –
- Velocity type:
0: phase velocity,
1: group velocity.
ifunc (int, optional, default 2) –
- Select wave type and algorithm for period equation:
1: Love-wave (Thomson-Haskell method),
2: Rayleigh-wave (Dunkin’s matrix),
3: Rayleigh-wave (fast delta matrix).
dc (scalar, optional, default 0.005) – Phase velocity increment for root finding.
dt (scalar, optional, default 0.025) – Frequency increment (%) for calculating group velocity.
- Returns
Phase or group dispersion velocity.
- Return type
array_like
swegn96¶
- disba.swegn96(t, d, a, b, rho, mode, ifunc=2, dc=0.005)[source]¶
Get eigenfunctions for a given period and mode.
- Parameters
t (scalar) – Period (in s).
d (array_like) – Layer thickness (in km).
a (array_like) – Layer P-wave velocity (in km/s).
b (array_like) – Layer S-wave velocity (in km/s).
rho (array_like) – Layer density (in g/cm3).
mode (int, optional, default 0) – Mode number (0 if fundamental).
ifunc (int, optional, default 2) –
- Select wave type and algorithm for period equation:
1: Love-wave (Thomson-Haskell method),
2: Rayleigh-wave (Dunkin’s matrix),
3: Rayleigh-wave (fast delta matrix).
dc (scalar, optional, default 0.005) – Phase velocity increment for root finding.
- Returns
Eigenfunctions.
- Return type
array_like
srfker96¶
- disba.srfker96(t, d, a, b, rho, mode=0, itype=0, ifunc=2, ipar=2, dc=0.005, dt=0.025, dp=0.025)[source]¶
Get phase or group velocity sensitivity kernel.
- Parameters
t (scalar) – Period (in s).
d (array_like) – Layer thickness (in km).
a (array_like) – Layer P-wave velocity (in km/s).
b (array_like) – Layer S-wave velocity (in km/s).
rho (array_like) – Layer density (in g/cm3).
mode (int, optional, default 0) – Mode number (0 if fundamental).
itype (int, optional, default 0) –
- Velocity type:
0: phase velocity,
1: group velocity.
ifunc (int, optional, default 2) –
- Select wave type and algorithm for period equation:
1: Love-wave (Thomson-Haskell method),
2: Rayleigh-wave (Dunkin’s matrix),
3: Rayleigh-wave (fast delta matrix).
ipar (int, optional, default 2) –
- Select parameter with respect to which sensitivity kernel is calculated:
0: Layer thickness,
1: P-wave velocity,
2: S-wave velocity,
3: Layer density.
dc (scalar, optional, default 0.005) – Phase velocity increment for root finding.
dt (scalar, optional, default 0.025) – Frequency increment (%) for calculating group velocity.
dp (scalar, optional, default 0.025) – Parameter increment (%) for numerical partial derivatives.
- Returns
Phase or group velocity sensitivity kernel.
- Return type
array_like