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

trunk/SophyaExt/XephemAstroLib/eq_gal.c

@@ -4 +4 @@
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 
-static void galeq_aux P_((int sw, double x, double y, double *p, double *q));
-static void galeq_init P_((void));
+static void galeq_aux (int sw, double x, double y, double *p, double *q);
+static void galeq_init (void);
 
 #define EQtoGAL 1
@@ -18 +17 @@
 static double gpd = degrad(27.12825);   /* Dec of  " */
 static double cgpd, sgpd;               /* cos() and sin() of gpd */
-static double mjd2000;                  /* mjd of 2000 */
+static double mj2000;                   /* mj of 2000 */
 static int before;                      /* whether these have been set yet */
 
 /* given ra and dec, each in radians, for the given epoch, find the
- * corresponding galactic latitude, *lat, and longititude, *lng, also each in
+ * corresponding galactic latitude, *lt, and longititude, *lg, also each in
  * radians.
  */
 void
-eq_gal (mjd, ra, dec, lat, lng)
-double mjd, ra, dec;
-double *lat, *lng;
+eq_gal (double mj, double ra, double dec, double *lt, double *lg)
 {
         galeq_init();
-        precess (mjd, mjd2000, &ra, &dec);
-        galeq_aux (EQtoGAL, ra, dec, lng, lat);
+        precess (mj, mj2000, &ra, &dec);
+        galeq_aux (EQtoGAL, ra, dec, lg, lt);
 }
 
-/* given galactic latitude, lat, and longititude, lng, each in radians, find
+/* given galactic latitude, lt, and longititude, lg, each in radians, find
  * the corresponding equitorial ra and dec, also each in radians, at the 
  * given epoch.
  */
 void
-gal_eq (mjd, lat, lng, ra, dec)
-double mjd, lat, lng;
-double *ra, *dec;
+gal_eq (double mj, double lt, double lg, double *ra, double *dec)
 {
         galeq_init();
-        galeq_aux (GALtoEQ, lng, lat, ra, dec);
-        precess (mjd2000, mjd, ra, dec);
+        galeq_aux (GALtoEQ, lg, lt, ra, dec);
+        precess (mj2000, mj, ra, dec);
 }
 
 static void
-galeq_aux (sw, x, y, p, q)
-int sw;                 /* +1 for eq to gal, -1 for vv. */
-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. */
+galeq_aux (
+int sw,                 /* +1 for eq to gal, -1 for vv. */
+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. */
 {
         double sy, cy, a, ca, sa, b, sq, c, d; 
@@ -96 +91 @@
             cgpd = cos (gpd);
             sgpd = sin (gpd);
-            mjd2000 = J2000;
+            mj2000 = J2000;
             before = 1;
         }
@@ -102 +97 @@
 
 /* For RCS Only -- Do Not Edit */
-static 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 $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: eq_gal.c,v $ $Date: 2004-06-15 16:52:38 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};