diff --git a/CHANGES.rst b/CHANGES.rst index e6d8e67a4..611ccc84c 100644 --- a/CHANGES.rst +++ b/CHANGES.rst @@ -16,6 +16,10 @@ API Changes - Deprecated `sbpy.ginga_plugins` in favor of using `sbpy-ginga` at https://github.com/NASA-Planetary-Science/sbpy-ginga [#413] +sbpy.data.ephem +^^^^^^^^^^^^^^^ + +- New command-line script: ``sbpy-ephem``. [#396] Other Changes and Additions --------------------------- diff --git a/docs/sbpy/data/ephem.rst b/docs/sbpy/data/ephem.rst index 01eb52bc6..9f50801d3 100644 --- a/docs/sbpy/data/ephem.rst +++ b/docs/sbpy/data/ephem.rst @@ -309,4 +309,43 @@ Telescope: The properties computed by pyoorb and listed in the resulting table are defined in the `pyoorb documentation `_. Note that this function -requires pyoorb to be installed, which is not a requirement for `sbpy`. \ No newline at end of file +requires pyoorb to be installed, which is not a requirement for `sbpy`. + + +Executable scripts +------------------ + +Ephemerides may be generated from the command line using the ``sbpy-ephem`` +script installed during sbpy's installation. The script supports ephemeris +generation from Horizons, the Minor Planet Center, and Miriade. + +For example, to display the ephemeris of comet 2P/Encke from the Minor Planet +Center: + +.. code:: bash + + $ sbpy-ephem mpc 2P + # requested target: 2P + # returned target: 2P + # location: 500 + date RA Dec rh delta phase solarelong V Proper motion Direction + deg deg AU AU deg deg mag arcsec / h deg + ----------------------- ----------- ----------- ----- ----- ----- ---------- ---- ------------- --------- + 2024-08-17 00:00:00.000 21:55:57.10 -15:12:47.0 3.368 2.356 0.9 177.1 21.1 47.72 253.7 + 2024-08-18 00:00:00.000 21:54:41.20 -15:18:06.0 3.373 2.361 0.7 177.5 21.1 47.54 253.9 + 2024-08-19 00:00:00.000 21:53:25.60 -15:23:21.0 3.378 2.367 0.8 177.4 21.1 47.32 254.0 + 2024-08-20 00:00:00.000 21:52:10.30 -15:28:32.0 3.383 2.372 1.0 176.7 21.2 47.07 254.1 + 2024-08-21 00:00:00.000 21:50:55.30 -15:33:38.0 3.388 2.378 1.3 175.8 21.2 46.79 254.3 + 2024-08-22 00:00:00.000 21:49:40.70 -15:38:40.0 3.393 2.385 1.6 174.7 21.2 46.47 254.4 + 2024-08-23 00:00:00.000 21:48:26.50 -15:43:37.0 3.398 2.391 1.9 173.5 21.3 46.13 254.6 + 2024-08-24 00:00:00.000 21:47:12.90 -15:48:29.0 3.403 2.398 2.3 172.3 21.3 45.75 254.7 + 2024-08-25 00:00:00.000 21:45:59.80 -15:53:16.0 3.408 2.406 2.6 171.1 21.3 45.34 254.9 + 2024-08-26 00:00:00.000 21:44:47.30 -15:57:57.0 3.413 2.413 3.0 169.9 21.4 44.9 255.0 + 2024-08-27 00:00:00.000 21:43:35.40 -16:02:33.0 3.418 2.421 3.3 168.7 21.4 44.44 255.2 + 2024-08-28 00:00:00.000 21:42:24.30 -16:07:04.0 3.423 2.43 3.7 167.4 21.4 43.94 255.3 + 2024-08-29 00:00:00.000 21:41:13.90 -16:11:28.0 3.427 2.438 4.0 166.2 21.5 43.41 255.5 + 2024-08-30 00:00:00.000 21:40:04.30 -16:15:47.0 3.432 2.447 4.4 165.0 21.5 42.86 255.6 + ... + +A limited number of input parameters are supported. Run with the ``--help`` +option for more details. \ No newline at end of file diff --git a/sbpy/data/__init__.py b/sbpy/data/__init__.py index 68aa6d44e..0a5a471b2 100644 --- a/sbpy/data/__init__.py +++ b/sbpy/data/__init__.py @@ -48,8 +48,19 @@ class Conf: fieldnames_info = [ # General { - "description": "Target Identifier", - "fieldnames": ["targetname", "id", "Object"], + "description": "Target Name", + "fieldnames": [ + "targetname", + "Targetname", + "target", + "Target", + "name", + "Name", + "id", + "ID", + "Object", + "object", + ], "provenance": ["orbit", "ephem", "obs", "phys"], "dimension": None, }, @@ -61,13 +72,7 @@ class Conf: }, { "description": "Target Number", - "fieldnames": ["number"], - "provenance": ["orbit", "ephem", "obs", "phys"], - "dimension": None, - }, - { - "description": "Target Name", - "fieldnames": ["name"], + "fieldnames": ["number", "target number"], "provenance": ["orbit", "ephem", "obs", "phys"], "dimension": None, }, diff --git a/sbpy/data/ephem/__init__.py b/sbpy/data/ephem/__init__.py new file mode 100644 index 000000000..6d0ee0cbe --- /dev/null +++ b/sbpy/data/ephem/__init__.py @@ -0,0 +1 @@ +from .core import Ephem # noqa: F401 diff --git a/sbpy/data/ephem/cli.py b/sbpy/data/ephem/cli.py new file mode 100644 index 000000000..ba92306c4 --- /dev/null +++ b/sbpy/data/ephem/cli.py @@ -0,0 +1,351 @@ +# Licensed under a 3-clause BSD style license - see LICENSE.rst +""" +=================== +sbpy.data.ephem.cli +=================== + +Ephemerides may be generated by this module via the ``sbpy-ephem`` command line script. + +""" + +import numpy as np +import sys +import argparse +from typing import Any, Dict, List, Optional, Set, Union + +import astropy.units as u +from astropy.time import Time +from astropy.coordinates import Angle + +from .core import Ephem, EphemerisService +from ..names import Names, TargetNameParseError +from ...utils.decorators import requires + + +class StoreEphemerisService(argparse.Action): + """Helper class to support ``EphemerisService`` objects in argparse + arguments.""" + + def __call__(self, parser, namespace, values, option_string=None): + value = EphemerisService(values) + setattr(namespace, self.dest, value) + + +class EphemerisCLI: + """Command-line interface for ephemeris generation. + + Use this class via the command line-line script ``sbpy-ephem``. + + + Parameters + ---------- + argv : list of str, optional + Parse command line parameters from this list, rather than sys.argv. + + """ + + @requires("astroquery") + def __init__(self, argv: Optional[List[str]] = None) -> None: + self.args: argparse.Namespace = self.parse_args(argv) + self.epochs: Dict[str, Any] = self._format_epochs(self.args) + + self.eph: Ephem + if self.args.service == EphemerisService.mpc: + self.eph = self.mpc() + elif self.args.service == EphemerisService.horizons: + self.eph = self.horizons() + elif self.args.service == EphemerisService.miriade: + self.eph = self.miriade() + + self.eph: Ephem = self._format_eph(self.eph) + + @classmethod + def run(cls) -> None: # pragma: no cover + """Command-line script entry point.""" + try: + cli = cls() + except Exception as e: + # check for debug flag + if any([x == "--debug" for x in sys.argv]): + raise e + print(str(e), file=sys.stderr) + sys.exit(1) + + print( + f"""# requested target: {cli.args.target} +# returned target: {cli.eph.meta['target']} +# location: {cli.args.location}""" + ) + cli.eph.table.pprint_all() + + @staticmethod + def parse_args(argv: Optional[List[str]] = None) -> argparse.Namespace: + """Parse command-line arguments. + + Parameters + ---------- + argv : list of str, optional + Parse command line parameters from this list, rather than sys.argv. + + """ + + today: Time = Time(Time.now().iso[:10]) + + parser: argparse.ArgumentParser = argparse.ArgumentParser() + # names in choices must have entries in `EphemerisService` + parser.add_argument( + "service", + choices=[ + str(EphemerisService.horizons), + str(EphemerisService.miriade), + str(EphemerisService.mpc), + ], + action=StoreEphemerisService, + help="online service to use for ephemeris", + ) + parser.add_argument( + "target", + ) + parser.add_argument( + "--location", + "-l", + default="500", + help="observer location as a Minor Planet Center observatory code, e.g., G37", + ) + parser.add_argument( + "--start", + default=today, + type=Time, + help="start time (e.g., YYYY-MM-DD), default today [UTC]", + ) + parser.add_argument( + "--step", + default="1 d", + type=u.Quantity, + help="ephemeris step size with units (d, h, min)", + ) + + stop_number_mutex = parser.add_mutually_exclusive_group() + stop_number_mutex.add_argument( + "--stop", + "--end", + type=Time, + help="stop time (e.g., YYYY-MM-DD) [UTC]", + ) + stop_number_mutex.add_argument( + "--number", + "-n", + type=int, + help="number of time steps", + ) + + parser.add_argument( + "--radec", + "-r", + default="hmsdms", + choices=("hmsdms", "deg"), + help="RA, Dec coordinate output format", + ) + + horizons: argparse._ArgumentGroup = parser.add_argument_group( + "JPL Horizons options" + ) + horizons.add_argument( + "--id-type", + "-i", + default=None, + choices=( + "smallbody", + "designation", + "name", + "asteroid_name", + "comet_name", + ), + help=( + "target identifer type; automatically set to 'designation' for " + "target names that appear to be a cometary designation" + ), + ) + horizons.add_argument( + "--quantities", + "-q", + default="1,3,9,19,20,23,24,27,33", + help="table quantities to return as a comma-separated list of integers (default: 1,3,9,19,20,23,24,27,33)", + ) + + miriade: argparse._ArgumentGroup = parser.add_argument_group( + "IMCCE Miriade options" + ) + miriade.add_argument( + "--type", + "-t", + default="asteroid", + choices=( + "asteroid", + "comet", + "dwarf planet", + "planet", + "satellite", + ), + help="target identifer type", + ) + + parser.add_argument( + "--debug", action="store_true", help="print debugging information" + ) + + args = parser.parse_args(argv) + + # parameter checks: + if args.stop is not None and args.start > args.stop: + raise ValueError("start cannot be after stop") + + return args + + @staticmethod + def _format_epochs(args) -> Dict[str, Any]: + """Format the epochs dictionary.""" + + epochs: dict = { + "start": args.start, + "step": args.step, + "stop": args.stop, + } + defined_attributes: Set[str] = { + k + for k in ["start", "stop", "step", "number"] + if getattr(args, k) is not None + } + + def given(*attributes): + return set(attributes) == defined_attributes + + # start and step are always defined by parse_args + if given("start", "step"): + # default is 60 steps + epochs["stop"] = epochs["start"] + 60 * epochs["step"] + elif given("start", "step", "number"): + epochs["stop"] = epochs["start"] + args.number * epochs["step"] + + return epochs + + def mpc(self) -> Ephem: + """Get ephemeris from the Minor Planet Center.""" + + eph = Ephem.from_mpc( + self.args.target, + epochs=self.epochs, + location=self.args.location, + ) + return eph + + def horizons(self) -> Ephem: + """Get ephemeris from Horizons.""" + + # comet specific options: always avoid multiple fragment matches, always + # return closest apparition + closest_apparition: bool = False + no_fragments: bool = False + id_type: Union[str, None] = self.args.id_type + try: + Names.parse_comet(self.args.target) + jd_mid = np.mean((self.epochs["start"].jd, self.epochs["stop"].jd)) + closest_apparition = f"<{jd_mid:.0f}" + no_fragments = True + id_type = "designation" if id_type is None else id_type + except TargetNameParseError: + pass + + eph = Ephem.from_horizons( + self.args.target, + id_type=id_type, + closest_apparition=closest_apparition, + no_fragments=no_fragments, + epochs=self.epochs, + location=self.args.location, + quantities=self.args.quantities, + ) + + # simplify output + eph["epoch"].format = "iso" + eph["ra"] = eph["ra"].unmasked + eph["dec"] = eph["dec"].unmasked + for k in ["M1", "k1", "solar_presence", "lunar_presence"]: + if k in eph.table.colnames: + del eph.table[k] + + return eph + + def miriade(self) -> Ephem: + """Get ephemeris from Miriade.""" + + eph = Ephem.from_miriade( + self.args.target, + objtype=self.args.type, + epochs=self.epochs, + location=self.args.location, + coordtype=1, + ) + + # simplify output + eph["epoch"].format = "iso" + + for k in eph.field_names: + if hasattr(eph[k], "unmasked"): + eph[k] = eph[k].unmasked + + return eph + + def _format_eph(self, eph: Ephem) -> Ephem: + # convert RA and Dec to Angle in units of degree + eph["ra"] = Angle(eph["ra"], eph["ra"].unit).to("deg") + eph["dec"] = Angle(eph["dec"], eph["dec"].unit).to("deg") + + if self.args.radec == "hmsdms": + eph["ra"].info.format = lambda x: x.to_string( + sep=":", precision=2, unit="hourangle" + ) + eph["dec"].info.format = lambda x: x.to_string(sep=":", precision=1) + else: + eph["ra"].info.format = lambda x: x.to_string( + unit="deg", decimal=True, precision=5 + ) + eph["dec"].info.format = lambda x: x.to_string( + unit="deg", decimal=True, precision=5 + ) + + # normalize output for most common columns + first_fields: List[str] = [ + "target", + "date", + "RA", + "Dec", + "rh", + "delta", + "phase", + "solarelong", + "DEC_rate", + "RA*cos(Dec)_rate", + ] + + k: str + translated: str + for k in first_fields: + if k not in eph: + continue + translated = eph._translate_columns(k)[0] + eph.table.rename_column(translated, k) + + # re-order + fields: List[str] = eph.field_names + for k in reversed(first_fields): + if k not in fields: + continue + fields.insert(0, fields.pop(fields.index(k))) + + eph = eph[fields] + target = eph["target"][0] + del eph.table["target"] + eph.table.meta["target"] = target + + return eph diff --git a/sbpy/data/ephem.py b/sbpy/data/ephem/core.py similarity index 66% rename from sbpy/data/ephem.py rename to sbpy/data/ephem/core.py index 2328cd3db..37a2fa5f5 100644 --- a/sbpy/data/ephem.py +++ b/sbpy/data/ephem/core.py @@ -1,14 +1,17 @@ # Licensed under a 3-clause BSD style license - see LICENSE.rst """ -====================== -sbpy data.Ephem Module -====================== +==================== +sbpy.data.ephem.core +==================== -Class for storing and querying ephemerides +Defines the Ephem class for storing and querying ephemerides. created on June 04, 2017 + """ + import os +import enum from warnings import warn import numpy as np @@ -24,21 +27,21 @@ from astroquery.mpc import MPC from astroquery.imcce import Miriade from astroquery.exceptions import InvalidQueryError -except ImportError: +except ImportError: # pragma: no cover pass try: import pyoorb -except ImportError: +except ImportError: # pragma: no cover pass -from ..bib import cite -from .core import DataClass, Conf, QueryError, TimeScaleWarning -from ..exceptions import RequiredPackageUnavailable -from .orbit import Orbit, OpenOrbError -from ..utils.decorators import requires +from ...bib import cite +from ..core import DataClass, Conf, QueryError, TimeScaleWarning +from ...exceptions import RequiredPackageUnavailable +from ..orbit import Orbit, OpenOrbError +from ...utils.decorators import requires -__all__ = ['Ephem'] +__all__ = ["Ephem"] __doctest_requires__ = { @@ -47,15 +50,34 @@ } +# include one entry for every supported ephemeris service +class EphemerisService(enum.Enum): + horizons = "horizons" + mpc = "mpc" + miriade = "miriade" + oorb = "oorb" + + def __str__(self): + return self.value + + def __repr__(self): + return self.value + + class Ephem(DataClass): """Class for querying, manipulating, and calculating ephemerides""" @classmethod @requires("astroquery") - @cite({'data source': '1996DPS....28.2504G', - 'software: astroquery': '2019AJ....157...98G'}) - def from_horizons(cls, targetids, id_type='smallbody', - epochs=None, location='500', **kwargs): + @cite( + { + "data source": "1996DPS....28.2504G", + "software: astroquery": "2019AJ....157...98G", + } + ) + def from_horizons( + cls, targetids, id_type="smallbody", epochs=None, location="500", **kwargs + ): """Load target ephemerides from `JPL Horizons `_ using `astroquery.jplhorizons.HorizonsClass.ephemerides` @@ -140,10 +162,14 @@ def from_horizons(cls, targetids, id_type='smallbody', if epochs is None: _epochs = [Time.now().utc.jd] elif isinstance(epochs, Time): - if epochs.scale != 'utc': - warn(('converting {} epochs to utc for use in ' - 'astroquery.jplhorizons').format(epochs.scale), - TimeScaleWarning) + if epochs.scale != "utc": + warn( + ( + "converting {} epochs to utc for use in " + "astroquery.jplhorizons" + ).format(epochs.scale), + TimeScaleWarning, + ) _epochs = epochs.utc.jd # ~1 ms precision for modern dates # 2022 Feb 04: JPL Horizons API v1.1 will return sorted dates, but @@ -158,22 +184,23 @@ def from_horizons(cls, targetids, id_type='smallbody', elif isinstance(epochs, dict): _epochs = epochs.copy() - if 'start' in _epochs and 'stop' in _epochs and 'number' in epochs: - _epochs['step'] = _epochs['number']*u.dimensionless_unscaled + if "start" in _epochs and "stop" in _epochs and "number" in epochs: + _epochs["step"] = _epochs["number"] * u.dimensionless_unscaled # convert to utc and iso for astroquery.jplhorizons - _epochs['start'] = _epochs['start'].utc.iso - _epochs['stop'] = _epochs['stop'].utc.iso - if 'step' in _epochs: - if _epochs['step'].unit is not u.dimensionless_unscaled: - _epochs['step'] = '{:d}{:s}'.format( - int(_epochs['step'].value), - {u.minute: 'm', u.hour: 'h', u.d: 'd', - u.year: 'y'}[_epochs['step'].unit]) + _epochs["start"] = _epochs["start"].utc.iso + _epochs["stop"] = _epochs["stop"].utc.iso + if "step" in _epochs: + if _epochs["step"].unit is not u.dimensionless_unscaled: + _epochs["step"] = "{:d}{:s}".format( + int(_epochs["step"].value), + {u.minute: "m", u.hour: "h", u.d: "d", u.year: "y"}[ + _epochs["step"].unit + ], + ) else: - _epochs['step'] = '{:d}'.format( - int(_epochs['step'].value-1)) + _epochs["step"] = "{:d}".format(int(_epochs["step"].value - 1)) else: - raise ValueError('Invalid `epochs` parameter') + raise ValueError("Invalid `epochs` parameter") # if targetids is a list, run separate Horizons queries and append if not isinstance(targetids, (list, ndarray, tuple)): @@ -182,24 +209,32 @@ def from_horizons(cls, targetids, id_type='smallbody', # turn EarthLocation into dictionary of strings as used by # astroquery.jplhorizons if isinstance(location, EarthLocation): - location = {'lon': location.lon.deg, - 'lat': location.lat.deg, - 'elevation': location.height.to('km')} + location = { + "lon": location.lon.deg, + "lat": location.lat.deg, + "elevation": location.height.to("km"), + } - # append ephemerides table for each targetid + # append ephemerides table for each targetid, collect all URLs all_eph = None + urls = [] for targetid in targetids: # load ephemerides using astroquery.jplhorizons - obj = Horizons(id=targetid, id_type=id_type, - location=location, epochs=_epochs) + obj = Horizons( + id=targetid, id_type=id_type, location=location, epochs=_epochs + ) try: eph = obj.ephemerides(**kwargs) except ValueError as e: raise QueryError( - ('Error raised by astroquery.jplhorizons: {:s}\n' - 'The following query was attempted: {:s}').format( - str(e), obj.uri)) + ( + "Error raised by astroquery.jplhorizons: {:s}\n" + "The following query was attempted: {:s}" + ).format(str(e), obj.uri) + ) + + urls.append(obj.uri) if _sorted: # restore user's original order @@ -210,14 +245,14 @@ def from_horizons(cls, targetids, id_type='smallbody', # should really change '---' units to None in # astroquery.jplhorizons.__init__.py for column_name in eph.columns: - if eph[column_name].unit == '---': + if eph[column_name].unit == "---": eph[column_name].unit = None # workaround for astroquery 0.3.9.dev5056 and earlier, # Horizons column named RA_rate always includes the # cos(Dec) term: - if 'RA_rate' in eph.colnames: - eph['RA_rate'].name = 'RA*cos(Dec)_rate' + if "RA_rate" in eph.colnames: + eph["RA_rate"].name = "RA*cos(Dec)_rate" if all_eph is None: all_eph = eph @@ -227,38 +262,47 @@ def from_horizons(cls, targetids, id_type='smallbody', # turn epochs into astropy.time.Time and apply timescale # convert ut1 epochs to utc # https://ssd.jpl.nasa.gov/horizons/manual.html - if any(all_eph['datetime_jd'] < 2437665.5): - all_eph['datetime_jd'][all_eph['datetime_jd'] < - 2437665.5] = Time( - all_eph['datetime_jd'][all_eph['datetime_jd'] < - 2437665.5], - scale='ut1', format='jd').utc.jd - all_eph['epoch'] = Time(all_eph['datetime_jd'], format='jd', - scale='utc') - if 'siderealtime' in all_eph.colnames: - all_eph['siderealtime'].unit = u.Unit('hour') - - all_eph.remove_column('datetime_jd') - all_eph.remove_column('datetime_str') + if any(all_eph["datetime_jd"] < 2437665.5): + all_eph["datetime_jd"][all_eph["datetime_jd"] < 2437665.5] = Time( + all_eph["datetime_jd"][all_eph["datetime_jd"] < 2437665.5], + scale="ut1", + format="jd", + ).utc.jd + all_eph["epoch"] = Time(all_eph["datetime_jd"], format="jd", scale="utc") + if "siderealtime" in all_eph.colnames: + all_eph["siderealtime"].unit = u.Unit("hour") + + all_eph.remove_column("datetime_jd") + all_eph.remove_column("datetime_str") # convert solartime and hour angle into Quantity - if 'solartime' in all_eph.colnames: - all_eph['solartime'] = Angle( - all_eph['solartime'], 'hr').hourangle * u.hr + if "solartime" in all_eph.colnames: + all_eph["solartime"] = Angle(all_eph["solartime"], "hr").hourangle * u.hr - if 'locapp_hourangle' in all_eph.colnames: - all_eph['locapp_hourangle'] = u.Quantity(Angle( - all_eph['locapp_hourangle'], 'hr'), u.hourangle) + if "locapp_hourangle" in all_eph.colnames: + all_eph["locapp_hourangle"] = u.Quantity( + Angle(all_eph["locapp_hourangle"], "hr"), u.hourangle + ) return cls.from_table(all_eph) @classmethod @requires("astroquery") - @cite({'data source': - 'https://minorplanetcenter.net/iau/MPEph/MPEph.html', - 'software: astroquery': '2019AJ....157...98G'}) - def from_mpc(cls, targetids, epochs=None, location='500', ra_format=None, - dec_format=None, **kwargs): + @cite( + { + "data source": "https://minorplanetcenter.net/iau/MPEph/MPEph.html", + "software: astroquery": "2019AJ....157...98G", + } + ) + def from_mpc( + cls, + targetids, + epochs=None, + location="500", + ra_format=None, + dec_format=None, + **kwargs, + ): """Load ephemerides from the `Minor Planet Center `_. @@ -382,49 +426,61 @@ def from_mpc(cls, targetids, epochs=None, location='500', ra_format=None, else: _epochs = epochs - if _epochs.scale != 'utc': - warn(('converting {} epochs to utc for use in ' - 'astroquery.mpc').format(_epochs.scale), - TimeScaleWarning) + if _epochs.scale != "utc": + warn( + ("converting {} epochs to utc for use in " "astroquery.mpc").format( + _epochs.scale + ), + TimeScaleWarning, + ) _epochs = _epochs.utc elif isinstance(epochs, dict): _epochs = epochs.copy() - start = _epochs['start'] # required - if start.scale != 'utc': - warn(('converting {} start epoch to utc for use in ' - 'astroquery.mpc').format(start.scale), - TimeScaleWarning) + start = _epochs["start"] # required + if start.scale != "utc": + warn( + ( + "converting {} start epoch to utc for use in " "astroquery.mpc" + ).format(start.scale), + TimeScaleWarning, + ) start = start.utc - step = _epochs.get('step') - stop = _epochs.get('stop') - if stop is not None and stop.scale != 'utc': - warn(('converting {} stop epoch to utc for use in ' - 'astroquery.mpc').format(stop.scale), - TimeScaleWarning) + step = _epochs.get("step") + stop = _epochs.get("stop") + if stop is not None and stop.scale != "utc": + warn( + ( + "converting {} stop epoch to utc for use in " "astroquery.mpc" + ).format(stop.scale), + TimeScaleWarning, + ) stop = stop.utc - number = _epochs.get('number') + number = _epochs.get("number") if step is not None and stop is None: step = u.Quantity(step) if step.unit not in (u.d, u.h, u.min, u.s): raise QueryError( - 'step must have units of days, hours, minutes,' - ' or seconds') + "step must have units of days, hours, minutes," " or seconds" + ) if stop is not None: if step is not None and number is None: # start and stop both defined, estimate number of steps dt = (Time(stop).jd - Time(start).jd) * u.d number = int((dt / step).decompose()) + 1 elif step is None and number is not None: - step = int((stop-start).jd*1440/(number-1))*u.minute + step = int((stop - start).jd * 1440 / (number - 1)) * u.minute else: raise QueryError( - ('epoch definition unclear; step xor number ' - 'must be provided with start and stop')) + ( + "epoch definition unclear; step xor number " + "must be provided with start and stop" + ) + ) else: - raise ValueError('Invalid `epochs` parameter') + raise ValueError("Invalid `epochs` parameter") - # append ephemerides table for each targetid + # append ephemerides table for each targetid and collect URLs all_eph = None for targetid in targetids: try: @@ -432,25 +488,31 @@ def from_mpc(cls, targetids, epochs=None, location='500', ra_format=None, if start is None: eph = [] for i in range(len(_epochs)): - e = MPC.get_ephemeris(targetid, location=location, - start=Time(_epochs[i], - scale='utc'), - number=1, **kwargs) - e['Date'] = e['Date'].iso # for vstack to work + e = MPC.get_ephemeris( + targetid, + location=location, + start=Time(_epochs[i], scale="utc"), + number=1, + **kwargs, + ) + e["Date"] = e["Date"].iso # for vstack to work eph.append(e) eph = QTable(vstack(eph)) - eph['Date'] = Time(eph['Date'], scale='utc') + eph["Date"] = Time(eph["Date"], scale="utc") else: - eph = MPC.get_ephemeris(targetid, location=location, - start=start, step=step, - number=number, **kwargs) + eph = MPC.get_ephemeris( + targetid, + location=location, + start=start, + step=step, + number=number, + **kwargs, + ) except InvalidQueryError as e: - raise QueryError( - 'Error raised by astroquery.mpc: {:s}'.format(str(e))) + raise QueryError("Error raised by astroquery.mpc: {:s}".format(str(e))) # add targetname column - eph.add_column(Column([targetid]*len(eph), - name='Targetname'), index=0) + eph.add_column(Column([targetid] * len(eph), name="Targetname"), index=0) if all_eph is None: all_eph = eph @@ -460,25 +522,30 @@ def from_mpc(cls, targetids, epochs=None, location='500', ra_format=None, all_eph = QTable(all_eph) # convert RA and Dec to Angle - all_eph['RA'] = Angle(all_eph['RA'], all_eph['RA'].unit) - all_eph['Dec'] = Angle(all_eph['Dec'], all_eph['Dec'].unit) + all_eph["RA"] = Angle(all_eph["RA"], all_eph["RA"].unit) + all_eph["Dec"] = Angle(all_eph["Dec"], all_eph["Dec"].unit) if ra_format is not None: - all_eph['RA'].info.format = lambda x: x.to_string(**ra_format) + all_eph["RA"].info.format = lambda x: x.to_string(**ra_format) if dec_format is not None: - all_eph['Dec'].info.format = lambda x: x.to_string(**dec_format) + all_eph["Dec"].info.format = lambda x: x.to_string(**dec_format) return cls.from_table(all_eph) @classmethod @requires("astroquery") - @cite({'data source': 'http://vo.imcce.fr/webservices/miriade/', - 'software: astroquery': '2019AJ....157...98G'}) - def from_miriade(cls, targetids, objtype='asteroid', - epochs=None, location='500', **kwargs): + @cite( + { + "data source": "https://vo.imcce.fr/webservices/miriade/", + "software: astroquery": "2019AJ....157...98G", + } + ) + def from_miriade( + cls, targetids, objtype="asteroid", epochs=None, location="500", **kwargs + ): """Load target ephemerides from - `IMCCE Miriade `_ using + `IMCCE Miriade `_ using `astroquery.imcce.MiriadeClass.get_ephemerides` Parameters @@ -525,7 +592,7 @@ def from_miriade(cls, targetids, objtype='asteroid', * For detailed explanations of the queried fields, refer to `astroquery.imcce.MiriadeClass.get_ephemerides` and the `Miriade documentation - `_. + `_. * By default, all properties are provided in the J2000.0 reference system. Different settings can be chosen using additional keyword arguments as used by @@ -553,41 +620,60 @@ def from_miriade(cls, targetids, objtype='asteroid', # format _epoch for astroquery.imcce.Miriade if epochs is None: - _epochs = {'start': Time.now().utc.jd} + _epochs = {"start": Time.now().utc.jd} elif isinstance(epochs, Time): - if epochs.scale != 'utc': - warn(('converting {} epochs to utc for use in ' - 'astroquery.imcce').format(epochs.scale), - TimeScaleWarning) - _epochs = {'start': epochs.utc} + if epochs.scale != "utc": + warn( + ( + "converting {} epochs to utc for use in " "astroquery.imcce" + ).format(epochs.scale), + TimeScaleWarning, + ) + _epochs = {"start": epochs.utc} elif isinstance(epochs, dict): _epochs = epochs.copy() - if _epochs['start'].scale != 'utc': - warn(('converting {} start epoch to utc for use in ' - 'astroquery.imcce').format(_epochs['start'].scale), - TimeScaleWarning) - _epochs['start'] = _epochs['start'].utc - if 'stop' in _epochs and _epochs['stop'].scale != 'utc': - warn(('converting {} stop epoch to utc for use in ' - 'astroquery.imcce').format(_epochs['stop'].scale), - TimeScaleWarning) - _epochs['stop'] = _epochs['stop'].utc - if 'number' in _epochs: + if _epochs["start"].scale != "utc": + warn( + ( + "converting {} start epoch to utc for use in " + "astroquery.imcce" + ).format(_epochs["start"].scale), + TimeScaleWarning, + ) + _epochs["start"] = _epochs["start"].utc + if "stop" in _epochs and _epochs["stop"].scale != "utc": + warn( + ( + "converting {} stop epoch to utc for use in " "astroquery.imcce" + ).format(_epochs["stop"].scale), + TimeScaleWarning, + ) + _epochs["stop"] = _epochs["stop"].utc + if "number" in _epochs: # turn interval/number into step size based on full minutes - _epochs['step'] = int((Time(_epochs['stop']) - - Time(_epochs['start'])).jd * - 86400 / (_epochs['number']-1)) * u.s - elif 'step' in _epochs: - _epochs['number'] = ((Time(_epochs['stop']) - - Time(_epochs['start'])).jd * - 86400 / _epochs['step'].to('s').value) + 1 - if 'step' in _epochs: - _epochs['step'] = '{:f}{:s}'.format( - _epochs['step'].value, - {u.s: 's', u.minute: 'm', u.hour: 'h', - u.day: 'd'}[_epochs['step'].unit]) + _epochs["step"] = ( + int( + (Time(_epochs["stop"]) - Time(_epochs["start"])).jd + * 86400 + / (_epochs["number"] - 1) + ) + * u.s + ) + elif "step" in _epochs: + _epochs["number"] = ( + (Time(_epochs["stop"]) - Time(_epochs["start"])).jd + * 86400 + / _epochs["step"].to("s").value + ) + 1 + if "step" in _epochs: + _epochs["step"] = "{:f}{:s}".format( + _epochs["step"].value, + {u.s: "s", u.minute: "m", u.hour: "h", u.day: "d"}[ + _epochs["step"].unit + ], + ) else: - raise ValueError('Invalid `epochs` parameter') + raise ValueError("Invalid `epochs` parameter") # if targetids is a list, run separate Horizons queries and append if not isinstance(targetids, (list, ndarray, tuple)): @@ -597,50 +683,56 @@ def from_miriade(cls, targetids, objtype='asteroid', # astroquery.jplhorizons if isinstance(location, EarthLocation): - location = '{:+f} {:+f} {:.1f}'.format( - location.lon.deg, - location.lat.deg, - location.height.to('m').value) + location = "{:+f} {:+f} {:.1f}".format( + location.lon.deg, location.lat.deg, location.height.to("m").value + ) # append ephemerides table for each targetid all_eph = None for targetid in targetids: query = Miriade() try: - if 'step' not in _epochs and 'number' not in _epochs: - if not iterable(_epochs['start']): + if "step" not in _epochs and "number" not in _epochs: + if not iterable(_epochs["start"]): # single epoch - eph = query.get_ephemerides(targetname=targetid, - objtype=objtype, - location=location, - epoch=_epochs['start'], - **kwargs) + eph = query.get_ephemerides( + targetname=targetid, + objtype=objtype, + location=location, + epoch=_epochs["start"], + **kwargs, + ) else: # multiple epochs eph = [] - for i in range(len(_epochs['start'])): + for i in range(len(_epochs["start"])): e = query.get_ephemerides( targetname=targetid, objtype=objtype, location=location, - epoch=_epochs['start'][i], - **kwargs + epoch=_epochs["start"][i], + **kwargs, ) eph.append(e) eph = vstack(eph) else: # dictionary eph = query.get_ephemerides( - targetname=targetid, objtype=objtype, - location=location, epoch=_epochs['start'], - epoch_step=_epochs['step'], - epoch_nsteps=_epochs['number'], - **kwargs) + targetname=targetid, + objtype=objtype, + location=location, + epoch=_epochs["start"], + epoch_step=_epochs["step"], + epoch_nsteps=_epochs["number"], + **kwargs, + ) except RuntimeError as e: raise QueryError( - ('Error raised by astroquery.imcce: {:s}\n' - 'The following query was attempted: {:s}').format( - str(e), query.uri)) + ( + "Error raised by astroquery.imcce: {:s}\n" + "The following query was attempted: {:s}" + ).format(str(e), query.uri) + ) if all_eph is None: all_eph = eph @@ -648,15 +740,21 @@ def from_miriade(cls, targetids, objtype='asteroid', all_eph = vstack([all_eph, eph]) all_eph = QTable(all_eph) - all_eph['epoch'] = Time(all_eph['epoch'], scale='utc', format='jd') + all_eph["epoch"] = Time(all_eph["epoch"], scale="utc", format="jd") return cls.from_table(all_eph) @classmethod @requires("pyoorb") - @cite({'method': '2009M&PS...44.1853G', - 'software': 'https://github.com/oorb/oorb'}) - def from_oo(cls, orbit, epochs=None, location='500', scope='full', - dynmodel='N', ephfile='de430'): + @cite({"method": "2009M&PS...44.1853G", "software": "https://github.com/oorb/oorb"}) + def from_oo( + cls, + orbit, + epochs=None, + location="500", + scope="full", + dynmodel="N", + ephfile="de430", + ): """Uses pyoorb to derive ephemerides from an `~Orbit` object. For a list of output parameters, please read the `pyoorb documentation `_. @@ -746,10 +844,10 @@ def from_oo(cls, orbit, epochs=None, location='500', scope='full', 1 Ceres 2458529.1502996054 ... 1.4771304981564537e-05 70.3246138990413 1 Ceres 2458529.191966272 ... 1.448582020449618e-05 70.33407221340468 1 Ceres 2458529.2336329385 ... 1.326517587380005e-05 70.34353031031534 - 1 Ceres 2458529.2752996054 ... 1.1193369555934085e-05 70.35298818987367 """ + 1 Ceres 2458529.2752996054 ... 1.1193369555934085e-05 70.35298818987367""" if pyoorb is None: - raise RequiredPackageUnavailable('pyoorb') + raise RequiredPackageUnavailable("pyoorb") # create a copy of orbit orb = Orbit.from_table(orbit.table) @@ -761,19 +859,17 @@ def from_oo(cls, orbit, epochs=None, location='500', scope='full', timescale = epochs.scale.upper() # initialize pyoorb - if os.getenv('OORB_DATA') is None: + if os.getenv("OORB_DATA") is None: # oorb installed using conda pyoorb.pyoorb.oorb_init() else: - ephfile = os.path.join(os.getenv('OORB_DATA'), ephfile+'.dat') + ephfile = os.path.join(os.getenv("OORB_DATA"), ephfile + ".dat") pyoorb.pyoorb.oorb_init(ephfile) # identify orbit type based on available table columns orbittype = None - for testtype in ['KEP', 'COM', 'CART']: - field_names = [ - field[0] for field in Conf.oorb_orbit_fields[testtype] - ] + for testtype in ["KEP", "COM", "CART"]: + field_names = [field[0] for field in Conf.oorb_orbit_fields[testtype]] try: orb._translate_columns(field_names[1:6]) orbittype = testtype @@ -782,11 +878,10 @@ def from_oo(cls, orbit, epochs=None, location='500', scope='full', pass if orbittype is None: - raise OpenOrbError( - 'orbit type cannot be determined from elements') + raise OpenOrbError("orbit type cannot be determined from elements") # add/update orbittype column - orb['orbittype'] = [orbittype] * len(orb) + orb["orbittype"] = [orbittype] * len(orb) # derive and apply default units default_units = {} @@ -800,66 +895,68 @@ def from_oo(cls, orbit, epochs=None, location='500', scope='full', default_units[primary_field_name] = field_unit for colname in orb.field_names: - if (colname in default_units.keys() and - not isinstance(orb[colname], - (u.Quantity, u.CompositeUnit, Time))): + if colname in default_units.keys() and not isinstance( + orb[colname], (u.Quantity, u.CompositeUnit, Time) + ): orb[colname].unit = default_units[colname] # convert epochs to TT - orb['epoch'] = orb['epoch'].tt + orb["epoch"] = orb["epoch"].tt epochs = epochs.tt try: - epochs = list(zip(epochs.mjd, - [Conf.oorb_timeScales['TT']]*len(epochs))) + epochs = list(zip(epochs.mjd, [Conf.oorb_timeScales["TT"]] * len(epochs))) except TypeError: - epochs = [(epochs.mjd, Conf.oorb_timeScales['TT'])] + epochs = [(epochs.mjd, Conf.oorb_timeScales["TT"])] - if scope == 'full': + if scope == "full": oo_eph, err = pyoorb.pyoorb.oorb_ephemeris_full( - orb._to_oo(), - location, - epochs, - dynmodel) - elif scope == 'basic': + orb._to_oo(), location, epochs, dynmodel + ) + elif scope == "basic": oo_eph, err = pyoorb.pyoorb.oorb_ephemeris_basic( - orb._to_oo(), - location, - epochs, - dynmodel) + orb._to_oo(), location, epochs, dynmodel + ) if err != 0: - OpenOrbError('pyoorb failed with error code {:d}'.format(err)) + OpenOrbError("pyoorb failed with error code {:d}".format(err)) # reorder data on per-column basis and apply units - oo_eph_col = hstack([oo_eph.transpose()[:, :, i] - for i in range(oo_eph.shape[0])]).tolist() + oo_eph_col = hstack( + [oo_eph.transpose()[:, :, i] for i in range(oo_eph.shape[0])] + ).tolist() oo_eph_col_u = [] - if scope == 'full': + if scope == "full": for i, col in enumerate(oo_eph_col): - oo_eph_col_u.append(Ephem._unit_apply( - col, Conf.oorb_ephem_full_units[i])) - ephem = cls.from_columns(oo_eph_col_u, - names=Conf.oorb_ephem_full_fields) - elif scope == 'basic': + oo_eph_col_u.append( + Ephem._unit_apply(col, Conf.oorb_ephem_full_units[i]) + ) + ephem = cls.from_columns(oo_eph_col_u, names=Conf.oorb_ephem_full_fields) + elif scope == "basic": for i, col in enumerate(oo_eph_col): - oo_eph_col_u.append(Ephem._unit_apply( - col, Conf.oorb_ephem_basic_units[i])) - ephem = cls.from_columns(oo_eph_col_u, - names=Conf.oorb_ephem_basic_fields) + oo_eph_col_u.append( + Ephem._unit_apply(col, Conf.oorb_ephem_basic_units[i]) + ) + ephem = cls.from_columns(oo_eph_col_u, names=Conf.oorb_ephem_basic_fields) # add targetname column - ephem.table.add_column(Column(data=sum([[orb['targetname'][i]] * - len(epochs) for i in - range(len(orb.table))], - []), - name='targetname'), - index=0) + ephem.table.add_column( + Column( + data=sum( + [ + [orb["targetname"][i]] * len(epochs) + for i in range(len(orb.table)) + ], + [], + ), + name="targetname", + ), + index=0, + ) # convert MJD to astropy.time.TimeJulian Date - ephem.table['epoch'] = Time(ephem['MJD'], format='mjd', - scale=timescale.lower()) - ephem.table['epoch'].format = 'jd' - ephem.table.remove_column('MJD') + ephem.table["epoch"] = Time(ephem["MJD"], format="mjd", scale=timescale.lower()) + ephem.table["epoch"].format = "jd" + ephem.table.remove_column("MJD") return ephem diff --git a/sbpy/data/tests/data/TestEphemerisCLI-2p-horizons.txt b/sbpy/data/tests/data/TestEphemerisCLI-2p-horizons.txt new file mode 100644 index 000000000..84661e236 --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-2p-horizons.txt @@ -0,0 +1,294 @@ +API VERSION: 1.2 +API SOURCE: NASA/JPL Horizons API + +******************************************************************************* +JPL/HORIZONS 2P/Encke 2024-Sep-02 08:00:31 +Rec #:90000091 (+COV) Soln.date: 2024-Aug-26_10:56:41 # obs: 556 (2018-2024) + +IAU76/J2000 helio. ecliptic osc. elements (au, days, deg., period=Julian yrs): + + EPOCH= 2459752.5 ! 2022-Jun-22.0000000 (TDB) RMSW= n.a. + EC= .8485141889848308 QR= .3362300806790429 TP= 2460239.0189482248 + OM= 334.3120522286535 W= 187.0124965530834 IN= 11.50170416921873 + A= 2.219548342025076 MA= 214.9870056150526 ADIST= 4.10286660337111 + PER= 3.3067785736152 N= .298062377 ANGMOM= .013561606 + DAN= 3.93787 DDN= .33739 L= 161.1850899 + B= -1.3949208 MOID= .168138 TP= 2023-Oct-21.5189482248 + +Comet physical (GM= km^3/s^2; RAD= km): + GM= n.a. RAD= 2.4 + M1= 15.6 M2= n.a. k1= 4.5 k2= n.a. PHCOF= n.a. +Comet non-gravitational force model (AMRAT=m^2/kg;A1-A3=au/d^2;DT=days;R0=au): + AMRAT= 0. DT= 0. + A1= 3.368033468723E-10 A2= 8.548115729354E-12 A3= 0. + Standard model: + ALN= .1112620426 NK= 4.6142 NM= 2.15 NN= 5.093 R0= 2.808 + +COMET comments +1: soln ref.= JPL#K232/4, data arc: 2018-11-05 to 2024-08-18 +2: k1=4.5; +******************************************************************************* + + +******************************************************************************* +Ephemeris / API_USER Mon Sep 2 08:00:31 2024 Pasadena, USA / Horizons +******************************************************************************* +Target body name: 2P/Encke {source: JPL#K232/4} +Center body name: Earth (399) {source: DE441} +Center-site name: GEOCENTRIC +******************************************************************************* +Start time : A.D. 2024-Aug-16 00:00:00.0000 UT +Stop time : A.D. 2024-Oct-15 00:00:00.0000 UT +Step-size : 1440 minutes +******************************************************************************* +Target pole/equ : undefined +Target radii : 2.4 km +Center geodetic : 0.0, 0.0, -6378.137 {E-lon(deg),Lat(deg),Alt(km)} +Center cylindric: 0.0, 0.0, 0.0 {E-lon(deg),Dxy(km),Dz(km)} +Center pole/equ : ITRF93 {East-longitude positive} +Center radii : 6378.137, 6378.137, 6356.752 km {Equator_a, b, pole_c} +Target primary : Sun +Vis. interferer : MOON (R_eq= 1737.400) km {source: DE441} +Rel. light bend : Sun {source: DE441} +Rel. lght bnd GM: 1.3271E+11 km^3/s^2 +Small-body perts: Yes {source: SB441-N16} +Atmos refraction: NO (AIRLESS) +RA format : DEG +Time format : BOTH +Calendar mode : Mixed Julian/Gregorian +EOP file : eop.240829.p241122 +EOP coverage : DATA-BASED 1962-JAN-20 TO 2024-AUG-29. PREDICTS-> 2024-NOV-21 +Units conversion: 1 au= 149597870.700 km, c= 299792.458 km/s, 1 day= 86400.0 s +Table cut-offs 1: Elevation (-90.0deg=NO ),Airmass (>38.000=NO), Daylight (NO ) +Table cut-offs 2: Solar elongation ( 0.0,180.0=NO ),Local Hour Angle( 0.0=NO ) +Table cut-offs 3: RA/DEC angular rate ( 0.0=NO ) +Table format : Comma Separated Values (spreadsheet) +******************************************************************************* +Initial IAU76/J2000 heliocentric ecliptic osculating elements (au, days, deg.): + EPOCH= 2459752.5 ! 2022-Jun-22.0000000 (TDB) RMSW= n.a. + EC= .8485141889848308 QR= .3362300806790429 TP= 2460239.0189482248 + OM= 334.3120522286535 W= 187.0124965530834 IN= 11.50170416921873 + Equivalent ICRF heliocentric cartesian coordinates (au, au/d): + X= 3.886668467170212E+00 Y=-9.188393246574216E-01 Z=-2.098903569670719E-01 + VX=-9.846074938312395E-04 VY= 3.120416928338697E-03 VZ= 1.988497527345202E-03 +Comet physical (GM= km^3/s^2; RAD= km): + GM= n.a. RAD= 2.4 + M1= 15.6 M2= n.a. k1= 4.5 k2= n.a. PHCOF= n.a. +Comet non-gravitational force model (AMRAT=m^2/kg;A1-A3=au/d^2;DT=days;R0=au): + AMRAT= 0. DT= 0. + A1= 3.368033468723E-10 A2= 8.548115729354E-12 A3= 0. + Standard model: + ALN= .1112620426 NK= 4.6142 NM= 2.15 NN= 5.093 R0= 2.808 +******************************************************************************************************************************************************************************************************************************************** + Date__(UT)__HR:MN, Date_________JDUT, , , R.A._(ICRF), DEC_(ICRF), dRA*cosD, d(DEC)/dt, T-mag, N-mag, r, rdot, delta, deldot, S-O-T,/r, S-T-O, PsAng, PsAMV, GlxLon, GlxLat, +******************************************************************************************************************************************************************************************************************************************** +$$SOE + 2024-Aug-16 00:00, 2460538.500000000, , , 329.30474, -15.12338, -45.9107, -13.5742, 19.827, n.a., 3.362689434345, 8.8625509, 2.35165022935209, 7.7914783, 176.1737,/L, 1.1540, 290.572, 237.575, 40.373744,-47.873291, + 2024-Aug-17 00:00, 2460539.500000000, , , 328.98787, -15.21298, -45.7974, -13.4254, 19.834, n.a., 3.367797212379, 8.8257359, 2.35631688615317, 8.3691715, 177.0443,/L, 0.8899, 306.748, 237.776, 40.067975,-47.629799, + 2024-Aug-18 00:00, 2460540.500000000, , , 328.67178, -15.30155, -45.6512, -13.2664, 19.842, n.a., 3.372883747218, 8.7890194, 2.36131704455975, 8.9463448, 177.5296,/L, 0.7414, 332.189, 237.975, 39.766426,-47.386178, + 2024-Aug-19 00:00, 2460541.500000000, , , 328.35667, -15.38902, -45.4728, -13.0977, 19.850, n.a., 3.377949095643, 8.7524007, 2.36665032393581, 9.5227309, 177.4051,/T, 0.7757, 1.325, 238.170, 39.469241,-47.142623, + 2024-Aug-20 00:00, 2460542.500000000, , , 328.04277, -15.47533, -45.2628, -12.9198, 19.858, n.a., 3.382993313982, 8.7158789, 2.37231620908038, 10.0981297, 176.7412,/T, 0.9716, 22.924, 238.363, 39.176553,-46.899328, + 2024-Aug-21 00:00, 2460543.500000000, , , 327.73031, -15.56043, -45.0217, -12.7332, 19.866, n.a., 3.388016458114, 8.6794533, 2.37831408460100, 10.6723897, 175.7868,/T, 1.2539, 36.043, 238.552, 38.888486,-46.656486, + 2024-Aug-22 00:00, 2460544.500000000, , , 327.41949, -15.64425, -44.7499, -12.5382, 19.875, n.a., 3.393018583468, 8.6431232, 2.38464325373495, 11.2453738, 174.6977,/T, 1.5752, 44.092, 238.738, 38.605153,-46.414286, + 2024-Aug-23 00:00, 2460545.500000000, , , 327.11054, -15.72675, -44.4476, -12.3351, 19.884, n.a., 3.397999745033, 8.6068877, 2.39130293467752, 11.8169183, 173.5427,/T, 1.9145, 49.363, 238.921, 38.326659,-46.172919, + 2024-Aug-24 00:00, 2460546.500000000, , , 326.80367, -15.80786, -44.1151, -12.1243, 19.893, n.a., 3.402959997362, 8.5707461, 2.39829223520000, 12.3868005, 172.3527,/T, 2.2624, 53.039, 239.101, 38.053104,-45.932573, + 2024-Aug-25 00:00, 2460547.500000000, , , 326.49908, -15.88755, -43.7527, -11.9058, 19.902, n.a., 3.407899394576, 8.5346977, 2.40561011380958, 12.9547255, 171.1428,/T, 2.6141, 55.740, 239.278, 37.784581,-45.693437, + 2024-Aug-26 00:00, 2460548.500000000, , , 326.19699, -15.96575, -43.3605, -11.6801, 19.912, n.a., 3.412817990372, 8.4987416, 2.41325533858522, 13.5203308, 169.9211,/T, 2.9670, 57.807, 239.452, 37.521178,-45.455699, + 2024-Aug-27 00:00, 2460549.500000000, , , 325.89761, -16.04242, -42.9389, -11.4474, 19.922, n.a., 3.417715838026, 8.4628773, 2.42122645250697, 14.0832026, 168.6925,/T, 3.3195, 59.443, 239.623, 37.262978,-45.219546, + 2024-Aug-28 00:00, 2460550.500000000, , , 325.60114, -16.11753, -42.4885, -11.2081, 19.932, n.a., 3.422592990400, 8.4271039, 2.42952174934040, 14.6428948, 167.4599,/T, 3.6706, 60.773, 239.790, 37.010058,-44.985161, + 2024-Aug-29 00:00, 2460551.500000000, , , 325.30778, -16.19101, -42.0098, -10.9625, 19.943, n.a., 3.427449499947, 8.3914207, 2.43813926003474, 15.1989457, 166.2253,/T, 4.0196, 61.878, 239.954, 36.762487,-44.752727, + 2024-Aug-30 00:00, 2460552.500000000, , , 325.01774, -16.26285, -41.5036, -10.7109, 19.953, n.a., 3.432285418714, 8.3558271, 2.44707674750675, 15.7508906, 164.9902,/T, 4.3660, 62.813, 240.115, 36.520330,-44.522421, + 2024-Aug-31 00:00, 2460553.500000000, , , 324.73121, -16.33298, -40.9707, -10.4539, 19.964, n.a., 3.437100798347, 8.3203222, 2.45633170736155, 16.2982697, 163.7555,/T, 4.7091, 63.617, 240.272, 36.283644,-44.294419, + 2024-Sep-01 00:00, 2460554.500000000, , , 324.44837, -16.40139, -40.4119, -10.1918, 19.975, n.a., 3.441895690095, 8.2849056, 2.46590137268890, 16.8406345, 162.5221,/T, 5.0489, 64.317, 240.426, 36.052478,-44.068889, + 2024-Sep-02 00:00, 2460555.500000000, , , 324.16941, -16.46804, -39.8283, -9.92505, 19.987, n.a., 3.446670144816, 8.2495763, 2.47578272183622, 17.3775520, 161.2906,/T, 5.3848, 64.934, 240.577, 35.826876,-43.845996, + 2024-Sep-03 00:00, 2460556.500000000, , , 323.89449, -16.53290, -39.2209, -9.65417, 19.998, n.a., 3.451424212974, 8.2143338, 2.48597248859938, 17.9086090, 160.0614,/T, 5.7167, 65.483, 240.725, 35.606872,-43.625899, + 2024-Sep-04 00:00, 2460557.500000000, , , 323.62379, -16.59595, -38.5908, -9.37960, 20.010, n.a., 3.456157944652, 8.1791775, 2.49646717448255, 18.4334148, 158.8351,/T, 6.0443, 65.976, 240.869, 35.392496,-43.408751, + 2024-Sep-05 00:00, 2460558.500000000, , , 323.35746, -16.65716, -37.9391, -9.10179, 20.022, n.a., 3.460871389545, 8.1441065, 2.50726306262885, 18.9516041, 157.6119,/T, 6.3674, 66.422, 241.010, 35.183770,-43.194698, + 2024-Sep-06 00:00, 2460559.500000000, , , 323.09567, -16.71651, -37.2670, -8.82122, 20.035, n.a., 3.465564596972, 8.1091203, 2.51835623286801, 19.4628384, 156.3921,/T, 6.6859, 66.828, 241.148, 34.980708,-42.983881, + 2024-Sep-07 00:00, 2460560.500000000, , , 322.83853, -16.77399, -36.5757, -8.53834, 20.047, n.a., 3.470237615869, 8.0742182, 2.52974257722419, 19.9668059, 155.1761,/T, 6.9995, 67.200, 241.282, 34.783319,-42.776431, + 2024-Sep-08 00:00, 2460561.500000000, , , 322.58621, -16.82958, -35.8665, -8.25361, 20.060, n.a., 3.474890494800, 8.0393995, 2.54141781528446, 20.4632210, 153.9640,/T, 7.3082, 67.542, 241.413, 34.591605,-42.572476, + 2024-Sep-09 00:00, 2460562.500000000, , , 322.33881, -16.88327, -35.1405, -7.96746, 20.072, n.a., 3.479523281952, 8.0046637, 2.55337750908490, 20.9518236, 152.7561,/T, 7.6117, 67.859, 241.540, 34.405561,-42.372133, + 2024-Sep-10 00:00, 2460563.500000000, , , 322.09646, -16.93506, -34.3990, -7.68032, 20.085, n.a., 3.484136025138, 7.9700099, 2.56561707760337, 21.4323785, 151.5525,/T, 7.9100, 68.153, 241.664, 34.225178,-42.175515, + 2024-Sep-11 00:00, 2460564.500000000, , , 321.85927, -16.98493, -33.6432, -7.39261, 20.099, n.a., 3.488728771801, 7.9354377, 2.57813181150516, 21.9046770, 150.3533,/T, 8.2029, 68.427, 241.785, 34.050440,-41.982725, + 2024-Sep-12 00:00, 2460565.500000000, , , 321.62733, -17.03290, -32.8744, -7.10473, 20.112, n.a., 3.493301569013, 7.9009463, 2.59091688939235, 22.3685407, 149.1588,/T, 8.4904, 68.683, 241.902, 33.881326,-41.793859, + 2024-Sep-13 00:00, 2460566.500000000, , , 321.40074, -17.07895, -32.0938, -6.81710, 20.125, n.a., 3.497854463477, 7.8665352, 2.60396739733102, 22.8238290, 147.9691,/T, 8.7723, 68.924, 242.016, 33.717809,-41.609008, + 2024-Sep-14 00:00, 2460567.500000000, , , 321.17958, -17.12309, -31.3027, -6.53008, 20.139, n.a., 3.502387501528, 7.8322036, 2.61727835361104, 23.2704498, 146.7841,/T, 9.0486, 69.150, 242.127, 33.559859,-41.428253, + 2024-Sep-15 00:00, 2460568.500000000, , , 320.96392, -17.16533, -30.5023, -6.24406, 20.153, n.a., 3.506900729137, 7.7979510, 2.63084474006346, 23.7083707, 145.6041,/T, 9.3192, 69.363, 242.234, 33.407438,-41.251668, + 2024-Sep-16 00:00, 2460569.500000000, , , 320.75383, -17.20568, -29.6937, -5.95936, 20.167, n.a., 3.511394191910, 7.7637767, 2.64466153923797, 24.1376273, 144.4291,/T, 9.5841, 69.564, 242.338, 33.260508,-41.079319, + 2024-Sep-17 00:00, 2460570.500000000, , , 320.54937, -17.24414, -28.8781, -5.67626, 20.181, n.a., 3.515867935091, 7.7296802, 2.65872377309445, 24.5583186, 143.2591,/T, 9.8431, 69.755, 242.438, 33.119023,-40.911268, + 2024-Sep-18 00:00, 2460571.500000000, , , 320.35059, -17.28073, -28.0562, -5.39500, 20.195, n.a., 3.520322003567, 7.6956607, 2.67302653461664, 24.9705862, 142.0941,/T, 10.0963, 69.935, 242.536, 32.982937,-40.747568, + 2024-Sep-19 00:00, 2460572.500000000, , , 320.15753, -17.31545, -27.2288, -5.11570, 20.209, n.a., 3.524756441867, 7.6617178, 2.68756500170856, 25.3745752, 140.9343,/T, 10.3437, 70.106, 242.630, 32.852201,-40.588269, + 2024-Sep-20 00:00, 2460573.500000000, , , 319.97023, -17.34833, -26.3964, -4.83843, 20.223, n.a., 3.529171294169, 7.6278507, 2.70233442604059, 25.7703888, 139.7795,/T, 10.5851, 70.269, 242.720, 32.726768,-40.433415, + 2024-Sep-21 00:00, 2460574.500000000, , , 319.78873, -17.37937, -25.5593, -4.56322, 20.238, n.a., 3.533566604305, 7.5940590, 2.71733009824313, 26.1580548, 138.6299,/T, 10.8205, 70.424, 242.808, 32.606588,-40.283048, + 2024-Sep-22 00:00, 2460575.500000000, , , 319.61307, -17.40860, -24.7181, -4.29006, 20.252, n.a., 3.537942415765, 7.5603419, 2.73254730016307, 26.5375157, 137.4854,/T, 11.0500, 70.571, 242.892, 32.491615,-40.137207, + 2024-Sep-23 00:00, 2460576.500000000, , , 319.44328, -17.43602, -23.8731, -4.01895, 20.267, n.a., 3.542298771701, 7.5266990, 2.74798125861797, 26.9086423, 136.3460,/T, 11.2734, 70.711, 242.973, 32.381804,-39.995928, + 2024-Sep-24 00:00, 2460577.500000000, , , 319.27938, -17.46165, -23.0247, -3.74992, 20.282, n.a., 3.546635714933, 7.4931297, 2.76362711148269, 27.2712604, 135.2119,/T, 11.4909, 70.845, 243.051, 32.277108,-39.859242, + 2024-Sep-25 00:00, 2460578.500000000, , , 319.12140, -17.48550, -22.1737, -3.48301, 20.296, n.a., 3.550953287952, 7.4596334, 2.77947988976684, 27.6251783, 134.0829,/T, 11.7023, 70.972, 243.126, 32.177483,-39.727178, + 2024-Sep-26 00:00, 2460579.500000000, , , 318.96936, -17.50759, -21.3205, -3.21827, 20.311, n.a., 3.555251532925, 7.4262095, 2.79553451343772, 27.9702058, 132.9591,/T, 11.9076, 71.093, 243.197, 32.082884,-39.599762, + 2024-Sep-27 00:00, 2460580.500000000, , , 318.82327, -17.52793, -20.4659, -2.95580, 20.326, n.a., 3.559530491695, 7.3928575, 2.81178579623575, 28.3061663, 131.8406,/T, 12.1069, 71.209, 243.266, 31.993266,-39.477013, + 2024-Sep-28 00:00, 2460581.500000000, , , 318.68315, -17.54654, -19.6105, -2.69566, 20.341, n.a., 3.563790205787, 7.3595768, 2.82822845498139, 28.6329014, 130.7273,/T, 12.3002, 71.319, 243.331, 31.908581,-39.358948, + 2024-Sep-29 00:00, 2460582.500000000, , , 318.54900, -17.56343, -18.7551, -2.43796, 20.356, n.a., 3.568030716411, 7.3263669, 2.84485712033785, 28.9502711, 129.6192,/T, 12.4873, 71.424, 243.393, 31.828784,-39.245578, + 2024-Sep-30 00:00, 2460583.500000000, , , 318.42082, -17.57861, -17.9002, -2.18277, 20.371, n.a., 3.572252064459, 7.2932272, 2.86166634752825, 29.2581538, 128.5164,/T, 12.6684, 71.525, 243.452, 31.753825,-39.136910, + 2024-Oct-01 00:00, 2460584.500000000, , , 318.29861, -17.59212, -17.0467, -1.93018, 20.386, n.a., 3.576454290517, 7.2601573, 2.87865062660956, 29.5564455, 127.4189,/T, 12.8435, 71.620, 243.508, 31.683656,-39.032948, + 2024-Oct-02 00:00, 2460585.500000000, , , 318.18236, -17.60396, -16.1951, -1.68027, 20.402, n.a., 3.580637434858, 7.2271565, 2.89580439250246, 29.8450594, 126.3266,/T, 13.0125, 71.711, 243.562, 31.618228,-38.933688, + 2024-Oct-03 00:00, 2460586.500000000, , , 318.07205, -17.61416, -15.3461, -1.43313, 20.417, n.a., 3.584801537451, 7.1942244, 2.91312203515841, 30.1239267, 125.2396,/T, 13.1755, 71.798, 243.612, 31.557488,-38.839127, + 2024-Oct-04 00:00, 2460587.500000000, , , 317.96766, -17.62272, -14.5004, -1.18883, 20.432, n.a., 3.588946637960, 7.1613604, 2.93059791015383, 30.3929972, 124.1579,/T, 13.3326, 71.881, 243.659, 31.501386,-38.749253, + 2024-Oct-05 00:00, 2460588.500000000, , , 317.86917, -17.62968, -13.6587, -0.94745, 20.447, n.a., 3.593072775750, 7.1285640, 2.94822634978532, 30.6522407, 123.0814,/T, 13.4836, 71.959, 243.704, 31.449869,-38.664052, + 2024-Oct-06 00:00, 2460589.500000000, , , 317.77655, -17.63504, -12.8214, -0.70905, 20.463, n.a., 3.597179989884, 7.0958347, 2.96600167454008, 30.9016466, 122.0101,/T, 13.6287, 72.034, 243.745, 31.402884,-38.583507, + 2024-Oct-07 00:00, 2460590.500000000, , , 317.68978, -17.63884, -11.9892, -0.47371, 20.478, n.a., 3.601268319129, 7.0631720, 2.98391820474236, 31.1412249, 120.9441,/T, 13.7680, 72.105, 243.784, 31.360376,-38.507595, + 2024-Oct-08 00:00, 2460591.500000000, , , 317.60880, -17.64109, -11.1627, -0.24149, 20.493, n.a., 3.605337801956, 7.0305754, 3.00197027229463, 31.3710058, 119.8832,/T, 13.9013, 72.172, 243.820, 31.322290,-38.436292, + 2024-Oct-09 00:00, 2460592.500000000, , , 317.53359, -17.64181, -10.3425, -0.01243, 20.509, n.a., 3.609388476546, 6.9980444, 3.02015223273410, 31.5910412, 118.8275,/T, 14.0289, 72.236, 243.853, 31.288571,-38.369567, + 2024-Oct-10 00:00, 2460593.500000000, , , 317.46410, -17.64102, -9.52890, 0.213405, 20.524, n.a., 3.613420380787, 6.9655784, 3.03845847824449, 31.8014068, 117.7768,/T, 14.1507, 72.296, 243.883, 31.259162,-38.307390, + 2024-Oct-11 00:00, 2460594.500000000, , , 317.40027, -17.63875, -8.72255, 0.435966, 20.539, n.a., 3.617433552280, 6.9331771, 3.05688345266062, 32.0022066, 116.7313,/T, 14.2668, 72.353, 243.911, 31.234007,-38.249723, + 2024-Oct-12 00:00, 2460595.500000000, , , 317.34207, -17.63501, -7.92389, 0.655211, 20.554, n.a., 3.621428028337, 6.9008398, 3.07542166965788, 32.1935796, 115.6908,/T, 14.3772, 72.407, 243.936, 31.213048,-38.196529, + 2024-Oct-13 00:00, 2460596.500000000, , , 317.28942, -17.62983, -7.13332, 0.871104, 20.570, n.a., 3.625403845987, 6.8685662, 3.09406773488878, 32.3757064, 114.6552,/T, 14.4820, 72.458, 243.959, 31.196226,-38.147766, + 2024-Oct-14 00:00, 2460597.500000000, , , 317.24228, -17.62324, -6.35119, 1.083629, 20.585, n.a., 3.629361041977, 6.8363557, 3.11281637140458, 32.5488131, 113.6246,/T, 14.5812, 72.505, 243.979, 31.183485,-38.103390, + 2024-Oct-15 00:00, 2460598.500000000, , , 317.20058, -17.61525, -5.57775, 1.292795, 20.600, n.a., 3.633299652770, 6.8042079, 3.13166244504940, 32.7131663, 112.5988,/T, 14.6750, 72.550, 243.997, 31.174766,-38.063355, +$$EOE +******************************************************************************************************************************************************************************************************************************************** +Column meaning: + +TIME + + Times PRIOR to 1962 are UT1, a mean-solar time closely related to the +prior but now-deprecated GMT. Times AFTER 1962 are in UTC, the current +civil or "wall-clock" time-scale. UTC is kept within 0.9 seconds of UT1 +using integer leap-seconds for 1972 and later years. + + Conversion from the internal Barycentric Dynamical Time (TDB) of solar +system dynamics to the non-uniform civil UT time-scale requested for output +has not been determined for UTC times after the next July or January 1st. +Therefore, the last known leap-second is used as a constant over future +intervals. + + Time tags refer to the UT time-scale conversion from TDB on Earth +regardless of observer location within the solar system, although clock +rates may differ due to the local gravity field and no analog to "UT" +may be defined for that location. + + Any 'b' symbol in the 1st-column denotes a B.C. date. First-column blank +(" ") denotes an A.D. date. + +CALENDAR SYSTEM + + Mixed calendar mode was active such that calendar dates after AD 1582-Oct-15 +(if any) are in the modern Gregorian system. Dates prior to 1582-Oct-5 (if any) +are in the Julian calendar system, which is automatically extended for dates +prior to its adoption on 45-Jan-1 BC. The Julian calendar is useful for +matching historical dates. The Gregorian calendar more accurately corresponds +to the Earth's orbital motion and seasons. A "Gregorian-only" calendar mode is +available if such physical events are the primary interest. + + NOTE: "n.a." in output means quantity "not available" at the print-time. + + 'R.A._(ICRF), DEC_(ICRF),' = + Astrometric right ascension and declination of the target center with +respect to the observing site (coordinate origin) in the reference frame of +the planetary ephemeris (ICRF). Compensated for down-leg light-time delay +aberration. + + Units: RA in decimal degrees, ddd.fffff{ffff} + DEC in decimal degrees, sdd.fffff{ffff} + + 'dRA*cosD, d(DEC)/dt,' = + The angular rate of change in aparent RA and DEC of the target. This is +with respect to the non-inertial IAU76/80 Earth true equator and equinox +of-date reference frame. d(RA)/dt is multiplied by the cosine of declination +to provide a linear rate in the plane-of-sky. Units: ARCSECONDS PER HOUR + + 'T-mag, N-mag,' = + Comets' apparent visual total magnitude ("T-mag") and nuclear magnitude +("N-mag") using the standard IAU model: + + T-mag= M1 + 5*log10(delta) + k1*log10(r) + N-mag= M2 + 5*log10(delta) + k2*log10(r) + phcof*beta + + Units: MAGNITUDES + + 'r, rdot,' = + The Sun's apparent range ("r", light-time aberrated) and range-rate ("rdot") +relative to the target center, as seen by the observer. A positive "rdot" means +the target center was moving away from the Sun, negative means moving toward +the Sun. Units: AU and KM/S + + 'delta, deldot,' = + Apparent range ("delta", light-time aberrated) and range-rate ("delta-dot") +of the target center relative to the observer. A positive "deldot" means the +target center is moving away from the observer, negative indicates movement +toward the observer. Units: AU and KM/S + + 'S-O-T,/r,' = + Sun-Observer-Target apparent SOLAR ELONGATION ANGLE seen from the observers' +location at print-time. + + The '/r' column provides a code indicating the targets' apparent position +relative to the Sun in the observers' sky, as described below: + + Case A: For an observing location on the surface of a rotating body, that +body rotational sense is considered: + + /T indicates target TRAILS Sun (evening sky: rises and sets AFTER Sun) + /L indicates target LEADS Sun (morning sky: rises and sets BEFORE Sun) + + Case B: For an observing point that does not have a rotational model (such +as a spacecraft), the "leading" and "trailing" condition is defined by the +observers' heliocentric ORBITAL motion: + + * If continuing in the observers' current direction of heliocentric + motion would encounter the targets' apparent longitude first, followed + by the Sun's, the target LEADS the Sun as seen by the observer. + + * If the Sun's apparent longitude would be encountered first, followed + by the targets', the target TRAILS the Sun. + + Two other codes can be output: + /* indicates observer is Sun-centered (undefined) + /? Target is aligned with Sun center (no lead or trail) + + The S-O-T solar elongation angle is numerically the minimum separation +angle of the Sun and target in the sky in any direction. It does NOT indicate +the amount of separation in the leading or trailing directions, which would +be defined along the equator of a spherical coordinate system. + + Units: DEGREES + + 'S-T-O,' = + The Sun-Target-Observer angle; the interior vertex angle at target center +formed by a vector from the target to the apparent center of the Sun (at +reflection time on the target) and the apparent vector from target to the +observer at print-time. Slightly different from true PHASE ANGLE (requestable +separately) at the few arcsecond level in that it includes stellar aberration +on the down-leg from target to observer. Units: DEGREES + + 'PsAng, PsAMV,' = + The position angles of the extended Sun-to-target radius vector ("PsAng") +and the negative of the targets' heliocentric velocity vector ("PsAMV"), as +seen in the observers' plane-of-sky, measured counter-clockwise (east) from +reference-frame north-pole. Primarily intended for ACTIVE COMETS, "PsAng" +is an indicator of the comets' gas-tail orientation in the sky (being in the +anti-sunward direction) while "PsAMV" is an indicator of dust-tail orientation. +Units: DEGREES + + 'GlxLon, GlxLat,' = + Observer-centered Galactic System II (post WW II) longitude and latitude +of the target centers' apparent position, with light-time, gravitational +deflection of light, and stellar aberrations. Units: DEGREES + +Computations by ... + + Solar System Dynamics Group, Horizons On-Line Ephemeris System + 4800 Oak Grove Drive, Jet Propulsion Laboratory + Pasadena, CA 91109 USA + + General site: https://ssd.jpl.nasa.gov/ + Mailing list: https://ssd.jpl.nasa.gov/email_list.html + System news : https://ssd.jpl.nasa.gov/horizons/news.html + User Guide : https://ssd.jpl.nasa.gov/horizons/manual.html + Connect : browser https://ssd.jpl.nasa.gov/horizons/app.html#/x + API https://ssd-api.jpl.nasa.gov/doc/horizons.html + command-line telnet ssd.jpl.nasa.gov 6775 + e-mail/batch https://ssd.jpl.nasa.gov/ftp/ssd/hrzn_batch.txt + scripts https://ssd.jpl.nasa.gov/ftp/ssd/SCRIPTS + Author : Jon.D.Giorgini@jpl.nasa.gov + +******************************************************************************************************************************************************************************************************************************************** diff --git a/sbpy/data/tests/data/TestEphemerisCLI-2p-miriade.txt b/sbpy/data/tests/data/TestEphemerisCLI-2p-miriade.txt new file mode 100644 index 000000000..88da4d40f --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-2p-miriade.txt @@ -0,0 +1,123 @@ + + + + + + +Miriade - Positional ephemeris + +Positional ephemeris + + + + + + + + + + + + + +Name or designation of the target + + +Julian day (UTC) + + +Right ascension of target centre + + +Declination of target centre + + +Distance from observer + + +Apparent visual magnitude + + +Phase angle + + +Solar elongation + + +Rate of change of target center RA (multiplied by the cosine of the declination) + + +Rate of change of target center DE + + +Rate of change of target distance from observer + + + +2P 2460538.5000000000021 57 13.24921-15 07 24.0954 2.351652518 19.85 1.15176.18-0.76545E+00-0.22522E+00 7.7912170 +2P 2460539.5000000000021 55 57.20239-15 12 46.6437 2.356319023 19.86 0.89177.05-0.76356E+00-0.22273E+00 8.3689102 +2P 2460540.5000000000021 54 41.33882-15 18 5.4862 2.361319030 19.86 0.74177.53-0.76112E+00-0.22007E+00 8.9460836 +2P 2460541.5000000000021 53 25.71289-15 23 20.3861 2.366652158 19.87 0.78177.40-0.75815E+00-0.21725E+00 9.5224699 +2P 2460542.5000000000021 52 10.37805-15 28 31.1187 2.372317893 19.88 0.98176.73-0.75464E+00-0.21429E+00 10.0978690 +2P 2460543.5000000000021 50 55.38688-15 33 37.4701 2.378315617 19.89 1.26175.78-0.75062E+00-0.21117E+00 10.6721296 +2P 2460544.5000000000021 49 40.79134-15 38 39.2364 2.384644636 19.90 1.58174.69-0.74609E+00-0.20792E+00 11.2451144 +2P 2460545.5000000000021 48 26.64285-15 43 36.2216 2.391304167 19.91 1.92173.53-0.74105E+00-0.20453E+00 11.8166596 +2P 2460546.5000000000021 47 12.99256-15 48 28.2368 2.398293318 19.92 2.27172.34-0.73550E+00-0.20101E+00 12.3865427 +2P 2460547.5000000000021 45 59.89134-15 53 15.0990 2.405611047 19.93 2.62171.13-0.72945E+00-0.19737E+00 12.9544687 +2P 2460548.5000000000021 44 47.38972-15 57 56.6318 2.413256124 19.94 2.97169.91-0.72291E+00-0.19361E+00 13.5200752 +2P 2460549.5000000000021 43 35.53778-16 02 32.6654 2.421227090 19.95 3.33168.68-0.71588E+00-0.18974E+00 14.0829483 +2P 2460550.5000000000021 42 24.38497-16 07 3.0371 2.429522240 19.96 3.68167.45-0.70837E+00-0.18575E+00 14.6426419 +2P 2460551.5000000000021 41 13.97986-16 11 27.5930 2.438139604 19.97 4.03166.21-0.70039E+00-0.18166E+00 15.1986944 +2P 2460552.5000000000021 40 4.36999-16 15 46.1879 2.447076947 19.98 4.37164.98-0.69195E+00-0.17747E+00 15.7506409 +2P 2460553.5000000000021 38 55.60170-16 19 58.6863 2.456331763 19.99 4.72163.74-0.68306E+00-0.17319E+00 16.2980217 +2P 2460554.5000000000021 37 47.71991-16 24 4.9628 2.465901285 20.00 5.06162.51-0.67374E+00-0.16883E+00 16.8403884 +2P 2460555.5000000000021 36 40.76807-16 28 4.9027 2.475782492 20.01 5.39161.28-0.66401E+00-0.16439E+00 17.3773079 +2P 2460556.5000000000021 35 34.78795-16 31 58.4023 2.485972118 20.03 5.72160.05-0.65388E+00-0.15989E+00 17.9083669 +2P 2460557.5000000000021 34 29.81959-16 35 45.3685 2.496466664 20.04 6.05158.82-0.64337E+00-0.15532E+00 18.4331748 +2P 2460558.5000000000021 33 25.90116-16 39 25.7199 2.507262414 20.05 6.37157.60-0.63250E+00-0.15070E+00 18.9513665 +2P 2460559.5000000000021 32 23.06895-16 42 59.3858 2.518355447 20.06 6.69156.38-0.62129E+00-0.14604E+00 19.4626031 +2P 2460560.5000000000021 31 21.35726-16 46 26.3067 2.529741656 20.08 7.01155.17-0.60977E+00-0.14133E+00 19.9665730 +2P 2460561.5000000000021 30 20.79839-16 49 46.4339 2.541416760 20.09 7.32153.95-0.59794E+00-0.13660E+00 20.4629906 +2P 2460562.5000000000021 29 21.42260-16 52 59.7294 2.553376321 20.10 7.62152.75-0.58583E+00-0.13185E+00 20.9515957 +2P 2460563.5000000000021 28 23.25810-16 56 6.1654 2.565615758 20.11 7.92151.54-0.57347E+00-0.12707E+00 21.4321532 +2P 2460564.5000000000021 27 26.33102-16 59 5.7242 2.578130362 20.13 8.21150.34-0.56086E+00-0.12229E+00 21.9044544 +2P 2460565.5000000000021 26 30.66540-17 01 58.3980 2.590915312 20.14 8.50149.15-0.54804E+00-0.11751E+00 22.3683207 +2P 2460566.5000000000021 25 36.28317-17 04 44.1886 2.603965693 20.16 8.78147.96-0.53503E+00-0.11274E+00 22.8236118 +2P 2460567.5000000000021 24 43.20413-17 07 23.1072 2.617276525 20.17 9.06146.77-0.52183E+00-0.10797E+00 23.2702355 +2P 2460568.5000000000021 23 51.44601-17 09 55.1739 2.630842788 20.18 9.33145.59-0.50849E+00-0.10322E+00 23.7081593 +2P 2460569.5000000000021 23 1.02449-17 12 20.4176 2.644659465 20.20 9.59144.42-0.49501E+00-0.98493E-01 24.1374188 +2P 2460570.5000000000021 22 11.95339-17 14 38.8743 2.658721579 20.21 9.85143.25-0.48140E+00-0.93794E-01 24.5581131 +2P 2460571.5000000000021 21 24.24489-17 16 50.5863 2.673024223 20.23 10.10142.08-0.46770E+00-0.89126E-01 24.9703837 +2P 2460572.5000000000021 20 37.90983-17 18 55.5999 2.687562573 20.24 10.35140.92-0.45390E+00-0.84491E-01 25.3743757 +2P 2460573.5000000000021 19 52.95799-17 20 53.9641 2.702331883 20.26 10.59139.77-0.44002E+00-0.79890E-01 25.7701923 +2P 2460574.5000000000021 19 9.39836-17 22 45.7286 2.717327442 20.27 10.83138.62-0.42606E+00-0.75325E-01 26.1578614 +2P 2460575.5000000000021 18 27.23922-17 24 30.9433 2.732544533 20.28 11.06137.48-0.41204E+00-0.70793E-01 26.5373254 +2P 2460576.5000000000021 17 46.48812-17 26 9.6579 2.747978382 20.30 11.28136.34-0.39795E+00-0.66296E-01 26.9084551 +2P 2460577.5000000000021 17 7.15185-17 27 41.9222 2.763624127 20.31 11.50135.20-0.38380E+00-0.61835E-01 27.2710764 +2P 2460578.5000000000021 16 29.23622-17 29 7.7870 2.779476800 20.33 11.71134.07-0.36961E+00-0.57408E-01 27.6249973 +2P 2460579.5000000000021 15 52.74605-17 30 27.3043 2.795531320 20.34 11.91132.95-0.35539E+00-0.53019E-01 27.9700279 +2P 2460580.5000000000021 15 17.68502-17 31 40.5276 2.811782500 20.36 12.11131.83-0.34114E+00-0.48667E-01 28.3059916 +2P 2460581.5000000000021 14 44.05561-17 32 47.5126 2.828225059 20.38 12.31130.72-0.32687E+00-0.44354E-01 28.6327298 +2P 2460582.5000000000021 14 11.85916-17 33 48.3169 2.844853625 20.39 12.49129.61-0.31261E+00-0.40083E-01 28.9501026 +2P 2460583.5000000000021 13 41.09576-17 34 43.0000 2.861662756 20.41 12.67128.51-0.29836E+00-0.35853E-01 29.2579885 +2P 2460584.5000000000021 13 11.76434-17 35 31.6237 2.878646940 20.42 12.85127.41-0.28412E+00-0.31667E-01 29.5562833 +2P 2460585.5000000000021 12 43.86264-17 36 14.2516 2.895800613 20.44 13.02126.32-0.26992E+00-0.27526E-01 29.8449003 +2P 2460586.5000000000021 12 17.38724-17 36 50.9492 2.913118164 20.45 13.18125.23-0.25577E+00-0.23431E-01 30.1237707 +2P 2460587.5000000000021 11 52.33358-17 37 21.7840 2.930593949 20.47 13.34124.15-0.24167E+00-0.19383E-01 30.3928443 +2P 2460588.5000000000021 11 28.69597-17 37 46.8251 2.948222301 20.48 13.49123.07-0.22763E+00-0.15384E-01 30.6520908 +2P 2460589.5000000000021 11 6.46766-17 38 6.1434 2.965997540 20.50 13.63122.00-0.21367E+00-0.11435E-01 30.9014998 +2P 2460590.5000000000021 10 45.64084-17 38 19.8113 2.983913986 20.51 13.77120.94-0.19980E+00-0.75365E-02 31.1410811 +2P 2460591.5000000000021 10 26.20670-17 38 27.9029 3.001965971 20.53 13.91119.88-0.18602E+00-0.36900E-02 31.3708650 +2P 2460592.5000000000021 10 8.15546-17 38 30.4933 3.020147851 20.55 14.03118.82-0.17234E+00 0.10375E-03 31.5909034 +2P 2460593.5000000000021 09 51.47641-17 38 27.6592 3.038454017 20.56 14.16117.77-0.15878E+00 0.38439E-02 31.8012719 +2P 2460594.5000000000021 09 36.15795-17 38 19.4783 3.056878914 20.58 14.27116.72-0.14533E+00 0.75297E-02 32.0020747 +2P 2460595.5000000000021 09 22.18766-17 38 6.0294 3.075417055 20.59 14.38115.68-0.13202E+00 0.11160E-01 32.1934506 +2P 2460596.5000000000021 09 9.55234-17 37 47.3925 3.094063047 20.61 14.49114.65-0.11883E+00 0.14735E-01 32.3755803 +2P 2460597.5000000000021 08 58.23811-17 37 23.6479 3.112811611 20.62 14.59113.62-0.10579E+00 0.18254E-01 32.5486899 +2P 2460598.5000000000021 08 48.23056-17 36 54.8761 3.131657614 20.64 14.68112.59-0.92899E-01 0.21717E-01 32.7130459 + + + + + + diff --git a/sbpy/data/tests/data/TestEphemerisCLI-2p-mpc.txt b/sbpy/data/tests/data/TestEphemerisCLI-2p-mpc.txt new file mode 100644 index 000000000..f645a1d47 --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-2p-mpc.txt @@ -0,0 +1,111 @@ + + + + + Minor Planet Ephemeris Service: Query Results + + +

