Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/deltat.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/deltat.c (modified) (6 diffs)

trunk/SophyaExt/XephemAstroLib/deltat.c

@@ -29 +29 @@
  * in 2130.
  *
- * Input is mjd (modified julian date from MJD0 on). [stern]
- *  Note that xephem uses a different epoch for this "mjd" than the
+ * Input is mj (modified julian date from MJD0 on). [stern]
+ *  Note that xephem uses a different epoch for this "mj" than the
  *  normal value of JD=240000.5.
  * See AA page B4.
  *
- * Output double deltat(mjd) is ET-UT1 in seconds.
+ * Output double deltat(mj) is ET-UT1 in seconds.
  *
  *
@@ -67 +67 @@
  *   - adopted #include's for xephem
  *   - made dt[] static
- *   - made mjd the time argument [was: year Y].
+ *   - made mj the time argument [was: year Y].
  *   - updated observed and extrapolated data from tables at
  *      ftp://maia.usno.navy.mil/ser7/ -- data deviated by up to 0.8 s
@@ -74 +74 @@
  *   - replaced treatment after TABEND by linear extrapolation instead
  *      of second order version
- *   - installed lastmjd cache (made ans static)
+ *   - installed lastmj cache (made ans static)
  *
  *   - no changes to table interpolation scheme and past extrapolations */
 
-#include "P_.h"
+#include <math.h>
+
 #include "astro.h"
 
@@ -161 +162 @@
  * of the Earth rotation rate in the ET time scale.
  */
-double deltat(mjd)
-double mjd;
+double deltat(double mj)
 {
         double Y;
@@ -168 +168 @@
         int d[6];
         int i, iy, k;
-        double floor();
         static double ans;
-        static double lastmjd = -10000;
-
-        if (mjd == lastmjd) {
+        static double lastmj = -10000;
+
+        if (mj == lastmj) {
             return(ans);
         }
-        lastmjd = mjd;
-
-        Y = 2000.0 + (mjd - J2000)/365.25;
+        lastmj = mj;
+
+        Y = 2000.0 + (mj - J2000)/365.25;
 
         if( Y > TABEND  &&  Y < 2130.0 ) {
@@ -305 +304 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2004-06-15 16:52:38 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};