From e8c25e1ae98faeca40e3dcb53f604da1c638aadf Mon Sep 17 00:00:00 2001
From: Brandon Rhodes <brandon@rhodesmill.org>
Date: Thu, 12 Dec 2019 21:17:22 -0500
Subject: [PATCH] Fix PCK rate computation

---
 jplephem/__init__.py | 2 +-
 jplephem/pck.py      | 1 +
 2 files changed, 2 insertions(+), 1 deletion(-)

diff --git a/jplephem/__init__.py b/jplephem/__init__.py
index 06270c5..8b7c3e7 100644
--- 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
 -------------------------
diff --git a/jplephem/pck.py b/jplephem/pck.py
index 295b156..0efc3ce 100644
--- 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]