Initial commit

1 files changed, 196 insertions, 0 deletions

diff --git a/vendor/voclient/libvoclient/vocSkybot_f77.c b/vendor/voclient/libvoclient/vocSkybot_f77.c
new file mode 100644
index 00000000..e95bc29c
--- /dev/null
+++ b/vendor/voclient/libvoclient/vocSkybot_f77.c
@@ -0,0 +1,196 @@
+/***************************************************************************
+ *
+ *  SkyBoT Minor Planet Service Interface:
+ *  --------------------------------------
+ *
+ *  A utility class for accessing the IMCCE Skybot ephemerides service.  
+ *  The service is called as part of the constructor for the object 
+ *  allowing easy access to the result data.  Class methods also exist so
+ *  a VO Client interface the same access to the data.
+ *
+ *  Reference:      http://www.imcce.fr/webservices/skybot/
+ *
+ *  Class Methods:      
+ *  --------------
+ *
+ *              sb = voc_skybot  (ra, dec, rsz, dsz, epoch)
+ *
+ *      count = voc_skybotNObjs  (sb)
+ *      str = voc_skybotStrAttr  (sb, attrname, index)
+ *     dval = voc_skybotDblAttr  (sb, attrname, index)
+ *
+ *  Available Attributes:
+ *  ---------------------
+ *
+ *      number  string      Asteroid number
+ *      name    string      Asteroid name
+ *      ra      double      J2000 Equatorial RA
+ *      dec     double      J2000 Equatorial Dec
+ *      class   string      Object classification
+ *      vmag    double      Visual magnitude
+ *      poserr  double      Error on position (arcsec)
+ *      cdist   double      Body-to-center angular distance
+ *      dra     double      RA motion (arcsec/hr)
+ *      ddec    double      Dec motion (arcsec/hr)
+ *      dgeo    double      Geocentric distance (AU)
+ *      dhelio  double      Heliocentric distance (AU)
+ *      px      double      Mean J2000 heliocentric position vector (AU)
+ *      py      double          "       "       "       "       "
+ *      pz      double          "       "       "       "       "
+ *      vx      double      Mean J2000 heliocentric position vector (AU/day)
+ *      vy      double          "       "       "       "       "
+ *      vz      double          "       "       "       "       "       
+ *      JD0     double      T0, epoch of position vector (JD)
+ * 
+ *
+ *  Example Usage:      
+ *  ---------------
+ *
+ *  1) Find bodies in a 900" radius around (0.0,0.0) on JD 2453939.123
+ *
+ *      VOCSkybot skybot = new VOCSkybot (0.0,0.0,900,2453939.123);
+ *      for (int i=0; i < skybot.getNObjs(); i++)
+ *          System.out.print   ("Name: '" + sb.getStrAttr("name",i)+"' " +
+ *                              "RA: " + sb.getDblAttr("ra",i) + " " +
+ *                              "Dec: " + sb.getDblAttr("dec",i) + " " +
+ *                              "Mv: " + sb.getDblAttr("vmag",i) + " "); 
+ * 
+ *
+ *
+ *  @file       vocSkybot_f77.c
+ *  @author     Michael Fitzpatrick
+ *  @version    June 2006
+ *
+ **************************************************************************/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <unistd.h>
+#include <signal.h>
+#include <errno.h>
+
+#define _VOCLIENT_LIB_
+#include "VOClient.h"
+
+
+#ifdef _NO_US_
+
+#define VF_SKYBOT 		vfskybot
+#define VF_SKYBOTNOBJS 		vfskybotnobjs
+#define VF_SKYBOTSTR 		vfskybotstr
+#define VF_SKYBOTDBL 		vfskybotdbl
+
+#else
+
+#define VF_SKYBOT 		vfskybot_
+#define VF_SKYBOTNOBJS 		vfskybotnobjs_
+#define VF_SKYBOTSTR 		vfskybotstr_
+#define VF_SKYBOTDBL 		vfskybotdbl_
+
+#endif
+
+
+/*  Private interface procedures.
+*/
+extern char *sstrip (char *instr, int len);
+extern void  spad (char *outstr, int len);
+
+
+
+/**
+ *  VFSKYBOT --  Call the SkyBoT (Sky Bodies Tracker) service from IMCCE.
+ *  This service returns a list of minor planets withing the specified
+ *  search radius/box about the poition at the epoch given.  Epoch is 
+ *  assumed to be JD, ra/dec are J2000 decimal degrees, search size is 
+ *  given in arcsec.
+ *
+ *  @brief      Call the SkyBoT (Sky Bodies Tracker) service from IMCCE.
+ *  @fn         call vfskybot (double *ra, double *dec, double *rsz,
+ *                      double *dsz, double *epoch, int *sb)
+ *
+ *  @param  ra          RA position of query (decimal degrees)
+ *  @param  dec         Dec position of query (decimal degrees)
+ *  @param  rsz         RA size of query (arcsec)
+ *  @param  dsz         Dec size of query (arcsec)
+ *  @param  epoch       epoch of query (JD)
+ *  @param  sb          handle to Skybot return object
+ *  @returns            nothing
+ */
+void      
+VF_SKYBOT (double *ra, double *dec, double *rsz, double *dsz, double *epoch,
+	int *sb)
+{
+    *sb = (int) voc_skybot (*ra, *dec, *rsz, *dsz, *epoch);
+}
+
+
+/**
+ *  VFSKYBOTSTR -- Return a real-valued field as a string value. 
+ *  Notice the index begins at 1 and not 0.
+ *
+ *  @brief      Return a real-valued field as a string value.
+ *  @fn         call vfskybotstr (int *sb, char *attr, int *index, char *out)
+ *
+ *  @param  sb          handle to skybot query
+ *  @param  attr        attribute name
+ *  @param  index       object index
+ *  @param  out         output string value
+ *  @returns            nothing
+ */
+void
+VF_SKYBOTSTR (int *sb, char *attrname, int *index, char *outstr, int *len,
+	int alen, int olen)
+{
+    char *_attrname = sstrip (attrname, alen);
+
+    char *_result = voc_skybotStrAttr ((Skybot) *sb, _attrname, *index-1);
+    memset (outstr, 0, olen);
+    if ((*len = strlen(_result)) > olen)
+        fprintf (stderr, "Warning: truncating attr string: len=%d maxch=%d\n",
+            *len, olen);
+    spad (strncpy (outstr, _result, *len), olen);
+
+    free ((char *) _result);
+    free ((char *) _attrname);
+}
+
+
+/**
+ *  VFSKYBOTDBL -- Return a real-valued field as a double-precision value.
+ *  Notice the index begins at 1 and not 0.
+ *
+ *  @brief      Return a real-valued field as a double-precision value.
+ *  @fn         call vfskybotdbl (int *sb, char *attr, int *index, double *dval)
+ *
+ *  @param  sb          handle to skybot query
+ *  @param  attr        attribute name
+ *  @param  index       object index
+ *  @param  out         output string value
+ *  @returns            nothing
+ */
+void      
+VF_SKYBOTDBL (int *sb, char *attrname, int *index, double *dval, int alen)
+{
+    char *_attrname = sstrip (attrname, alen);
+
+    *dval = voc_skybotDblAttr (*sb, _attrname, *index-1);
+
+    free ((char *) _attrname);
+}
+
+
+/**
+ *  VFSKYBOTNOBJS -- Return the number of objects found from the query.
+ *
+ *  @brief      Return the number of objects found from the query.
+ *  @fn         call vfskybotnobjs (int *sb, int *nobjs)
+ *
+ *  @param  sb          handle to skybot query
+ *  @param  nobjs       number of objects in the query
+ *  @returns            nothing
+ */
+void      
+VF_SKYBOTNOBJS (int *sb, int *nobjs)
+{
+    *nobjs = voc_skybotNObjs (*sb);
+}