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

trunk/SophyaExt/XephemAstroLib/eq_ecl.c

@@ -2 +2 @@
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 
-static void ecleq_aux P_((int sw, double mjd, double x, double y,
-    double *p, double *q));
+static void ecleq_aux (int sw, double mj, double x, double y,
+    double *p, double *q);
 
 #define EQtoECL 1
@@ -12 +11 @@
 
 
-/* given the modified Julian date, mjd, and an equitorial ra and dec, each in
- * radians, find the corresponding geocentric ecliptic latitude, *lat, and
- * longititude, *lng, also each in radians.
+/* given the modified Julian date, mj, and an equitorial ra and dec, each in
+ * radians, find the corresponding geocentric ecliptic latitude, *lt, and
+ * longititude, *lg, also each in radians.
  * correction for the effect on the angle of the obliquity due to nutation is
  * not included.
  */
 void
-eq_ecl (mjd, ra, dec, lat, lng)
-double mjd;
-double ra, dec;
-double *lat, *lng;
+eq_ecl (double mj, double ra, double dec, double *lt, double *lg)
 {
-        ecleq_aux (EQtoECL, mjd, ra, dec, lng, lat);
+        ecleq_aux (EQtoECL, mj, ra, dec, lg, lt);
 }
 
-/* given the modified Julian date, mjd, and a geocentric ecliptic latitude,
- * *lat, and longititude, *lng, each in radians, find the corresponding
+/* given the modified Julian date, mj, and a geocentric ecliptic latitude,
+ * *lt, and longititude, *lg, each in radians, find the corresponding
  * equitorial ra and dec, also each in radians.
  * correction for the effect on the angle of the obliquity due to nutation is
@@ -34 +30 @@
  */
 void
-ecl_eq (mjd, lat, lng, ra, dec)
-double mjd;
-double lat, lng;
-double *ra, *dec;
+ecl_eq (double mj, double lt, double lg, double *ra, double *dec)
 {
-        ecleq_aux (ECLtoEQ, mjd, lng, lat, ra, dec);
+        ecleq_aux (ECLtoEQ, mj, lg, lt, ra, dec);
 }
 
 static void
-ecleq_aux (sw, mjd, x, y, p, q)
-int sw;                 /* +1 for eq to ecliptic, -1 for vv. */
-double mjd;
-double x, y;            /* sw==1: x==ra, y==dec.  sw==-1: x==lng, y==lat. */
-double *p, *q;          /* sw==1: p==lng, q==lat. sw==-1: p==ra, q==dec. */
+ecleq_aux (
+int sw,                 /* +1 for eq to ecliptic, -1 for vv. */
+double mj,
+double x, double y,     /* sw==1: x==ra, y==dec.  sw==-1: x==lg, y==lt. */
+double *p, double *q)   /* sw==1: p==lg, q==lt. sw==-1: p==ra, q==dec. */
 {
-        static double lastmjd = -10000; /* last mjd calculated */
+        static double lastmj = -10000;  /* last mj calculated */
         static double seps, ceps;       /* sin and cos of mean obliquity */
         double sx, cx, sy, cy, ty, sq;
 
-        if (mjd != lastmjd) {
+        if (mj != lastmj) {
             double eps;
-            obliquity (mjd, &eps);              /* mean obliquity for date */
+            obliquity (mj, &eps);               /* mean obliquity for date */
             seps = sin(eps);
             ceps = cos(eps);
-            lastmjd = mjd;
+            lastmj = mj;
         }
 
@@ -77 +70 @@
 
 /* For RCS Only -- Do Not Edit */
-static 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 $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_ecl.c,v $ $Date: 2004-06-15 16:52:38 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};