Skip to content

Commit 5cfafa7

Browse files
committed
Revert "Fix seis/utilmeca (#8821)"
This reverts commit bf5ba5b.
1 parent 41423e5 commit 5cfafa7

File tree

1 file changed

+7
-15
lines changed

1 file changed

+7
-15
lines changed

src/seis/utilmeca.c

Lines changed: 7 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -910,21 +910,13 @@ void meca_axe2dc (struct SEIS_AXIS T, struct SEIS_AXIS P, struct SEIS_NODAL_PLAN
910910
}
911911
if (p2 < 0.0) p2 += 360.0;
912912

913-
NP1->dip = d1; NP1->str = p1;
914-
NP2->dip = d2; NP2->str = p2;
915-
916-
/* Epsilon-aware tie-break: avoid quadrant swap when P.dip is almost equal T.dip */
917-
if (P.dip > T.dip + SEIS_EPSILON) {
918-
im = -1;
919-
} else if (P.dip < T.dip - SEIS_EPSILON) {
920-
im = 1;
921-
} else {
922-
/* |P.dip - T.dip| <= SEIS_EPSILON: deterministic choice */
923-
im = -1;
924-
}
925-
926-
NP1->rake = meca_computed_rake2 (NP2->str, NP2->dip, NP1->str, NP1->dip, im);
927-
NP2->rake = meca_computed_rake2 (NP1->str, NP1->dip, NP2->str, NP2->dip, im);
913+
NP1->dip = d1; NP1->str = p1;
914+
NP2->dip = d2; NP2->str = p2;
915+
916+
im = 1;
917+
if (P.dip > T.dip) im = -1;
918+
NP1->rake = meca_computed_rake2 (NP2->str, NP2->dip, NP1->str, NP1->dip, im);
919+
NP2->rake = meca_computed_rake2 (NP1->str, NP1->dip, NP2->str, NP2->dip, im);
928920
}
929921

930922
void meca_dc2axe (st_me meca, struct SEIS_AXIS *T, struct SEIS_AXIS *N, struct SEIS_AXIS *P) {

0 commit comments

Comments
 (0)