Context Navigation

obliq.c@ 1558

Visit:

Last change on this file since 1558 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

Rev	Line
[1457]	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	*/
	11	void
	12	obliquity (mjd, eps)
	13	double mjd;
	14	double *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 */
	30	static 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.