Changeset 3111 in Sophya for trunk/SophyaExt/XephemAstroLib/deltat.c

Timestamp:

Nov 22, 2006, 2:53:31 PM (19 years ago)

Author:

cmv

Message:

mise en conformite xephem 3.7.2 cmv 22/11/2006

File:

• trunk/SophyaExt/XephemAstroLib/deltat.c (modified) (3 diffs)

trunk/SophyaExt/XephemAstroLib/deltat.c

@@ -144 +144 @@
     /* 1999 USNO data 1998.0 thru 2000.0 (McBurnett) */
      6297, 6347, 6383, 
-    /* 1999 extrapolation (McBurnett), 2001.0 thru 2006.0
-       Ramp up to 1.6 s/yr to transition smoothly to Stevenson formula
-       in 2130.0 */
-     6440, 6510, 6600, 6750, 6900, 7060
+    /* 1999 extrapolation (McBurnett), 2001.0 thru 2006.0 */
+     /* 6440, 6510, 6600, 6750, 6900, 7060 */
+     6409, 6430, 6447, 6507, 6578, 6610 /* ECD */
 
     /* original 1997 USNO extrapolation (stern), 1998.0 thru 2004.0
@@ -178 +177 @@
         Y = 2000.0 + (mj - J2000)/365.25;
 
-        if( Y > TABEND  &&  Y < 2130.0 ) {
+        if( Y > TABEND) {
             /* linear interpolation from table end; stern */
             B = Y - TABEND;
-            ans = dt[TABSIZ-1] + B * (dt[TABSIZ-1]  - dt[TABSIZ-2]);
+            ans = dt[TABSIZ-1] + B * (dt[TABSIZ-1]  - dt[TABSIZ-11])/10;
             ans *= 0.01;
             return(ans);
         }
 
-        if( Y < TABSTART  ||  Y >= 2130.0 ) {
-            if( (Y >= 948.0 - 15.0  &&  Y < TABSTART) ||  Y >= 2130.0 ) {
+        if( Y < TABSTART) {
+            if( Y >= 948.0 - 15.0 ) {
                 /* Stephenson and Morrison, stated domain is 948 to 1600:
                  * 25.5(centuries from 1800)^2 - 1.9159(centuries from 1955)^2
@@ -304 +303 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2005-08-21 10:02:37 $ $Revision: 1.5 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2006-11-22 13:53:29 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};