Minor Planet Ephemeris Service: Query Results

+ Below are the results of your request from the Minor Planet Center's + Minor Planet Ephemeris Service. +

+ Newly designated objects may take up to 1 day to show up in this service. +

+

+ Orbits and ephemerides of unnumbered NEOs are up-to-date. Other orbits are in + the process of being refreshed. +

+

+ The current system is not completely reliable in the case of objects + with very-close approaches with the Earth. +

We are working on a completely new system, + but for the time being we encourage the users to double check + the results with other ephemeris generators when the object is very close to + Earth. +

+ Ephemerides are for + the geocenter. +

+ 2P/Encke +

Perturbed ephemeris below is based on + unperturbed elements from +MPEC 2024-R10. +

+ 

+0002P
+Date       UT      R.A. (J2000) Decl.    Delta     r     El.    Ph.   m1     Sky Motion
+            h m s                                                            "/hr     P.A.
+2024 08 16 000000 21 57 13.1 -15 07 25   2.352   3.363  176.2   1.1  21.2   47.87    253.6
+2024 08 17 000000 21 55 57.1 -15 12 47   2.356   3.368  177.1   0.9  21.1   47.72    253.7
+2024 08 18 000000 21 54 41.2 -15 18 06   2.361   3.373  177.5   0.7  21.1   47.54    253.9
+2024 08 19 000000 21 53 25.6 -15 23 21   2.367   3.378  177.4   0.8  21.1   47.32    254.0
+2024 08 20 000000 21 52 10.3 -15 28 32   2.372   3.383  176.7   1.0  21.2   47.07    254.1
+2024 08 21 000000 21 50 55.3 -15 33 38   2.378   3.388  175.8   1.3  21.2   46.79    254.3
+2024 08 22 000000 21 49 40.7 -15 38 40   2.385   3.393  174.7   1.6  21.2   46.47    254.4
+2024 08 23 000000 21 48 26.5 -15 43 37   2.391   3.398  173.5   1.9  21.3   46.13    254.6
+2024 08 24 000000 21 47 12.9 -15 48 29   2.398   3.403  172.3   2.3  21.3   45.75    254.7
+2024 08 25 000000 21 45 59.8 -15 53 16   2.406   3.408  171.1   2.6  21.3   45.34    254.9
+2024 08 26 000000 21 44 47.3 -15 57 57   2.413   3.413  169.9   3.0  21.4   44.90    255.0
+2024 08 27 000000 21 43 35.4 -16 02 33   2.421   3.418  168.7   3.3  21.4   44.44    255.2
+2024 08 28 000000 21 42 24.3 -16 07 04   2.430   3.423  167.4   3.7  21.4   43.94    255.3
+2024 08 29 000000 21 41 13.9 -16 11 28   2.438   3.427  166.2   4.0  21.5   43.41    255.5
+2024 08 30 000000 21 40 04.3 -16 15 47   2.447   3.432  165.0   4.4  21.5   42.86    255.6
+2024 08 31 000000 21 38 55.5 -16 19 59   2.456   3.437  163.7   4.7  21.5   42.28    255.8
+2024 09 01 000000 21 37 47.6 -16 24 06   2.466   3.442  162.5   5.1  21.5   41.67    255.9
+2024 09 02 000000 21 36 40.7 -16 28 05   2.476   3.447  161.3   5.4  21.5   41.04    256.1
+2024 09 03 000000 21 35 34.7 -16 31 59   2.486   3.451  160.0   5.7  21.5   40.39    256.3
+2024 09 04 000000 21 34 29.7 -16 35 46   2.496   3.456  158.8   6.1  21.6   39.71    256.4
+2024 09 05 000000 21 33 25.8 -16 39 26   2.507   3.461  157.6   6.4  21.6   39.01    256.6
+2024 09 06 000000 21 32 23.0 -16 43 00   2.518   3.466  156.4   6.7  21.6   38.29    256.8
+2024 09 07 000000 21 31 21.3 -16 46 27   2.530   3.470  155.2   7.0  21.6   37.56    257.0
+2024 09 08 000000 21 30 20.7 -16 49 47   2.541   3.475  154.0   7.3  21.6   36.80    257.1
+2024 09 09 000000 21 29 21.3 -16 53 00   2.553   3.480  152.7   7.6  21.7   36.03    257.3
+2024 09 10 000000 21 28 23.2 -16 56 07   2.566   3.484  151.5   7.9  21.7   35.24    257.5
+2024 09 11 000000 21 27 26.2 -16 59 06   2.578   3.489  150.3   8.2  21.7   34.44    257.7
+2024 09 12 000000 21 26 30.6 -17 01 59   2.591   3.493  149.1   8.5  21.7   33.63    257.9
+2024 09 13 000000 21 25 36.2 -17 04 45   2.604   3.498  148.0   8.8  21.7   32.81    258.1
+2024 09 14 000000 21 24 43.1 -17 07 24   2.617   3.502  146.8   9.1  21.8   31.97    258.3
+2024 09 15 000000 21 23 51.4 -17 09 56   2.631   3.507  145.6   9.3  21.8   31.13    258.5
+2024 09 16 000000 21 23 00.9 -17 12 21   2.645   3.511  144.4   9.6  21.8   30.28    258.7
+2024 09 17 000000 21 22 11.9 -17 14 39   2.659   3.516  143.2   9.8  21.8   29.43    259.0
+2024 09 18 000000 21 21 24.2 -17 16 51   2.673   3.520  142.1  10.1  21.8   28.57    259.2
+2024 09 19 000000 21 20 37.8 -17 18 56   2.688   3.525  140.9  10.3  21.9   27.70    259.5
+2024 09 20 000000 21 19 52.9 -17 20 54   2.702   3.529  139.8  10.6  21.9   26.83    259.7
+2024 09 21 000000 21 19 09.3 -17 22 46   2.717   3.534  138.6  10.8  21.9   25.96    260.0
+2024 09 22 000000 21 18 27.2 -17 24 31   2.733   3.538  137.5  11.1  21.9   25.08    260.3
+2024 09 23 000000 21 17 46.4 -17 26 10   2.748   3.542  136.3  11.3  21.9   24.21    260.5
+2024 09 24 000000 21 17 07.1 -17 27 42   2.764   3.547  135.2  11.5  22.0   23.33    260.8
+2024 09 25 000000 21 16 29.2 -17 29 08   2.779   3.551  134.1  11.7  22.0   22.44    261.2
+2024 09 26 000000 21 15 52.7 -17 30 28   2.796   3.555  132.9  11.9  22.0   21.56    261.5
+2024 09 27 000000 21 15 17.6 -17 31 41   2.812   3.560  131.8  12.1  22.0   20.68    261.9
+2024 09 28 000000 21 14 44.0 -17 32 48   2.828   3.564  130.7  12.3  22.0   19.79    262.3
+2024 09 29 000000 21 14 11.8 -17 33 49   2.845   3.568  129.6  12.5  22.1   18.91    262.7
+2024 09 30 000000 21 13 41.0 -17 34 43   2.862   3.572  128.5  12.7  22.1   18.03    263.1
+2024 10 01 000000 21 13 11.7 -17 35 32   2.879   3.576  127.4  12.8  22.1   17.15    263.6
+2024 10 02 000000 21 12 43.8 -17 36 15   2.896   3.581  126.3  13.0  22.1   16.28    264.2
+2024 10 03 000000 21 12 17.3 -17 36 51   2.913   3.585  125.2  13.2  22.1   15.41    264.8
+2024 10 04 000000 21 11 52.3 -17 37 22   2.931   3.589  124.1  13.3  22.2   14.55    265.4
+2024 10 05 000000 21 11 28.6 -17 37 47   2.948   3.593  123.1  13.5  22.2   13.69    266.1
+2024 10 06 000000 21 11 06.4 -17 38 06   2.966   3.597  122.0  13.6  22.2   12.84    266.9
+2024 10 07 000000 21 10 45.6 -17 38 20   2.984   3.601  120.9  13.8  22.2   12.00    267.8
+2024 10 08 000000 21 10 26.1 -17 38 28   3.002   3.605  119.9  13.9  22.2   11.16    268.9
+2024 10 09 000000 21 10 08.1 -17 38 31   3.020   3.609  118.8  14.0  22.3   10.34    270.0
+2024 10 10 000000 21 09 51.4 -17 38 28   3.038   3.613  117.8  14.2  22.3    9.53    271.4
+2024 10 11 000000 21 09 36.1 -17 38 20   3.057   3.617  116.7  14.3  22.3    8.73    273.0
+2024 10 12 000000 21 09 22.1 -17 38 06   3.075   3.621  115.7  14.4  22.3    7.95    274.8
+2024 10 13 000000 21 09 09.5 -17 37 48   3.094   3.625  114.6  14.5  22.3    7.19    277.1
+2024 10 14 000000 21 08 58.2 -17 37 24   3.113   3.629  113.6  14.6  22.4    6.44    279.8
+2024 10 15 000000 21 08 48.2 -17 36 55   3.132   3.633  112.6  14.7  22.4    5.72    283.2
+
+

