1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
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);
}
|
