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

trunk/SophyaExt/XephemAstroLib/sun.c

@@ -2 +2 @@
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 #include "vsop87.h"
 
-/* given the modified JD, mjd, return the true geocentric ecliptic longitude
+/* given the modified JD, mj, return the true geocentric ecliptic longitude
  *   of the sun for the mean equinox of the date, *lsn, in radians, the
  *   sun-earth distance, *rsn, in AU, and the latitude *bsn, in radians
@@ -17 +16 @@
  */
 void
-sunpos (mjd, lsn, rsn, bsn)
-double mjd;
-double *lsn, *rsn, *bsn;
+sunpos (double mj, double *lsn, double *rsn, double *bsn)
 {
-        static double last_mjd = -3691, last_lsn, last_rsn, last_bsn;
+        static double last_mj = -3691, last_lsn, last_rsn, last_bsn;
         double ret[6];
 
-        if (mjd == last_mjd) {
+        if (mj == last_mj) {
             *lsn = last_lsn;
             *rsn = last_rsn;
@@ -31 +28 @@
         }
 
-        vsop87(mjd, SUN, 0.0, ret);     /* full precision earth pos */
+        vsop87(mj, SUN, 0.0, ret);      /* full precision earth pos */
 
         *lsn = ret[0] - PI;             /* revert to sun pos */
@@ -39 +36 @@
         last_rsn = *rsn = ret[2];
         last_bsn = -ret[1];
-        last_mjd = mjd;
+        last_mj = mj;
 
         if (bsn) *bsn = last_bsn;       /* assign only if non-NULL pointer */
@@ -45 +42 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2001-10-22 12:08:28 $ $Revision: 1.2 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sun.c,v $ $Date: 2004-06-15 16:52:40 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};