@@ -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
930922void meca_dc2axe (st_me meca , struct SEIS_AXIS * T , struct SEIS_AXIS * N , struct SEIS_AXIS * P ) {
0 commit comments