ssb_to_psb_xyz_ICRS

nanograv · Sep 24, 2023 · c22a006 · c22a006
1 parent d3c7be9
commit c22a006
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/pint/models/astrometry.py b/src/pint/models/astrometry.py
@@ -453,9 +453,11 @@ def ssb_to_psb_xyz_ICRS(self, epoch=None):
             ra = ra0 + pmra * dt / np.cos(dec0)
             dec = dec0 + pmdec * dt
 
-        return u.Quantity(
-            [np.cos(dec) * np.cos(ra), np.cos(dec) * np.sin(ra), np.sin(dec)]
-        ).T
+        ra = ra.to_value(u.rad)
+        dec = dec.to_value(u.rad)
+
+        cos_dec = np.cos(dec)
+        return u.Quantity([cos_dec * np.cos(ra), cos_dec * np.sin(ra), np.sin(dec)]).T
 
     def get_params_as_ICRS(self):
         return {