view release on metacpan or  search on metacpan

palsrc/palEtrms.c  view on Meta::CPAN

*     constants the latest value is used here.
*
*  See also:
*     - Smith, C.A. et al., 1989.  Astr.J. 97, 265.
*     - Yallop, B.D. et al., 1989.  Astr.J. 97, 274.

*  Authors:
*     PTW: Pat Wallace (STFC)
*     TIMJ: Tim Jenness (JAC, Hawaii)
*     {enter_new_authors_here}

*  History:
*     2012-02-12 (TIMJ):
*        Initial version with documentation taken from Fortran SLA
*        Adapted with permission from the Fortran SLALIB library.
*     {enter_further_changes_here}

*  Copyright:
*     Copyright (C) 1996 Rutherford Appleton Laboratory
*     Copyright (C) 2012 Science and Technology Facilities Council.
*     All Rights Reserved.

*  Licence:
*     This program is free software: you can redistribute it and/or
*     modify it under the terms of the GNU Lesser General Public
*     License as published by the Free Software Foundation, either
*     version 3 of the License, or (at your option) any later
*     version.
*
*     This program is distributed in the hope that it will be useful,
*     but WITHOUT ANY WARRANTY; without even the implied warranty of
*     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*     GNU Lesser General Public License for more details.
*
*     You should have received a copy of the GNU Lesser General
*     License along with this program.  If not, see
*     <http://www.gnu.org/licenses/>.

*  Bugs:
*     {note_any_bugs_here}
*-
*/

#include "pal.h"
#include "palmac.h"

void palEtrms ( double ep, double ev[3] ) {

  /* Use the J2000 aberration constant */
  const double ABCONST = 20.49552;

  double t, e, e0, p, ek, cp;

  /*  Julian centuries since B1950 */
  t = (ep - 1950.) * .0100002135903;

  /*  Eccentricity */
  e = .01673011 - (t * 1.26e-7 + 4.193e-5) * t;

  /*  Mean obliquity */
  e0 = (84404.836 - ((t * .00181 + .00319) * t + 46.8495) * t) *
    PAL__DAS2R;

  /*  Mean longitude of perihelion */
  p = (((t * .012 + 1.65) * t + 6190.67) * t + 1015489.951) *
    PAL__DAS2R;

  /*  E-terms */
  ek = e * ABCONST * PAL__DAS2R;
  cp = cos(p);
  ev[0] = ek * sin(p);
  ev[1] = -ek * cp * cos(e0);
  ev[2] = -ek * cp * sin(e0);

}