FAILURE MAP
← Case archive

FA-71661 / Seismic magnitude estimation / Open access

Estimated Mercalli intensity: lower intensity clamp · case 01

Distant small events show negative intensities and wrap around the numeral table.

Verified by executionVariant 1 · 8 checks per implementationDownload source bundle ↓JSON ↗

ROOT CAUSE

Intensity is not clamped at I.

VERIFIED REPAIR

Clamp to 1..12.

Unsuccessful approach: Clamping at 0 produces intensity 0 and the XII label.

Case contract

R = max(hypocentral distance, 1 km). I = 1.45M - 2.46 log10(R) + 1.5, rounded half up to an integer and clamped to 1..12. Return [I, Roman numeral].

Why this case matters

ShakeMap-style intensity summaries are published as Roman numerals.

1 / The failure

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(mag, epi_km, depth_km):
    roman = ["I", "II", "III", "IV", "V", "VI", "VII", "VIII", "IX", "X", "XI", "XII"]
    r = max(math.hypot(epi_km, depth_km), 1.0)
    i = 1.45 * mag - 2.46 * math.log10(r) + 1.5
    v = int(math.floor(i + 0.5))
    v = min(v, 12)
    return [v, roman[v - 1]]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['M2.0 at 0 km depth 0.5', [2.0, 0, 0.5], [4, 'IV']], ['M2.0 at 0 km depth 5', [2.0, 0, 5], [3, 'III']], ['M2.0 at 3 km depth 2', [2.0, 3, 2], [3, 'III']], ['M2.0 at 20 km depth 10', [2.0, 20, 10], [1, 'I']], ['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']], ['M4.5 at 0 km depth 5', [4.5, 0, 5], [6, 'VI']], ['M4.5 at 3 km depth 2', [4.5, 3, 2], [7, 'VII']], ['M4.5 at 20 km depth 10', [4.5, 20, 10], [5, 'V']], ['M4.5 at 80 km depth 15', [4.5, 80, 15], [3, 'III']], ['M4.5 at 300 km depth 30', [4.5, 300, 30], [2, 'II']], ['M4.5 at 900 km depth 10', [4.5, 900, 10], [1, 'I']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M6.0 at 0 km depth 0.5', [6.0, 0, 0.5], [10, 'X']], ['M6.0 at 0 km depth 5', [6.0, 0, 5], [8, 'VIII']], ['M6.0 at 3 km depth 2', [6.0, 3, 2], [9, 'IX']], ['M6.0 at 20 km depth 10', [6.0, 20, 10], [7, 'VII']], ['M6.0 at 80 km depth 15', [6.0, 80, 15], [5, 'V']], ['M6.0 at 300 km depth 30', [6.0, 300, 30], [4, 'IV']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M7.2 at 0 km depth 0.5', [7.2, 0, 0.5], [12, 'XII']], ['M7.2 at 0 km depth 5', [7.2, 0, 5], [10, 'X']], ['M7.2 at 3 km depth 2', [7.2, 3, 2], [11, 'XI']], ['M7.2 at 20 km depth 10', [7.2, 20, 10], [9, 'IX']], ['M7.2 at 80 km depth 15', [7.2, 80, 15], [7, 'VII']], ['M7.2 at 300 km depth 30', [7.2, 300, 30], [6, 'VI']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M8.5 at 0 km depth 0.5', [8.5, 0, 0.5], [12, 'XII']], ['M8.5 at 0 km depth 5', [8.5, 0, 5], [12, 'XII']], ['M8.5 at 3 km depth 2', [8.5, 3, 2], [12, 'XII']], ['M8.5 at 20 km depth 10', [8.5, 20, 10], [11, 'XI']], ['M8.5 at 80 km depth 15', [8.5, 80, 15], [9, 'IX']], ['M8.5 at 300 km depth 30', [8.5, 300, 30], [8, 'VIII']], ['M8.5 at 900 km depth 10', [8.5, 900, 10], [7, 'VII']]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
M2.0 at 0 km depth 0.5[4, 'IV'][4, 'IV']Passed
M2.0 at 0 km depth 5[3, 'III'][3, 'III']Passed
M2.0 at 3 km depth 2[3, 'III'][3, 'III']Passed
M2.0 at 20 km depth 10[1, 'I'][1, 'I']Passed
M2.0 at 80 km depth 15[0, 'XII'][1, 'I']Failed
M2.0 at 300 km depth 30[-2, 'X'][1, 'I']Failed
M2.0 at 900 km depth 10[-3, 'IX'][1, 'I']Failed
M4.5 at 0 km depth 0.5[8, 'VIII'][8, 'VIII']Passed

SHA-256 / 6eaba1194635bff7e9b61d4425161f0e1272b7b5d7e6b3ad3026b0f894d66a80

2 / The unsuccessful fix

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(mag, epi_km, depth_km):
    roman = ["I", "II", "III", "IV", "V", "VI", "VII", "VIII", "IX", "X", "XI", "XII"]
    r = max(math.hypot(epi_km, depth_km), 1.0)
    i = 1.45 * mag - 2.46 * math.log10(r) + 1.5
    v = int(math.floor(i + 0.5))
    v = min(max(v, 0), 12)
    return [v, roman[v - 1]]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['M2.0 at 0 km depth 0.5', [2.0, 0, 0.5], [4, 'IV']], ['M2.0 at 0 km depth 5', [2.0, 0, 5], [3, 'III']], ['M2.0 at 3 km depth 2', [2.0, 3, 2], [3, 'III']], ['M2.0 at 20 km depth 10', [2.0, 20, 10], [1, 'I']], ['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']], ['M4.5 at 0 km depth 5', [4.5, 0, 5], [6, 'VI']], ['M4.5 at 3 km depth 2', [4.5, 3, 2], [7, 'VII']], ['M4.5 at 20 km depth 10', [4.5, 20, 10], [5, 'V']], ['M4.5 at 80 km depth 15', [4.5, 80, 15], [3, 'III']], ['M4.5 at 300 km depth 30', [4.5, 300, 30], [2, 'II']], ['M4.5 at 900 km depth 10', [4.5, 900, 10], [1, 'I']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M6.0 at 0 km depth 0.5', [6.0, 0, 0.5], [10, 'X']], ['M6.0 at 0 km depth 5', [6.0, 0, 5], [8, 'VIII']], ['M6.0 at 3 km depth 2', [6.0, 3, 2], [9, 'IX']], ['M6.0 at 20 km depth 10', [6.0, 20, 10], [7, 'VII']], ['M6.0 at 80 km depth 15', [6.0, 80, 15], [5, 'V']], ['M6.0 at 300 km depth 30', [6.0, 300, 30], [4, 'IV']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M7.2 at 0 km depth 0.5', [7.2, 0, 0.5], [12, 'XII']], ['M7.2 at 0 km depth 5', [7.2, 0, 5], [10, 'X']], ['M7.2 at 3 km depth 2', [7.2, 3, 2], [11, 'XI']], ['M7.2 at 20 km depth 10', [7.2, 20, 10], [9, 'IX']], ['M7.2 at 80 km depth 15', [7.2, 80, 15], [7, 'VII']], ['M7.2 at 300 km depth 30', [7.2, 300, 30], [6, 'VI']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M8.5 at 0 km depth 0.5', [8.5, 0, 0.5], [12, 'XII']], ['M8.5 at 0 km depth 5', [8.5, 0, 5], [12, 'XII']], ['M8.5 at 3 km depth 2', [8.5, 3, 2], [12, 'XII']], ['M8.5 at 20 km depth 10', [8.5, 20, 10], [11, 'XI']], ['M8.5 at 80 km depth 15', [8.5, 80, 15], [9, 'IX']], ['M8.5 at 300 km depth 30', [8.5, 300, 30], [8, 'VIII']], ['M8.5 at 900 km depth 10', [8.5, 900, 10], [7, 'VII']]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
M2.0 at 0 km depth 0.5[4, 'IV'][4, 'IV']Passed
M2.0 at 0 km depth 5[3, 'III'][3, 'III']Passed
M2.0 at 3 km depth 2[3, 'III'][3, 'III']Passed
M2.0 at 20 km depth 10[1, 'I'][1, 'I']Passed
M2.0 at 80 km depth 15[0, 'XII'][1, 'I']Failed
M2.0 at 300 km depth 30[0, 'XII'][1, 'I']Failed
M2.0 at 900 km depth 10[0, 'XII'][1, 'I']Failed
M4.5 at 0 km depth 0.5[8, 'VIII'][8, 'VIII']Passed

SHA-256 / 2875ce704cbfedc68c2b00688b47093ae52dd905cdb20764a3718deee91ed51f

3 / The verified repair

Exit 0
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(mag, epi_km, depth_km):
    roman = ["I", "II", "III", "IV", "V", "VI", "VII", "VIII", "IX", "X", "XI", "XII"]
    r = max(math.hypot(epi_km, depth_km), 1.0)
    i = 1.45 * mag - 2.46 * math.log10(r) + 1.5
    v = int(math.floor(i + 0.5))
    v = min(max(v, 1), 12)
    return [v, roman[v - 1]]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['M2.0 at 0 km depth 0.5', [2.0, 0, 0.5], [4, 'IV']], ['M2.0 at 0 km depth 5', [2.0, 0, 5], [3, 'III']], ['M2.0 at 3 km depth 2', [2.0, 3, 2], [3, 'III']], ['M2.0 at 20 km depth 10', [2.0, 20, 10], [1, 'I']], ['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M4.5 at 0 km depth 0.5', [4.5, 0, 0.5], [8, 'VIII']], ['M4.5 at 0 km depth 5', [4.5, 0, 5], [6, 'VI']], ['M4.5 at 3 km depth 2', [4.5, 3, 2], [7, 'VII']], ['M4.5 at 20 km depth 10', [4.5, 20, 10], [5, 'V']], ['M4.5 at 80 km depth 15', [4.5, 80, 15], [3, 'III']], ['M4.5 at 300 km depth 30', [4.5, 300, 30], [2, 'II']], ['M4.5 at 900 km depth 10', [4.5, 900, 10], [1, 'I']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 900 km depth 10', [2.0, 900, 10], [1, 'I']], ['M6.0 at 0 km depth 0.5', [6.0, 0, 0.5], [10, 'X']], ['M6.0 at 0 km depth 5', [6.0, 0, 5], [8, 'VIII']], ['M6.0 at 3 km depth 2', [6.0, 3, 2], [9, 'IX']], ['M6.0 at 20 km depth 10', [6.0, 20, 10], [7, 'VII']], ['M6.0 at 80 km depth 15', [6.0, 80, 15], [5, 'V']], ['M6.0 at 300 km depth 30', [6.0, 300, 30], [4, 'IV']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M2.0 at 300 km depth 30', [2.0, 300, 30], [1, 'I']], ['M7.2 at 0 km depth 0.5', [7.2, 0, 0.5], [12, 'XII']], ['M7.2 at 0 km depth 5', [7.2, 0, 5], [10, 'X']], ['M7.2 at 3 km depth 2', [7.2, 3, 2], [11, 'XI']], ['M7.2 at 20 km depth 10', [7.2, 20, 10], [9, 'IX']], ['M7.2 at 80 km depth 15', [7.2, 80, 15], [7, 'VII']], ['M7.2 at 300 km depth 30', [7.2, 300, 30], [6, 'VI']]], [['M2.0 at 80 km depth 15', [2.0, 80, 15], [1, 'I']], ['M8.5 at 0 km depth 0.5', [8.5, 0, 0.5], [12, 'XII']], ['M8.5 at 0 km depth 5', [8.5, 0, 5], [12, 'XII']], ['M8.5 at 3 km depth 2', [8.5, 3, 2], [12, 'XII']], ['M8.5 at 20 km depth 10', [8.5, 20, 10], [11, 'XI']], ['M8.5 at 80 km depth 15', [8.5, 80, 15], [9, 'IX']], ['M8.5 at 300 km depth 30', [8.5, 300, 30], [8, 'VIII']], ['M8.5 at 900 km depth 10', [8.5, 900, 10], [7, 'VII']]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
M2.0 at 0 km depth 0.5[4, 'IV'][4, 'IV']Passed
M2.0 at 0 km depth 5[3, 'III'][3, 'III']Passed
M2.0 at 3 km depth 2[3, 'III'][3, 'III']Passed
M2.0 at 20 km depth 10[1, 'I'][1, 'I']Passed
M2.0 at 80 km depth 15[1, 'I'][1, 'I']Passed
M2.0 at 300 km depth 30[1, 'I'][1, 'I']Passed
M2.0 at 900 km depth 10[1, 'I'][1, 'I']Passed
M4.5 at 0 km depth 0.5[8, 'VIII'][8, 'VIII']Passed

SHA-256 / 8840399c775b03ffba3b29bed11d42e427a3d8bbf86d992f8289eafaea147c0c

Verification & scope

Stipulated deterministic teaching model of a seismological magnitude procedure; constants and tables are fixed by the contract and no claim of agency or standards conformance is made. This reproducer isolates one failure mechanism. Results cover the supplied fixtures. Variants within a family share a test contract and should remain grouped when constructing evaluation splits. Related mechanisms with a shared evaluation_group must also remain together; these controlled models are not independent production incidents.

Observations recorded using Python 3.12.14 at 2026-09-29T14:48:31.637851+00:00.

Case digest / 825ca47f98bc8a7eb3ec0057517bb4f7e6a76d0686733b3440390cb0d56550d9