+ These calculations have been performed on the + Tamkin + Foundation Computing Network. +

+

+ Valid HTML 4.01! +

+ + diff --git a/sbpy/data/tests/data/TestEphemerisCLI-c1995o1-horizons.txt b/sbpy/data/tests/data/TestEphemerisCLI-c1995o1-horizons.txt new file mode 100644 index 000000000..884bd5108 --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-c1995o1-horizons.txt @@ -0,0 +1,283 @@ +API VERSION: 1.2 +API SOURCE: NASA/JPL Horizons API + +******************************************************************************* +JPL/HORIZONS Hale-Bopp (C/1995 O1) 2024-Aug-16 13:11:45 +Rec #:90002212 (+COV) Soln.date: 2022-Aug-01_01:56:46 # obs: 66 (1993-2022) + +IAU76/J2000 helio. ecliptic osc. elements (au, days, deg., period=Julian yrs): + + EPOCH= 2459837.5 ! 2022-Sep-15.0000000 (TDB) RMSW= n.a. + EC= .9949810027633206 QR= .890537663547794 TP= 2450537.1349071441 + OM= 282.7334213961641 W= 130.4146670659176 IN= 89.28759424740302 + A= 177.4333839117583 MA= 3.878386339423163 ADIST= 353.9762301599687 + PER= 2363.5304681429 N= .000417014 ANGMOM= .02292857 + DAN= 5.00538 DDN= 1.07996 L= 101.8968625 + B= 49.580132 MOID= .0878151 TP= 1997-Mar-29.6349071441 + +Comet physical (GM= km^3/s^2; RAD= km): + GM= n.a. RAD= 30. + M1= 4.8 M2= n.a. k1= 4. k2= n.a. PHCOF= n.a. +COMET comments +1: soln ref.= JPL#226, data arc: 1993-04-27 to 2022-07-09 +2: k1=4.; +******************************************************************************* + + +******************************************************************************* +Ephemeris / API_USER Fri Aug 16 13:11:45 2024 Pasadena, USA / Horizons +******************************************************************************* +Target body name: Hale-Bopp (C/1995 O1) {source: JPL#226} +Center body name: Earth (399) {source: DE441} +Center-site name: GEOCENTRIC +******************************************************************************* +Start time : A.D. 2024-Aug-16 00:00:00.0000 UT +Stop time : A.D. 2024-Oct-15 00:00:00.0000 UT +Step-size : 1440 minutes +******************************************************************************* +Target pole/equ : undefined +Target radii : 30. km +Center geodetic : 0.0, 0.0, -6378.137 {E-lon(deg),Lat(deg),Alt(km)} +Center cylindric: 0.0, 0.0, 0.0 {E-lon(deg),Dxy(km),Dz(km)} +Center pole/equ : ITRF93 {East-longitude positive} +Center radii : 6378.137, 6378.137, 6356.752 km {Equator_a, b, pole_c} +Target primary : Sun +Vis. interferer : MOON (R_eq= 1737.400) km {source: DE441} +Rel. light bend : Sun {source: DE441} +Rel. lght bnd GM: 1.3271E+11 km^3/s^2 +Small-body perts: Yes {source: SB441-N16} +Atmos refraction: NO (AIRLESS) +RA format : DEG +Time format : BOTH +Calendar mode : Mixed Julian/Gregorian +EOP file : eop.240815.p241108 +EOP coverage : DATA-BASED 1962-JAN-20 TO 2024-AUG-15. PREDICTS-> 2024-NOV-07 +Units conversion: 1 au= 149597870.700 km, c= 299792.458 km/s, 1 day= 86400.0 s +Table cut-offs 1: Elevation (-90.0deg=NO ),Airmass (>38.000=NO), Daylight (NO ) +Table cut-offs 2: Solar elongation ( 0.0,180.0=NO ),Local Hour Angle( 0.0=NO ) +Table cut-offs 3: RA/DEC angular rate ( 0.0=NO ) +Table format : Comma Separated Values (spreadsheet) +******************************************************************************* +Initial IAU76/J2000 heliocentric ecliptic osculating elements (au, days, deg.): + EPOCH= 2459837.5 ! 2022-Sep-15.0000000 (TDB) RMSW= n.a. + EC= .9949810027633206 QR= .890537663547794 TP= 2450537.1349071441 + OM= 282.7334213961641 W= 130.4146670659176 IN= 89.28759424740302 + Equivalent ICRF heliocentric cartesian coordinates (au, au/d): + X= 3.907631452214869E+00 Y=-1.373895334060347E+00 Z=-4.624358508575312E+01 + VX= 3.778244409519935E-04 VY=-5.803173067116371E-04 VZ=-3.255716412104052E-03 +Comet physical (GM= km^3/s^2; RAD= km): + GM= n.a. RAD= 30. + M1= 4.8 M2= n.a. k1= 4. k2= n.a. PHCOF= n.a. +******************************************************************************************************************************************************************************************************************************************** + Date__(UT)__HR:MN, Date_________JDUT, , , R.A._(ICRF), DEC_(ICRF), dRA*cosD, d(DEC)/dt, T-mag, N-mag, r, rdot, delta, deldot, S-O-T,/r, S-T-O, PsAng, PsAMV, GlxLon, GlxLat, +******************************************************************************************************************************************************************************************************************************************** +$$SOE + 2024-Aug-16 00:00, 2460538.500000000, , , 339.94076, -85.76646, -2.77385, -0.95007, 19.974, n.a., 48.70398705476, 5.5481506, 48.3840175329577, 14.3142929, 107.8545,/L, 1.1334, 345.385, 125.938, 305.601288,-30.660941, + 2024-Aug-17 00:00, 2460539.500000000, , , 339.68838, -85.77237, -2.78898, -0.89086, 19.974, n.a., 48.70719116755, 5.5479480, 48.3923308961369, 14.4742785, 107.5528,/L, 1.1350, 346.601, 126.285, 305.615209,-30.645503, + 2024-Aug-18 00:00, 2460540.500000000, , , 339.43441, -85.77787, -2.80292, -0.83127, 19.975, n.a., 48.71039516167, 5.5477455, 48.4007359344561, 14.6317701, 107.2462,/L, 1.1366, 347.815, 126.633, 305.629399,-30.630285, + 2024-Aug-19 00:00, 2460541.500000000, , , 339.17897, -85.78297, -2.81568, -0.77129, 19.975, n.a., 48.71359903713, 5.5475429, 48.4092312233878, 14.7868250, 106.9346,/L, 1.1382, 349.028, 126.982, 305.643857,-30.615293, + 2024-Aug-20 00:00, 2460542.500000000, , , 338.92218, -85.78767, -2.82727, -0.71097, 19.976, n.a., 48.71680279399, 5.5473404, 48.4178153742781, 14.9395076, 106.6182,/L, 1.1398, 350.239, 127.330, 305.658577,-30.600529, + 2024-Aug-21 00:00, 2460543.500000000, , , 338.66419, -85.79196, -2.83768, -0.65029, 19.976, n.a., 48.72000643227, 5.5471379, 48.4264870339100, 15.0898732, 106.2971,/L, 1.1414, 351.449, 127.679, 305.673556,-30.585998, + 2024-Aug-22 00:00, 2460544.500000000, , , 338.40511, -85.79584, -2.84691, -0.58926, 19.977, n.a., 48.72320995201, 5.5469354, 48.4352448737881, 15.2379494, 105.9713,/L, 1.1430, 352.656, 128.027, 305.688790,-30.571703, + 2024-Aug-23 00:00, 2460545.500000000, , , 338.14510, -85.79931, -2.85492, -0.52790, 19.977, n.a., 48.72641335325, 5.5467329, 48.4440875702588, 15.3837240, 105.6409,/L, 1.1446, 353.862, 128.376, 305.704276,-30.557649, + 2024-Aug-24 00:00, 2460546.500000000, , , 337.88427, -85.80236, -2.86171, -0.46621, 19.978, n.a., 48.72961663603, 5.5465305, 48.4530137800802, 15.5271418, 105.3061,/L, 1.1461, 355.065, 128.723, 305.720009,-30.543838, + 2024-Aug-25 00:00, 2460547.500000000, , , 337.62278, -85.80500, -2.86724, -0.40423, 19.978, n.a., 48.73281980038, 5.5463281, 48.4620221173881, 15.6681122, 104.9669,/L, 1.1477, 356.267, 129.071, 305.735987,-30.530276, + 2024-Aug-26 00:00, 2460548.500000000, , , 337.36077, -85.80723, -2.87149, -0.34197, 19.979, n.a., 48.73602284634, 5.5461257, 48.4711111365470, 15.8065210, 104.6235,/L, 1.1492, 357.466, 129.417, 305.752205,-30.516965, + 2024-Aug-27 00:00, 2460549.500000000, , , 337.09838, -85.80903, -2.87443, -0.27948, 19.979, n.a., 48.73922577397, 5.5459233, 48.4802793224628, 15.9422438, 104.2758,/L, 1.1507, 358.662, 129.763, 305.768658,-30.503911, + 2024-Aug-28 00:00, 2460550.500000000, , , 336.83575, -85.81041, -2.87604, -0.21679, 19.980, n.a., 48.74242858331, 5.5457210, 48.4895250874486, 16.0751550, 103.9241,/L, 1.1521, 359.857, 130.107, 305.785344,-30.491118, + 2024-Aug-29 00:00, 2460551.500000000, , , 336.57305, -85.81138, -2.87631, -0.15395, 19.980, n.a., 48.74563127441, 5.5455187, 48.4988467726150, 16.2051343, 103.5684,/T, 1.1536, 1.049, 130.450, 305.802257,-30.478588, + 2024-Aug-30 00:00, 2460552.500000000, , , 336.31042, -85.81192, -2.87522, -0.09099, 19.981, n.a., 48.74883384734, 5.5453165, 48.5082426518212, 16.3320694, 103.2088,/T, 1.1550, 2.238, 130.792, 305.819393,-30.466328, + 2024-Aug-31 00:00, 2460553.500000000, , , 336.04801, -85.81203, -2.87277, -0.02797, 19.981, n.a., 48.75203630214, 5.5451143, 48.5177109368424, 16.4558577, 102.8454,/T, 1.1563, 3.424, 131.132, 305.836747,-30.454340, + 2024-Sep-01 00:00, 2460554.500000000, , , 335.78598, -85.81173, -2.86895, 0.035077, 19.982, n.a., 48.75523863888, 5.5449122, 48.5272497830553, 16.5764066, 102.4784,/T, 1.1577, 4.608, 131.470, 305.854314,-30.442628, + 2024-Sep-02 00:00, 2460555.500000000, , , 335.52448, -85.81100, -2.86376, 0.098099, 19.983, n.a., 48.75844085763, 5.5447101, 48.5368572953757, 16.6936341, 102.1078,/T, 1.1590, 5.789, 131.806, 305.872090,-30.431197, + 2024-Sep-03 00:00, 2460556.500000000, , , 335.26366, -85.80985, -2.85720, 0.161051, 19.983, n.a., 48.76164295846, 5.5445081, 48.5465315343617, 16.8074686, 101.7337,/T, 1.1602, 6.967, 132.140, 305.890069,-30.420049, + 2024-Sep-04 00:00, 2460557.500000000, , , 335.00367, -85.80828, -2.84928, 0.223886, 19.984, n.a., 48.76484494143, 5.5443061, 48.5562705223968, 16.9178494, 101.3563,/T, 1.1614, 8.141, 132.471, 305.908245,-30.409188, + 2024-Sep-05 00:00, 2460558.500000000, , , 334.74467, -85.80629, -2.84001, 0.286561, 19.984, n.a., 48.76804680662, 5.5441042, 48.5660722497966, 17.0247260, 100.9757,/T, 1.1626, 9.312, 132.800, 305.926615,-30.398617, + 2024-Sep-06 00:00, 2460559.500000000, , , 334.48680, -85.80389, -2.82939, 0.349028, 19.985, n.a., 48.77124855410, 5.5439024, 48.5759346806325, 17.1280580, 100.5919,/T, 1.1638, 10.480, 133.127, 305.945172,-30.388340, + 2024-Sep-07 00:00, 2460560.500000000, , , 334.23021, -85.80106, -2.81745, 0.411245, 19.985, n.a., 48.77445018395, 5.5437006, 48.5858557580899, 17.2278134, 100.2051,/T, 1.1648, 11.645, 133.450, 305.963910,-30.378360, + 2024-Sep-08 00:00, 2460561.500000000, , , 333.97504, -85.79782, -2.80420, 0.473167, 19.986, n.a., 48.77765169624, 5.5434989, 48.5958334092954, 17.3239683, 99.8154,/T, 1.1659, 12.806, 133.770, 305.982826,-30.368679, + 2024-Sep-09 00:00, 2460562.500000000, , , 333.72143, -85.79417, -2.78965, 0.534752, 19.986, n.a., 48.78085309104, 5.5432972, 48.6058655497306, 17.4165062, 99.4229,/T, 1.1669, 13.964, 134.088, 306.001912,-30.359301, + 2024-Sep-10 00:00, 2460563.500000000, , , 333.46954, -85.79011, -2.77383, 0.595961, 19.987, n.a., 48.78405436842, 5.5430956, 48.6159500875566, 17.5054185, 99.0277,/T, 1.1678, 15.117, 134.401, 306.021163,-30.350227, + 2024-Sep-11 00:00, 2460564.500000000, , , 333.21948, -85.78564, -2.75675, 0.656754, 19.988, n.a., 48.78725552845, 5.5428940, 48.6260849283667, 17.5907052, 98.6298,/T, 1.1687, 16.268, 134.712, 306.040574,-30.341461, + 2024-Sep-12 00:00, 2460565.500000000, , , 332.97139, -85.78077, -2.73846, 0.717094, 19.988, n.a., 48.79045657120, 5.5426925, 48.6362679809981, 17.6723780, 98.2295,/T, 1.1695, 17.414, 135.018, 306.060139,-30.333004, + 2024-Sep-13 00:00, 2460566.500000000, , , 332.72541, -85.77550, -2.71896, 0.776949, 19.989, n.a., 48.79365749673, 5.5424910, 48.6464971649951, 17.7504629, 97.8267,/T, 1.1703, 18.557, 135.321, 306.079853,-30.324858, + 2024-Sep-14 00:00, 2460567.500000000, , , 332.48166, -85.76983, -2.69830, 0.836289, 19.989, n.a., 48.79685830508, 5.5422895, 48.6567704199900, 17.8250038, 97.4216,/T, 1.1711, 19.696, 135.620, 306.099710,-30.317026, + 2024-Sep-15 00:00, 2460568.500000000, , , 332.24025, -85.76377, -2.67651, 0.895091, 19.990, n.a., 48.80005899632, 5.5420881, 48.6670857164922, 17.8960637, 97.0144,/T, 1.1717, 20.831, 135.914, 306.119704,-30.309510, + 2024-Sep-16 00:00, 2460569.500000000, , , 332.00131, -85.75732, -2.65361, 0.953341, 19.990, n.a., 48.80325957049, 5.5418867, 48.6774410662775, 17.9637215, 96.6050,/T, 1.1724, 21.962, 136.204, 306.139830,-30.302310, + 2024-Sep-17 00:00, 2460570.500000000, , , 331.76494, -85.75049, -2.62963, 1.011031, 19.991, n.a., 48.80646002763, 5.5416853, 48.6878345290444, 18.0280619, 96.1936,/T, 1.1729, 23.089, 136.490, 306.160083,-30.295430, + 2024-Sep-18 00:00, 2460571.500000000, , , 331.53126, -85.74327, -2.60460, 1.068160, 19.992, n.a., 48.80966036778, 5.5414839, 48.6982642111620, 18.0891585, 95.7803,/T, 1.1734, 24.213, 136.772, 306.180458,-30.288869, + 2024-Sep-19 00:00, 2460572.500000000, , , 331.30037, -85.73567, -2.57852, 1.124726, 19.992, n.a., 48.81286059098, 5.5412824, 48.7087282534616, 18.1470540, 95.3651,/T, 1.1739, 25.332, 137.048, 306.200948,-30.282631, + 2024-Sep-20 00:00, 2460573.500000000, , , 331.07237, -85.72771, -2.55142, 1.180728, 19.993, n.a., 48.81606069726, 5.5410810, 48.7192248086170, 18.2017451, 94.9482,/T, 1.1743, 26.448, 137.320, 306.221551,-30.276716, + 2024-Sep-21 00:00, 2460574.500000000, , , 330.84737, -85.71937, -2.52328, 1.236155, 19.993, n.a., 48.81926068666, 5.5408796, 48.7297520131740, 18.2531799, 94.5296,/T, 1.1746, 27.560, 137.587, 306.242259,-30.271126, + 2024-Sep-22 00:00, 2460575.500000000, , , 330.62545, -85.71067, -2.49411, 1.290990, 19.994, n.a., 48.82246055921, 5.5406782, 48.7403079617142, 18.3012674, 94.1095,/T, 1.1749, 28.667, 137.849, 306.263069,-30.265864, + 2024-Sep-23 00:00, 2460576.500000000, , , 330.40673, -85.70160, -2.46391, 1.345210, 19.995, n.a., 48.82566031494, 5.5404768, 48.7508906892172, 18.3458952, 93.6878,/T, 1.1751, 29.772, 138.106, 306.283974,-30.260930, + 2024-Sep-24 00:00, 2460577.500000000, , , 330.19129, -85.69218, -2.43268, 1.398785, 19.995, n.a., 48.82885995389, 5.5402754, 48.7614981636654, 18.3869464, 93.2648,/T, 1.1752, 30.872, 138.358, 306.304970,-30.256328, + 2024-Sep-25 00:00, 2460578.500000000, , , 329.97923, -85.68241, -2.40044, 1.451686, 19.996, n.a., 48.83205947610, 5.5400739, 48.7721282871823, 18.4243112, 92.8404,/T, 1.1753, 31.968, 138.604, 306.326051,-30.252058, + 2024-Sep-26 00:00, 2460579.500000000, , , 329.77063, -85.67229, -2.36720, 1.503882, 19.996, n.a., 48.83525888161, 5.5398725, 48.7827789023065, 18.4578935, 92.4149,/T, 1.1753, 33.061, 138.844, 306.347212,-30.248123, + 2024-Sep-27 00:00, 2460580.500000000, , , 329.56559, -85.66183, -2.33298, 1.555342, 19.997, n.a., 48.83845817045, 5.5396711, 48.7934478002502, 18.4876118, 91.9883,/T, 1.1753, 34.150, 139.079, 306.368446,-30.244523, + 2024-Sep-28 00:00, 2460581.500000000, , , 329.36419, -85.65103, -2.29779, 1.606036, 19.997, n.a., 48.84165734269, 5.5394697, 48.8041327291393, 18.5133981, 91.5607,/T, 1.1752, 35.235, 139.309, 306.389747,-30.241262, + 2024-Sep-29 00:00, 2460582.500000000, , , 329.16651, -85.63990, -2.26167, 1.655936, 19.998, n.a., 48.84485639835, 5.5392683, 48.8148314013761, 18.5351963, 91.1323,/T, 1.1750, 36.316, 139.532, 306.411110,-30.238339, + 2024-Sep-30 00:00, 2460583.500000000, , , 328.97261, -85.62844, -2.22463, 1.705015, 19.999, n.a., 48.84805533749, 5.5390669, 48.8255415000383, 18.5529607, 90.7030,/T, 1.1748, 37.393, 139.750, 306.432529,-30.235756, + 2024-Oct-01 00:00, 2460584.500000000, , , 328.78259, -85.61666, -2.18671, 1.753245, 19.999, n.a., 48.85125416015, 5.5388655, 48.8362606845982, 18.5666553, 90.2731,/T, 1.1745, 38.467, 139.961, 306.453998,-30.233515, + 2024-Oct-02 00:00, 2460585.500000000, , , 328.59649, -85.60457, -2.14793, 1.800601, 20.000, n.a., 48.85445286639, 5.5386641, 48.8469865963240, 18.5762531, 89.8427,/T, 1.1741, 39.537, 140.166, 306.475509,-30.231616, + 2024-Oct-03 00:00, 2460586.500000000, , , 328.41439, -85.59217, -2.10833, 1.847060, 20.000, n.a., 48.85765145626, 5.5384628, 48.8577168636373, 18.5817365, 89.4118,/T, 1.1736, 40.603, 140.366, 306.497058,-30.230059, + 2024-Oct-04 00:00, 2460587.500000000, , , 328.23634, -85.57946, -2.06794, 1.892598, 20.001, n.a., 48.86084992982, 5.5382614, 48.8684491075486, 18.5830974, 88.9806,/T, 1.1731, 41.666, 140.558, 306.518637,-30.228846, + 2024-Oct-05 00:00, 2460588.500000000, , , 328.06240, -85.56646, -2.02679, 1.937196, 20.002, n.a., 48.86404828711, 5.5380601, 48.8791809471527, 18.5803374, 88.5491,/T, 1.1725, 42.725, 140.745, 306.540240,-30.227977, + 2024-Oct-06 00:00, 2460589.500000000, , , 327.89261, -85.55318, -1.98493, 1.980833, 20.002, n.a., 48.86724652819, 5.5378588, 48.8899100050949, 18.5734672, 88.1176,/T, 1.1719, 43.780, 140.925, 306.561862,-30.227453, + 2024-Oct-07 00:00, 2460590.500000000, , , 327.72702, -85.53960, -1.94238, 2.023495, 20.003, n.a., 48.87044465313, 5.5376576, 48.9006339129425, 18.5625068, 87.6860,/T, 1.1711, 44.832, 141.099, 306.583496,-30.227272, + 2024-Oct-08 00:00, 2460591.500000000, , , 327.56567, -85.52576, -1.89918, 2.065166, 20.003, n.a., 48.87364266197, 5.5374563, 48.9113503165065, 18.5474857, 87.2545,/T, 1.1704, 45.880, 141.266, 306.605135,-30.227436, + 2024-Oct-09 00:00, 2460592.500000000, , , 327.40859, -85.51164, -1.85539, 2.105837, 20.004, n.a., 48.87684055478, 5.5372551, 48.9220568813085, 18.5284428, 86.8233,/T, 1.1695, 46.925, 141.426, 306.626773,-30.227943, + 2024-Oct-10 00:00, 2460593.500000000, , , 327.25582, -85.49726, -1.81102, 2.145497, 20.005, n.a., 48.88003833162, 5.5370539, 48.9327512985250, 18.5054282, 86.3923,/T, 1.1686, 47.966, 141.580, 306.648404,-30.228794, + 2024-Oct-11 00:00, 2460594.500000000, , , 327.10737, -85.48262, -1.76612, 2.184144, 20.005, n.a., 48.88323599254, 5.5368527, 48.9434312917552, 18.4785044, 85.9618,/T, 1.1676, 49.004, 141.727, 306.670023,-30.229989, + 2024-Oct-12 00:00, 2460595.500000000, , , 326.96328, -85.46774, -1.72073, 2.221776, 20.006, n.a., 48.88643353760, 5.5366516, 48.9540946247641, 18.4477482, 85.5319,/T, 1.1666, 50.039, 141.868, 306.691622,-30.231525, + 2024-Oct-13 00:00, 2460596.500000000, , , 326.82356, -85.45261, -1.67489, 2.258399, 20.006, n.a., 48.88963096687, 5.5364505, 48.9647391097999, 18.4132506, 85.1026,/T, 1.1654, 51.071, 142.002, 306.713196,-30.233403, + 2024-Oct-14 00:00, 2460597.500000000, , , 326.68822, -85.43725, -1.62863, 2.294021, 20.007, n.a., 48.89282828039, 5.5362494, 48.9753626151231, 18.3751135, 84.6740,/T, 1.1643, 52.099, 142.129, 306.734739,-30.235622, + 2024-Oct-15 00:00, 2460598.500000000, , , 326.55727, -85.42166, -1.58197, 2.328660, 20.007, n.a., 48.89602547824, 5.5360484, 48.9859630691605, 18.3334414, 84.2462,/T, 1.1630, 53.125, 142.249, 306.756245,-30.238180, +$$EOE +******************************************************************************************************************************************************************************************************************************************** +Column meaning: + +TIME + + Times PRIOR to 1962 are UT1, a mean-solar time closely related to the +prior but now-deprecated GMT. Times AFTER 1962 are in UTC, the current +civil or "wall-clock" time-scale. UTC is kept within 0.9 seconds of UT1 +using integer leap-seconds for 1972 and later years. + + Conversion from the internal Barycentric Dynamical Time (TDB) of solar +system dynamics to the non-uniform civil UT time-scale requested for output +has not been determined for UTC times after the next July or January 1st. +Therefore, the last known leap-second is used as a constant over future +intervals. + + Time tags refer to the UT time-scale conversion from TDB on Earth +regardless of observer location within the solar system, although clock +rates may differ due to the local gravity field and no analog to "UT" +may be defined for that location. + + Any 'b' symbol in the 1st-column denotes a B.C. date. First-column blank +(" ") denotes an A.D. date. + +CALENDAR SYSTEM + + Mixed calendar mode was active such that calendar dates after AD 1582-Oct-15 +(if any) are in the modern Gregorian system. Dates prior to 1582-Oct-5 (if any) +are in the Julian calendar system, which is automatically extended for dates +prior to its adoption on 45-Jan-1 BC. The Julian calendar is useful for +matching historical dates. The Gregorian calendar more accurately corresponds +to the Earth's orbital motion and seasons. A "Gregorian-only" calendar mode is +available if such physical events are the primary interest. + + NOTE: "n.a." in output means quantity "not available" at the print-time. + + 'R.A._(ICRF), DEC_(ICRF),' = + Astrometric right ascension and declination of the target center with +respect to the observing site (coordinate origin) in the reference frame of +the planetary ephemeris (ICRF). Compensated for down-leg light-time delay +aberration. + + Units: RA in decimal degrees, ddd.fffff{ffff} + DEC in decimal degrees, sdd.fffff{ffff} + + 'dRA*cosD, d(DEC)/dt,' = + The angular rate of change in aparent RA and DEC of the target. This is +with respect to the non-inertial IAU76/80 Earth true equator and equinox +of-date reference frame. d(RA)/dt is multiplied by the cosine of declination +to provide a linear rate in the plane-of-sky. Units: ARCSECONDS PER HOUR + + 'T-mag, N-mag,' = + Comets' apparent visual total magnitude ("T-mag") and nuclear magnitude +("N-mag") using the standard IAU model: + + T-mag= M1 + 5*log10(delta) + k1*log10(r) + N-mag= M2 + 5*log10(delta) + k2*log10(r) + phcof*beta + + Units: MAGNITUDES + + 'r, rdot,' = + The Sun's apparent range ("r", light-time aberrated) and range-rate ("rdot") +relative to the target center, as seen by the observer. A positive "rdot" means +the target center was moving away from the Sun, negative means moving toward +the Sun. Units: AU and KM/S + + 'delta, deldot,' = + Apparent range ("delta", light-time aberrated) and range-rate ("delta-dot") +of the target center relative to the observer. A positive "deldot" means the +target center is moving away from the observer, negative indicates movement +toward the observer. Units: AU and KM/S + + 'S-O-T,/r,' = + Sun-Observer-Target apparent SOLAR ELONGATION ANGLE seen from the observers' +location at print-time. + + The '/r' column provides a code indicating the targets' apparent position +relative to the Sun in the observers' sky, as described below: + + Case A: For an observing location on the surface of a rotating body, that +body rotational sense is considered: + + /T indicates target TRAILS Sun (evening sky: rises and sets AFTER Sun) + /L indicates target LEADS Sun (morning sky: rises and sets BEFORE Sun) + + Case B: For an observing point that does not have a rotational model (such +as a spacecraft), the "leading" and "trailing" condition is defined by the +observers' heliocentric ORBITAL motion: + + * If continuing in the observers' current direction of heliocentric + motion would encounter the targets' apparent longitude first, followed + by the Sun's, the target LEADS the Sun as seen by the observer. + + * If the Sun's apparent longitude would be encountered first, followed + by the targets', the target TRAILS the Sun. + + Two other codes can be output: + /* indicates observer is Sun-centered (undefined) + /? Target is aligned with Sun center (no lead or trail) + + The S-O-T solar elongation angle is numerically the minimum separation +angle of the Sun and target in the sky in any direction. It does NOT indicate +the amount of separation in the leading or trailing directions, which would +be defined along the equator of a spherical coordinate system. + + Units: DEGREES + + 'S-T-O,' = + The Sun-Target-Observer angle; the interior vertex angle at target center +formed by a vector from the target to the apparent center of the Sun (at +reflection time on the target) and the apparent vector from target to the +observer at print-time. Slightly different from true PHASE ANGLE (requestable +separately) at the few arcsecond level in that it includes stellar aberration +on the down-leg from target to observer. Units: DEGREES + + 'PsAng, PsAMV,' = + The position angles of the extended Sun-to-target radius vector ("PsAng") +and the negative of the targets' heliocentric velocity vector ("PsAMV"), as +seen in the observers' plane-of-sky, measured counter-clockwise (east) from +reference-frame north-pole. Primarily intended for ACTIVE COMETS, "PsAng" +is an indicator of the comets' gas-tail orientation in the sky (being in the +anti-sunward direction) while "PsAMV" is an indicator of dust-tail orientation. +Units: DEGREES + + 'GlxLon, GlxLat,' = + Observer-centered Galactic System II (post WW II) longitude and latitude +of the target centers' apparent position, with light-time, gravitational +deflection of light, and stellar aberrations. Units: DEGREES + +Computations by ... + + Solar System Dynamics Group, Horizons On-Line Ephemeris System + 4800 Oak Grove Drive, Jet Propulsion Laboratory + Pasadena, CA 91109 USA + + General site: https://ssd.jpl.nasa.gov/ + Mailing list: https://ssd.jpl.nasa.gov/email_list.html + System news : https://ssd.jpl.nasa.gov/horizons/news.html + User Guide : https://ssd.jpl.nasa.gov/horizons/manual.html + Connect : browser https://ssd.jpl.nasa.gov/horizons/app.html#/x + API https://ssd-api.jpl.nasa.gov/doc/horizons.html + command-line telnet ssd.jpl.nasa.gov 6775 + e-mail/batch https://ssd.jpl.nasa.gov/ftp/ssd/hrzn_batch.txt + scripts https://ssd.jpl.nasa.gov/ftp/ssd/SCRIPTS + Author : Jon.D.Giorgini@jpl.nasa.gov + +******************************************************************************************************************************************************************************************************************************************** diff --git a/sbpy/data/tests/data/TestEphemerisCLI-ceres-horizons.txt b/sbpy/data/tests/data/TestEphemerisCLI-ceres-horizons.txt new file mode 100644 index 000000000..606789c82 --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-ceres-horizons.txt @@ -0,0 +1,295 @@ +API VERSION: 1.2 +API SOURCE: NASA/JPL Horizons API + +******************************************************************************* +JPL/HORIZONS 1 Ceres (A801 AA) 2024-Sep-02 16:04:00 +Rec #: 1 (+COV) Soln.date: 2021-Apr-13_11:04:44 # obs: 1075 (1995-2021) + +IAU76/J2000 helio. ecliptic osc. elements (au, days, deg., period=Julian yrs): + + EPOCH= 2458849.5 ! 2020-Jan-01.00 (TDB) Residual RMS= .24563 + EC= .07687465013145245 QR= 2.556401146697176 TP= 2458240.1791309435 + OM= 80.3011901917491 W= 73.80896808746482 IN= 10.59127767086216 + A= 2.769289292143484 MA= 130.3159688200986 ADIST= 2.982177437589792 + PER= 4.60851 N= .213870844 ANGMOM= .028541613 + DAN= 2.69515 DDN= 2.81323 L= 153.8445987 + B= 10.1666387 MOID= 1.59231997 TP= 2018-May-01.6791309435 + +Asteroid physical parameters (km, seconds, rotational period in hours): + GM= 62.6284 RAD= 469.7 ROTPER= 9.07417 + H= 3.34 G= .120 B-V= .713 + ALBEDO= .090 STYP= C + +ASTEROID comments: +1: soln ref.= JPL#48, OCC=0 radar(60 delay, 0 Dop.) +2: source=ORB +******************************************************************************* + + +******************************************************************************* +Ephemeris / API_USER Mon Sep 2 16:04:00 2024 Pasadena, USA / Horizons +******************************************************************************* +Target body name: 1 Ceres (A801 AA) {source: JPL#48} +Center body name: Earth (399) {source: DE441} +Center-site name: GEOCENTRIC +******************************************************************************* +Start time : A.D. 2024-Aug-16 00:00:00.0000 UT +Stop time : A.D. 2024-Oct-15 00:00:00.0000 UT +Step-size : 1440 minutes +******************************************************************************* +Target pole/equ : IAU {East-longitude positive} +Target radii : 487.3, 487.3, 446.0 km {Equator_a, b, pole_c} +Center geodetic : 0.0, 0.0, -6378.137 {E-lon(deg),Lat(deg),Alt(km)} +Center cylindric: 0.0, 0.0, 0.0 {E-lon(deg),Dxy(km),Dz(km)} +Center pole/equ : ITRF93 {East-longitude positive} +Center radii : 6378.137, 6378.137, 6356.752 km {Equator_a, b, pole_c} +Target primary : Sun +Vis. interferer : MOON (R_eq= 1737.400) km {source: DE441} +Rel. light bend : Sun {source: DE441} +Rel. lght bnd GM: 1.3271E+11 km^3/s^2 +Small-body perts: Yes {source: SB441-N16} +Atmos refraction: NO (AIRLESS) +RA format : DEG +Time format : BOTH +Calendar mode : Mixed Julian/Gregorian +EOP file : eop.240902.p241126 +EOP coverage : DATA-BASED 1962-JAN-20 TO 2024-SEP-02. PREDICTS-> 2024-NOV-25 +Units conversion: 1 au= 149597870.700 km, c= 299792.458 km/s, 1 day= 86400.0 s +Table cut-offs 1: Elevation (-90.0deg=NO ),Airmass (>38.000=NO), Daylight (NO ) +Table cut-offs 2: Solar elongation ( 0.0,180.0=NO ),Local Hour Angle( 0.0=NO ) +Table cut-offs 3: RA/DEC angular rate ( 0.0=NO ) +Table format : Comma Separated Values (spreadsheet) +******************************************************************************* +Initial IAU76/J2000 heliocentric ecliptic osculating elements (au, days, deg.): + EPOCH= 2458849.5 ! 2020-Jan-01.00 (TDB) Residual RMS= .24563 + EC= .07687465013145245 QR= 2.556401146697176 TP= 2458240.1791309435 + OM= 80.3011901917491 W= 73.80896808746482 IN= 10.59127767086216 + Equivalent ICRF heliocentric cartesian coordinates (au, au/d): + X= 1.007608869613381E+00 Y=-2.390064275223502E+00 Z=-1.332124522752402E+00 + VX= 9.201724467227128E-03 VY= 3.370381135398406E-03 VZ=-2.850337057661093E-04 +Asteroid physical parameters (km, seconds, rotational period in hours): + GM= 62.6284 RAD= 469.7 ROTPER= 9.07417 + H= 3.34 G= .120 B-V= .713 + ALBEDO= .090 STYP= C +******************************************************************************************************************************************************************************************************************************************** + Date__(UT)__HR:MN, Date_________JDUT, , , R.A._(ICRF), DEC_(ICRF), dRA*cosD, d(DEC)/dt, APmag, S-brt, r, rdot, delta, deldot, S-O-T,/r, S-T-O, PsAng, PsAMV, GlxLon, GlxLat, +******************************************************************************************************************************************************************************************************************************************** +$$SOE + 2024-Aug-16 00:00, 2460538.500000000, , , 278.62427, -30.91559, -9.07628, -1.29427, 8.161, 6.877, 2.923502437983, 0.9363365, 2.13108333791703, 18.0098740, 133.6759,/T, 14.5034, 79.132, 277.528, 3.328003,-10.250208, + 2024-Aug-17 00:00, 2460539.500000000, , , 278.55714, -30.92349, -8.19537, -1.12419, 8.181, 6.886, 2.924042154404, 0.9327665, 2.14156758454322, 18.2946213, 132.6732,/T, 14.7433, 79.412, 277.581, 3.295436,-10.201622, + 2024-Aug-18 00:00, 2460540.500000000, , , 278.49685, -30.93028, -7.31382, -0.95882, 8.200, 6.894, 2.924579805950, 0.9291873, 2.15221429009955, 18.5724724, 131.6757,/T, 14.9775, 79.680, 277.632, 3.266442,-10.157875, + 2024-Aug-19 00:00, 2460541.500000000, , , 278.44342, -30.93601, -6.43236, -0.79805, 8.220, 6.902, 2.925115387354, 0.9255990, 2.16301953036993, 18.8436353, 130.6834,/T, 15.2059, 79.937, 277.678, 3.240999,-10.118976, + 2024-Aug-20 00:00, 2460542.500000000, , , 278.39683, -30.94069, -5.55150, -0.64173, 8.239, 6.910, 2.925648893374, 0.9220016, 2.17397950692828, 19.1083328, 129.6962,/T, 15.4284, 80.185, 277.720, 3.219082,-10.084932, + 2024-Aug-21 00:00, 2460543.500000000, , , 278.35709, -30.94437, -4.67158, -0.48971, 8.258, 6.918, 2.926180318796, 0.9183952, 2.18509054651981, 19.3667695, 128.7142,/T, 15.6451, 80.422, 277.759, 3.200667,-10.055747, + 2024-Aug-22 00:00, 2460544.500000000, , , 278.32419, -30.94706, -3.79275, -0.34182, 8.278, 6.925, 2.926709658434, 0.9147799, 2.19634908007312, 19.6190963, 127.7373,/T, 15.8561, 80.650, 277.794, 3.185729,-10.031425, + 2024-Aug-23 00:00, 2460545.500000000, , , 278.29811, -30.94879, -2.91506, -0.19790, 8.296, 6.932, 2.927236907127, 0.9111555, 2.20775160402326, 19.8653867, 126.7656,/T, 16.0612, 80.869, 277.825, 3.174243,-10.011968, + 2024-Aug-24 00:00, 2460546.500000000, , , 278.27886, -30.94959, -2.03858, -0.05779, 8.315, 6.939, 2.927762059745, 0.9075223, 2.21929463337660, 20.1056327, 125.7989,/T, 16.2606, 81.079, 277.853, 3.166183, -9.997377, + 2024-Aug-25 00:00, 2460547.500000000, , , 278.26642, -30.94949, -1.16342, 0.078656, 8.334, 6.946, 2.928285111184, 0.9038803, 2.23097465829085, 20.3397599, 124.8374,/T, 16.4541, 81.281, 277.876, 3.161525, -9.987656, + 2024-Aug-26 00:00, 2460548.500000000, , , 278.26077, -30.94851, -0.28982, 0.211593, 8.352, 6.952, 2.928806056370, 0.9002295, 2.24278811278660, 20.5676527, 123.8809,/T, 16.6419, 81.476, 277.896, 3.160244, -9.982804, + 2024-Aug-27 00:00, 2460549.500000000, , , 278.26192, -30.94667, 0.581904, 0.341159, 8.371, 6.958, 2.929324890258, 0.8965699, 2.25473135834962, 20.7891793, 122.9296,/T, 16.8239, 81.662, 277.912, 3.162315, -9.982821, + 2024-Aug-28 00:00, 2460550.500000000, , , 278.26983, -30.94399, 1.451330, 0.467489, 8.389, 6.964, 2.929841607831, 0.8929017, 2.26680068032992, 21.0042105, 121.9833,/T, 17.0002, 81.841, 277.924, 3.167711, -9.987703, + 2024-Aug-29 00:00, 2460551.500000000, , , 278.28449, -30.94050, 2.317979, 0.590720, 8.407, 6.970, 2.930356204103, 0.8892248, 2.27899229295528, 21.2126304, 121.0421,/T, 17.1706, 82.013, 277.932, 3.176406, -9.997447, + 2024-Aug-30 00:00, 2460552.500000000, , , 278.30588, -30.93621, 3.181329, 0.710988, 8.424, 6.976, 2.930868674117, 0.8855393, 2.29130234901791, 21.4143426, 120.1059,/T, 17.3354, 82.178, 277.936, 3.188370,-10.012042, + 2024-Aug-31 00:00, 2460553.500000000, , , 278.33396, -30.93114, 4.040837, 0.828430, 8.442, 6.981, 2.931379012947, 0.8818453, 2.30372695157877, 21.6092717, 119.1749,/T, 17.4944, 82.336, 277.936, 3.203574,-10.031478, + 2024-Sep-01 00:00, 2460554.500000000, , , 278.36871, -30.92532, 4.895950, 0.943183, 8.459, 6.986, 2.931887215695, 0.8781429, 2.31626216633511, 21.7973637, 118.2488,/T, 17.6478, 82.488, 277.933, 3.221987,-10.055740, + 2024-Sep-02 00:00, 2460555.500000000, , , 278.41008, -30.91876, 5.746117, 1.055385, 8.477, 6.991, 2.932393277495, 0.8744320, 2.32890403414578, 21.9785865, 117.3279,/T, 17.7955, 82.633, 277.926, 3.243575,-10.084810, + 2024-Sep-03 00:00, 2460556.500000000, , , 278.45804, -30.91147, 6.590794, 1.165176, 8.494, 6.996, 2.932897193508, 0.8707127, 2.34164858355722, 22.1529295, 116.4119,/T, 17.9376, 82.772, 277.914, 3.268305,-10.118667, + 2024-Sep-04 00:00, 2460557.500000000, , , 278.51253, -30.90347, 7.429456, 1.272696, 8.511, 7.001, 2.933398958925, 0.8669851, 2.35449184317403, 22.3204035, 115.5010,/T, 18.0741, 82.906, 277.900, 3.296143,-10.157285, + 2024-Sep-05 00:00, 2460558.500000000, , , 278.57351, -30.89478, 8.261600, 1.378084, 8.527, 7.005, 2.933898568965, 0.8632492, 2.36742985358321, 22.4810390, 114.5950,/T, 18.2050, 83.033, 277.881, 3.327051,-10.200638, + 2024-Sep-06 00:00, 2460559.500000000, , , 278.64093, -30.88541, 9.086753, 1.481480, 8.544, 7.010, 2.934396018874, 0.8595051, 2.38045867844869, 22.6348850, 113.6940,/T, 18.3305, 83.155, 277.858, 3.360992,-10.248695, + 2024-Sep-07 00:00, 2460560.500000000, , , 278.71472, -30.87537, 9.904472, 1.583023, 8.560, 7.014, 2.934891303924, 0.8557528, 2.39357441445142, 22.7820065, 112.7979,/T, 18.4505, 83.272, 277.832, 3.397928,-10.301422, + 2024-Sep-08 00:00, 2460561.500000000, , , 278.79483, -30.86468, 10.71435, 1.682850, 8.576, 7.018, 2.935384419410, 0.8519923, 2.40677319997890, 22.9224826, 111.9066,/T, 18.5651, 83.383, 277.803, 3.437819,-10.358784, + 2024-Sep-09 00:00, 2460562.500000000, , , 278.88119, -30.85334, 11.51599, 1.781097, 8.592, 7.021, 2.935875360652, 0.8482238, 2.42005122282099, 23.0564052, 111.0202,/T, 18.6743, 83.489, 277.769, 3.480626,-10.420742, + 2024-Sep-10 00:00, 2460563.500000000, , , 278.97373, -30.84136, 12.30906, 1.877897, 8.608, 7.025, 2.936364122989, 0.8444471, 2.43340472752630, 23.1838792, 110.1386,/T, 18.7783, 83.590, 277.732, 3.526309,-10.487255, + 2024-Sep-11 00:00, 2460564.500000000, , , 279.07238, -30.82876, 13.09321, 1.973378, 8.624, 7.028, 2.936850701784, 0.8406625, 2.44683002342992, 23.3050245, 109.2618,/T, 18.8770, 83.687, 277.691, 3.574825,-10.558281, + 2024-Sep-12 00:00, 2460565.500000000, , , 279.17707, -30.81554, 13.86814, 2.067668, 8.639, 7.031, 2.937335092414, 0.8368698, 2.46032349457724, 23.4199808, 108.3896,/T, 18.9705, 83.779, 277.647, 3.626134,-10.633775, + 2024-Sep-13 00:00, 2460566.500000000, , , 279.28773, -30.80171, 14.63358, 2.160891, 8.655, 7.034, 2.937817290280, 0.8330692, 2.47388161268740, 23.5289132, 107.5220,/T, 19.0589, 83.866, 277.599, 3.680193,-10.713689, + 2024-Sep-14 00:00, 2460567.500000000, , , 279.40428, -30.78727, 15.38928, 2.253169, 8.670, 7.037, 2.938297290798, 0.8292606, 2.48750095367901, 23.6320182, 106.6590,/T, 19.1422, 83.949, 277.548, 3.736960,-10.797975, + 2024-Sep-15 00:00, 2460568.500000000, , , 279.52663, -30.77224, 16.13509, 2.344626, 8.685, 7.040, 2.938775089401, 0.8254441, 2.50117821681100, 23.7295262, 105.8005,/T, 19.2206, 84.027, 277.493, 3.796393,-10.886584, + 2024-Sep-16 00:00, 2460569.500000000, , , 279.65472, -30.75661, 16.87098, 2.435390, 8.699, 7.043, 2.939250681543, 0.8216197, 2.51491024304133, 23.8216956, 104.9464,/T, 19.2940, 84.101, 277.435, 3.858450,-10.979465, + 2024-Sep-17 00:00, 2460570.500000000, , , 279.78846, -30.74040, 17.59704, 2.525595, 8.714, 7.045, 2.939724062696, 0.8177875, 2.52869402632583, 23.9087934, 104.0967,/T, 19.3625, 84.172, 277.374, 3.923089,-11.076566, + 2024-Sep-18 00:00, 2460571.500000000, , , 279.92776, -30.72359, 18.31352, 2.615379, 8.728, 7.047, 2.940195228352, 0.8139475, 2.54252670999696, 23.9910637, 103.2512,/T, 19.4262, 84.238, 277.310, 3.990269,-11.177839, + 2024-Sep-19 00:00, 2460572.500000000, , , 280.07257, -30.70620, 19.02082, 2.704885, 8.742, 7.050, 2.940664174022, 0.8100997, 2.55640556250219, 24.0686898, 102.4099,/T, 19.4852, 84.301, 277.242, 4.059954,-11.283234, + 2024-Sep-20 00:00, 2460573.500000000, , , 280.22280, -30.68823, 19.71939, 2.794249, 8.756, 7.051, 2.941130895242, 0.8062442, 2.57032793357556, 24.1417668, 101.5726,/T, 19.5396, 84.359, 277.171, 4.132105,-11.392706, + 2024-Sep-21 00:00, 2460574.500000000, , , 280.37838, -30.66968, 20.40964, 2.883597, 8.770, 7.053, 2.941595387569, 0.8023810, 2.58429120043038, 24.2102961, 100.7394,/T, 19.5892, 84.415, 277.097, 4.206688,-11.506212, + 2024-Sep-22 00:00, 2460575.500000000, , , 280.53926, -30.65054, 21.09193, 2.973042, 8.783, 7.055, 2.942057646582, 0.7985102, 2.59829271811011, 24.2742042, 99.9101,/T, 19.6344, 84.466, 277.019, 4.283671,-11.623710, + 2024-Sep-23 00:00, 2460576.500000000, , , 280.70536, -30.63081, 21.76644, 3.062677, 8.797, 7.057, 2.942517667889, 0.7946317, 2.61232978539736, 24.3333756, 99.0847,/T, 19.6749, 84.515, 276.939, 4.363022,-11.745160, + 2024-Sep-24 00:00, 2460577.500000000, , , 280.87662, -30.61050, 22.43327, 3.152584, 8.810, 7.058, 2.942975447121, 0.7907457, 2.62639963006079, 24.3876846, 98.2631,/T, 19.7111, 84.560, 276.856, 4.444711,-11.870523, + 2024-Sep-25 00:00, 2460578.500000000, , , 281.05298, -30.58961, 23.09240, 3.242833, 8.823, 7.059, 2.943430979933, 0.7868522, 2.64049941013030, 24.4370182, 97.4453,/T, 19.7427, 84.602, 276.770, 4.528706,-11.999761, + 2024-Sep-26 00:00, 2460579.500000000, , , 281.23438, -30.56812, 23.74375, 3.333485, 8.836, 7.061, 2.943884262009, 0.7829512, 2.65462622471897, 24.4812863, 96.6312,/T, 19.7701, 84.640, 276.680, 4.614979,-12.132832, + 2024-Sep-27 00:00, 2460580.500000000, , , 281.42075, -30.54603, 24.38721, 3.424598, 8.849, 7.062, 2.944335289058, 0.7790428, 2.66877712841490, 24.5204247, 95.8208,/T, 19.7930, 84.676, 276.588, 4.703500,-12.269698, + 2024-Sep-28 00:00, 2460581.500000000, , , 281.61203, -30.52335, 25.02265, 3.516225, 8.861, 7.062, 2.944784056816, 0.7751271, 2.68294914544770, 24.5543919, 95.0140,/T, 19.8118, 84.709, 276.493, 4.794239,-12.410317, + 2024-Sep-29 00:00, 2460582.500000000, , , 281.80815, -30.50007, 25.64993, 3.608414, 8.873, 7.063, 2.945230561046, 0.7712040, 2.69713928201387, 24.5831662, 94.2108,/T, 19.8263, 84.739, 276.396, 4.887167,-12.554647, + 2024-Sep-30 00:00, 2460583.500000000, , , 282.00905, -30.47618, 26.26893, 3.701211, 8.885, 7.064, 2.945674797537, 0.7672736, 2.71134453660760, 24.6067424, 93.4111,/T, 19.8366, 84.766, 276.295, 4.982254,-12.702646, + 2024-Oct-01 00:00, 2460584.500000000, , , 282.21465, -30.45168, 26.87950, 3.794658, 8.897, 7.064, 2.946116762105, 0.7633360, 2.72556190890574, 24.6251296, 92.6150,/T, 19.8428, 84.790, 276.192, 5.079471,-12.854269, + 2024-Oct-02 00:00, 2460585.500000000, , , 282.42489, -30.42656, 27.48153, 3.888794, 8.909, 7.065, 2.946556450592, 0.7593912, 2.73978840791017, 24.6383508, 91.8223,/T, 19.8450, 84.812, 276.086, 5.178787,-13.009471, + 2024-Oct-03 00:00, 2460586.500000000, , , 282.63970, -30.40083, 28.07489, 3.983651, 8.920, 7.065, 2.946993858867, 0.7554392, 2.75402105988979, 24.6464429, 91.0330,/T, 19.8431, 84.831, 275.977, 5.280175,-13.168207, + 2024-Oct-04 00:00, 2460587.500000000, , , 282.85900, -30.37447, 28.65945, 4.079263, 8.932, 7.065, 2.947428982827, 0.7514802, 2.76825691637979, 24.6494565, 90.2471,/T, 19.8373, 84.847, 275.866, 5.383604,-13.330431, + 2024-Oct-05 00:00, 2460588.500000000, , , 283.08273, -30.34748, 29.23511, 4.175654, 8.943, 7.065, 2.947861818392, 0.7475141, 2.78249306223278, 24.6474564, 89.4646,/T, 19.8277, 84.862, 275.753, 5.489045,-13.496094, + 2024-Oct-06 00:00, 2460589.500000000, , , 283.31081, -30.31986, 29.80178, 4.272850, 8.954, 7.065, 2.948292361513, 0.7435409, 2.79672662357929, 24.6405206, 88.6853,/T, 19.8141, 84.874, 275.636, 5.596469,-13.665150, + 2024-Oct-07 00:00, 2460590.500000000, , , 283.54317, -30.29159, 30.35937, 4.370871, 8.965, 7.065, 2.948720608162, 0.7395609, 2.81095477559160, 24.6287408, 87.9094,/T, 19.7969, 84.883, 275.518, 5.705848,-13.837549, + 2024-Oct-08 00:00, 2460591.500000000, , , 283.77972, -30.26268, 30.90780, 4.469734, 8.975, 7.065, 2.949146554341, 0.7355739, 2.82517475014633, 24.6122219, 87.1366,/T, 19.7759, 84.891, 275.397, 5.817151,-14.013241, + 2024-Oct-09 00:00, 2460592.500000000, , , 284.02041, -30.23312, 31.44701, 4.569454, 8.986, 7.064, 2.949570196076, 0.7315801, 2.83938384377555, 24.5910831, 86.3670,/T, 19.7512, 84.897, 275.274, 5.930352,-14.192177, + 2024-Oct-10 00:00, 2460593.500000000, , , 284.26515, -30.20289, 31.97697, 4.670041, 8.996, 7.064, 2.949991529420, 0.7275794, 2.85357942656845, 24.5654604, 85.6004,/T, 19.7230, 84.900, 275.148, 6.045421,-14.374306, + 2024-Oct-11 00:00, 2460594.500000000, , , 284.51387, -30.17200, 32.49765, 4.771505, 9.006, 7.063, 2.950410550451, 0.7235719, 2.86775895277410, 24.5355099, 84.8370,/T, 19.6912, 84.902, 275.021, 6.162331,-14.559577, + 2024-Oct-12 00:00, 2460595.500000000, , , 284.76649, -30.14044, 33.00906, 4.873856, 9.016, 7.063, 2.950827255274, 0.7195577, 2.88191997355603, 24.5014116, 84.0766,/T, 19.6560, 84.902, 274.891, 6.281054,-14.747940, + 2024-Oct-13 00:00, 2460596.500000000, , , 285.02293, -30.10821, 33.51125, 4.977107, 9.025, 7.062, 2.951241640017, 0.7155369, 2.89606015141041, 24.4633709, 83.3191,/T, 19.6173, 84.900, 274.759, 6.401564,-14.939343, + 2024-Oct-14 00:00, 2460597.500000000, , , 285.28313, -30.07529, 34.00436, 5.081278, 9.035, 7.061, 2.951653700836, 0.7115094, 2.91017727401272, 24.4216145, 82.5645,/T, 19.5753, 84.897, 274.625, 6.523834,-15.133735, + 2024-Oct-15 00:00, 2460598.500000000, , , 285.54700, -30.04168, 34.48860, 5.186400, 9.044, 7.060, 2.952063433911, 0.7074752, 2.92426926290955, 24.3763759, 81.8127,/T, 19.5300, 84.892, 274.489, 6.647840,-15.331068, +$$EOE +******************************************************************************************************************************************************************************************************************************************** +Column meaning: + +TIME + + Times PRIOR to 1962 are UT1, a mean-solar time closely related to the +prior but now-deprecated GMT. Times AFTER 1962 are in UTC, the current +civil or "wall-clock" time-scale. UTC is kept within 0.9 seconds of UT1 +using integer leap-seconds for 1972 and later years. + + Conversion from the internal Barycentric Dynamical Time (TDB) of solar +system dynamics to the non-uniform civil UT time-scale requested for output +has not been determined for UTC times after the next July or January 1st. +Therefore, the last known leap-second is used as a constant over future +intervals. + + Time tags refer to the UT time-scale conversion from TDB on Earth +regardless of observer location within the solar system, although clock +rates may differ due to the local gravity field and no analog to "UT" +may be defined for that location. + + Any 'b' symbol in the 1st-column denotes a B.C. date. First-column blank +(" ") denotes an A.D. date. + +CALENDAR SYSTEM + + Mixed calendar mode was active such that calendar dates after AD 1582-Oct-15 +(if any) are in the modern Gregorian system. Dates prior to 1582-Oct-5 (if any) +are in the Julian calendar system, which is automatically extended for dates +prior to its adoption on 45-Jan-1 BC. The Julian calendar is useful for +matching historical dates. The Gregorian calendar more accurately corresponds +to the Earth's orbital motion and seasons. A "Gregorian-only" calendar mode is +available if such physical events are the primary interest. + + NOTE: "n.a." in output means quantity "not available" at the print-time. + + 'R.A._(ICRF), DEC_(ICRF),' = + Astrometric right ascension and declination of the target center with +respect to the observing site (coordinate origin) in the reference frame of +the planetary ephemeris (ICRF). Compensated for down-leg light-time delay +aberration. + + Units: RA in decimal degrees, ddd.fffff{ffff} + DEC in decimal degrees, sdd.fffff{ffff} + + 'dRA*cosD, d(DEC)/dt,' = + The angular rate of change in aparent RA and DEC of the target. This is +with respect to the non-inertial IAU76/80 Earth true equator and equinox +of-date reference frame. d(RA)/dt is multiplied by the cosine of declination +to provide a linear rate in the plane-of-sky. Units: ARCSECONDS PER HOUR + + 'APmag, S-brt,' = + The asteroids' approximate apparent airless visual magnitude and surface +brightness using the standard IAU H-G system magnitude model: + + APmag = H + 5*log10(delta) + 5*log10(r) - 2.5*log10((1-G)*phi_1 + G*phi_2) + + For solar phase angles >90 deg, the error could exceed 1 magnitude. For +phase angles >120 degrees, output values are rounded to the nearest integer to +indicate error could be large and unknown. For Earth-based observers, the +estimated dimming due to atmospheric absorption (extinction) is available as +a separate, requestable quantity. + + Surface brightness is the average airless visual magnitude of a +square-arcsecond of the illuminated portion of the apparent disk. It is +computed only if the target radius is known. + + Units: MAGNITUDES & MAGNITUDES PER SQUARE ARCSECOND + + 'r, rdot,' = + The Sun's apparent range ("r", light-time aberrated) and range-rate ("rdot") +relative to the target center, as seen by the observer. A positive "rdot" means +the target center was moving away from the Sun, negative means moving toward +the Sun. Units: AU and KM/S + + 'delta, deldot,' = + Apparent range ("delta", light-time aberrated) and range-rate ("delta-dot") +of the target center relative to the observer. A positive "deldot" means the +target center is moving away from the observer, negative indicates movement +toward the observer. Units: AU and KM/S + + 'S-O-T,/r,' = + Sun-Observer-Target apparent SOLAR ELONGATION ANGLE seen from the observers' +location at print-time. + + The '/r' column provides a code indicating the targets' apparent position +relative to the Sun in the observers' sky, as described below: + + Case A: For an observing location on the surface of a rotating body, that +body rotational sense is considered: + + /T indicates target TRAILS Sun (evening sky: rises and sets AFTER Sun) + /L indicates target LEADS Sun (morning sky: rises and sets BEFORE Sun) + + Case B: For an observing point that does not have a rotational model (such +as a spacecraft), the "leading" and "trailing" condition is defined by the +observers' heliocentric ORBITAL motion: + + * If continuing in the observers' current direction of heliocentric + motion would encounter the targets' apparent longitude first, followed + by the Sun's, the target LEADS the Sun as seen by the observer. + + * If the Sun's apparent longitude would be encountered first, followed + by the targets', the target TRAILS the Sun. + + Two other codes can be output: + /* indicates observer is Sun-centered (undefined) + /? Target is aligned with Sun center (no lead or trail) + + The S-O-T solar elongation angle is numerically the minimum separation +angle of the Sun and target in the sky in any direction. It does NOT indicate +the amount of separation in the leading or trailing directions, which would +be defined along the equator of a spherical coordinate system. + + Units: DEGREES + + 'S-T-O,' = + The Sun-Target-Observer angle; the interior vertex angle at target center +formed by a vector from the target to the apparent center of the Sun (at +reflection time on the target) and the apparent vector from target to the +observer at print-time. Slightly different from true PHASE ANGLE (requestable +separately) at the few arcsecond level in that it includes stellar aberration +on the down-leg from target to observer. Units: DEGREES + + 'PsAng, PsAMV,' = + The position angles of the extended Sun-to-target radius vector ("PsAng") +and the negative of the targets' heliocentric velocity vector ("PsAMV"), as +seen in the observers' plane-of-sky, measured counter-clockwise (east) from +reference-frame north-pole. Primarily intended for ACTIVE COMETS, "PsAng" +is an indicator of the comets' gas-tail orientation in the sky (being in the +anti-sunward direction) while "PsAMV" is an indicator of dust-tail orientation. +Units: DEGREES + + 'GlxLon, GlxLat,' = + Observer-centered Galactic System II (post WW II) longitude and latitude +of the target centers' apparent position, with light-time, gravitational +deflection of light, and stellar aberrations. Units: DEGREES + +Computations by ... + + Solar System Dynamics Group, Horizons On-Line Ephemeris System + 4800 Oak Grove Drive, Jet Propulsion Laboratory + Pasadena, CA 91109 USA + + General site: https://ssd.jpl.nasa.gov/ + Mailing list: https://ssd.jpl.nasa.gov/email_list.html + System news : https://ssd.jpl.nasa.gov/horizons/news.html + User Guide : https://ssd.jpl.nasa.gov/horizons/manual.html + Connect : browser https://ssd.jpl.nasa.gov/horizons/app.html#/x + API https://ssd-api.jpl.nasa.gov/doc/horizons.html + command-line telnet ssd.jpl.nasa.gov 6775 + e-mail/batch https://ssd.jpl.nasa.gov/ftp/ssd/hrzn_batch.txt + scripts https://ssd.jpl.nasa.gov/ftp/ssd/SCRIPTS + Author : Jon.D.Giorgini@jpl.nasa.gov + +******************************************************************************************************************************************************************************************************************************************** diff --git a/sbpy/data/tests/data/TestEphemerisCLI-ceres-miriade.txt b/sbpy/data/tests/data/TestEphemerisCLI-ceres-miriade.txt new file mode 100644 index 000000000..efb5e91a1 --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-ceres-miriade.txt @@ -0,0 +1,128 @@ + + + + + + +Miriade - Positional ephemeris + +Positional ephemeris + + + + + + + + + + + + + + + + + + +Name or designation of the target + + +Julian day (UTC) + + +Right ascension of target centre + + +Declination of target centre + + +Distance from observer + + +Apparent visual magnitude + + +Phase angle + + +Solar elongation + + +Rate of change of target center RA (multiplied by the cosine of the declination) + + +Rate of change of target center DE + + +Rate of change of target distance from observer + + + +Ceres 2460538.5000000000018 34 29.78990-30 54 55.9994 2.131085077 8.13 14.51133.67-0.15128E+00-0.21151E-01 18.0099122 +Ceres 2460539.5000000000018 34 13.67853-30 55 24.4291 2.141569346 8.15 14.75132.66-0.13659E+00-0.18358E-01 18.2946586 +Ceres 2460540.5000000000018 33 59.21049-30 55 48.8930 2.152216073 8.17 14.99131.67-0.12189E+00-0.15643E-01 18.5725087 +Ceres 2460541.5000000000018 33 46.38663-30 56 9.5031 2.163021335 8.19 15.21130.67-0.10719E+00-0.13005E-01 18.8436706 +Ceres 2460542.5000000000018 33 35.20658-30 56 26.3679 2.173981332 8.21 15.44129.69-0.92507E-01-0.10441E-01 19.1083672 +Ceres 2460543.5000000000018 33 25.66914-30 56 39.5923 2.185092391 8.23 15.65128.71-0.77835E-01-0.79488E-02 19.3668030 +Ceres 2460544.5000000000018 33 17.77260-30 56 49.2777 2.196350944 8.24 15.86127.73-0.63182E-01-0.55253E-02 19.6191289 +Ceres 2460545.5000000000018 33 11.51502-30 56 55.5210 2.207753487 8.26 16.07126.76-0.48548E-01-0.31679E-02 19.8654184 +Ceres 2460546.5000000000018 33 6.89430-30 56 58.4156 2.219296535 8.28 16.27125.79-0.33934E-01-0.87394E-03 20.1056634 +Ceres 2460547.5000000000018 33 3.90815-30 56 58.0509 2.230976578 8.30 16.46124.83-0.19342E-01 0.13591E-02 20.3397896 +Ceres 2460548.5000000000018 33 2.55391-30 56 54.5130 2.242790050 8.32 16.65123.87-0.47768E-02 0.35338E-02 20.5676815 +Ceres 2460549.5000000000018 33 2.82834-30 56 47.8841 2.254733312 8.34 16.83122.92 0.97571E-02 0.56523E-02 20.7892072 +Ceres 2460550.5000000000018 33 4.72748-30 56 38.2435 2.266802650 8.35 17.01121.97 0.24252E-01 0.77171E-02 21.0042374 +Ceres 2460551.5000000000018 33 8.24648-30 56 25.6669 2.278994279 8.37 17.18121.03 0.38701E-01 0.97303E-02 21.2126564 +Ceres 2460552.5000000000018 33 13.37950-30 56 10.2270 2.291304350 8.39 17.34120.10 0.53095E-01 0.11694E-01 21.4143677 +Ceres 2460553.5000000000018 33 20.11969-30 55 51.9928 2.303728967 8.41 17.50119.17 0.67425E-01 0.13612E-01 21.6092958 +Ceres 2460554.5000000000018 33 28.45912-30 55 31.0303 2.316264196 8.42 17.66118.24 0.81681E-01 0.15484E-01 21.7973868 +Ceres 2460555.5000000000018 33 38.38881-30 55 7.4019 2.328906077 8.44 17.80117.32 0.95855E-01 0.17314E-01 21.9786087 +Ceres 2460556.5000000000018 33 49.89871-30 54 41.1668 2.341650639 8.46 17.95116.40 0.10994E+00 0.19105E-01 22.1529508 +Ceres 2460557.5000000000018 34 2.97775-30 54 12.3807 2.354493911 8.47 18.08115.49 0.12392E+00 0.20858E-01 22.3204238 +Ceres 2460558.5000000000018 34 17.61388-30 53 41.0959 2.367431934 8.49 18.21114.59 0.13779E+00 0.22576E-01 22.4810584 +Ceres 2460559.5000000000018 34 33.79410-30 53 7.3614 2.380460770 8.51 18.34113.69 0.15155E+00 0.24261E-01 22.6349035 +Ceres 2460560.5000000000018 34 51.50452-30 52 31.2228 2.393576517 8.52 18.46112.79 0.16518E+00 0.25915E-01 22.7820241 +Ceres 2460561.5000000000018 35 10.73046-30 51 52.7223 2.406775312 8.54 18.57111.90 0.17868E+00 0.27542E-01 22.9224993 +Ceres 2460562.5000000000018 35 31.45645-30 51 11.8988 2.420053345 8.56 18.68111.01 0.19204E+00 0.29142E-01 23.0564210 +Ceres 2460563.5000000000018 35 53.66632-30 50 28.7881 2.433406859 8.57 18.79110.13 0.20527E+00 0.30719E-01 23.1838941 +Ceres 2460564.5000000000018 36 17.34326-30 49 43.4227 2.446832163 8.59 18.88109.25 0.21834E+00 0.32274E-01 23.3050386 +Ceres 2460565.5000000000018 36 42.46981-30 48 55.8321 2.460325643 8.60 18.98108.38 0.23126E+00 0.33810E-01 23.4199940 +Ceres 2460566.5000000000018 37 9.02798-30 48 6.0426 2.473883768 8.62 19.07107.51 0.24402E+00 0.35328E-01 23.5289255 +Ceres 2460567.5000000000018 37 36.99926-30 47 14.0778 2.487503117 8.63 19.15106.65 0.25661E+00 0.36832E-01 23.6320297 +Ceres 2460568.5000000000018 38 6.36477-30 46 19.9580 2.501180386 8.65 19.23105.79 0.26905E+00 0.38322E-01 23.7295369 +Ceres 2460569.5000000000018 38 37.10541-30 45 23.7007 2.514912419 8.66 19.30104.94 0.28132E+00 0.39801E-01 23.8217054 +Ceres 2460570.5000000000018 39 9.20215-30 44 25.3200 2.528696208 8.68 19.37104.09 0.29342E+00 0.41271E-01 23.9088025 +Ceres 2460571.5000000000018 39 42.63624-30 43 24.8269 2.542528897 8.69 19.43103.24 0.30536E+00 0.42735E-01 23.9910720 +Ceres 2460572.5000000000018 40 17.38954-30 42 22.2291 2.556407754 8.70 19.49102.40 0.31715E+00 0.44195E-01 24.0686973 +Ceres 2460573.5000000000018 40 53.44470-30 41 17.5308 2.570330130 8.72 19.55101.57 0.32880E+00 0.45653E-01 24.1417736 +Ceres 2460574.5000000000018 41 30.78517-30 40 10.7332 2.584293400 8.73 19.60100.73 0.34031E+00 0.47111E-01 24.2103021 +Ceres 2460575.5000000000018 42 9.39512-30 39 1.8348 2.598294922 8.75 19.64 99.90 0.35168E+00 0.48571E-01 24.2742095 +Ceres 2460576.5000000000018 42 49.25920-30 37 50.8315 2.612331992 8.76 19.68 99.08 0.36293E+00 0.50035E-01 24.3333801 +Ceres 2460577.5000000000018 43 30.36232-30 36 37.7171 2.626401839 8.77 19.72 98.26 0.37404E+00 0.51504E-01 24.3876884 +Ceres 2460578.5000000000018 44 12.68942-30 35 22.4837 2.640501622 8.78 19.75 97.44 0.38503E+00 0.52978E-01 24.4370212 +Ceres 2460579.5000000000018 44 56.22534-30 34 5.1217 2.654628438 8.80 19.78 96.62 0.39589E+00 0.54460E-01 24.4812886 +Ceres 2460580.5000000000018 45 40.95475-30 32 45.6200 2.668779343 8.81 19.80 95.81 0.40662E+00 0.55951E-01 24.5204262 +Ceres 2460581.5000000000018 46 26.86205-30 31 23.9663 2.682951361 8.82 19.82 95.01 0.41721E+00 0.57450E-01 24.5543927 +Ceres 2460582.5000000000018 47 13.93140-30 30 0.1469 2.697141498 8.83 19.83 94.20 0.42767E+00 0.58959E-01 24.5831663 +Ceres 2460583.5000000000018 48 2.14671-30 28 34.1470 2.711346753 8.85 19.84 93.41 0.43799E+00 0.60478E-01 24.6067418 +Ceres 2460584.5000000000018 48 51.49160-30 27 5.9507 2.725564124 8.86 19.85 92.61 0.44817E+00 0.62010E-01 24.6251282 +Ceres 2460585.5000000000018 49 41.94947-30 25 35.5412 2.739790623 8.87 19.85 91.82 0.45820E+00 0.63553E-01 24.6383488 +Ceres 2460586.5000000000018 50 33.50347-30 24 2.9009 2.754023273 8.88 19.85 91.03 0.46810E+00 0.65108E-01 24.6464401 +Ceres 2460587.5000000000018 51 26.13650-30 22 28.0112 2.768259128 8.89 19.84 90.24 0.47784E+00 0.66677E-01 24.6494531 +Ceres 2460588.5000000000018 52 19.83126-30 20 50.8531 2.782495272 8.90 19.83 89.46 0.48744E+00 0.68259E-01 24.6474522 +Ceres 2460589.5000000000018 53 14.57024-30 19 11.4067 2.796728831 8.92 19.82 88.68 0.49689E+00 0.69855E-01 24.6405158 +Ceres 2460590.5000000000018 54 10.33576-30 17 29.6519 2.810956980 8.93 19.80 87.90 0.50618E+00 0.71466E-01 24.6287353 +Ceres 2460591.5000000000018 55 7.10998-30 15 45.5679 2.825176951 8.94 19.78 87.13 0.51533E+00 0.73091E-01 24.6122157 +Ceres 2460592.5000000000018 56 4.87494-30 13 59.1336 2.839386041 8.95 19.76 86.36 0.52432E+00 0.74730E-01 24.5910763 +Ceres 2460593.5000000000018 57 3.61258-30 12 10.3277 2.853581620 8.96 19.73 85.60 0.53315E+00 0.76385E-01 24.5654529 +Ceres 2460594.5000000000018 58 3.30476-30 10 19.1285 2.867761142 8.97 19.70 84.83 0.54183E+00 0.78054E-01 24.5355018 +Ceres 2460595.5000000000018 59 3.93333-30 08 25.5139 2.881922158 8.98 19.66 84.07 0.55036E+00 0.79740E-01 24.5014029 +Ceres 2460596.5000000000019 00 5.48019-30 06 29.4618 2.896062331 8.99 19.62 83.31 0.55873E+00 0.81440E-01 24.4633616 +Ceres 2460597.5000000000019 01 7.92741-30 04 30.9495 2.910179448 9.00 19.58 82.56 0.56696E+00 0.83157E-01 24.4216046 +Ceres 2460598.5000000000019 02 11.25738-30 02 29.9536 2.924271431 9.01 19.54 81.81 0.57503E+00 0.84890E-01 24.3763655 + + + + + + diff --git a/sbpy/data/tests/data/TestEphemerisCLI-ceres-mpc.txt b/sbpy/data/tests/data/TestEphemerisCLI-ceres-mpc.txt new file mode 100644 index 000000000..755414dfc --- /dev/null +++ b/sbpy/data/tests/data/TestEphemerisCLI-ceres-mpc.txt @@ -0,0 +1,117 @@ + + + + + Minor Planet Ephemeris Service: Query Results + + +

