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
|
include <fset.h>
include "../lib/daophotdef.h"
# T_GRPSELECT -- Select groups from a GROUP file on the basis of the size
# of the group. Only groups of the sizes specified are copied into the
# output table.
procedure t_grpselect ()
pointer ingroup # the input GROUP file
pointer outgroup # the output GROUP file
int min_group # the minimum group size
int max_group # the maximum group size
bool gr_text
int ilist, lilist, olist, lolist, verbose
pointer sp, tp_in, tp_out, dao
bool clgetb(), itob()
int open(), tbtopn(), clgeti(), fstati(), btoi(), access()
int fntopnb(), fntlenb(), fntgfnb()
begin
# Set the standard output to flush on newline.
if (fstati (STDOUT, F_REDIR) == NO)
call fseti (STDOUT, F_FLUSHNL, YES)
# Get some memory.
call smark (sp)
call salloc (ingroup, SZ_FNAME, TY_CHAR)
call salloc (outgroup, SZ_FNAME, TY_CHAR)
# Get the various task parameters.
call clgstr ("ingroupfile", Memc[ingroup], SZ_FNAME)
call clgstr ("outgroupfile", Memc[outgroup], SZ_FNAME)
min_group = clgeti ("min_group")
max_group = clgeti ("max_group")
verbose = btoi (clgetb ("verbose"))
# Open the daophot structure.
call dp_init (dao)
# Open the photometry structure.
call dp_apsetup (dao)
# Set some parameters.
call dp_seti (dao, VERBOSE, verbose)
# Get the lists.
ilist = fntopnb (Memc[ingroup], NO)
lilist = fntlenb (ilist)
olist = fntopnb (Memc[outgroup], NO)
lolist = fntlenb (olist)
# Check the list lengths.
if (lilist != lolist) {
call fntclsb (ilist)
call fntclsb (olist)
call sfree (sp)
call error (0,
"The input and output list lengths are not compatible")
}
# Loop over the files.
while ((fntgfnb (ilist, Memc[ingroup], SZ_FNAME) != EOF) &&
(fntgfnb (olist, Memc[outgroup], SZ_FNAME) != EOF)) {
# Open the input file.
gr_text = itob (access (Memc[ingroup], 0, TEXT_FILE))
if (gr_text)
tp_in = open (Memc[ingroup], READ_ONLY, TEXT_FILE)
else
tp_in = tbtopn (Memc[ingroup], READ_ONLY, 0)
# Open an output file of the same type as the input file.
if (gr_text)
tp_out = open (Memc[outgroup], NEW_FILE, TEXT_FILE)
else
tp_out = tbtopn (Memc[outgroup], NEW_COPY, tp_in)
# Read in the groups and select by group size.
call dp_sgroup (dao, tp_in, tp_out, gr_text, min_group, max_group)
# Close the input and output files.
if (gr_text) {
call close (tp_in)
call close (tp_out)
} else {
call tbtclo (tp_in)
call tbtclo (tp_out)
}
}
# Close the lists.
call fntclsb (ilist)
call fntclsb (olist)
# Free the photometry structure.
call dp_apclose (dao)
# Free the daophot structure.
call dp_free (dao)
call sfree (sp)
end
|
