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

trunk/SophyaExt/XephemAstroLib/jupmoon.c

@@ -4 +4 @@
 #include <stdlib.h>
 #include <string.h>
+#include <errno.h>
 #include <math.h>
 
@@ -31 +32 @@
 static double cmlIImjd; /*    "                      II      " */
 
-/* file containing BDL coefficients */
-static char jbdlfn[] = "jupiter.9910";
-
 /* These values are from the Explanatory Supplement.
  * Precession degrades them gradually over time.
@@ -134 +132 @@
         char buf[1024];
         FILE *fp;
-        int i;
-
-        /* only valid 1999 through 2010 */
-        if (JD < 2451179.50000 || JD >= 2455562.5)
+        char *fn;
+        int i;
+
+        /* check ranges and appropriate data file */
+        if (JD < 2451179.50000)         /* Jan 1 1999 UTC */
+            return (-1);
+        if (JD < 2455562.5)             /* Jan 1 2011 UTC */
+            fn = "jupiter.9910";
+        else if (JD < 2459215.5)        /* Jan 1 2021 UTC */
+            fn = "jupiter.1020";
+        else
             return (-1);
 
         /* open */
-        (void) sprintf (buf, "%s/%s", dir, jbdlfn);
+        (void) sprintf (buf, "%s/%s", dir, fn);
         fp = fopen (buf, "r");
-        if (!fp)
-            return (-1);
+        if (!fp) {
+            fprintf (stderr, "%s: %s\n", fn, strerror(errno));
+            return (-1);
+        }
 
         /* use it */
         if ((i = read_bdl (fp, JD, x, y, z, buf)) < 0) {
-            fprintf (stderr, "%s: %s\n", jbdlfn, buf);
+            fprintf (stderr, "%s: %s\n", fn, buf);
             fclose (fp);
             return (-1);
         }
         if (i != J_NMOONS-1) {
-            fprintf (stderr, "%s: BDL says %d moons, code expects %d", jbdlfn, 
+            fprintf (stderr, "%s: BDL says %d moons, code expects %d", fn, 
                                                                 i, J_NMOONS-1);
             fclose (fp);
@@ -382 +389 @@
 
 /* For RCS Only -- Do Not Edit */
-static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: jupmoon.c,v $ $Date: 2005-08-21 10:02:37 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};
+static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: jupmoon.c,v $ $Date: 2006-11-22 13:53:29 $ $Revision: 1.4 $ $Name: not supported by cvs2svn $"};