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
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
|
#
# SBQUERY -- Call the IMCCE SkyBoT (Sky Bodies Tracker) service to locate
# minor planets within a radius at a given position and epoch.
#
include <time.h>
# Possible fields.
define SB_FIELDS "|num|name|class|ra|dec|vmag|poserr|cdist|dra|ddec|\
|dgeo|dhelio|px|py|pz|vx|vy|vz|jd0|"
define SB_NFIELDS 20 # accounts for the newline
define SB_KNUM "NUM" # keywords
define SB_KNAME "NAME"
define SB_KCLASS "CLASS"
define SB_KRA "RA"
define SB_KDEC "DEC"
define SB_KVMAG "VMAG"
define SB_KPOSERR "POSERR"
define SB_KCDIST "CDIST"
define SB_KDRA "DRA"
define SB_KDDEC "DDEC"
define SB_KDGEO "DGEO"
define SB_KDHELIO "DHELIO"
define SB_KPX "PX"
define SB_KPY "PY"
define SB_KPZ "PZ"
define SB_KVX "VX"
define SB_KVY "VY"
define SB_KVZ "VZ"
define SB_KJD0 "JD0"
define SB_LNUM "number" # attribute keywords
define SB_LNAME "name"
define SB_LCLASS "class"
define SB_LRA "ra"
define SB_LDEC "dec"
define SB_LVMAG "vmag"
define SB_LPOSERR "poserr"
define SB_LCDIST "cdist"
define SB_LDRA "dra"
define SB_LDDEC "ddec"
define SB_LDGEO "dgeo"
define SB_LDHELIO "dhelio"
define SB_LPX "px"
define SB_LPY "py"
define SB_LPZ "pz"
define SB_LVX "vx"
define SB_LVY "vy"
define SB_LVZ "vz"
define SB_LJD0 "jd0"
define SB_FNUM 1 # field numbers
define SB_FNAME 2
define SB_FCLASS 3
define SB_FRA 4
define SB_FDEC 5
define SB_FVMAG 6
define SB_FPOSERR 7
define SB_FCDIST 8
define SB_FDRA 9
define SB_FDDEC 10
# NL
define SB_FDGEO 12
define SB_FDHELIO 13
define SB_FPX 14
define SB_FPY 15
define SB_FPZ 16
define SB_FVX 17
define SB_FVY 18
define SB_FVZ 19
define SB_FJD0 20
define SB_FREAL "%-10.4f"
define SB_FWSTRING "%-12s"
define SB_FSTRING "%-10s"
define SB_FINTEGER "%10d"
define SB_FSEX "%-012h"
procedure t_sbquery ()
pointer sp, fields, fieldstr, epoch
double ra, dec, sr, ut, jd_epoch
int sb, nobjs, nfields, i, gmtco
int tm[LEN_TMSTRUCT]
bool fmt
long otime, ntime
real clgetr()
bool clgetb(), streq()
int vx_initVOClient(), vx_skybotnobjs(), vx_skybot()
int sb_fields()
long clktime()
double ast_date_to_julday()
begin
call smark (sp)
call salloc (epoch, SZ_FNAME, TY_CHAR)
call salloc (fields, SZ_LINE, TY_INT)
call salloc (fieldstr, SZ_LINE, TY_CHAR)
call aclrs (Memc[epoch], SZ_FNAME)
call aclri (Memc[fields], SZ_LINE)
call aclrs (Memc[fieldstr], SZ_LINE)
ra = clgetr ("ra") # Get the parameters
dec = clgetr ("dec")
sr = clgetr ("sr")
call clgstr ("epoch", Memc[epoch], SZ_LINE)
if (streq (Memc[epoch], "now")) {
# Handle the special case of 'now', i.e. find out what time it is.
call zgmtco (gmtco)
call brktime (clktime(0)+gmtco, tm)
ut = real (TM_HOUR(tm)) +
(real(TM_MIN(tm))/60.) +
(real(TM_SEC(tm))/3600.)
jd_epoch = ast_date_to_julday (TM_YEAR(tm), TM_MONTH(tm),
TM_MDAY(tm), ut)
} else {
# Convert the user-supplied epoch.
call sscan (Memc[epoch])
call gargd (jd_epoch)
}
call clgstr ("fields", Memc[fieldstr], SZ_LINE)
fmt = clgetb ("fmt")
if (vx_initVOClient("") == ERR) { # Initialize
call clputi ("status", ERR)
return
}
if (fmt) {
call printf ("#\n# Search Terms: ")
call printf ("Position=(%H,%h) sr=%.1f epoch=%.4f\n")
call pargd (ra)
call pargd (dec)
call pargd (sr)
call pargd (jd_epoch)
call printf ("#\n")
call flush (STDOUT)
}
# Call the service
otime = clktime (0)
sb = vx_skybot (ra, dec, sr, sr, jd_epoch)
ntime = clktime (otime)
if (sb != NULL) {
nobjs = vx_skybotnobjs (sb)
if (fmt) {
call printf ("# Found %d objects in %d seconds\n#\n")
call pargi (nobjs)
call pargl (ntime)
}
} else {
call eprintf ("Error in calling SkyBoT service\n")
call sfree (sp)
call clputi ("status", ERR)
call vx_closeVOClient (0)
return
}
# Get the selected fields.
nfields = sb_fields (Memc[fieldstr], Memi[fields], SB_NFIELDS)
if (nfields <= 0) {
call vx_closeVOClient (0)
call sfree (sp)
return
}
# Finally, print the table of results.
if (fmt)
call sb_pheader (Memi[fields], nfields)
for (i=0; i < nobjs; i=i+1)
call sb_print (sb, i, Memi[fields], nfields)
call vx_closeVOClient (0)
call clputi ("status", OK)
end
# SB_FIELDS -- Procedure to decode the fields string into a list of the
# fields to be printed.
int procedure sb_fields (fieldstr, fields, max_nfields)
char fieldstr[ARB] #I string containing the list of fields
int fields[ARB] #O fields array
int max_nfields #I maximum number of fields
int nfields, flist, field
pointer sp, fname
int fntopenb(), fntgfnb(), strdic()
begin
nfields = 0
call smark (sp)
call salloc (fname, SZ_FNAME, TY_CHAR)
flist = fntopenb (fieldstr, NO)
while (fntgfnb (flist, Memc[fname], SZ_FNAME) != EOF &&
(nfields < max_nfields)) {
field = strdic (Memc[fname], Memc[fname], SZ_FNAME, SB_FIELDS)
if (field == 0)
next
nfields = nfields + 1
fields[nfields] = field
}
call fntclsb (flist)
call sfree (sp)
return (nfields)
end
# SB_PHEADER -- Print the banner fields.
procedure sb_pheader (fields, nfields)
int fields[ARB] # fields to be printed
int nfields # number of fields
int i
begin
call printf ("#")
do i = 1, nfields {
switch (fields[i]) {
case SB_FNUM:
call printf (SB_FSTRING)
call pargstr (SB_KNUM)
case SB_FNAME:
call printf (SB_FWSTRING)
call pargstr (SB_KNAME)
case SB_FCLASS:
call printf (SB_FSTRING)
call pargstr (SB_KCLASS)
case SB_FRA:
call printf (SB_FWSTRING)
call pargstr (SB_KRA)
case SB_FDEC:
call printf (SB_FWSTRING)
call pargstr (SB_KDEC)
case SB_FVMAG:
call printf (SB_FSTRING)
call pargstr (SB_KVMAG)
case SB_FPOSERR:
call printf (SB_FSTRING)
call pargstr (SB_KPOSERR)
case SB_FCDIST:
call printf (SB_FSTRING)
call pargstr (SB_KCDIST)
case SB_FDRA:
call printf (SB_FSTRING)
call pargstr (SB_KDRA)
case SB_FDDEC:
call printf (SB_FSTRING)
call pargstr (SB_KDDEC)
case SB_FDGEO:
call printf (SB_FSTRING)
call pargstr (SB_KDGEO)
case SB_FDHELIO:
call printf (SB_FSTRING)
call pargstr (SB_KDHELIO)
case SB_FPX:
call printf (SB_FSTRING)
call pargstr (SB_KPX)
case SB_FPY:
call printf (SB_FSTRING)
call pargstr (SB_KPY)
case SB_FPZ:
call printf (SB_FSTRING)
call pargstr (SB_KPZ)
case SB_FVX:
call printf (SB_FSTRING)
call pargstr (SB_KVX)
case SB_FVY:
call printf (SB_FSTRING)
call pargstr (SB_KVY)
case SB_FVZ:
call printf (SB_FSTRING)
call pargstr (SB_KVZ)
case SB_FJD0:
call printf (SB_FSTRING)
call pargstr (SB_KJD0)
}
}
call printf ("\n")
call flush (STDOUT)
end
# SB_PRINT -- Print the banner fields.
procedure sb_print (sb, objnum, fields, nfields)
int sb # skybot pointer
int objnum # skybot result number
int fields[ARB] # fields to be printed
int nfields # number of fields
int i, len
char buf[SZ_LINE]
double dval
int vx_skybotstr()
double vx_skybotdbl()
begin
do i = 1, nfields {
call aclrs (buf, SZ_LINE)
switch (fields[i]) {
case SB_FNUM:
len = vx_skybotstr (sb, SB_LNUM, objnum, buf, SZ_FNAME)
call printf (SB_FSTRING)
call pargstr (buf)
case SB_FNAME:
len = vx_skybotstr (sb, SB_LNAME, objnum, buf, SZ_FNAME)
call printf (SB_FWSTRING)
call pargstr (buf)
case SB_FCLASS:
len = vx_skybotstr (sb, SB_LCLASS, objnum, buf, SZ_FNAME)
call printf (SB_FSTRING)
call pargstr (buf)
case SB_FRA:
dval = vx_skybotdbl (sb, SB_LRA, objnum)
call printf (SB_FSEX)
call pargd (dval)
case SB_FDEC:
dval = vx_skybotdbl (sb, SB_LDEC, objnum)
call printf (SB_FSEX)
call pargd (dval)
case SB_FVMAG:
dval = vx_skybotdbl (sb, SB_LVMAG, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FPOSERR:
dval = vx_skybotdbl (sb, SB_LPOSERR, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FCDIST:
dval = vx_skybotdbl (sb, SB_LCDIST, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FDRA:
dval = vx_skybotdbl (sb, SB_LDRA, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FDDEC:
dval = vx_skybotdbl (sb, SB_LDDEC, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FDGEO:
dval = vx_skybotdbl (sb, SB_LDGEO, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FDHELIO:
dval = vx_skybotdbl (sb, SB_LDHELIO, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FPX:
dval = vx_skybotdbl (sb, SB_LPX, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FPY:
dval = vx_skybotdbl (sb, SB_LPY, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FPZ:
dval = vx_skybotdbl (sb, SB_LPZ, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FVX:
dval = vx_skybotdbl (sb, SB_LVX, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FVY:
dval = vx_skybotdbl (sb, SB_LVY, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FVZ:
dval = vx_skybotdbl (sb, SB_LVZ, objnum)
call printf (SB_FREAL)
call pargd (dval)
case SB_FJD0:
dval = vx_skybotdbl (sb, SB_LJD0, objnum)
call printf (SB_FREAL)
call pargd (dval)
}
}
call printf ("\n")
call flush (STDOUT)
end
|
