Changeset 3111 in Sophya for trunk/SophyaExt/XephemAstroLib/ap_as.c

Timestamp:

Nov 22, 2006, 2:53:31 PM (19 years ago)

Author:

cmv

Message:

mise en conformite xephem 3.7.2 cmv 22/11/2006

File:

• trunk/SophyaExt/XephemAstroLib/ap_as.c (modified) (4 diffs)

trunk/SophyaExt/XephemAstroLib/ap_as.c

@@ -1 +1 @@
 #include <string.h>
+#include <math.h>
 
 #include "astro.h"
@@ -11 +12 @@
 ap_as (Now *np, double Mjd, double *rap, double *decp)
 {
+        double r0 = *rap, d0 = *decp;
         Obj o;
         Now n;
 
+        /* as -> ap */
         zero_mem ((void *)&o, sizeof(o));
         o.o_type = FIXED;
@@ -24 +27 @@
         *rap -= o.s_ra - *rap;
         *decp -= o.s_dec - *decp;
+
+        /* then back to start for second order correction */
+        o.o_type = FIXED;
+        o.f_RA = (float)*rap;
+        o.f_dec = (float)*decp;
+        o.f_epoch = (float)mjd;
+        memcpy ((void *)&n, (void *)np, sizeof(Now));
+        n.n_epoch = EOD;
+        obj_cir (&n, &o);
+        *rap -= o.s_ra - r0;
+        *decp -= o.s_dec - d0;
+
         radecrange (rap, decp);
         precess (mjd, Mjd, rap, decp);
@@ -51 +66 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2005-08-21 10:02:36 $ $Revision: 1.5 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: ap_as.c,v $ $Date: 2006-11-22 13:53:28 $ $Revision: 1.6 $ $Name: not supported by cvs2svn $"};