Changeset 4017 in Sophya


Ignore:
Timestamp:
Sep 21, 2011, 6:17:52 PM (14 years ago)
Author:
cmv
Message:

fichiers de Xephem 3.7.5 update, cmv 21/09/2011

Location:
trunk/SophyaExt/XephemAstroLib
Files:
1 deleted
65 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaExt/XephemAstroLib/aa_hadec.c

    r3654 r4017  
    7575
    7676/* For RCS Only -- Do Not Edit */
    77 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aa_hadec.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     77static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aa_hadec.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/aberration.c

    r3654 r4017  
    159159
    160160/* For RCS Only -- Do Not Edit */
    161 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aberration.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     161static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aberration.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/actan.c

    r3654 r4017  
    11#include <math.h>
    22
    3 /* @(#) $Id: actan.c,v 1.8 2009-07-16 10:34:35 cmv Exp $ */
     3/* @(#) $Id: actan.c,v 1.9 2011-09-21 16:17:47 cmv Exp $ */
    44
    55/* commonly in math.h, but not in strict ANSI C */
    … …  
    6565
    6666/* For RCS Only -- Do Not Edit */
    67 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: actan.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     67static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: actan.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/airmass.c

    r3654 r4017  
    2424
    2525/* For RCS Only -- Do Not Edit */
    26 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: airmass.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     26static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: airmass.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/anomaly.c

    r3654 r4017  
    6161
    6262/* For RCS Only -- Do Not Edit */
    63 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: anomaly.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     63static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: anomaly.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/ap_as.c

    r3654 r4017  
    6666
    6767/* For RCS Only -- Do Not Edit */
    68 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     68static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/astro.h

    r3654 r4017  
    808808
    809809/* For RCS Only -- Do Not Edit
    810  * @(#) $RCSfile: astro.h,v $ $Date: 2009-07-16 10:34:35 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     810 * @(#) $RCSfile: astro.h,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    811811 */

  • trunk/SophyaExt/XephemAstroLib/atlas.c

    r3654 r4017  
    194194
    195195/* For RCS Only -- Do Not Edit */
    196 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: atlas.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     196static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: atlas.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/auxil.c

    r3654 r4017  
    4040
    4141/* For RCS Only -- Do Not Edit */
    42 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: auxil.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     42static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: auxil.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/bdl.c

    r3654 r4017  
    236236
    237237/* For RCS Only -- Do Not Edit */
    238 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: bdl.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     238static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: bdl.c,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/bdl.h

    r3654 r4017  
    33
    44/* For RCS Only -- Do Not Edit
    5  * @(#) $RCSfile: bdl.h,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $
     5 * @(#) $RCSfile: bdl.h,v $ $Date: 2011-09-21 16:17:47 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $
    66 */

  • trunk/SophyaExt/XephemAstroLib/chap95.c

    r3654 r4017  
    172172
    173173/* For RCS Only -- Do Not Edit */
    174 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     174static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/chap95.h

    r3654 r4017  
    6565
    6666/* For RCS Only -- Do Not Edit
    67  * @(#) $RCSfile: chap95.h,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     67 * @(#) $RCSfile: chap95.h,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    6868 */

  • trunk/SophyaExt/XephemAstroLib/chap95_data.c

    r3654 r4017  
    781781
    782782/* For RCS Only -- Do Not Edit */
    783 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95_data.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     783static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95_data.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/circum.c

    r3654 r4017  
    861861
    862862/* For RCS Only -- Do Not Edit */
    863 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: circum.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     863static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: circum.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/comet.c

    r3654 r4017  
    8181
    8282/* For RCS Only -- Do Not Edit */
    83 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     83static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/compare_with_xephem.csh

    r2818 r4017  
    99
    1010set dir = $SOPHYASOURCE/XephemAstroLib/
    11 set log = `pwd`/compare.log
     11set log = /tmp/compare.log
    1212
    1313rm -f $log; touch $log

  • trunk/SophyaExt/XephemAstroLib/constel.c

    r3654 r4017  
    18011801
    18021802/* For RCS Only -- Do Not Edit */
    1803 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: constel.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     1803static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: constel.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/dbfmt.c

    r3654 r4017  
    1313int get_fields (char *s, int delim, char *fields[]);
    1414
    15 #define MAXDBLINE       256     /* longest allowed db line */
     15#define MAXDBLINE       512     /* longest allowed db line */
    1616
    1717#define FLDSEP          ','     /* major field separator */
    … …  
    7474
    7575        /* do all the parsing on a copy */
    76         (void) strcpy (copy, s);
     76        (void) strncpy (copy, s, MAXDBLINE-1);
     77        copy[MAXDBLINE-1] = '\0';
    7778        i = strlen(copy);
    7879        if (copy[i-1] == '\n')
    … …  
    10261027
    10271028/* For RCS Only -- Do Not Edit */
    1028 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     1029static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/deep.c

    r3654 r4017  
    788788
    789789/* For RCS Only -- Do Not Edit */
    790 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deep.c,v $ $Date: 2009-07-16 10:34:36 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     790static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deep.c,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/deepconst.h

    r3654 r4017  
    22#define _CONST_H
    33
    4 /* $Id: deepconst.h,v 1.8 2009-07-16 10:34:37 cmv Exp $ */
     4/* $Id: deepconst.h,v 1.9 2011-09-21 16:17:48 cmv Exp $ */
    55
    66
    … …  
    3131
    3232/* For RCS Only -- Do Not Edit
    33  * @(#) $RCSfile: deepconst.h,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     33 * @(#) $RCSfile: deepconst.h,v $ $Date: 2011-09-21 16:17:48 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    3434 */

  • trunk/SophyaExt/XephemAstroLib/deltat.c

    r3654 r4017  
    11/* DeltaT = Ephemeris Time - Universal Time
    22 *
    3  * original code by Stephen Moshier <moshier@world.std.com>,
    4  * adapted to xephem by Michael Sternberg <sternberg@physik.tu-chemnitz.de>
    5  * smoothed transitions and updated by Neal McBurnett <nealmcb@bell-labs.com>
    6  *
    7  **********************************************************************
    8  *
    9  * The tabulated values of deltaT, in hundredths of a second,
    10  * were taken from The Astronomical Almanac, page K8.  The program
    11  * adjusts for a value of secular tidal acceleration ndot = -25.8
    12  * arcsec per century squared, the value used in JPL's DE403 ephemeris.
    13  * ELP2000 (and DE200) used the value -23.8946.
    14  *
    15  * The tabulated range is 1620.0 through 1998.0.  Bessel's interpolation
    16  * formula is implemented to obtain fourth order interpolated values at
    17  * intermediate times.
    18  *
    19  * For dates earlier than the tabulated range, the program
    20  * calculates approximate formulae of Stephenson and Morrison
    21  * or K. M. Borkowski.  These approximations have an estimated
    22  * error of 15 minutes at 1500 B.C.  They are not adjusted for small
    23  * improvements in the current estimate of ndot because the formulas
    24  * were derived from studies of ancient eclipses and other historical
    25  * information, whose interpretation depends only partly on ndot.
    26  *
    27  * For future vaues of deltaT, the function smoothly transitions with
    28  * a linear segment back to Stephenson & Morrison's quadratic formula
    29  * in 2130.
    30  *
    31  * Input is mj (modified julian date from MJD0 on). [stern]
    32  *  Note that xephem uses a different epoch for this "mj" than the
    33  *  normal value of JD=240000.5.
    34  * See AA page B4.
    35  *
    36  * Output double deltat(mj) is ET-UT1 in seconds.
     3 * Adapted 2011/4/14 from Stephen Moshier <moshier@world.std.com>,
     4 * cosmetic changes only.
     5 *
     6 * Compile as follows to create stand-alone test program:
     7 *   cc -DTEST_MAIN deltat.c libastro.a
     8 *
     9 * Tabulated values of deltaT, in hundredths of a second, are
     10 * from The Astronomical Almanac and current IERS reports.
     11 * A table of values for the pre-telescopic period was taken from
     12 * Morrison and Stephenson (2004).  The overall tabulated range is
     13 * -1000.0 through 2011.0.  Values at intermediate times are interpolated
     14 * from the tables.
     15 *
     16 * For dates earlier and later than the tabulated range, the program
     17 * calculates a polynomial extrapolation formula.
     18 *
     19 * Updated deltaT predictions can be obtained from this network archive,
     20 *    http://maia.usno.navy.mil
     21 * then appended to the dt[] table and update TABEND.
     22 *
     23 * Input is XEphem's MJD, output is ET-UT in seconds.
    3724 *
    3825 *
    3926 * References:
     27 *
     28 * Morrison, L. V., and F. R. Stephenson, Historical values of the Earth's
     29 * clock error deltat T and the calculation of eclipses. Journal for the
     30 * History of Astronomy 35, 327-336 (2004)
    4031 *
    4132 * Stephenson, F. R., and L. V. Morrison, "Long-term changes
    … …  
    4435 * Series A 313, 47-70 (1984)
    4536 *
    46  * Borkowski, K. M., "ELP2000-85 and the Dynamical Time
    47  * - Universal Time relation," Astronomy and Astrophysics
    48  * 205, L8-L10 (1988)
    49  * Borkowski's formula is derived from eclipses going back to 2137 BC
    50  * and uses lunar position based on tidal coefficient of -23.9 arcsec/cy^2.
    51  *
    5237 * Chapront-Touze, Michelle, and Jean Chapront, _Lunar Tables
    5338 * and Programs from 4000 B.C. to A.D. 8000_, Willmann-Bell 1991
    54  * Their table agrees with the one here, but the entries are
    55  * rounded to the nearest whole second.
    5639 *
    5740 * Stephenson, F. R., and M. A. Houlden, _Atlas of Historical
    5841 * Eclipse Maps_, Cambridge U. Press (1986)
    5942 *
    60  * from obsolete extrapolation code [stern]:
    61  * Morrison, L. V. and F. R. Stephenson, "Sun and Planetary System"
    62  * vol 96,73 eds. W. Fricke, G. Teleki, Reidel, Dordrecht (1982)
    63  *
    64  **********************************************************************
    65  *
    66  * changes by stern:
    67  *   - adopted #include's for xephem
    68  *   - made dt[] static
    69  *   - made mj the time argument [was: year Y].
    70  *   - updated observed and extrapolated data from tables at
    71  *      ftp://maia.usno.navy.mil/ser7/ -- data deviated by up to 0.8 s
    72  *   - removed references to "extern double dtgiven"
    73  *   - removed DEMO #define and its references
    74  *   - replaced treatment after TABEND by linear extrapolation instead
    75  *      of second order version
    76  *   - installed lastmj cache (made ans static)
    77  *
    78  *   - no changes to table interpolation scheme and past extrapolations */
     43 */
    7944
    8045#include <math.h>
    … …  
    8247#include "astro.h"
    8348
    84 #define TABSTART 1620.0
    85 #define TABEND 2017.0
    86 #define TABSIZ 398
    87 
    88 /* Note, Stephenson and Morrison's table starts at the year 1630.
    89  * The Chapronts' table does not agree with the Almanac prior to 1630.
    90  * The actual accuracy decreases rapidly prior to 1780.
    91  */
    92 static short dt[TABSIZ] = {
     49#define TABSTART 1620
     50#define TABEND 2011
     51#define TABSIZ (TABEND - TABSTART + 1)
     52
     53/* Morrison and Stephenson (2004)
     54 * This table covers -1000 through 1700 in 100-year steps.
     55 * Values are in whole seconds.
     56 * Estimated standard error at -1000 is 640 seconds; at 1600, 20 seconds.
     57 * The first value in the table has been adjusted 28 sec for
     58 * continuity with their long-term quadratic extrapolation formula.
     59 * The last value in this table agrees with the AA table at 1700,
     60 * so there is no discontinuity at either endpoint.
     61 */
     62#define MS_SIZ 28
     63short m_s[MS_SIZ] = {
     64    /* -1000 to -100 */
     65    25428, 23700, 22000, 21000, 19040, 17190, 15530, 14080, 12790, 11640,
     66
     67    /* 0 to 900 */
     68    10580, 9600, 8640, 7680, 6700, 5710, 4740, 3810, 2960, 2200,
     69
     70    /* 1000 to 1700 */
     71    1570, 1090, 740, 490, 320, 200, 120, 9,
     72};
     73
     74
     75/* Entries prior to 1955 in the following table are from
     76 * the 1984 Astronomical Almanac and assume ndot = -26.0.
     77 * For dates prior to 1700, the above table is used instead of this one.
     78 */
     79short dt[TABSIZ] = {
    9380    /* 1620.0 thru 1659.0 */
    9481    12400, 11900, 11500, 11000, 10600, 10200, 9800, 9500, 9100, 8800,
    … …  
    9683    6200, 6000, 5800, 5700, 5500, 5400, 5300, 5100, 5000, 4900,
    9784    4800, 4700, 4600, 4500, 4400, 4300, 4200, 4100, 4000, 3800,
     85
    9886    /* 1660.0 thru 1699.0 */
    9987    3700, 3600, 3500, 3400, 3300, 3200, 3100, 3000, 2800, 2700,
    … …  
    10189    1600, 1500, 1400, 1400, 1300, 1200, 1200, 1100, 1100, 1000,
    10290    1000, 1000, 900, 900, 900, 900, 900, 900, 900, 900,
     91
    10392    /* 1700.0 thru 1739.0 */
    10493    900, 900, 900, 900, 900, 900, 900, 900, 1000, 1000,
    … …  
    10695    1100, 1100, 1100, 1100, 1100, 1100, 1100, 1100, 1100, 1100,
    10796    1100, 1100, 1100, 1100, 1200, 1200, 1200, 1200, 1200, 1200,
     97
    10898    /* 1740.0 thru 1779.0 */
    10999    1200, 1200, 1200, 1200, 1300, 1300, 1300, 1300, 1300, 1300,
    … …  
    111101    1500, 1500, 1500, 1500, 1500, 1600, 1600, 1600, 1600, 1600,
    112102    1600, 1600, 1600, 1600, 1600, 1700, 1700, 1700, 1700, 1700,
     103
    113104    /* 1780.0 thru 1799.0 */
    114105    1700, 1700, 1700, 1700, 1700, 1700, 1700, 1700, 1700, 1700,
    115106    1700, 1700, 1600, 1600, 1600, 1600, 1500, 1500, 1400, 1400,
     107
    116108    /* 1800.0 thru 1819.0 */
    117109    1370, 1340, 1310, 1290, 1270, 1260, 1250, 1250, 1250, 1250,
    118110    1250, 1250, 1250, 1250, 1250, 1250, 1250, 1240, 1230, 1220,
     111
    119112    /* 1820.0 thru 1859.0 */
    120113    1200, 1170, 1140, 1110, 1060, 1020, 960, 910, 860, 800,
    … …  
    122115    570, 580, 590, 610, 620, 630, 650, 660, 680, 690,
    123116    710, 720, 730, 740, 750, 760, 770, 770, 780, 780,
     117
    124118    /* 1860.0 thru 1899.0 */
    125119    788, 782, 754, 697, 640, 602, 541, 410, 292, 182,
    … …  
    127121    -540, -542, -520, -546, -546, -579, -563, -564, -580, -566,
    128122    -587, -601, -619, -664, -644, -647, -609, -576, -466, -374,
     123
    129124    /* 1900.0 thru 1939.0 */
    130125    -272, -154, -2, 124, 264, 386, 537, 614, 775, 913,
    … …  
    132127    2116, 2225, 2241, 2303, 2349, 2362, 2386, 2449, 2434, 2408,
    133128    2402, 2400, 2387, 2395, 2386, 2393, 2373, 2392, 2396, 2402,
     129
    134130    /* 1940.0 thru 1979.0 */
    135      2433, 2483, 2530, 2570, 2624, 2677, 2728, 2778, 2825, 2871,
    136      2915, 2957, 2997, 3036, 3072, 3107, 3135, 3168, 3218, 3268,
    137      3315, 3359, 3400, 3447, 3503, 3573, 3654, 3743, 3829, 3920,
    138      4018, 4117, 4223, 4337, 4449, 4548, 4646, 4752, 4853, 4959,
    139     /* 1980.0 thru 1995.0 */
    140      5054, 5138, 5217, 5296, 5379, 5434, 5487, 5532, 5582, 5630,
    141      5686, 5757, 5831, 5912, 5998, 6078,
    142     /* new USNO data (stern) */
    143      6163, 6230,
    144     /* 1999 USNO data 1998.0 thru 2000.0 (McBurnett) */
    145      6297, 6347, 6383,
    146 
    147      6409, 6430, 6447, 6507, 6546, 6546, 6546,
    148 
    149     /* 2008+ added 5/19/2009 from http://maia.usno.navy.mil/ser7/deltat.preds */
    150 
    151     6546,
    152     6570,
    153     6650,
    154     6710,
    155     6800,
    156     6800,
    157     6900,
    158     6900,
    159     7000,
    160     7000,
     131    2433, 2483, 2530, 2570, 2624, 2677, 2728, 2778, 2825, 2871,
     132    2915, 2957, 2997, 3036, 3072, 3107, 3135, 3168, 3218, 3268,
     133    3315, 3359, 3400, 3447, 3503, 3573, 3654, 3743, 3829, 3920,
     134    4018, 4117, 4223, 4337, 4449, 4548, 4646, 4752, 4853, 4959,
     135
     136    /* 1980.0 thru 2011.0 */
     137    5054, 5138, 5217, 5296, 5379, 5434, 5487, 5532, 5582, 5630,
     138    5686, 5757, 5831, 5912, 5998, 6078, 6163, 6230, 6297, 6347,
     139    6383, 6409, 6430, 6447, 6457, 6469, 6485, 6515, 6546, 6578,
     140    6607, 6632,
    161141};
    162142
    163 /* calculate  DeltaT = ET - UT1 in seconds.  Describes the irregularities
     143
     144/* Given MJD return DeltaT = ET - UT1 in seconds.  Describes the irregularities
    164145 * of the Earth rotation rate in the ET time scale.
    165146 */
    166 double deltat(double mj)
     147double
     148deltat(double mj)
    167149{
    168         double Y;
    169         double p, B;
     150        static double ans, lastmj;
     151        double Y, p, B;
    170152        int d[6];
    171153        int i, iy, k;
    172         static double ans;
    173         static double lastmj = -10000;
    174 
    175         if (mj == lastmj) {
    176             return(ans);
     154
     155        if (mj == lastmj)
     156            return (ans);
     157        lastmj = mj;
     158
     159        mjd_year (mj, &Y);
     160
     161        if( Y > TABEND ) {
     162            /* Extrapolate future values beyond the lookup table.  */
     163            if (Y > (TABEND + 100.0)) {
     164                /* Morrison & Stephenson (2004) long-term curve fit.  */
     165                B = 0.01 * (Y - 1820.0);
     166                ans = 32.0 * B * B - 20.0;
     167
     168            } else {
     169
     170                double a, b, c, d, m0, m1;
     171
     172                /* Cubic interpolation between last tabulated value
     173                 * and long-term curve evaluated at 100 years later.
     174                 */
     175
     176                /* Last tabulated delta T value. */
     177                a = 0.01 * dt[TABSIZ-1];
     178                /* Approximate slope in past 10 years. */
     179                b = 0.001 * (dt[TABSIZ-1] - dt[TABSIZ - 11]);
     180
     181                /* Long-term curve 100 years hence. */
     182                B = 0.01 * (TABEND + 100.0 - 1820.0);
     183                m0 = 32.0 * B*B - 20.0;
     184                /* Its slope. */
     185                m1 = 0.64 * B;
     186
     187                /* Solve for remaining coefficients of an interpolation polynomial
     188                 * that agrees in value and slope at both ends of the 100-year
     189                 * interval.
     190                 */
     191                d = 2.0e-6 * (50.0 * (m1 + b) - m0 + a);
     192                c = 1.0e-4 * (m0 - a - 100.0 * b - 1.0e6 * d);
     193
     194                /* Note, the polynomial coefficients do not depend on Y.
     195                 * A given tabulation and long-term formula
     196                 * determine the polynomial.
     197                 * Thus, for the IERS table ending at 2011.0, the coefficients are
     198                 * a = 66.32
     199                 * b = 0.223
     200                 * c = 0.03231376
     201                 * d = -0.0001607784
     202                 */
     203
     204                /* Compute polynomial value at desired time. */
     205                p = Y - TABEND;
     206                ans = a + p * (b  + p * (c + p * d));
     207            }
     208
     209            return (ans);
    177210        }
    178         lastmj = mj;
    179 
    180         Y = 2000.0 + (mj - J2000)/365.25;
    181 
    182         if( Y > TABEND) {
    183             /* linear interpolation from table end; stern */
    184             B = Y - TABEND;
    185             ans = dt[TABSIZ-1] + B * (dt[TABSIZ-1]  - dt[TABSIZ-11])/10;
    186             ans *= 0.01;
    187             return(ans);
     211
     212
     213        /* Use Morrison and Stephenson (2004) prior to the year 1700.  */
     214        if( Y < 1700.0 ) {
     215            if (Y <= -1000.0) {
     216                /* Morrison and Stephenson long-term fit.  */
     217                B = 0.01 * (Y - 1820.0);
     218                ans = 32.0 * B * B - 20.0;
     219
     220            } else {
     221
     222                /* Morrison and Stephenson recommend linear interpolation
     223                 * between tabulations.
     224                 */
     225                iy = Y;
     226                iy = (iy + 1000) / 100;  /* Integer index into the table. */
     227                B = -1000 + 100 * iy;    /* Starting year of tabulated interval.  */
     228                p = m_s[iy];
     229                ans = p + 0.01 * (Y - B) * (m_s[iy + 1] - p);
     230            }
     231
     232            return (ans);
    188233        }
    189234
    190         if( Y < TABSTART) {
    191             if( Y >= 948.0 - 15.0 ) {
    192                 /* Stephenson and Morrison, stated domain is 948 to 1600:
    193                  * 25.5(centuries from 1800)^2 - 1.9159(centuries from 1955)^2
    194                  * Here we offset by -15 y to minimize the discontinuity,
    195                  * thus we use it from 933.0 to 1620.0,
    196                  * and from the end of the table to 2130.0.
    197                  * f(1620.0) = 60.955200, slope -0.079 s/y
    198                  * f(2004.0) = 105.649728, slope 1.02 s/y
    199                  * f(2048.0) = 155.176, slope 1.23 s/y
    200                  * f(2084.0) = 202.49, slope 1.4 s/y
    201                  * f(2130.0) = 272, slope .1616
    202                  * f(2150.0) = 305, slope .17
    203                  */
    204                 B = 0.01*(Y - 2000.0);
    205                 ans = (23.58 * B + 100.3)*B + 101.6;
    206             } else {
    207                 /* Borkowski */
    208                 /* f(2004.0) = 542.7435, slope 2.65 s/y */
    209                 B = 0.01*(Y - 2000.0)  +  3.75;
    210                 ans = 35.0 * B * B  +  40.;
    211             }
    212             return(ans);
    213         }
    214 
    215         /* Besselian interpolation from tabulated values.
     235        /* Besselian interpolation between tabulated values
     236         * in the telescopic era.
    216237         * See AA page K11.
    217238         */
    218239
    219         /* value for 1620.1 is 121.96 or so, not 124.0 */
    220 
    221         /* Index into the table.
    222          */
     240        /* Index into the table.  */
    223241        p = floor(Y);
    224         iy = (int)(p - TABSTART);
    225         /* Zeroth order estimate is value at start of year
    226          */
     242        iy = (int) (p - TABSTART);
     243        /* Zeroth order estimate is value at start of year */
    227244        ans = dt[iy];
    228245        k = iy + 1;
    … …  
    230247            goto done; /* No data, can't go on. */
    231248
    232         /* The fraction of tabulation interval
    233          */
     249        /* The fraction of tabulation interval */
    234250        p = Y - p;
    235251
    236         /* First order interpolated value
    237          */
     252        /* First order interpolated value */
    238253        ans += p*(dt[k] - dt[iy]);
    239254        if( (iy-1 < 0) || (iy+2 >= TABSIZ) )
    240255            goto done; /* can't do second differences */
    241256
    242         /* Make table of first differences
    243          */
     257        /* Make table of first differences */
    244258        k = iy - 2;
    245         for( i=0; i<5; i++ ) {
     259        for (i=0; i<5; i++) {
    246260            if( (k < 0) || (k+1 >= TABSIZ) )
    247261                d[i] = 0;
    248             else d[i] = dt[k+1] - dt[k];
    249                 k += 1;
     262            else
     263                d[i] = dt[k+1] - dt[k];
     264            k += 1;
    250265        }
    251266
    252         /* Compute second differences
    253          */
     267        /* Compute second differences */
    254268        for( i=0; i<4; i++ )
    255269            d[i] = d[i+1] - d[i];
    256270        B = 0.25*p*(p-1.0);
    257271        ans += B*(d[1] + d[2]);
    258         if( iy+2 >= TABSIZ )
     272        if (iy+2 >= TABSIZ)
    259273            goto done;
    260274
    261         /* Compute third differences
    262          */
     275        /* Compute third differences */
    263276        for( i=0; i<3; i++ )
    264277            d[i] = d[i+1] - d[i];
    265278        B = 2.0*B/3.0;
    266279        ans += (p-0.5)*B*d[1];
    267         if( (iy-2 < 0) || (iy+3 > TABSIZ) )
     280        if ((iy-2 < 0) || (iy+3 > TABSIZ) )
    268281            goto done;
    269282
    270         /* Compute fourth differences
    271          */
     283        /* Compute fourth differences */
    272284        for( i=0; i<2; i++ )
    273285            d[i] = d[i+1] - d[i];
    … …  
    275287        ans += B*(d[0] + d[1]);
    276288
    277         done:
     289    done:
     290
     291        ans *= 0.01;
     292
     293#if 0 /* ndot = -26.0 assumed; no correction. */
     294
    278295        /* Astronomical Almanac table is corrected by adding the expression
    279296         *     -0.000091 (ndot + 26)(year-1955)^2  seconds
    … …  
    284301         * are not affected by errors in Lunar or planetary theory.
    285302         */
    286         ans *= 0.01;
    287         if( Y < 1955.0 ) {
    288             B = (Y - 1955.0);
    289             ans += -0.000091 * (-25.8 + 26.0) * B * B;
    290         }
     303        if( Y < 1955.0 )
     304                {
     305                B = (Y - 1955.0);
     306        #if 1
     307                ans += -0.000091 * (-25.8 + 26.0) * B * B;
     308        #else
     309                ans += -0.000091 * (-23.8946 + 26.0) * B * B;
     310        #endif
     311                }
     312
     313#endif /* 0 */
     314
    291315        return( ans );
    292316}
    293317
    294318
    295 #ifdef TEST_DT
    296 main()
     319#ifdef TEST_MAIN
     320
     321/* Exercise program.
     322 */
     323#include <stdio.h>
     324#include <stdlib.h>
     325
     326int main(int ac, char *av[])
    297327{
    298         double ans, y;
    299 
    300         while (scanf("%lf", &y) == 1) {
    301                 ans = deltat((y - 2000.0)*365.25 + J2000);
    302                 printf("%.4lf %.4lf\n", y, ans);
    303         }
     328        double ans, mj, y = atof(av[1]);
     329        year_mjd (y, &mj);
     330        ans = deltat(mj);
     331        printf( "%.4lf\n", ans );
     332        return (0);
    304333}
    305334#endif
    306 
    307 /* For RCS Only -- Do Not Edit */
    308 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/earthsat.c

    r3654 r4017  
    784784
    785785/* For RCS Only -- Do Not Edit */
    786 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: earthsat.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     786static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: earthsat.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/eq_ecl.c

    r3654 r4017  
    7070
    7171/* For RCS Only -- Do Not Edit */
    72 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     72static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/eq_gal.c

    r3654 r4017  
    9797
    9898/* For RCS Only -- Do Not Edit */
    99 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_gal.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     99static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_gal.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/formats.c

    r3654 r4017  
    210210
    211211/* For RCS Only -- Do Not Edit */
    212 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: formats.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     212static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: formats.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/helio.c

    r3654 r4017  
    4949
    5050/* For RCS Only -- Do Not Edit */
    51 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: helio.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     51static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: helio.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/jupmoon.c

    r3654 r4017  
    389389
    390390/* For RCS Only -- Do Not Edit */
    391 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: jupmoon.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     391static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: jupmoon.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/libration.c

    r3654 r4017  
    22162216
    22172217/* For RCS Only -- Do Not Edit */
    2218 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: libration.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     2218static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: libration.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/magdecl.c

    r3654 r4017  
    379379
    380380/* For RCS Only -- Do Not Edit */
    381 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: magdecl.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};
     381static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: magdecl.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/marsmoon.c

    r3654 r4017  
    263263
    264264/* For RCS Only -- Do Not Edit */
    265 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: marsmoon.c,v $ $Date: 2009-07-16 10:34:37 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     265static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: marsmoon.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/misc.c

    r3654 r4017  
    501501
    502502/* For RCS Only -- Do Not Edit */
    503 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     503static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/mjd.c

    r3654 r4017  
    254254
    255255/* For RCS Only -- Do Not Edit */
    256 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mjd.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     256static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mjd.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/moon.c

    r3654 r4017  
    34653465
    34663466/* For RCS Only -- Do Not Edit */
    3467 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moon.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     3467static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moon.c,v $ $Date: 2011-09-21 16:17:49 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/mooncolong.c

    r3654 r4017  
    234234
    235235/* For RCS Only -- Do Not Edit */
    236 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     236static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/moonnf.c

    r3654 r4017  
    6767
    6868/* For RCS Only -- Do Not Edit */
    69 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moonnf.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     69static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moonnf.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/nutation.c

    r3654 r4017  
    439439
    440440/* For RCS Only -- Do Not Edit */
    441 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     441static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/obliq.c

    r3654 r4017  
    2525
    2626/* For RCS Only -- Do Not Edit */
    27 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: obliq.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     27static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: obliq.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/parallactic.c

    r3654 r4017  
    5252
    5353/* For RCS Only -- Do Not Edit */
    54 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallactic.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     54static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallactic.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/parallax.c

    r3654 r4017  
    4040
    4141/* For RCS Only -- Do Not Edit */
    42 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallax.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     42static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallax.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/plans.c

    r3654 r4017  
    225225
    226226/* For RCS Only -- Do Not Edit */
    227 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plans.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     227static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plans.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/plshadow.c

    r3654 r4017  
    4848
    4949/* For RCS Only -- Do Not Edit */
    50 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plshadow.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     50static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plshadow.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/precess.c

    r3654 r4017  
    144144
    145145/* For RCS Only -- Do Not Edit */
    146 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: precess.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     146static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: precess.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/preferences.h

    r3654 r4017  
    2929
    3030/* For RCS Only -- Do Not Edit
    31  * @(#) $RCSfile: preferences.h,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     31 * @(#) $RCSfile: preferences.h,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    3232 */

  • trunk/SophyaExt/XephemAstroLib/reduce.c

    r3654 r4017  
    7676
    7777/* For RCS Only -- Do Not Edit */
    78 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     78static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/refract.c

    r3654 r4017  
    8989
    9090/* For RCS Only -- Do Not Edit */
    91 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: refract.c,v $ $Date: 2009-07-16 10:34:38 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     91static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: refract.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/rings.c

    r3654 r4017  
    4343
    4444/* For RCS Only -- Do Not Edit */
    45 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: rings.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};
     45static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: rings.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/riset.c

    r3654 r4017  
    9898
    9999/* For RCS Only -- Do Not Edit */
    100 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     100static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/riset_cir.c

    r3654 r4017  
    383383
    384384/* For RCS Only -- Do Not Edit */
    385 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset_cir.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     385static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset_cir.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/satlib.h

    r3654 r4017  
    22#define __SATLIB_H
    33
    4 /* $Id: satlib.h,v 1.8 2009-07-16 10:34:39 cmv Exp $ */
     4/* $Id: satlib.h,v 1.9 2011-09-21 16:17:50 cmv Exp $ */
    55
    66typedef struct _SatElem {
    … …  
    203203
    204204/* For RCS Only -- Do Not Edit
    205  * @(#) $RCSfile: satlib.h,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     205 * @(#) $RCSfile: satlib.h,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    206206 */

  • trunk/SophyaExt/XephemAstroLib/satmoon.c

    r3654 r4017  
    508508
    509509/* For RCS Only -- Do Not Edit */
    510 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: satmoon.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     510static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: satmoon.c,v $ $Date: 2011-09-21 16:17:50 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/satspec.h

    r3654 r4017  
    22#define __SATSPEC_H
    33
    4 /* $Id: satspec.h,v 1.8 2009-07-16 10:34:39 cmv Exp $ */
     4/* $Id: satspec.h,v 1.9 2011-09-21 16:17:51 cmv Exp $ */
    55
    66#include "sattypes.h"
    … …  
    4040
    4141/* For RCS Only -- Do Not Edit
    42  * @(#) $RCSfile: satspec.h,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     42 * @(#) $RCSfile: satspec.h,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    4343 */

  • trunk/SophyaExt/XephemAstroLib/sattypes.h

    r3654 r4017  
    22#define __SATTYPES_H
    33
    4 /* $Id: sattypes.h,v 1.8 2009-07-16 10:34:39 cmv Exp $ */
     4/* $Id: sattypes.h,v 1.9 2011-09-21 16:17:51 cmv Exp $ */
    55
    66typedef struct _Vec3 {
    … …  
    2525
    2626/* For RCS Only -- Do Not Edit
    27  * @(#) $RCSfile: sattypes.h,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     27 * @(#) $RCSfile: sattypes.h,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    2828 */

  • trunk/SophyaExt/XephemAstroLib/sdp4.c

    r3654 r4017  
    428428
    429429/* For RCS Only -- Do Not Edit */
    430 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sdp4.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     430static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sdp4.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sgp4.c

    r3654 r4017  
    399399
    400400/* For RCS Only -- Do Not Edit */
    401 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sgp4.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     401static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sgp4.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sphcart.c

    r3654 r4017  
    4141
    4242/* For RCS Only -- Do Not Edit */
    43 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sphcart.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     43static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sphcart.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sun.c

    r3654 r4017  
    4242
    4343/* For RCS Only -- Do Not Edit */
    44 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     44static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/thetag.c

    r3654 r4017  
    33#include "deepconst.h"
    44
    5 /* @(#) $Id: thetag.c,v 1.8 2009-07-16 10:34:39 cmv Exp $ */
     5/* @(#) $Id: thetag.c,v 1.9 2011-09-21 16:17:51 cmv Exp $ */
    66
    77
    … …  
    8888
    8989/* For RCS Only -- Do Not Edit */
    90 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: thetag.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     90static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: thetag.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/twobody.c

    r3654 r4017  
    241241
    242242/* For RCS Only -- Do Not Edit */
    243 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: twobody.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};
     243static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: twobody.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/umoon.c

    r3654 r4017  
    268268
    269269/* For RCS Only -- Do Not Edit */
    270 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: umoon.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};
     270static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: umoon.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.7 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/utc_gst.c

    r3654 r4017  
    9393
    9494/* For RCS Only -- Do Not Edit */
    95 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: utc_gst.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     95static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: utc_gst.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/vector.h

    r3654 r4017  
    22#define __SATVECTOR_H
    33
    4 /* $Id: vector.h,v 1.8 2009-07-16 10:34:39 cmv Exp $ */
     4/* $Id: vector.h,v 1.9 2011-09-21 16:17:51 cmv Exp $ */
    55
    66#define dotp(A,B) ((A).x*(B).x+(A).y*(B).y+(A).z*(B).z)
    … …  
    1616
    1717/* For RCS Only -- Do Not Edit
    18  * @(#) $RCSfile: vector.h,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     18 * @(#) $RCSfile: vector.h,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    1919 */

  • trunk/SophyaExt/XephemAstroLib/vsop87.c

    r3654 r4017  
    207207
    208208/* For RCS Only -- Do Not Edit */
    209 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
     209static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2011-09-21 16:17:51 $ $Revision: 1.10 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/vsop87.h

    r3654 r4017  
    9090
    9191/* For RCS Only -- Do Not Edit
    92  * @(#) $RCSfile: vsop87.h,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $
     92 * @(#) $RCSfile: vsop87.h,v $ $Date: 2011-09-21 16:17:52 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $
    9393 */

  • trunk/SophyaExt/XephemAstroLib/vsop87_data.c

    r3654 r4017  
    69866986
    69876987/* For RCS Only -- Do Not Edit */
    6988 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87_data.c,v $ $Date: 2009-07-16 10:34:39 $ $Revision: 1.8 $ $Name: not supported by cvs2svn $"};
     6988static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87_data.c,v $ $Date: 2011-09-21 16:17:52 $ $Revision: 1.9 $ $Name: not supported by cvs2svn $"};
Note: See TracChangeset for help on using the changeset viewer.