Changeset 2818 in Sophya for trunk/SophyaExt/XephemAstroLib/misc.c

Timestamp:

Aug 21, 2005, 12:02:40 PM (20 years ago)

Author:

cmv

Message:

Update de Xephem 3.7 cmv 21/08/2005

File:

• trunk/SophyaExt/XephemAstroLib/misc.c (modified) (3 diffs)

trunk/SophyaExt/XephemAstroLib/misc.c

@@ -180 +180 @@
 }
 
-/* convert ra to ha, in range -PI .. PI.
+/* convert ra to ha, in range 0..2*PI
  * need dec too if not already apparent.
  */
@@ -192 +192 @@
         now_lst (np, &lst);
         ha = hrrad(lst) - ra;
-        if (ha < -PI) ha += 2*PI;
-        if (ha >  PI) ha -= 2*PI;
+        if (ha < 0)
+            ha += 2*PI;
         *hap = ha;
+}
+
+/* find Greenwich Hour Angle of the given object at the given time, 0..2*PI.
+ */
+void
+gha (Now *np, Obj *op, double *ghap)
+{
+        Now n = *np;
+        Obj o = *op;
+        double tmp;
+
+        n.n_epoch = EOD;
+        n.n_lng = 0.0;
+        n.n_lat = 0.0;
+        obj_cir (&n, &o);
+        now_lst (&n, &tmp);
+        tmp = hrrad(tmp) - o.s_ra;
+        if (tmp < 0)
+            tmp += 2*PI;
+        *ghap = tmp;
 }
 
@@ -481 +501 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2005-03-07 16:47:18 $ $Revision: 1.5 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: misc.c,v $ $Date: 2005-08-21 10:02:38 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};