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
|
C
C F77CONE -- Fortran example task of VOClient interface
C
C M.Fitzpatrick, NOAO, Jul 2006
program f77cone
double precision ra, dec, sr
character service*80
C Initialize values. Note the VOClient library assumes all floating
C point values are double precision.
ra = 12.0d0
dec = 12.0d0
sr = 0.1d0
service = "http://www.nofs.navy.mil/cgi-bin/vo_cone.cgi?CAT=USNO-B1&"
call coneservice (service, ra, dec, sr)
stop
end
C Simple test routine to call a Cone search service and summarize results.
subroutine coneservice (service, ra, dec, sr)
character*80 service, fmt, sid, sra, sdec
character*160 attrlist, qstring
integer nrec, nattr, len, ier, index
integer cone, query, qr, rec
double precision ra, dec, sr
C Initialize the VOClient
call vfinitvoclient ("", ier)
C Get a new connection to the named service and form a query.
call vfopenconeconnection (service, cone, ier)
call vfgetconequery (cone, ra, dec, sr, query)
C Print the query string we're about to execute.
index = 1
call vfgetquerystring (query, 3, index, qstring, len)
print *, "Executing Query: "
print *, qstring
print *, ""
C Execute the query, get back a response handle 'qr'.
call vfexecquery (query, qr)
C Now summarize the response.
call vfgetrecordcount (qr, nrec)
if (nrec <= 0) then
print *, "No records matched"
return
endif
C Get a sample record so we can summarize columns.
call vfgetrecord (qr, 1, rec)
if (rec > 0) then
call vfgetattrcount (rec, nattr)
endif
print *, "# returns ", nrec, " records of ", nattr, " attrs"
print *, "#"
print *, "# --- Summary output ---"
print *, "#"
C Get the list of attributes. Note that when receiving a string value
C the interface requires you specify the max length of the string and
C will return the actual length of the string.
call vfgetattrlist (rec, attrlist, len)
print *, "# ", attrlist
print *, "#"
C Summarize and print selected query results. Note the results table
C is one-indexed in this interface. We also show how to catch an invalid
C record descriptor. To try this, change the starting index to 0 instead
C of 1.
do i = 1, nrec, 1
C Get the row in the table and then select values to print
call vfgetrecord (qr, i, rec)
if (rec > 0) then
call vfgetstringattr (rec, "ID_MAIN", sid, len)
call vfgetstringattr (rec, "POS_EQ_RA_MAIN", sra, len)
call vfgetstringattr (rec, "POS_EQ_DEC_MAIN", sdec, len)
else
print *, "Error getting record number: ", i
goto 100
endif
fmt = '(i3,a6,a18,a6,a18,a6,a18)'
write (*, fmt) i, " id= ", sid, " ra= ", sra, " dec= ", sdec
enddo
100 continue
C Close the connection and shut down the VOClient server.
call vfcloseconnection (cone)
call vfclosevoclient (1, ier)
return
end
|
