From 0c324cbe4b9f329f8b77e6ae1f4a894e4134ec6c Mon Sep 17 00:00:00 2001 From: luca-s Date: Thu, 5 Feb 2026 18:28:37 +0100 Subject: [PATCH] [homogeneous] fix TravelTimeList.delta unit TravelTimeList.delta was in km but it should be degree --- libs/seiscomp/seismology/ttt/homogeneous.cpp | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/libs/seiscomp/seismology/ttt/homogeneous.cpp b/libs/seiscomp/seismology/ttt/homogeneous.cpp index 5edc3a213..d5718766d 100644 --- a/libs/seiscomp/seismology/ttt/homogeneous.cpp +++ b/libs/seiscomp/seismology/ttt/homogeneous.cpp @@ -111,7 +111,7 @@ double computeDistance(double lat1, double lon1, double *backAzimuth = nullptr) { double dist; Math::Geo::delazi(lat1, lon1, lat2, lon2, &dist, azimuth, backAzimuth); - return Math::Geo::deg2km(dist); + return dist; } @@ -189,7 +189,8 @@ Homogeneous::isInside(double lat, double lon, double dep) if ( dep < _minDepth || dep > _maxDepth ) { return false; } - double dist = computeDistance(lat, lon, _centerLat, _centerLon); + double dist = + Math::Geo::deg2km(computeDistance(lat, lon, _centerLat, _centerLon)); if ( dist > _radius ) { return false; } @@ -220,7 +221,7 @@ Homogeneous::compute(const char *phase, } // straight ray path since we are in a homogeneous media - double Hdist = computeDistance(lat1, lon1, lat2, lon2); + double Hdist = Math::Geo::deg2km(computeDistance(lat1, lon1, lat2, lon2)); double Vdist = dep1 + alt2/1000.; double distance = sqrt(Hdist*Hdist + Vdist*Vdist); // [km] @@ -262,7 +263,7 @@ Homogeneous::computeTime(const char *phase, } // straight ray path since we are in a homogeneous media - double Hdist = computeDistance(lat1, lon1, lat2, lon2); + double Hdist = Math::Geo::deg2km(computeDistance(lat1, lon1, lat2, lon2)); double Vdist = dep1 + alt2/1000.; double distance = sqrt(Hdist*Hdist + Vdist*Vdist); // [km]