Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/nutation.c

Timestamp:

Jun 15, 2004, 6:54:12 PM (21 years ago)

Author:

cmv

Message:

nouvelle version de xephem/libastro (3.6) cmv 15/6/04

File:

• trunk/SophyaExt/XephemAstroLib/nutation.c (modified) (7 diffs)

trunk/SophyaExt/XephemAstroLib/nutation.c

@@ -3 +3 @@
  * (given to 0.001") EXACTLY over 750 days (1995 and 1996)
  */
+#include <stdio.h>
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 
@@ -268 +268 @@
 };
 
-/* given the modified JD, mjd, find the nutation in obliquity, *deps, and
+/* given the modified JD, mj, find the nutation in obliquity, *deps, and
  * the nutation in longitude, *dpsi, each in radians.
  */
 void
-nutation (mjd, deps, dpsi)
-double mjd;
-double *deps;   /* on input:  precision parameter in arc seconds */
-double *dpsi;
+nutation (
+double mj,
+double *deps,   /* on input:  precision parameter in arc seconds */
+double *dpsi)
 {
-        static double lastmjd = -10000, lastdeps, lastdpsi;
+        static double lastmj = -10000, lastdeps, lastdpsi;
         double T, T2, T3, T10;                  /* jul cent since J2000 */
         double prec;                            /* series precis in arc sec */
@@ -287 +287 @@
                          */
 
-        if (mjd == lastmjd) {
+        if (mj == lastmj) {
             *deps = lastdeps;
             *dpsi = lastdpsi;
@@ -304 +304 @@
         prec *= NUT_SCALE/10;
 
-        T = (mjd - J2000)/36525.;
+        T = (mj - J2000)/36525.;
         T2 = T * T;
         T3 = T2 * T;
@@ -362 +362 @@
         lastdeps = degrad(lastdeps/3600./NUT_SCALE);
 
-        lastmjd = mjd;
+        lastmj = mj;
         *deps = lastdeps;
         *dpsi = lastdpsi;
 }
 
-/* given the modified JD, mjd, correct, IN PLACE, the right ascension *ra
+/* given the modified JD, mj, correct, IN PLACE, the right ascension *ra
  * and declination *dec (both in radians) for nutation.
  */
 void
-nut_eq (mjd, ra, dec)
-double mjd, *ra, *dec;
+nut_eq (double mj, double *ra, double *dec)
 {
-        static double lastmjd = -10000;
+        static double lastmj = -10000;
         static double a[3][3];          /* rotation matrix */
         double xold, yold, zold, x, y, z;
 
-        if (mjd != lastmjd) {
+        if (mj != lastmj) {
             double epsilon, dpsi, deps;
             double se, ce, sp, cp, sede, cede;
 
-            obliquity(mjd, &epsilon);
-            nutation(mjd, &deps, &dpsi);
+            obliquity(mj, &epsilon);
+            nutation(mj, &deps, &dpsi);
 
             /* the rotation matrix a applies the nutation correction to
@@ -428 +427 @@
             a[2][2] = sede*cp*se+cede*ce;
 
-            lastmjd = mjd;
+            lastmj = mj;
         }
 
@@ -440 +439 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2004-06-15 16:52:39 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};