5.7.1.4. 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:

float

Returns:

Magnitude at maximum curvature

Note

Should be used as a guide, often under-estimates Mc. Personally not fond of this method.

Example

>>> from obspy.clients.fdsn import Client
>>> from obspy import UTCDateTime
>>> from eqcorrscan.utils.mag_calc import calc_max_curv
>>> client = Client('IRIS')
>>> t1 = UTCDateTime('2012-03-26T00:00:00')
>>> t2 = t1 + (3 * 86400)
>>> catalog = client.get_events(starttime=t1, endtime=t2, minmagnitude=3)
>>> magnitudes = [event.magnitudes[0].mag for event in catalog]
>>> calc_max_curv(magnitudes, plotvar=False)
3.6000000000000001