diff options
Diffstat (limited to 'vendor/voclient/libvoclient/vocSkybot_f77.c')
| -rw-r--r-- | vendor/voclient/libvoclient/vocSkybot_f77.c | 196 |
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); +} |