Minor Planet Ephemeris Service: Query Results

+ Below are the results of your request from the Minor Planet Center's + Minor Planet Ephemeris Service. +

+ Newly designated objects may take up to 1 day to show up in this service. +

+

+ Orbits and ephemerides of unnumbered NEOs are up-to-date. Other orbits are in + the process of being refreshed. +

+

+ The current system is not completely reliable in the case of objects + with very-close approaches with the Earth. +

We are working on a completely new system, + but for the time being we encourage the users to double check + the results with other ephemeris generators when the object is very close to + Earth. +

+ Ephemerides are for + the geocenter. +

+

+ (1) Ceres +

Number of variant orbits available: 11

+

Perturbed ephemeris below is based on + 125-opp + unperturbed elements from + MPO 840811. + Last observed on 2022 Sept. 16. +

Discovery date : 1801 01 01 +
Discovery site : Palermo +
Discoverer(s) : Piazzi, G. + 

+00001              [H= 3.33]
+Date       UT      R.A. (J2000) Decl.    Delta     r     El.    Ph.   V      Sky Motion       Uncertainty info
+            h m s                                                            "/hr     P.A.    3-sig/" P.A.
+2024 08 16 000000 18 34 29.8 -30 54 56   2.131   2.924  133.7  14.5   8.1    9.17    262.0       N/A   N/A / Map / Offsets
+2024 08 17 000000 18 34 13.7 -30 55 25   2.142   2.924  132.7  14.8   8.1    8.27    262.3       N/A   N/A / Map / Offsets
+2024 08 18 000000 18 33 59.2 -30 55 49   2.152   2.925  131.7  15.0   8.2    7.38    262.7       N/A   N/A / Map / Offsets
+2024 08 19 000000 18 33 46.4 -30 56 10   2.163   2.925  130.7  15.2   8.2    6.48    263.1       N/A   N/A / Map / Offsets
+2024 08 20 000000 18 33 35.2 -30 56 26   2.174   2.926  129.7  15.4   8.2    5.59    263.6       N/A   N/A / Map / Offsets
+2024 08 21 000000 18 33 25.7 -30 56 40   2.185   2.926  128.7  15.7   8.2    4.70    264.2       N/A   N/A / Map / Offsets
+2024 08 22 000000 18 33 17.8 -30 56 49   2.196   2.927  127.7  15.9   8.2    3.81    265.0       N/A   N/A / Map / Offsets
+2024 08 23 000000 18 33 11.5 -30 56 56   2.208   2.927  126.8  16.1   8.3    2.92    266.3       N/A   N/A / Map / Offsets
+2024 08 24 000000 18 33 06.9 -30 56 59   2.219   2.928  125.8  16.3   8.3    2.04    268.5       N/A   N/A / Map / Offsets
+2024 08 25 000000 18 33 03.9 -30 56 58   2.231   2.928  124.8  16.5   8.3    1.17    274.0       N/A   N/A / Map / Offsets
+2024 08 26 000000 18 33 02.6 -30 56 55   2.243   2.929  123.9  16.6   8.3    0.36    306.4       N/A   N/A / Map / Offsets
+2024 08 27 000000 18 33 02.9 -30 56 48   2.255   2.929  122.9  16.8   8.3    0.68    059.8       N/A   N/A / Map / Offsets
+2024 08 28 000000 18 33 04.8 -30 56 38   2.267   2.930  122.0  17.0   8.3    1.53    072.3       N/A   N/A / Map / Offsets
+2024 08 29 000000 18 33 08.3 -30 56 26   2.279   2.930  121.0  17.2   8.4    2.39    075.9       N/A   N/A / Map / Offsets
+2024 08 30 000000 18 33 13.4 -30 56 10   2.291   2.931  120.1  17.3   8.4    3.26    077.6       N/A   N/A / Map / Offsets
+2024 08 31 000000 18 33 20.2 -30 55 52   2.304   2.931  119.2  17.5   8.4    4.13    078.6       N/A   N/A / Map / Offsets
+2024 09 01 000000 18 33 28.5 -30 55 31   2.316   2.932  118.2  17.7   8.4    4.99    079.3       N/A   N/A / Map / Offsets
+2024 09 02 000000 18 33 38.4 -30 55 07   2.329   2.932  117.3  17.8   8.4    5.84    079.8       N/A   N/A / Map / Offsets
+2024 09 03 000000 18 33 49.9 -30 54 41   2.342   2.933  116.4  17.9   8.4    6.69    080.1       N/A   N/A / Map / Offsets
+2024 09 04 000000 18 34 03.0 -30 54 12   2.355   2.933  115.5  18.1   8.5    7.54    080.4       N/A   N/A / Map / Offsets
+2024 09 05 000000 18 34 17.7 -30 53 41   2.367   2.934  114.6  18.2   8.5    8.38    080.7       N/A   N/A / Map / Offsets
+2024 09 06 000000 18 34 33.8 -30 53 07   2.380   2.934  113.7  18.3   8.5    9.21    080.9       N/A   N/A / Map / Offsets
+2024 09 07 000000 18 34 51.5 -30 52 31   2.394   2.935  112.8  18.5   8.5   10.03    081.1       N/A   N/A / Map / Offsets
+2024 09 08 000000 18 35 10.8 -30 51 53   2.407   2.935  111.9  18.6   8.5   10.85    081.2       N/A   N/A / Map / Offsets
+2024 09 09 000000 18 35 31.5 -30 51 12   2.420   2.936  111.0  18.7   8.5   11.65    081.4       N/A   N/A / Map / Offsets
+2024 09 10 000000 18 35 53.7 -30 50 29   2.433   2.936  110.1  18.8   8.6   12.45    081.5       N/A   N/A / Map / Offsets
+2024 09 11 000000 18 36 17.4 -30 49 43   2.447   2.937  109.3  18.9   8.6   13.24    081.6       N/A   N/A / Map / Offsets
+2024 09 12 000000 18 36 42.5 -30 48 56   2.460   2.937  108.4  19.0   8.6   14.02    081.7       N/A   N/A / Map / Offsets
+2024 09 13 000000 18 37 09.1 -30 48 06   2.474   2.938  107.5  19.1   8.6   14.79    081.8       N/A   N/A / Map / Offsets
+2024 09 14 000000 18 37 37.0 -30 47 14   2.488   2.938  106.7  19.1   8.6   15.55    081.8       N/A   N/A / Map / Offsets
+2024 09 15 000000 18 38 06.4 -30 46 20   2.501   2.939  105.8  19.2   8.6   16.30    081.9       N/A   N/A / Map / Offsets
+2024 09 16 000000 18 38 37.2 -30 45 24   2.515   2.939  104.9  19.3   8.7   17.05    081.9       N/A   N/A / Map / Offsets
+2024 09 17 000000 18 39 09.3 -30 44 25   2.529   2.940  104.1  19.4   8.7   17.78    082.0       N/A   N/A / Map / Offsets
+2024 09 18 000000 18 39 42.7 -30 43 25   2.543   2.940  103.2  19.4   8.7   18.50    082.0       N/A   N/A / Map / Offsets
+2024 09 19 000000 18 40 17.4 -30 42 22   2.556   2.941  102.4  19.5   8.7   19.21    082.1       N/A   N/A / Map / Offsets
+2024 09 20 000000 18 40 53.5 -30 41 18   2.570   2.941  101.6  19.5   8.7   19.92    082.1       N/A   N/A / Map / Offsets
+2024 09 21 000000 18 41 30.8 -30 40 11   2.584   2.942  100.7  19.6   8.7   20.61    082.1       N/A   N/A / Map / Offsets
+2024 09 22 000000 18 42 09.5 -30 39 02   2.598   2.942   99.9  19.6   8.7   21.30    082.1       N/A   N/A / Map / Offsets
+2024 09 23 000000 18 42 49.3 -30 37 51   2.612   2.943   99.1  19.7   8.7   21.98    082.1       N/A   N/A / Map / Offsets
+2024 09 24 000000 18 43 30.4 -30 36 38   2.626   2.943   98.3  19.7   8.8   22.65    082.2       N/A   N/A / Map / Offsets
+2024 09 25 000000 18 44 12.7 -30 35 23   2.641   2.943   97.4  19.7   8.8   23.32    082.2       N/A   N/A / Map / Offsets
+2024 09 26 000000 18 44 56.3 -30 34 05   2.655   2.944   96.6  19.8   8.8   23.98    082.2       N/A   N/A / Map / Offsets
+2024 09 27 000000 18 45 41.0 -30 32 46   2.669   2.944   95.8  19.8   8.8   24.62    082.2       N/A   N/A / Map / Offsets
+2024 09 28 000000 18 46 26.9 -30 31 24   2.683   2.945   95.0  19.8   8.8   25.27    082.2       N/A   N/A / Map / Offsets
+2024 09 29 000000 18 47 14.0 -30 30 00   2.697   2.945   94.2  19.8   8.8   25.90    082.1       N/A   N/A / Map / Offsets
+2024 09 30 000000 18 48 02.2 -30 28 34   2.711   2.946   93.4  19.8   8.8   26.53    082.1       N/A   N/A / Map / Offsets
+2024 10 01 000000 18 48 51.6 -30 27 06   2.726   2.946   92.6  19.8   8.8   27.14    082.1       N/A   N/A / Map / Offsets
+2024 10 02 000000 18 49 42.0 -30 25 36   2.740   2.947   91.8  19.8   8.9   27.75    082.1       N/A   N/A / Map / Offsets
+2024 10 03 000000 18 50 33.6 -30 24 03   2.754   2.947   91.0  19.8   8.9   28.35    082.1       N/A   N/A / Map / Offsets
+2024 10 04 000000 18 51 26.2 -30 22 28   2.768   2.947   90.2  19.8   8.9   28.95    082.1       N/A   N/A / Map / Offsets
+2024 10 05 000000 18 52 19.9 -30 20 51   2.783   2.948   89.5  19.8   8.9   29.53    082.0       N/A   N/A / Map / Offsets
+2024 10 06 000000 18 53 14.6 -30 19 11   2.797   2.948   88.7  19.8   8.9   30.10    082.0       N/A   N/A / Map / Offsets
+2024 10 07 000000 18 54 10.4 -30 17 30   2.811   2.949   87.9  19.8   8.9   30.67    082.0       N/A   N/A / Map / Offsets
+2024 10 08 000000 18 55 07.2 -30 15 46   2.825   2.949   87.1  19.8   8.9   31.23    081.9       N/A   N/A / Map / Offsets
+2024 10 09 000000 18 56 04.9 -30 13 59   2.839   2.950   86.4  19.8   8.9   31.77    081.9       N/A   N/A / Map / Offsets
+2024 10 10 000000 18 57 03.7 -30 12 10   2.854   2.950   85.6  19.7   8.9   32.31    081.8       N/A   N/A / Map / Offsets
+2024 10 11 000000 18 58 03.4 -30 10 19   2.868   2.950   84.8  19.7   9.0   32.84    081.8       N/A   N/A / Map / Offsets
+2024 10 12 000000 18 59 04.0 -30 08 26   2.882   2.951   84.1  19.7   9.0   33.36    081.8       N/A   N/A / Map / Offsets
+2024 10 13 000000 19 00 05.6 -30 06 29   2.896   2.951   83.3  19.6   9.0   33.88    081.7       N/A   N/A / Map / Offsets
+2024 10 14 000000 19 01 08.0 -30 04 31   2.910   2.952   82.6  19.6   9.0   34.38    081.7       N/A   N/A / Map / Offsets
+2024 10 15 000000 19 02 11.3 -30 02 30   2.924   2.952   81.8  19.5   9.0   34.87    081.6       N/A   N/A / Map / Offsets
+
+

