Fix bug in computing theta in _components.py

Cosmoglobe · Aug 23, 2021 · deb3e3c · deb3e3c
1 parent f76fafa
commit deb3e3c
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/zodipy/_components.py b/zodipy/_components.py
@@ -139,7 +139,11 @@ def get_coordinates(
 
         x_earth_prime = x_earth - self.x_0
         y_earth_prime = y_earth - self.y_0
-        θ_prime = np.arctan2(y_prime - y_earth_prime , x_prime - x_earth_prime)
+
+        θ_prime = (
+            np.arctan2(y_prime , x_prime) 
+            - np.arctan2(y_earth_prime , x_earth_prime)
+        )
 
         return (R_prime, Z_prime, θ_prime), R_helio