Eigenfunction¶

class disba.RayleighEigen(depth, ur, uz, tz, tr, period, mode)¶

Create new instance of RayleighEigen(depth, ur, uz, tz, tr, period, mode)

class disba.LoveEigen(depth, uu, tt, period, mode)¶

Create new instance of LoveEigen(depth, uu, tt, period, mode)

EigenFunction¶

class disba.EigenFunction(thickness, velocity_p, velocity_s, density, algorithm='dunkin', dc=0.005)[source]¶

Eigenfunction class.

Parameters

thickness (array_like) – Layer thickness (in km).
velocity_p (array_like) – Layer P-wave velocity (in km/s).
velocity_s (array_like) – Layer S-wave velocity (in km/s).
density (array_like) – Layer density (in g/cm3).
algorithm (str {'dunkin', 'fast-delta'}, optional, default 'dunkin') –
Algorithm to use for computation of Rayleigh-wave dispersion:
- ’dunkin’: Dunkin’s matrix (adapted from surf96),
- ’fast-delta’: fast delta matrix (after Buchen and Ben-Hador, 1996).
dc (scalar, optional, default 0.005) – Phase velocity increment for root finding.

__call__(t, mode=0, wave='rayleigh')[source]¶

Compute eigenfunctions for a given period and mode.

Parameters

t (scalar) – Period (in s).
mode (int, optional, default 0) – Mode number (0 if fundamental).
wave (str {'love', 'rayleigh'}, optional, default 'rayleigh') – Wave type.

Returns

Eigenfunction as a namedtuple:

If wave == 'love', (depth, uu, tt, period, mode),

If wave == 'rayleigh', (depth, ur, uz, tz, tr, period, mode).

Return type

disba.RayleighEigen or disba.LoveEigen

property algorithm¶: Return algorithm to use for computation of Rayleigh-wave dispersion.

resample(dz)¶

Resample velocity model.

Dispersion

Ellipticity