Fix PCK rate computation

2019-12-12 21:17:22 -05:00 · 2019-12-12 21:17:22 -05:00 · e8c25e1ae9
parent 6f792c5865
commit e8c25e1ae9
2 changed files with 2 additions and 1 deletions
--- a/jplephem/init.py
+++ b/jplephem/init.py
@ -292,7 +292,7 @@ You can ask for velocity as well.
 >>> print(r)
 [3.928e-02 3.878e-01 3.253e+03]
 >>> print(v)
-[2.318e-03 1.672e-04 9.177e-01]
+[6.707e-09 4.838e-10 2.655e-06]

 Legacy Ephemeris Packages
 -------------------------
--- a/jplephem/pck.py
+++ b/jplephem/pck.py
@ -216,6 +216,7 @@ class Segment(object):

        #rates = (dT.T * coefficients).sum(axis=2)
        rates = w0 + s * dw0 - dw1
+        rates /= intlen / 2.0
        if scalar:
            rates = rates[:,0]