Merge pull request #1576 from valdaarhun/fix_h_vec_units

Fix unit of h_vec in cross product calculation
poliastro · Sep 7, 2022 · 88c16fe · 88c16fe
2 parents a4f02b6 + 302b2ba
commit 88c16fe
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/poliastro/twobody/orbit/scalar.py b/src/poliastro/twobody/orbit/scalar.py
@@ -173,7 +173,7 @@ def e_vec(self):
     def h_vec(self):
         """Specific angular momentum vector."""
         h_vec = (
-            np.cross(self.r.to_value(u.km), self.v.to(u.km / u.s))
+            np.cross(self.r.to_value(u.km), self.v.to_value(u.km / u.s))
             * u.km**2
             / u.s
         )