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

Last change on this file since 3508 was 3477, checked in by cmv, 18 years ago

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