4.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: Return type: Returns: Magnitude at maximum curvature
Note
Should be used as a guide, often under-estimates Mc.
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.1000000000000001