From 5cfafa7ff9fa5f902e54f397a9658b70076ee160 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 14 Nov 2025 12:16:12 +0800 Subject: [PATCH] Revert "Fix seis/utilmeca (#8821)" This reverts commit bf5ba5b021f6a551688fd7e8b1f4133ac3a3bf92. --- src/seis/utilmeca.c | 22 +++++++--------------- 1 file changed, 7 insertions(+), 15 deletions(-) diff --git a/src/seis/utilmeca.c b/src/seis/utilmeca.c index c6210c05eff..31142a1bc8f 100644 --- a/src/seis/utilmeca.c +++ b/src/seis/utilmeca.c @@ -910,21 +910,13 @@ void meca_axe2dc (struct SEIS_AXIS T, struct SEIS_AXIS P, struct SEIS_NODAL_PLAN } if (p2 < 0.0) p2 += 360.0; - NP1->dip = d1; NP1->str = p1; - NP2->dip = d2; NP2->str = p2; - - /* Epsilon-aware tie-break: avoid quadrant swap when P.dip is almost equal T.dip */ - if (P.dip > T.dip + SEIS_EPSILON) { - im = -1; - } else if (P.dip < T.dip - SEIS_EPSILON) { - im = 1; - } else { - /* |P.dip - T.dip| <= SEIS_EPSILON: deterministic choice */ - im = -1; - } - - NP1->rake = meca_computed_rake2 (NP2->str, NP2->dip, NP1->str, NP1->dip, im); - NP2->rake = meca_computed_rake2 (NP1->str, NP1->dip, NP2->str, NP2->dip, im); + NP1->dip = d1; NP1->str = p1; + NP2->dip = d2; NP2->str = p2; + + im = 1; + if (P.dip > T.dip) im = -1; + NP1->rake = meca_computed_rake2 (NP2->str, NP2->dip, NP1->str, NP1->dip, im); + NP2->rake = meca_computed_rake2 (NP1->str, NP1->dip, NP2->str, NP2->dip, im); } void meca_dc2axe (st_me meca, struct SEIS_AXIS *T, struct SEIS_AXIS *N, struct SEIS_AXIS *P) {