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

Timestamp:

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

Author:

cmv

Message:

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

File:

• trunk/SophyaExt/XephemAstroLib/comet.c (modified) (5 diffs)

trunk/SophyaExt/XephemAstroLib/comet.c

@@ -1 +1 @@
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 
-/* given a modified Julian date, mjd, and a set of heliocentric parabolic
- * orbital elements referred to the epoch of date (mjd):
+/* given a modified Julian date, mj, and a set of heliocentric parabolic
+ * orbital elements referred to the epoch of date (mj):
  *   ep:   epoch of perihelion,
  *   inc:  inclination,
@@ -30 +29 @@
  */
 void
-comet (mjd, ep, inc, ap, qp, om, lpd, psi, rp, rho, lam, bet)
-double mjd;
-double ep, inc, ap, qp, om;
-double *lpd, *psi, *rp, *rho, *lam, *bet;
+comet (double mj, double ep, double inc, double ap, double qp, double om,
+double *lpd, double *psi, double *rp, double *rho, double *lam, double *bet)
 {
         double w, s, s2;
@@ -44 +41 @@
 
 #define ERRLMT  0.0001
-        w = ((mjd-ep)*3.649116e-02)/(qp*sqrt(qp));
+        w = ((mj-ep)*3.649116e-02)/(qp*sqrt(qp));
         s = w/3;
         for (;;) {
@@ -68 +65 @@
         range (lpd, 2*PI);
         rd = *rp * cpsi;
-        sunpos (mjd, &lsn, &rsn, 0);
+        sunpos (mj, &lsn, &rsn, 0);
         lg = lsn+PI;
         re = rsn;
@@ -84 +81 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2001-10-22 12:08:26 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: comet.c,v $ $Date: 2004-06-15 16:52:38 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};