+ These calculations have been performed on the + Tamkin + Foundation Computing Network. +

+

+ Valid HTML 4.01! +

+ + diff --git a/sbpy/data/tests/test_ephem.py b/sbpy/data/tests/test_ephem.py index 9c33dd72c..80bbc2d66 100644 --- a/sbpy/data/tests/test_ephem.py +++ b/sbpy/data/tests/test_ephem.py @@ -1,35 +1,171 @@ # Licensed under a 3-clause BSD style license - see LICENSE.rst +import os import pytest +from urllib.parse import urlencode +from types import SimpleNamespace import astropy.units as u from astropy.time import Time -from ... import exceptions as sbe +try: + import astroquery + from astroquery.utils.mocks import MockResponse + from astroquery.query import BaseQuery +except ImportError: # pragma: no cover + astroquery = None + + from ... import bib -from .. import Ephem, Orbit +from ..orbit import Orbit +from ..ephem import Ephem # retreived from Horizons on 23 Apr 2020 CERES = { - 'targetname': '1 Ceres', - 'H': u.Quantity(3.4, 'mag'), - 'G': 0.12, - 'e': 0.07741102040801928, - 'q': u.Quantity(2.55375156, 'au'), - 'incl': u.Quantity(10.58910839, 'deg'), - 'Omega': u.Quantity(80.29081558, 'deg'), - 'w': u.Quantity(73.7435117, 'deg'), - 'n': u.Quantity(0.21401711, 'deg / d'), - 'M': u.Quantity(154.70418799, 'deg'), - 'nu': u.Quantity(158.18663933, 'deg'), - 'a': u.Quantity(2.76802739, 'AU'), - 'Q': u.Quantity(2.98230321, 'AU'), - 'P': u.Quantity(1682.10848349, 'd'), - 'epoch': Time(2458963.26397076, scale='tdb', format='jd'), - 'Tp': Time(2458240.40500675, scale='tdb', format='jd') + "targetname": "1 Ceres", + "H": u.Quantity(3.4, "mag"), + "G": 0.12, + "e": 0.07741102040801928, + "q": u.Quantity(2.55375156, "au"), + "incl": u.Quantity(10.58910839, "deg"), + "Omega": u.Quantity(80.29081558, "deg"), + "w": u.Quantity(73.7435117, "deg"), + "n": u.Quantity(0.21401711, "deg / d"), + "M": u.Quantity(154.70418799, "deg"), + "nu": u.Quantity(158.18663933, "deg"), + "a": u.Quantity(2.76802739, "AU"), + "Q": u.Quantity(2.98230321, "AU"), + "P": u.Quantity(1682.10848349, "d"), + "epoch": Time(2458963.26397076, scale="tdb", format="jd"), + "Tp": Time(2458240.40500675, scale="tdb", format="jd"), } +def data_path(filename): + data_dir = os.path.join(os.path.dirname(__file__), "data") + return os.path.join(data_dir, filename) + + +def nonremote_request(self, request_type, url, **kwargs): + """monkeypatch replacement astroquery request function""" + + params = urlencode(kwargs.get("params", {})) + data = urlencode(kwargs.get("data", {})) + + key = (("URL", url), ("params", params), ("data", data)) + + # URL to name translation, files should be saved to the ./data directory. To + # add a new test URL-file pair, design the test as usual below and run the + # test with the --remote-data and --capture=tee-sys parameters. The error + # should give enough info to formulate the key. + # curl "URL" -o file.txt + # curl "URL" --data "DATA" -o file.txt + files = { + ( + ("URL", "https://ssd.jpl.nasa.gov/api/horizons.api"), + ( + "params", + "format=text&EPHEM_TYPE=OBSERVER&QUANTITIES=%271%2C3%2C9%2C19%2C20%2C23%2C24%2C27%2C33%27&" + "COMMAND=%22DES%3DC%2F1995+O1%3B+CAP%3C2460568%3B+NOFRAG%3B%22&SOLAR_ELONG=%220%2C180%22&" + "LHA_CUTOFF=0&CSV_FORMAT=YES&CAL_FORMAT=BOTH&ANG_FORMAT=DEG&APPARENT=AIRLESS&REF_SYSTEM=ICRF&" + "EXTRA_PREC=NO&CENTER=%27500%27&START_TIME=%222024-08-16+00%3A00%3A00.000%22&" + "STOP_TIME=%222024-10-15+00%3A00%3A00.000%22&STEP_SIZE=%221d%22&SKIP_DAYLT=NO", + ), + ("data", ""), + ): "TestEphemerisCLI-c1995o1-horizons.txt", + ( + ("URL", "https://ssd.jpl.nasa.gov/api/horizons.api"), + ( + "params", + "format=text&EPHEM_TYPE=OBSERVER&QUANTITIES=%271%2C3%2C9%2C19%2C20%2C23%2C24%2C27%2C33%27&" + "COMMAND=%22DES%3D2P%3B+CAP%3C2460568%3B+NOFRAG%3B%22&SOLAR_ELONG=%220%2C180%22&LHA_CUTOFF=0&" + "CSV_FORMAT=YES&CAL_FORMAT=BOTH&ANG_FORMAT=DEG&APPARENT=AIRLESS&REF_SYSTEM=ICRF&EXTRA_PREC=NO&" + "CENTER=%27500%27&START_TIME=%222024-08-16+00%3A00%3A00.000%22&" + "STOP_TIME=%222024-10-15+00%3A00%3A00.000%22&STEP_SIZE=%221d%22&SKIP_DAYLT=NO", + ), + ("data", ""), + ): "TestEphemerisCLI-2p-horizons.txt", + ( + ("URL", "https://cgi.minorplanetcenter.net/cgi-bin/mpeph2.cgi"), + ("params", ""), + ( + "data", + "ty=e&TextArea=2P&uto=0&igd=n&ibh=n&fp=y&adir=N&tit=&bu=&c=500&d=2024-08-16+000000&i=1&u=d&l=61&" + "raty=a&s=t&m=h", + ), + ): "TestEphemerisCLI-2p-mpc.txt", + ( + ("URL", "http://vo.imcce.fr/webservices/miriade/ephemcc_query.php"), + ( + "params", + "-name=2P&-type=Comet&-ep=2460538.5&-step=1.000000d&-nbd=61.0&-observer=500&-output=--jul&-tscale=UTC&" + "-theory=INPOP&-teph=1&-tcoor=1&-rplane=1&-oscelem=ASTORB&-mime=votable", + ), + ("data", ""), + ): "TestEphemerisCLI-2p-miriade.txt", + ( + ("URL", "http://vo.imcce.fr/webservices/miriade/ephemcc_query.php"), + ( + "params", + "-name=1&-type=Asteroid&-ep=2460538.5&-step=1.000000d&-nbd=61.0&-observer=500&-output=--jul&" + "-tscale=UTC&-theory=INPOP&-teph=1&-tcoor=1&-rplane=1&-oscelem=ASTORB&-mime=votable", + ), + ("data", ""), + ): "TestEphemerisCLI-ceres-miriade.txt", + ( + ("URL", "https://cgi.minorplanetcenter.net/cgi-bin/mpeph2.cgi"), + ("params", ""), + ( + "data", + "ty=e&TextArea=1&uto=0&igd=n&ibh=n&fp=y&adir=N&tit=&bu=&c=500&d=2024-08-16+000000&i=1&u=d&l=61&raty=a&" + "s=t&m=h", + ), + ): "TestEphemerisCLI-ceres-mpc.txt", + ( + ("URL", "https://ssd.jpl.nasa.gov/api/horizons.api"), + ( + "params", + "format=text&EPHEM_TYPE=OBSERVER&QUANTITIES=%271%2C3%2C9%2C19%2C20%2C23%2C24%2C27%2C33%27&" + "COMMAND=%221%3B%22&SOLAR_ELONG=%220%2C180%22&LHA_CUTOFF=0&CSV_FORMAT=YES&CAL_FORMAT=BOTH&" + "ANG_FORMAT=DEG&APPARENT=AIRLESS&REF_SYSTEM=ICRF&EXTRA_PREC=NO&CENTER=%27500%27&" + "START_TIME=%222024-08-16+00%3A00%3A00.000%22&STOP_TIME=%222024-10-15+00%3A00%3A00.000%22&" + "STEP_SIZE=%221d%22&SKIP_DAYLT=NO", + ), + ("data", ""), + ): "TestEphemerisCLI-ceres-horizons.txt", + } + + try: + fn = files[key] + except KeyError: # pragma: no cover + # printing this out helps with debugging + print("url:", url + (("?" + params) if len(params) > 0 else "")) + print("data:", data) + print("files key:", key) + raise KeyError("Request does not have corresponding mocked data") + + with open(data_path(fn), "rb") as inf: + response = MockResponse(content=inf.read(), url=url) + + # for MPC module + response.request = SimpleNamespace(body=data) + + return response + + +# use a pytest fixture to create a dummy 'requests.get' function, +# that mocks (monkeypatches) the actual 'requests.get' function: +@pytest.fixture +def patch_request(request): + if astroquery is None: + return None + + mp = request.getfixturevalue("monkeypatch") + + mp.setattr(BaseQuery, "_request", nonremote_request) + return mp + + class TestEphemFromOorb: def test_units(self): pytest.importorskip("pyoorb") @@ -37,44 +173,57 @@ def test_units(self): orbit1 = Orbit.from_dict(CERES) eph1 = Ephem.from_oo(orbit1) - orbit2 = Orbit.from_dict({ - 'targetname': orbit1['targetname'][0], - 'a': orbit1['a'].value[0] * u.au, - 'e': orbit1['e'][0], - 'i': orbit1['i'].value[0] * u.deg, - 'w': orbit1['w'].value[0] * u.deg, - 'Omega': orbit1['Omega'].value[0] * u.deg, - 'epoch': Time(orbit1['epoch'][0], format='jd'), - 'M': orbit1['M'].value[0] * u.deg, - 'H': orbit1['H'].value[0] * u.mag, - 'G': orbit1['G'][0] - }) + orbit2 = Orbit.from_dict( + { + "targetname": orbit1["targetname"][0], + "a": orbit1["a"].value[0] * u.au, + "e": orbit1["e"][0], + "i": orbit1["i"].value[0] * u.deg, + "w": orbit1["w"].value[0] * u.deg, + "Omega": orbit1["Omega"].value[0] * u.deg, + "epoch": Time(orbit1["epoch"][0], format="jd"), + "M": orbit1["M"].value[0] * u.deg, + "H": orbit1["H"].value[0] * u.mag, + "G": orbit1["G"][0], + } + ) eph2 = Ephem.from_oo(orbit2) - for k in ['ra', 'dec', 'RA*cos(Dec)_rate', 'dec_rate', 'alpha', 'r', - 'delta', 'V', 'hlon', 'hlat', 'EL']: + for k in [ + "ra", + "dec", + "RA*cos(Dec)_rate", + "dec_rate", + "alpha", + "r", + "delta", + "V", + "hlon", + "hlat", + "EL", + ]: assert u.isclose(eph1[k], eph2[k]) def test_basic(self): pytest.importorskip("pyoorb") orbit = Orbit.from_dict(CERES) - oo_ephem = Ephem.from_oo(orbit, scope='basic') - assert 'dec_rate' not in oo_ephem.field_names + oo_ephem = Ephem.from_oo(orbit, scope="basic") + assert "dec_rate" not in oo_ephem.field_names def test_timescale(self): pytest.importorskip("pyoorb") orbit = Orbit.from_dict(CERES) epoch = Time.now() - oo_ephem = Ephem.from_oo(orbit, epochs=epoch, scope='basic') - assert oo_ephem['epoch'].scale == epoch.scale + oo_ephem = Ephem.from_oo(orbit, epochs=epoch, scope="basic") + assert oo_ephem["epoch"].scale == epoch.scale def test_bib(self): pytest.importorskip("pyoorb") with bib.Tracking(): orbit = Orbit.from_dict(CERES) - oo_ephem = Ephem.from_oo(orbit, scope='basic') - assert 'sbpy.data.ephem.Ephem.from_oo' in bib.show() + Ephem.from_oo(orbit, scope="basic") + assert "sbpy.data.ephem.core.Ephem.from_oo" in bib.show() bib.reset() diff --git a/sbpy/data/tests/test_ephem_cli.py b/sbpy/data/tests/test_ephem_cli.py new file mode 100644 index 000000000..07c1bac1e --- /dev/null +++ b/sbpy/data/tests/test_ephem_cli.py @@ -0,0 +1,120 @@ +# Licensed under a 3-clause BSD style license - see LICENSE.rst + +import pytest +from types import SimpleNamespace + +import astropy.units as u +from astropy.time import Time + +from ..ephem.cli import EphemerisCLI +from .test_ephem import patch_request # noqa: F401 + + +class TestEphemCLI: + def test_format_epochs(self): + """EphemerisCLI.parse_args will always have start and step. Test + number, stop, and neither.""" + + epochs = { + "start": Time("2024-07-23"), + "step": u.Quantity(5, "day"), + } + epochs["stop"] = epochs["start"] + 10 * epochs["step"] + + # test stop + args = SimpleNamespace( + start=epochs["start"], step=epochs["step"], number=None, stop=epochs["stop"] + ) + + result = EphemerisCLI._format_epochs(args) + for k in epochs.keys(): + assert result[k] == epochs[k] + + # test number + args = SimpleNamespace( + start=epochs["start"], step=epochs["step"], number=10, stop=None + ) + + result = EphemerisCLI._format_epochs(args) + for k in epochs.keys(): + assert result[k] == epochs[k] + + # no number, no stop --> number = 60 + args = SimpleNamespace( + start=epochs["start"], step=epochs["step"], number=None, stop=None + ) + + result = EphemerisCLI._format_epochs(args) + epochs["stop"] = epochs["start"] + 60 * epochs["step"] + for k in epochs.keys(): + assert result[k] == epochs[k] + + def test_radec_format(self, patch_request): # noqa: F811 + """Test RA Dec formatting options.""" + + pytest.importorskip("astroquery") + + cli = EphemerisCLI( + ["horizons", "C/1995 O1", "--start=2024-08-16", "--radec=deg"] + ) + + row = str(cli.eph.table[0]).splitlines()[-1].split() + assert row[2] == "339.94076" + assert row[3] == "-85.76646" + + cli = EphemerisCLI( + ["horizons", "C/1995 O1", "--start=2024-08-16", "--radec=hmsdms"] + ) + row = str(cli.eph.table[0]).splitlines()[-1].split() + assert row[2] == "22:39:45.78" + assert row[3] == "-85:45:59.3" + + def test_start_stop_order(self, patch_request): # noqa: F811 + """Test start > stop.""" + + pytest.importorskip("astroquery") + + with pytest.raises(ValueError): + EphemerisCLI( + [ + "horizons", + "C/1995 O1", + "--start=2025-01-01", + "--stop=2024-01-01", + "--debug", + ] + ) + + @pytest.mark.parametrize( + "service, target", + [("horizons", "1 Ceres (A801 AA)"), ("mpc", "1"), ("miriade", "Ceres")], + ) + def test_asteroid(self, service, target, patch_request): # noqa: F811 + """Test services with a comet designation""" + + pytest.importorskip("astroquery") + + cmd = f"{service} 1 --start=2024-08-16" + if service == "horizons": + cmd += " --id-type=smallbody" + cli = EphemerisCLI(cmd.split()) + + assert cli.eph["date"][0].iso == "2024-08-16 00:00:00.000" + assert cli.eph.meta["target"] == target + + @pytest.mark.parametrize( + "service, target", + [("horizons", "2P/Encke"), ("mpc", "2P"), ("miriade", "2P")], + ) + def test_comet(self, service, target, patch_request): # noqa: F811 + """Test services with a comet designation""" + + pytest.importorskip("astroquery") + + cmd = f"{service} 2P --start=2024-08-16" + if service == "miriade": + cmd += " --type=comet" + cli = EphemerisCLI(cmd.split()) + + assert cli.eph["date"][0].iso == "2024-08-16 00:00:00.000" + assert cli.eph.meta["target"] == target diff --git a/sbpy/data/tests/test_ephem_cli_remote.py b/sbpy/data/tests/test_ephem_cli_remote.py new file mode 100644 index 000000000..f02a3a73d --- /dev/null +++ b/sbpy/data/tests/test_ephem_cli_remote.py @@ -0,0 +1,26 @@ +# Licensed under a 3-clause BSD style license - see LICENSE.rst + +import pytest +from ..ephem.cli import EphemerisCLI + +pytest.importorskip("astroquery") + + +@pytest.mark.remote_data() +class TestEphemCLI: + @pytest.mark.parametrize("service", ("horizons", "mpc", "miriade")) + def test_services(self, service): + """Spot check results from the services.""" + cli = EphemerisCLI( + f"{service} ceres --start=2024-01-01 --step=3d --stop=2024-01-31".split() + ) + assert len(cli.eph) == 11 + assert cli.eph["date"][0].iso == "2024-01-01 00:00:00.000" + assert cli.eph["date"][-1].iso == "2024-01-31 00:00:00.000" + + cli = EphemerisCLI( + f"{service} ceres --start=2024-01-01 --step=3d --number=10".split() + ) + assert len(cli.eph) == 11 + assert cli.eph["date"][0].iso == "2024-01-01 00:00:00.000" + assert cli.eph["date"][-1].iso == "2024-01-31 00:00:00.000" diff --git a/sbpy/data/tests/test_ephem_remote.py b/sbpy/data/tests/test_ephem_remote.py index 71a06766c..94a1ff136 100644 --- a/sbpy/data/tests/test_ephem_remote.py +++ b/sbpy/data/tests/test_ephem_remote.py @@ -1,7 +1,6 @@ # Licensed under a 3-clause BSD style license - see LICENSE.rst import pytest -from numpy import abs import warnings from numpy.testing import assert_allclose @@ -12,7 +11,9 @@ from astropy.tests.helper import assert_quantity_allclose from ... import bib -from .. import Ephem, Orbit, QueryError +from ..core import QueryError +from ..ephem import Ephem +from ..orbit import Orbit pytest.importorskip("astroquery") @@ -97,7 +98,7 @@ def test_Earthlocation(self): def test_bib(self): with bib.Tracking(): Ephem.from_horizons(['Ceres', 'Pallas']) - assert 'sbpy.data.ephem.Ephem.from_horizons' in bib.show() + assert 'sbpy.data.ephem.core.Ephem.from_horizons' in bib.show() bib.reset() def test_timescale(self): @@ -299,7 +300,7 @@ def test_queryfail(self): def test_bib(self): with bib.Tracking(): Ephem.from_mpc(['Ceres', 'Pallas']) - assert 'sbpy.data.ephem.Ephem.from_mpc' in bib.show() + assert 'sbpy.data.ephem.core.Ephem.from_mpc' in bib.show() bib.reset() @@ -398,8 +399,8 @@ def test_queryfail(self): def test_bib(self): with bib.Tracking(): - data = Ephem.from_miriade(['Ceres', 'Pallas']) - assert 'sbpy.data.ephem.Ephem.from_miriade' in bib.show() + Ephem.from_miriade(['Ceres', 'Pallas']) + assert 'sbpy.data.ephem.core.Ephem.from_miriade' in bib.show() bib.reset() diff --git a/sbpy/data/tests/test_obs_remote.py b/sbpy/data/tests/test_obs_remote.py index 8359b01ab..05db7513b 100644 --- a/sbpy/data/tests/test_obs_remote.py +++ b/sbpy/data/tests/test_obs_remote.py @@ -67,7 +67,7 @@ def test_jplhorizons(self): assert len(data.field_names) > len(obs.field_names) assert 'targetname_eph' in data.field_names - assert 'sbpy.data.ephem.Ephem.from_horizons' in bib.to_text() + assert 'sbpy.data.ephem.core.Ephem.from_horizons' in bib.to_text() def test_mpc(self): bib.track() @@ -77,7 +77,7 @@ def test_mpc(self): data = obs.supplement(service='mpc') assert len(data.field_names) > len(obs.field_names) - assert 'sbpy.data.ephem.Ephem.from_mpc' in bib.to_text() + assert 'sbpy.data.ephem.core.Ephem.from_mpc' in bib.to_text() def test_miriade(self): bib.track() @@ -87,7 +87,7 @@ def test_miriade(self): data = obs.supplement(service='miriade') assert len(data.field_names) > len(obs.field_names) - assert 'sbpy.data.ephem.Ephem.from_miriade' in bib.to_text() + assert 'sbpy.data.ephem.core.Ephem.from_miriade' in bib.to_text() def test_breaks(self): obs = Obs.from_dict({'epoch': Time([2451200, 2451201], format='jd'), diff --git a/sbpy/dynamics/state.py b/sbpy/dynamics/state.py index 2bae822e3..770f0c274 100644 --- a/sbpy/dynamics/state.py +++ b/sbpy/dynamics/state.py @@ -388,7 +388,7 @@ def from_ephem( Parameters ---------- - eph : ~sbpy.data.ephem.Ephem + eph : ~sbpy.data.ephem.core.Ephem Ephemeris object, must have time, position, and velocity. Position and velocity may be specified using ("x", "y", "z", "vx", "vy", and "vz"), or ("ra", "dec", "Delta", "RA*cos(Dec)_rate", "Dec_rate", diff --git a/sbpy/photometry/iau.py b/sbpy/photometry/iau.py index c8356f78b..6ad246cc0 100644 --- a/sbpy/photometry/iau.py +++ b/sbpy/photometry/iau.py @@ -249,7 +249,7 @@ def phasecoeff(self): @property def oe_amp(self): - """Opposition effect amplitude, :math:`\zeta-1` + r"""Opposition effect amplitude, :math:`\zeta-1` Based on Muinonen et al. (2010) Eq. 24) """ tmp = float(self._G1+self._G2) diff --git a/setup.cfg b/setup.cfg index 7e3cca17d..cea01f962 100644 --- a/setup.cfg +++ b/setup.cfg @@ -59,7 +59,7 @@ test = coverage docs = sbpy[all,test] - matplotlib>=3.1 + matplotlib>=3.3 sphinx-astropy>=1.3 astropy!=6.1.* scipy>=1.6 @@ -72,6 +72,9 @@ docs = ginga.rv.plugins = cometaryenhancements = sbpy.ginga_plugins:setup_cometaryenhancements +console_scripts = + sbpy-ephem = sbpy.data.ephem.cli:EphemerisCLI.run + [tool:pytest] minversion = 7.0 testpaths = "sbpy" "docs" diff --git a/tox.ini b/tox.ini index 465397dea..3283f2b8a 100644 --- a/tox.ini +++ b/tox.ini @@ -64,7 +64,7 @@ deps = # The oldestdeps factor is intended to be used to install the oldest versions of all # dependencies that have a minimum version. oldestdeps: numpy==1.21.* - oldestdeps: matplotlib==3.1.* + oldestdeps: matplotlib==3.3.* oldestdeps: scipy==1.6.* oldestdeps: synphot==1.1.* oldestdeps: astropy==5.3.*