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
|
include "../lib/apseldef.h"
# DP_TADINIT -- Initializie the column descriptors for the input ST photometry
# table.
procedure dp_tadinit (tp, column)
pointer tp # table descriptor
int column[ARB] # column pointer array
begin
# Find the column pointers
call tbcfnd (tp, ID, column[1], 1)
if (column[1] == NULL)
call tbcfnd (tp, "ID", column[1], 1)
if (column[1] == NULL)
call printf ("Error reading ID.\n")
call tbcfnd (tp, XCENTER, column[2], 1)
if (column[2] == NULL)
call tbcfnd (tp, "XCENTER", column[2], 1)
if (column[2] == NULL)
call printf ("Error reading XCENTER.\n")
call tbcfnd (tp, YCENTER, column[3], 1)
if (column[3] == NULL)
call tbcfnd (tp, "YCENTER", column[3], 1)
if (column[3] == NULL)
call printf ("Error reading YCENTER.\n")
call tbcfnd (tp, MAG, column[4], 1)
if (column[4] == NULL)
call tbcfnd (tp, APMAG, column[4], 1)
if (column[4] == NULL)
call printf ("Error reading MAG.\n")
end
# DP_TADREAD -- Read a record from the input ST photometry table.
procedure dp_tadread (tp, column, id, x, y, mag, row)
pointer tp # table descriptor
int column[ARB] # column pointer array
int id # output id
real x # output x value
real y # output y value
real mag # output magnitude
int row # integer row
bool nullflag
begin
call tbrgti (tp, column[1], id, nullflag, 1, row)
if (nullflag)
id = 0
call tbrgtr (tp, column[2], x, nullflag, 1, row)
call tbrgtr (tp, column[3], y, nullflag, 1, row)
call tbrgtr (tp, column[4], mag, nullflag, 1, row)
end
# DP_GCOORDS -- Read the coordinates and magnitudes from a simple text file.
int procedure dp_gcoords (cl, x, y, mag, id)
int cl # file descriptor
real x # x coordinate centers
real y # y coordinate centers
real mag # magnitudes
int id # id of the star
int fscan(), nscan()
begin
while (fscan (cl) != EOF) {
call gargr (x)
call gargr (y)
call gargr (mag)
if (nscan () < 3)
next
call gargi (id)
if (nscan() < 4)
id = 0
return (1)
}
return (EOF)
end
|
