Update phase angle

.github/workflows/flake8.yml

@@ -7,9 +7,9 @@ jobs:
     runs-on: ubuntu-latest
     name: flake8
     steps:
-    - uses: actions/checkout@v1
+    - uses: actions/checkout@v4
     - name: Set up Python 3.9
-      uses: actions/setup-python@v1
+      uses: actions/setup-python@v5
       with:
         python-version: 3.9
     - name: Install dependencies

.github/workflows/mypy.yml

@@ -7,9 +7,9 @@ jobs:
     runs-on: ubuntu-latest
     name: mypy
     steps:
-    - uses: actions/checkout@v1
+    - uses: actions/checkout@v4
     - name: Set up Python 3.9
-      uses: actions/setup-python@v1
+      uses: actions/setup-python@v5
       with:
         python-version: 3.9
     - name: Install dependencies

pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "tess-ephem"
-version = "0.6.3"
+version = "0.6.4"
 description = "Where are Solar System objects located in TESS FFI data?"
 license = "MIT"
 authors = ["Geert Barentsen <hello@geert.io>",

src/tess_ephem/__init__.py

@@ -6,5 +6,5 @@
 
 from .ephem import TessEphem, ephem  # noqa: E402
 
-__version__ = "0.6.3"
+__version__ = "0.6.4"
 __all__ = ["ephem", "TessEphem"]

src/tess_ephem/ephem.py

@@ -69,8 +69,8 @@ def __init__(
         self._deltaf = CubicSpline(eph["datetime_jd"], eph["delta"])
         # Delta-T (TDB - UT)
         self._deltat = CubicSpline(eph["datetime_jd"], eph["TDB-UT"])
-        # Phase angle
-        self._phif = CubicSpline(eph["datetime_jd"], eph["alpha_true"])
+        # Sun-Target-Observer (~phase) angle
+        self._alphaf = CubicSpline(eph["datetime_jd"], eph["alpha"])
         # Motion expressed in TESS pixels per hour
         motion = (
             np.hypot(eph["RA_rate"].quantity, eph["DEC_rate"].quantity)
@@ -105,7 +105,7 @@ def predict_sky(self, time: Time) -> DataFrame:
             H = np.full(len(time.jd), np.nan)
         r = self._rf(time.jd)
         delta = self._deltaf(time.jd)
-        phi = self._phif(time.jd)
+        sto_angle = self._alphaf(time.jd)
         motion = self._motionf(time.jd)
         deltat = self._deltat(time.jd)
         return DataFrame(
@@ -119,7 +119,7 @@ def predict_sky(self, time: Time) -> DataFrame:
                 "hmag": H,
                 "sun_distance": r,
                 "obs_distance": delta,
-                "phase_angle": phi,
+                "sto_angle": sto_angle,
             }
         )
 
@@ -288,7 +288,7 @@ def _get_horizons_ephem(
     step: str = "12H",
     id_type: str = "smallbody",
     location: str = "@TESS",
-    quantities: str = "1,3,9,19,20,30,43",
+    quantities: str = "1,3,9,19,20,24,30",
 ):
     """Returns JPL Horizons ephemeris.
 

tests/test_ephem.py

@@ -49,7 +49,7 @@ def test_predict():
     # 2. using tessrip to get average WCS from sector/camera/ccd
     # 3. converting RA,Dec to col,row with wcs_world2pix()
     assert np.round(pixel_locations.iloc[0]["row"], 1) == 1107.6
-    assert np.round(pixel_locations.iloc[0]["column"], 1) == 1087.8
+    assert np.round(pixel_locations.iloc[0]["column"], 1) == 1087.9
 
 
 def test_orbital_elements():