Round-offs in elements.

phoebe-project · Oct 22, 2024 · 214b0fb · 214b0fb
1 parent b5e61b1
commit 214b0fb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/phoebe/dynamics/orbel.py b/phoebe/dynamics/orbel.py
@@ -106,7 +106,7 @@ def orbel_xv2el(gm, r, v):
 
     fac = sqrt(h[0]*h[0] + h[1]*h[1])/absh
 
-    if fac < TINY:
+    if fac < TINY or inc == 0.0 or inc == pi:
         capom = 0.0
         u = arctan2(r[1], r[0])
         if abs(inc - pi) < 10.0*TINY: