From 56fd0590df812247ba1259414bc9de0fafd424cc Mon Sep 17 00:00:00 2001 From: Bill-Gray Date: Sun, 26 Nov 2023 12:14:37 -0500 Subject: [PATCH] The velocity due to topocentric offset was computed using the earth's angular rotation rate, resulting in incorrect velocity offsets for topocentric ephemerides on other planets. --- mpc_obs.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/mpc_obs.cpp b/mpc_obs.cpp index f3cef142..b5844570 100644 --- a/mpc_obs.cpp +++ b/mpc_obs.cpp @@ -1484,8 +1484,8 @@ static inline int compute_topocentric_offset( const double ut, spin_matrix( precess_matrix, precess_matrix + 3, lon); for( i = 0; i < 3; i++) { - const double omega = 2. * PI * 360.9856235 / 360.; - /* earth's rotation rate, in radians/day */ + const double omega = planet_rotation_rate( planet_no, 0) * PI / 180.; + /* planet rotation rate, in radians/day */ if( offset) offset[i] = (rho_cos_phi * precess_matrix[i]