Changeset 1719 in Sophya for trunk/SophyaExt/XephemAstroLib


Ignore:
Timestamp:
Oct 22, 2001, 2:08:28 PM (24 years ago)
Author:
cmv
Message:

Adapted to version 3.5 xephem cmv 22/10/2001

Location:
trunk/SophyaExt/XephemAstroLib
Files:
3 added
54 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaExt/XephemAstroLib/P_.h

    r1457 r1719  
    1111
    1212/* For RCS Only -- Do Not Edit
    13  * @(#) $RCSfile: P_.h,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     13 * @(#) $RCSfile: P_.h,v $ $Date: 2001-10-22 12:08:25 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    1414 */

  • trunk/SophyaExt/XephemAstroLib/README

    r1457 r1719  
    1919
    2020! For RCS Only -- Do Not Edit
    21 ! @(#) $RCSfile: README,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     21! @(#) $RCSfile: README,v $ $Date: 2001-10-22 12:08:25 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $

  • trunk/SophyaExt/XephemAstroLib/aa_hadec.c

    r1457 r1719  
    7676
    7777/* For RCS Only -- Do Not Edit */
    78 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aa_hadec.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     78static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aa_hadec.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/aberration.c

    r1457 r1719  
    134134
    135135/* For RCS Only -- Do Not Edit */
    136 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aberration.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     136static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aberration.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/actan.c

    r1457 r1719  
    11#include <math.h>
    22
    3 /* @(#) $Id: actan.c,v 1.1.1.1 2001-04-10 14:40:45 cmv Exp $ */
     3/* @(#) $Id: actan.c,v 1.2 2001-10-22 12:08:26 cmv Exp $ */
     4
     5/* commonly in math.h, but not in strict ANSI C */
     6#ifndef M_PI
     7#define M_PI            3.14159265358979323846
     8#define M_PI_2          1.57079632679489661923
     9#endif
    410
    511double
    … …  
    5965
    6066/* For RCS Only -- Do Not Edit */
    61 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: actan.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     67static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: actan.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/airmass.c

    r1457 r1719  
    2525
    2626/* For RCS Only -- Do Not Edit */
    27 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: airmass.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     27static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: airmass.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/anomaly.c

    r1457 r1719  
    6464
    6565/* For RCS Only -- Do Not Edit */
    66 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: anomaly.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     66static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: anomaly.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/ap_as.c

    r1457 r1719  
    2828        obj_cir (&n, &o);
    2929        *rap -= o.s_ra - *rap;
     30        *decp -= o.s_dec - *decp;
     31        if (*decp >  PI/2) {*decp =  PI - *decp; *rap += PI; }
     32        if (*decp < -PI/2) {*decp = -PI - *decp; *rap += PI; }
    3033        range (rap, 2*PI);
    31         *decp -= o.s_dec - *decp;
    32         if (*decp >  PI/2) *decp =  PI - *decp;
    33         if (*decp < -PI/2) *decp = -PI - *decp;
    3434        precess (mjd, Mjd, rap, decp);
    3535}
    … …  
    6060
    6161/* For RCS Only -- Do Not Edit */
    62 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     62static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/astro.h

    r1457 r1719  
    163163extern void refract P_((double pr, double tr, double ta, double *aa));
    164164
     165/* rings.c */
     166extern void satrings P_((double sb, double sl, double sr, double el, double er,
     167    double JD, double *etiltp, double *stiltp));
     168
    165169/* riset.c */
    166170extern void riset P_((double ra, double dec, double lat, double dis,
    … …  
    176180extern void sunpos P_((double mjd, double *lsn, double *rsn, double *bsn));
    177181
     182/* twobody.c */
     183extern void vrc P_((double *v, double *r, double tp, double e, double q));
     184
    178185/* utc_gst.c */
    179186extern void utc_gst P_((double mjd, double utc, double *gst));
    … …  
    184191
    185192/* For RCS Only -- Do Not Edit
    186  * @(#) $RCSfile: astro.h,v $ $Date: 2001-04-10 14:40:45 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     193 * @(#) $RCSfile: astro.h,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    187194 */

  • trunk/SophyaExt/XephemAstroLib/auxil.c

    r1457 r1719  
    1515static int prefs[NPREFS] = {
    1616    PREF_TOPO, PREF_METRIC, PREF_MDY, PREF_UTCTZ, PREF_HIPREC, PREF_NOMSGBELL,
    17     PREF_PREFILL, PREF_TIPSON, PREF_CONFIRMON
     17    PREF_PREFILL, PREF_TIPSON, PREF_CONFIRMON, PREF_WEEKSTART
    1818};
    1919
    … …  
    4848
    4949/* For RCS Only -- Do Not Edit */
    50 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: auxil.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     50static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: auxil.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/chap95.c

    r1457 r1719  
    178178
    179179/* For RCS Only -- Do Not Edit */
    180 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     180static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/chap95.h

    r1457 r1719  
    6666
    6767/* For RCS Only -- Do Not Edit
    68  * @(#) $RCSfile: chap95.h,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     68 * @(#) $RCSfile: chap95.h,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    6969 */

  • trunk/SophyaExt/XephemAstroLib/chap95_data.c

    r1457 r1719  
    781781
    782782/* For RCS Only -- Do Not Edit */
    783 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95_data.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     783static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: chap95_data.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/circum.c

    r1457 r1719  
    6767        double lam, bet;        /* geocentric ecliptic long and lat */
    6868        double dia, mag;        /* angular diameter at 1 AU and magnitude */
    69         double f;               /* fractional phase from earth */
    7069        int p;
    7170
    … …  
    8483        sunpos (mjed, &lsn, &rsn, 0);
    8584
    86         /* find helio long/lat; sun/planet and earth/plant dist; ecliptic
     85        /* find helio long/lat; sun/planet and earth/planet dist; ecliptic
    8786         * long/lat; diameter and mag.
    8887         */
    … …  
    9291        cir_sky (np, lpd, psi, rp, &rho, lam, bet, lsn, rsn, op);
    9392
    94         /* compute magnitude and angular size */
    95         f = op->s_phase ? 5*log10(rp*rho) - 5*log10(op->s_phase/100) : 100;
    96         set_smag (op, mag+f);
     93        /* set magnitude and angular size */
     94        set_smag (op, mag);
    9795        op->s_size = (float)(dia/rho);
    9896
    … …  
    129127        ra = op->f_RA;
    130128        dec = op->f_dec;
    131         precess (op->f_epoch, mjd, &ra, &dec);
     129        precess (op->f_epoch, mjed, &ra, &dec);
    132130
    133131        /* convert equatoreal ra/dec to mean geocentric ecliptic lat/long */
    134         eq_ecl (mjd, ra, dec, &bet, &lam);
     132        eq_ecl (mjed, ra, dec, &bet, &lam);
    135133
    136134        /* find solar ecliptical long.(mean equinox) and distance from earth */
    … …  
    138136
    139137        /* allow for relativistic light bending near the sun */
    140         deflect (mjd, lam, bet, lsn, rsn, 1e10, &ra, &dec);
     138        deflect (mjed, lam, bet, lsn, rsn, 1e10, &ra, &dec);
    141139
    142140        /* TODO: correction for annual parallax would go here */
    … …  
    145143         * geocentric
    146144         */
    147         nut_eq(mjd, &ra, &dec);
    148         ab_eq(mjd, lsn, &ra, &dec);
     145        nut_eq(mjed, &ra, &dec);
     146        ab_eq(mjed, lsn, &ra, &dec);
    149147        op->s_gaera = (float)ra;
    150148        op->s_gaedec = (float)dec;
    … …  
    193191        double dt;              /* light travel time to object */
    194192        double lg;              /* helio long of earth */
    195         double nu, ea;          /* true anomaly and eccentric anomaly */
    196         double ma;              /* mean anomaly */
     193        double nu;              /* true anomaly */
    197194        double rp=0;            /* distance from the sun */
    198195        double lo, slo, clo;    /* angle from ascending node */
    … …  
    206203        double lam;             /* geocentric ecliptic longitude */
    207204        double bet;             /* geocentric ecliptic latitude */
    208         double e;               /* fast eccentricity */
    209205        double ll=0, sll, cll;  /* helio angle between object and earth */
    210206        double mag;             /* magnitude */
    211207        double e_n;             /* mean daily motion */
     208        double tp;              /* time from perihelion (days) */
    212209        double rpd=0;
    213210        double y;
    … …  
    217214        sunpos (mjed, &lsn, &rsn, 0);
    218215        lg = lsn + PI;
    219 
    220         /* faster access to eccentricty */
    221         e = op->e_e;
    222216
    223217        /* mean daily motion is derived fro mean distance */
    … …  
    235229                                        &inc, &om, &Om);
    236230
    237             ma = degrad (op->e_M + (mjed - op->e_cepoch - dt) * e_n);
    238             anomaly (ma, e, &nu, &ea);
    239             rp = op->e_a * (1-e*e) / (1+e*cos(nu));
     231            tp = mjed - dt - (op->e_cepoch - op->e_M/e_n);
     232            vrc (&nu, &rp, tp, op->e_e, op->e_a*(1-op->e_e));
     233            nu = degrad(nu);
    240234            lo = nu + om;
    241235            slo = sin(lo);
    … …  
    301295        double dt;              /* light travel time to object */
    302296        double lg;              /* helio long of earth */
    303         double nu, ea;          /* true anomaly and eccentric anomaly */
    304         double ma;              /* mean anomaly */
     297        double nu;              /* true anomaly and eccentric anomaly */
    305298        double rp=0;            /* distance from the sun */
    306299        double lo, slo, clo;    /* angle from ascending node */
    … …  
    319312        double mag;             /* magnitude */
    320313        double a;               /* mean distance */
     314        double tp;              /* time from perihelion (days) */
    321315        double rpd=0;
    322316        double y;
    … …  
    342336                            &inc, &om, &Om);
    343337
    344             ma = degrad ((mjed - op->h_ep - dt) * n);
    345             anomaly (ma, e, &nu, &ea);
    346             rp = a * (e*e-1.0) / (1.0+e*cos(nu));
     338            tp = mjed - dt - op->h_ep;
     339            vrc (&nu, &rp, tp, op->h_e, op->h_qp);
     340            nu = degrad(nu);
    347341            lo = nu + om;
    348342            slo = sin(lo);
    … …  
    581575
    582576        /* convert to equatoreal [mean equator, with mean obliquity] */
    583         ecl_eq (mjd, bet, lam, &ra, &dec);
     577        ecl_eq (mjed, bet, lam, &ra, &dec);
    584578        tra = ra;       /* keep mean coordinates */
    585579        tdec = dec;
    … …  
    592586         */
    593587        if (!is_planet(op,SUN) && !is_planet(op,MOON))
    594             deflect (mjd, op->s_hlong, op->s_hlat, lsn, rsn, *rho, &ra, &dec);
     588            deflect (mjed, op->s_hlong, op->s_hlat, lsn, rsn, *rho, &ra, &dec);
    595589
    596590        /* correct ra/dec to form geocentric apparent */
    597         nut_eq (mjd, &ra, &dec);
     591        nut_eq (mjed, &ra, &dec);
    598592        if (!is_planet(op,MOON))
    599             ab_eq (mjd, lsn, &ra, &dec);
     593            ab_eq (mjed, lsn, &ra, &dec);
    600594        op->s_gaera = (float)ra;
    601595        op->s_gaedec = (float)dec;
    … …  
    635629            ra = tra + dra;
    636630            dec = tdec + ddec;
    637             precess (mjd, epoch, &ra, &dec);
     631            precess (mjed, epoch, &ra, &dec);
    638632        }
    639633        range(&ra, 2*PI);
    … …  
    755749
    756750/* For RCS Only -- Do Not Edit */
    757 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: circum.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     751static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: circum.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/circum.h

    r1457 r1719  
    6868
    6969/* longest object name, including trailing '\0' */
    70 #define MAXNM   14
     70#define MAXNM   22
    7171
    7272/* Obj is a massive union.
    … …  
    156156    float  eo_om;       /* argument of perihelion, degress */
    157157    float  eo_a;        /* mean distance, aka,semi-maj axis,AU */
    158     float  eo_e;        /* eccentricity */
    159158    float  eo_M;        /* mean anomaly, ie, degrees from perihelion at cepoch*/
    160159    float  eo_size;     /* angular size, in arc seconds at 1 AU */
     160    double eo_e;        /* eccentricity (double for when near 1 computing q) */
    161161    double eo_cepoch;   /* epoch date (M reference), as an mjd */
    162162    double eo_epoch;    /* equinox year (inc/Om/om reference), as an mjd. */
    … …  
    440440
    441441/* For RCS Only -- Do Not Edit
    442  * @(#) $RCSfile: circum.h,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
    443  */
     442 * @(#) $RCSfile: circum.h,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
     443 */

  • trunk/SophyaExt/XephemAstroLib/comet.c

    r1457 r1719  
    8484
    8585/* For RCS Only -- Do Not Edit */
    86 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     86static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/compare_with_xephem.csh

    r1680 r1719  
    44
    55set dir = $DPCSOURCE/XephemAstroLib/
     6set log = `pwd`/compare.log
    67
    7 rm -f compare.log; touch compare.log
     8rm -f $log; touch $log
    89
     10echo ======================= NOUVEAU REPERTOIRE >>& $log
     11ls *.h >>& $log
     12ls *.c >>& $log
     13
     14echo ======================= REPERTOIRE $dir >>& $log
     15pushd $dir
     16ls *.h  >>& $log
     17ls *.c >>& $log
     18popd
     19
     20echo ======================= INCLUDE >>& $log
    921foreach f ( *.h )
    10   echo ======================= $f >> compare.log
    11   diff $f $dir/$f >> compare.log
     22  echo ======================= $f >>& $log
     23  diff $f $dir/$f >>& $log
    1224end
    1325
     26echo ======================= FICHIERS >>& $log
    1427foreach f ( *.c )
    15   echo ======================= $f >> compare.log
    16   diff $f $dir/$f >> compare.log
     28  echo ======================= $f >>& $log
     29  diff $f $dir/$f >>& $log
    1730end
    1831

  • trunk/SophyaExt/XephemAstroLib/dbfmt.c

    r1457 r1719  
    510510            break;
    511511        case E_E:
    512             op->e_e = (float) atod (bp);
     512            op->e_e = atod (bp);
    513513            break;
    514514        case E_M:
    … …  
    10141014
    10151015/* For RCS Only -- Do Not Edit */
    1016 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     1016static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: dbfmt.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/deep.c

    r1457 r1719  
    788788
    789789/* For RCS Only -- Do Not Edit */
    790 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deep.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     790static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deep.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/deepconst.h

    r1457 r1719  
    22#define _CONST_H
    33
    4 /* $Id: deepconst.h,v 1.1.1.1 2001-04-10 14:40:46 cmv Exp $ */
     4/* $Id: deepconst.h,v 1.2 2001-10-22 12:08:26 cmv Exp $ */
    55
    66
    … …  
    3131
    3232/* For RCS Only -- Do Not Edit
    33  * @(#) $RCSfile: deepconst.h,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     33 * @(#) $RCSfile: deepconst.h,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    3434 */

  • trunk/SophyaExt/XephemAstroLib/deltat.c

    r1457 r1719  
    305305
    306306/* For RCS Only -- Do Not Edit */
    307 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     307static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: deltat.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/earthsat.c

    r1457 r1719  
    797797
    798798/* For RCS Only -- Do Not Edit */
    799 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: earthsat.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     799static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: earthsat.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/eq_ecl.c

    r1457 r1719  
    5151        static double lastmjd = -10000; /* last mjd calculated */
    5252        static double seps, ceps;       /* sin and cos of mean obliquity */
    53         double sx, cx, sy, cy, ty;
     53        double sx, cx, sy, cy, ty, sq;
    5454
    5555        if (mjd != lastmjd) {
    … …  
    6767        cx = cos(x);
    6868        sx = sin(x);
    69         /* modifs CMV *q = asin((sy*ceps)-(cy*seps*sx*sw)); */
    70         *q = (sy*ceps)-(cy*seps*sx*sw);
    71         if(*q<-1.) *q = -PI/2.; else if(*q>1.) *q = PI/2.; else *q = asin(*q);
     69        sq = (sy*ceps)-(cy*seps*sx*sw);
     70        if (sq < -1) sq = -1;
     71        if (sq >  1) sq =  1;
     72        *q = asin(sq);
    7273        *p = atan(((sx*ceps)+(ty*seps*sw))/cx);
    7374        if (cx<0) *p += PI;             /* account for atan quad ambiguity */
    … …  
    7677
    7778/* For RCS Only -- Do Not Edit */
    78 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     79static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/eq_gal.c

    r1457 r1719  
    102102
    103103/* For RCS Only -- Do Not Edit */
    104 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_gal.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     104static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_gal.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/formats.c

    r1457 r1719  
    305305
    306306/* For RCS Only -- Do Not Edit */
    307 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: formats.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     307static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: formats.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/helio.c

    r1457 r1719  
    5252
    5353/* For RCS Only -- Do Not Edit */
    54 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: helio.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     54static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: helio.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/libration.c

    r1457 r1719  
    22282228
    22292229/* For RCS Only -- Do Not Edit */
    2230 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: libration.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     2230static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: libration.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/misc.c

    r1457 r1719  
    473473
    474474/* For RCS Only -- Do Not Edit */
    475 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     475static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/mjd.c

    r1457 r1719  
    261261
    262262/* For RCS Only -- Do Not Edit */
    263 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mjd.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     263static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mjd.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/moon.c

    r1457 r1719  
    34853485
    34863486/* For RCS Only -- Do Not Edit */
    3487 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moon.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     3487static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: moon.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/mooncolong.c

    r1457 r1719  
    262262
    263263/* For RCS Only -- Do Not Edit */
    264 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     264static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/nutation.c

    r1457 r1719  
    440440
    441441/* For RCS Only -- Do Not Edit */
    442 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     442static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: nutation.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/obliq.c

    r1457 r1719  
    2828
    2929/* 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 $"};
     30static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: obliq.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/parallax.c

    r1457 r1719  
    4242
    4343/* For RCS Only -- Do Not Edit */
    44 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallax.c,v $ $Date: 2001-04-10 14:40:46 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     44static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: parallax.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/plans.c

    r1457 r1719  
    108108 * [planet][0] = angular size at 1 AU
    109109 * [planet][1] = magnitude at 1 AU from sun and earth and 0 deg phase angle
    110  */
    111 static double vis_elements[8][2] = {
    112         /* Mercury */   { 6.74, -0.42, },
    113         /* Venus */     { 16.92, -4.34, },
    114         /* Mars */      { 9.36, -1.20, },
    115         /* Jupiter */   { 196.74, -9.4, },
    116         /* Saturn */    { 165.6, -8.88, },
    117         /* Uranus */    { 65.8, -7.19, },
    118         /* Neptune */   { 62.2, -6.87, },
    119         /* Pluto */     { 8.2, -1.0, }
     110 * [planet][2] = A
     111 * [planet][3] = B
     112 * [planet][4] = C
     113 *   where mag correction = A*(i/100) + B*(i/100)^2 + C*(i/100)^3
     114 *      i = angle between sun and earth from planet, degrees
     115 * from Explanatory Supplement, 1992
     116 */
     117static double vis_elements[8][5] = {
     118        /* Mercury */   { 6.74, -0.36, 3.8, -2.73, 2.00},
     119        /* Venus */     { 16.92, -4.29, 0.09, 2.39, -.65},
     120        /* Mars */      { 9.36, -1.52, 1.60, 0., 0.},
     121        /* Jupiter */   { 196.74, -9.25, 0.50, 0., 0.},
     122        /* Saturn */    { 165.6, -8.88, 4.40, 0., 0.},
     123        /* Uranus */    { 65.8, -7.19, 0.28, 0., 0.},
     124        /* Neptune */   { 62.2, -6.87, 0., 0., 0.},
     125        /* Pluto */     { 8.2, -1.01, 4.1, 0., 0.}
    120126};
    121127
    … …  
    132138 *         seen from the center of the Earth for the given instant.
    133139 *   dia:  angular diameter in arcsec at 1 AU,
    134  *   mag:  visual magnitude when 1 AU from sun and earth at 0 phase angle.
     140 *   mag:  visual magnitude
    135141 *
    136142 * all angles are in radians, all distances in AU.
    … …  
    148154{
    149155        static double lastmjd = -10000;
    150         static double lsn, bsn, rsn;    /* geometric geocentric coords of sun */
    151         static double xsn, ysn, zsn;
     156        static double lsn, bsn, rsn;    /* geocentric coords of sun */
     157        static double xsn, ysn, zsn;    /* cartesian " */
    152158        double lp, bp, rp;              /* heliocentric coords of planet */
    153159        double xp, yp, zp, rho;         /* rect. coords and geocentric dist. */
    154160        double dt;                      /* light time */
     161        double *vp;                     /* vis_elements[p] */
     162        double ci, i;                   /* sun/earth angle: cos, degrees */
    155163        int pass;
    156164
    … …  
    202210        }
    203211
    204         *dia = vis_elements[p][0];
    205         *mag = vis_elements[p][1];
     212        vp = vis_elements[p];
     213        *dia = vp[0];
     214
     215        /* solve plane triangle, assume sun/earth dist == 1 */
     216        ci = (rp*rp + rho*rho - 1)/(2*rp*rho);
     217
     218        /* expl supp equation for mag */
     219        if (ci < -1) ci = -1;
     220        if (ci >  1) ci =  1;
     221        i = raddeg(acos(ci))/100.;
     222        *mag = vp[1] + 5*log10(rho*rp) + i*(vp[2] + i*(vp[3] + i*vp[4]));
     223
     224        /* rings contribution if SATURN */
     225        if (p == SATURN) {
     226            double et, st, set;
     227            satrings (bp, lp, rp, lsn+PI, rsn, mjd+MJD0, &et, &st);
     228            set = sin(fabs(et));
     229            *mag += (-2.60 + 1.25*set)*set;
     230        }
    206231}
    207232
    208233/* For RCS Only -- Do Not Edit */
    209 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plans.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     234static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: plans.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/precess.c

    r1457 r1719  
    146146
    147147/* For RCS Only -- Do Not Edit */
    148 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: precess.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     148static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: precess.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/preferences.h

    r1457 r1719  
    77
    88typedef enum {
    9     PREF_EQUATORIAL, PREF_UNITS, PREF_DATE_FORMAT, PREF_ZONE,
    10     PREF_DPYPREC, PREF_MSG_BELL, PREF_PRE_FILL, PREF_TIPS, PREF_CONFIRM,
     9    PREF_EQUATORIAL, PREF_UNITS, PREF_DATE_FORMAT, PREF_ZONE, PREF_DPYPREC,
     10    PREF_MSG_BELL, PREF_PRE_FILL, PREF_TIPS, PREF_CONFIRM, PREF_WEEKSTART,
    1111    NPREFS
    1212} Preferences;
    … …  
    2121typedef enum {PREF_TIPSON, PREF_NOTIPS} PrefTips;
    2222typedef enum {PREF_CONFIRMON, PREF_NOCONFIRM} PrefConfirm;
     23typedef enum {PREF_SAT, PREF_SUN, PREF_MON} PrefWeekStart;
    2324
    2425extern int pref_get P_((Preferences p));
    … …  
    2829
    2930/* For RCS Only -- Do Not Edit
    30  * @(#) $RCSfile: preferences.h,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     31 * @(#) $RCSfile: preferences.h,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    3132 */

  • trunk/SophyaExt/XephemAstroLib/reduce.c

    r1457 r1719  
    7575
    7676/* For RCS Only -- Do Not Edit */
    77 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     77static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: reduce.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/refract.c

    r1457 r1719  
    3232            *ta = taLT + (taGE - taLT)*p;
    3333        }
    34            
    3534}
    3635
    … …  
    6059        r = degrad(a/b);
    6160
    62         *ta  =  aa - r;
     61        *ta  =  (aa < 0 && r < 0) ? aa : aa - r;        /* 0 below ~5 degs */
    6362}
    6463
    … …  
    101100
    102101/* For RCS Only -- Do Not Edit */
    103 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: refract.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     102static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: refract.c,v $ $Date: 2001-10-22 12:08:27 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/riset.c

    r1457 r1719  
    103103
    104104/* For RCS Only -- Do Not Edit */
    105 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     105static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/riset_cir.c

    r1457 r1719  
    385385
    386386/* For RCS Only -- Do Not Edit */
    387 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset_cir.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     387static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: riset_cir.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/satlib.h

    r1457 r1719  
    22#define __SATLIB_H
    33
    4 /* $Id: satlib.h,v 1.1.1.1 2001-04-10 14:40:47 cmv Exp $ */
     4/* $Id: satlib.h,v 1.2 2001-10-22 12:08:28 cmv Exp $ */
    55
    66typedef struct _SatElem {
    … …  
    203203
    204204/* For RCS Only -- Do Not Edit
    205  * @(#) $RCSfile: satlib.h,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     205 * @(#) $RCSfile: satlib.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    206206 */

  • trunk/SophyaExt/XephemAstroLib/satspec.h

    r1457 r1719  
    22#define __SATSPEC_H
    33
    4 /* $Id: satspec.h,v 1.1.1.1 2001-04-10 14:40:47 cmv Exp $ */
     4/* $Id: satspec.h,v 1.2 2001-10-22 12:08:28 cmv Exp $ */
    55
    66#include "sattypes.h"
    … …  
    4040
    4141/* For RCS Only -- Do Not Edit
    42  * @(#) $RCSfile: satspec.h,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     42 * @(#) $RCSfile: satspec.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    4343 */

  • trunk/SophyaExt/XephemAstroLib/sattypes.h

    r1457 r1719  
    22#define __SATTYPES_H
    33
    4 /* $Id: sattypes.h,v 1.1.1.1 2001-04-10 14:40:47 cmv Exp $ */
     4/* $Id: sattypes.h,v 1.2 2001-10-22 12:08:28 cmv Exp $ */
    55
    66typedef struct _Vec3 {
    … …  
    2525
    2626/* For RCS Only -- Do Not Edit
    27  * @(#) $RCSfile: sattypes.h,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     27 * @(#) $RCSfile: sattypes.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    2828 */

  • trunk/SophyaExt/XephemAstroLib/sdp4.c

    r1457 r1719  
    110110    /* private temporary variables */
    111111    double A,AXN,AYN,AYNL,BETA,BETAL,CAPU,COS2U,COSEPW,
    112         COSIK,COSNOK,COSU,COSUK,E,ECOSE,ELSQ,EM,EPW,ESINE,OMGADF,PL,
     112        COSIK,COSNOK,COSU,COSUK,E,ECOSE,ELSQ,EM=0,EPW,ESINE,OMGADF,PL,
    113113        R,RDOT,RDOTK,RFDOT,RFDOTK,RK,SIN2U,SINEPW,SINIK,SINNOK,
    114114        SINU,SINUK,TEMP,TEMP1,TEMP2,TEMP3,TEMP4,TEMP5,TEMP6,TEMPA,
    115         TEMPE,TEMPL,TSQ,U,UK,UX,UY,UZ,VX,VY,VZ,XINC,XINCK,XL,XLL,XLT,
     115        TEMPE,TEMPL,TSQ,U,UK,UX,UY,UZ,VX,VY,VZ,XINC=0,XINCK,XL,XLL,XLT,
    116116        XMAM,XMDF,XMX,XMY,XN,XNODDF,XNODE,XNODEK;
    117117
    … …  
    432432
    433433/* For RCS Only -- Do Not Edit */
    434 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sdp4.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     434static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sdp4.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sgp4.c

    r1457 r1719  
    403403
    404404/* For RCS Only -- Do Not Edit */
    405 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sgp4.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     405static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sgp4.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sphcart.c

    r1457 r1719  
    4040
    4141/* For RCS Only -- Do Not Edit */
    42 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sphcart.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     42static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sphcart.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/sun.c

    r1457 r1719  
    4545
    4646/* For RCS Only -- Do Not Edit */
    47 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     47static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/thetag.c

    r1457 r1719  
    33#include "deepconst.h"
    44
    5 /* @(#) $Id: thetag.c,v 1.1.1.1 2001-04-10 14:40:47 cmv Exp $ */
     5/* @(#) $Id: thetag.c,v 1.2 2001-10-22 12:08:28 cmv Exp $ */
    66
    77
    … …  
    8888
    8989/* For RCS Only -- Do Not Edit */
    90 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: thetag.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     90static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: thetag.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/utc_gst.c

    r1457 r1719  
    100100
    101101/* For RCS Only -- Do Not Edit */
    102 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: utc_gst.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     102static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: utc_gst.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/vector.h

    r1457 r1719  
    22#define __SATVECTOR_H
    33
    4 /* $Id: vector.h,v 1.1.1.1 2001-04-10 14:40:47 cmv Exp $ */
     4/* $Id: vector.h,v 1.2 2001-10-22 12:08:28 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: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     18 * @(#) $RCSfile: vector.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    1919 */

  • trunk/SophyaExt/XephemAstroLib/vsop87.c

    r1681 r1719  
    152152
    153153            /* precision threshold */
    154             if(alpha>0) p = t_abs[alpha-1]; else p=0.;
    155             p = q/(t_abs[alpha] + alpha * p * 1e-4 + 1e-35);
     154            p= alpha ? q/(t_abs[alpha] + alpha*t_abs[alpha-1]*1e-4 + 1e-35) : q;
    156155#if VSOP_SPHERICAL
    157156            if (cooidx == 2)    /* scale by semimajor axis for radius */
    … …  
    212211
    213212/* For RCS Only -- Do Not Edit */
    214 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2001-10-11 12:40:41 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
     213static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};

  • trunk/SophyaExt/XephemAstroLib/vsop87.h

    r1457 r1719  
    9090
    9191/* For RCS Only -- Do Not Edit
    92  * @(#) $RCSfile: vsop87.h,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $
     92 * @(#) $RCSfile: vsop87.h,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $
    9393 */

  • trunk/SophyaExt/XephemAstroLib/vsop87_data.c

    r1457 r1719  
    69866986
    69876987/* For RCS Only -- Do Not Edit */
    6988 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87_data.c,v $ $Date: 2001-04-10 14:40:47 $ $Revision: 1.1.1.1 $ $Name: not supported by cvs2svn $"};
     6988static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: vsop87_data.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
Note: See TracChangeset for help on using the changeset viewer.