Fixed it so true anomaly is put on the range 0..2pi

MintonGroup · Dec 4, 2022 · eb34415 · eb34415
1 parent cfcb505
commit eb34415
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/orbel/orbel.f90 b/src/orbel/orbel.f90
@@ -1037,6 +1037,7 @@ pure module subroutine orbel_xv2el(mu, px, py, pz, vx, vy, vz, a, e, inc, capom,
          rdot = sign(sqrt(v2 - (h / r)**2),rdotv)
          sf = a * (1.0_DP - e**2) / (h * e) * rdot
          f = atan2(sf,cf)
+         if (f < 0.0_DP) f = f + TWOPI
       else
          f = u
       end if