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

trunk/SophyaExt/XephemAstroLib/aa_hadec.c

@@ -5 +5 @@
 #include <math.h>
 
-#include "P_.h"
 #include "astro.h"
 
-static void aaha_aux P_((double lat, double x, double y, double *p, double *q));
+static void aaha_aux (double lt, double x, double y, double *p, double *q);
 
-/* given geographical latitude (n+, radians), lat, altitude (up+, radians),
+/* given geographical latitude (n+, radians), lt, altitude (up+, radians),
  * alt, and azimuth (angle round to the east from north+, radians),
  * return hour angle (radians), ha, and declination (radians), dec.
  */
 void
-aa_hadec (lat, alt, az, ha, dec)
-double lat;
-double alt, az;
-double *ha, *dec;
+aa_hadec (
+double lt,
+double alt, double az,
+double *ha, double *dec)
 {
-        aaha_aux (lat, az, alt, ha, dec);
+        aaha_aux (lt, az, alt, ha, dec);
         if (*ha > PI)
             *ha -= 2*PI;
 }
 
-/* given geographical (n+, radians), lat, hour angle (radians), ha, and
+/* given geographical (n+, radians), lt, hour angle (radians), ha, and
  * declination (radians), dec, return altitude (up+, radians), alt, and
  * azimuth (angle round to the east from north+, radians),
  */
 void
-hadec_aa (lat, ha, dec, alt, az)
-double lat;
-double ha, dec;
-double *alt, *az;
+hadec_aa (
+double lt,
+double ha, double dec,
+double *alt, double *az)
 {
-        aaha_aux (lat, ha, dec, az, alt);
+        aaha_aux (lt, ha, dec, az, alt);
 }
 
@@ -43 +42 @@
  */
 double
-geoc_lat (phi)
-double phi;
+geoc_lat (
+double phi)
 {
 #define MAXLAT  degrad(89.9999) /* avoid tan() greater than this */
@@ -56 +55 @@
  */
 static void
-aaha_aux (lat, x, y, p, q)
-double lat;
-double x, y;
-double *p, *q;
+aaha_aux (
+double lt,
+double x, double y,
+double *p, double *q)
 {
-        static double last_lat = -3434, slat, clat;
+        static double last_lt = -3434, slt, clt;
         double cap, B;
 
-        if (lat != last_lat) {
-            slat = sin(lat);
-            clat = cos(lat);
-            last_lat = lat;
+        if (lt != last_lt) {
+            slt = sin(lt);
+            clt = cos(lt);
+            last_lt = lt;
         }
 
-        solve_sphere (-x, PI/2-y, slat, clat, &cap, &B);
+        solve_sphere (-x, PI/2-y, slt, clt, &cap, &B);
         *p = B;
         *q = PI/2 - acos(cap);
@@ -76 +75 @@
 
 /* For RCS Only -- Do Not Edit */
-static 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 $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: aa_hadec.c,v $ $Date: 2004-06-15 16:52:37 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};