astronomia/es/sidereal.js

/**
 * @copyright 2013 Sonia Keys
 * @copyright 2016 commenthol
 * @license MIT
 * @module sidereal
 */
/**
 * Sidereal: Chapter 12, Sidereal Time at Greenwich.
 */

import base from './base';
import nutation from './nutation';

/**
 * JDToCFrac returns values for use in computing sidereal time at Greenwich.
 *
 * Cen is centuries from J2000 of the JD at 0h UT of argument jd.  This is
 * the value to use for evaluating the IAU sidereal time polynomial.
 * DayFrac is the fraction of jd after 0h UT.  It is used to compute the
 * final value of sidereal time.
 *
 * @param {number} jd - Julian Days
 * @return {number[]} [century, fraction] century and fraction of jd after 0h UT
 */
export function JDToCFrac(jd) {
  var _base$modf = base.modf(jd + 0.5),
      j0 = _base$modf[0],
      f = _base$modf[1];

  return [base.J2000Century(j0 - 0.5), f]; // (cen, dayFrac /* float */)
}

/**
 * iau82 is a polynomial giving mean sidereal time at Greenwich at 0h UT.
 *
 * The polynomial is in centuries from J2000.0, as given by JDToCFrac.
 * Coefficients are those adopted in 1982 by the International Astronomical
 * Union and are given in (12.2) p. 87.
 */
export var iau82 = [24110.54841, 8640184.812866, 0.093104, 0.0000062];

/**
 * Mean returns mean sidereal time at Greenwich for a given JD.
 *
 * Computation is by IAU 1982 coefficients.  The result is in seconds of
 * time and is in the range [0,86400).
 *
 * @param {number} jd - Julian Days
 * @return {number}
 */
export function mean(jd) {
  return base.pmod(_mean(jd), 86400);
}

/**
 * @private
 */
function _mean(jd) {
  var _mean0UT2 = _mean0UT(jd),
      s = _mean0UT2[0],
      f = _mean0UT2[1];

  return s + f * 1.00273790935 * 86400;
}

/**
 * Mean0UT returns mean sidereal time at Greenwich at 0h UT on the given JD.
 *
 * The result is in seconds of time and is in the range [0,86400).
 *
 * @param {number} jd - Julian Days
 * @return {number}
 */
export function mean0UT(jd /* float */) {
  var _mean0UT3 = _mean0UT(jd),
      s = _mean0UT3[0],
      _ = _mean0UT3[1]; // eslint-disable-line


  return base.pmod(s, 86400);
}

/**
 * @private
 */
function _mean0UT(jd /* float */) {
  var _JDToCFrac = JDToCFrac(jd),
      cen = _JDToCFrac[0],
      f = _JDToCFrac[1];
  // (12.2) p. 87


  return [base.horner.apply(base, [cen].concat(iau82)), f]; // (sidereal, dayFrac /* float */)
}

/**
 * Apparent returns apparent sidereal time at Greenwich for the given JD.
 *
 * Apparent is mean plus the nutation in right ascension.
 *
 * The result is in seconds of time and is in the range [0,86400).
 *
 * @param {number} jd - Julian Days
 * @return {number}
 */
export function apparent(jd) {
  var s = _mean(jd); // seconds of time
  var n = nutation.nutationInRA(jd); // angle (radians) of RA
  var ns = n * 3600 * 180 / Math.PI / 15; // convert RA to time in seconds
  return base.pmod(s + ns, 86400);
}

/**
 * Apparent0UT returns apparent sidereal time at Greenwich at 0h UT
 * on the given JD.
 *
 * The result is in seconds of time and is in the range [0,86400).
 *
 * @param {number} jd - Julian Days
 * @return {number}
 */
export function apparent0UT(jd) {
  var _base$modf2 = base.modf(jd + 0.5),
      j0 = _base$modf2[0],
      f = _base$modf2[1];

  var cen = (j0 - 0.5 - base.J2000) / 36525;
  var s = base.horner.apply(base, [cen].concat(iau82)) + f * 1.00273790935 * 86400;
  var n = nutation.nutationInRA(j0); // angle (radians) of RA
  var ns = n * 3600 * 180 / Math.PI / 15; // convert RA to time in seconds
  return base.pmod(s + ns, 86400);
}

export default {
  JDToCFrac: JDToCFrac,
  iau82: iau82,
  mean: mean,
  mean0UT: mean0UT,
  apparent: apparent,
  apparent0UT: apparent0UT
};