5.8.1.3. eqcorrscan.utils.mag_calc.calc_max_curv¶
-
eqcorrscan.utils.mag_calc.
calc_max_curv
(magnitudes, plotvar=False)[source]¶ Calculate the magnitude of completeness using the maximum curvature method.
Parameters: - magnitudes (list) – List of magnitudes from which to compute the maximum curvature which will give an estimate of the magnitude of completeness given the assumption of a power-law scaling.
- plotvar (bool) – Turn plotting on and off
Return type: Returns: Magnitude at maximum curvature
Note
Should be used as a guide, often under-estimates Mc.
Example
>>> import numpy as np >>> mags = [] >>> for mag in np.arange(2.5,3, 0.1): ... mags.extend([mag] * int(20000 - 10 * mag)) >>> for mag in np.arange(3,7, 0.1): ... mags.extend([mag] * int(10 ** (7 - 1 * mag))) >>> calc_max_curv(mags, plotvar=False) 3.0