source: Sophya/trunk/SophyaExt/XephemAstroLib/obliq.c@ 1518

Last change on this file since 1518 was 1457, checked in by cmv, 25 years ago

import de la partie libastro de Xephem cmv+rz 10/4/2001
File size: 778 bytes
Line 
1#include <stdio.h>
2
3#include "P_.h"
4#include "astro.h"
5
6/* given the modified Julian date, mjd, find the mean obliquity of the
7 * ecliptic, *eps, in radians.
8 *
9 * IAU expression (see e.g. Astron. Almanac 1984); stern
10 */
11void
12obliquity (mjd, eps)
13double mjd;
14double *eps;
15{
16 static double lastmjd = -16347, lasteps;
17
18 if (mjd != lastmjd) {
19 double t = (mjd - J2000)/36525.; /* centuries from J2000 */
20 lasteps = degrad(23.4392911 + /* 23^ 26' 21".448 */
21 t * (-46.8150 +
22 t * ( -0.00059 +
23 t * ( 0.001813 )))/3600.0);
24 lastmjd = mjd;
25 }
26 *eps = lasteps;
27}
28
29/* For RCS Only -- Do Not Edit */
30static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: obliq.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
Note: See TracBrowser for help on using the repository browser.