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
|
include <fset.h>
# T_PFMERGE -- Merge photometry files written by PHOT, PSF, PEAK, GROUP,
# NSTAR or ALLSTAR into one photometry file. The fields ID, XCENTER,
# YCENTER, MAG, and MKSKY are read from each input file and written
# without change to the output file. The input files may be text files or
# binary ST tables files. The output file will have the same file type as the
# the first input file. The header of the output file will be the header
# of the first input file.
procedure t_pfmerge ()
pointer inphotfiles # the input photometry files
pointer outphotfile # the output photometry file
bool verbose # verbose output ?
int plist, lplist, infd, outfd, first_file, in_text, out_text
pointer sp, infname
bool clgetb()
int fstati(), fntopnb(), fntlenb(), fntgfnb(), access(), open()
pointer tbtopn()
begin
# Set the standard output to flush on newline.
if (fstati (STDOUT, F_REDIR) == NO)
call fseti (STDOUT, F_FLUSHNL, YES)
# Get some working memory.
call smark (sp)
call salloc (inphotfiles, SZ_FNAME, TY_CHAR)
call salloc (outphotfile, SZ_FNAME, TY_CHAR)
call salloc (infname, SZ_FNAME, TY_CHAR)
# Get the various task parameters.
call clgstr ("inphotfiles", Memc[inphotfiles], SZ_FNAME)
call clgstr ("outphotfile", Memc[outphotfile], SZ_FNAME)
verbose = clgetb ("verbose")
# Open the input file list and determine the file type of the first
# file.
plist = fntopnb (Memc[inphotfiles], NO)
lplist = fntlenb (plist)
if (lplist <= 0) {
call fntclsb (plist)
call sfree (sp)
return
}
if (fntgfnb (plist, Memc[infname], SZ_FNAME) == EOF) {
call fntclsb (plist)
call sfree (sp)
return
} else {
out_text = access (Memc[infname], 0, TEXT_FILE)
call fntrewb (plist)
}
# Open the output file.
if (out_text == YES)
outfd = open (Memc[outphotfile], NEW_FILE, TEXT_FILE)
else
outfd = tbtopn (Memc[outphotfile], NEW_FILE, 0)
# Loop over the list of input files
first_file = YES
while (fntgfnb (plist, Memc[infname], SZ_FNAME) != EOF) {
# Print message.
if (verbose) {
call printf ("Merging photometry file %s into %s\n")
call pargstr (Memc[infname])
call pargstr (Memc[outphotfile])
}
# Open the input file.
in_text = access (Memc[infname], 0, TEXT_FILE)
if (verbose && (in_text != out_text)) {
if (in_text == YES)
call eprintf ("File %s is not an ST table file.\n")
else
call eprintf ("File %s is not a text file.\n")
call pargstr (Memc[infname])
}
if (in_text == YES)
infd = open (Memc[infname], READ_ONLY, TEXT_FILE)
else
infd = tbtopn (Memc[infname], READ_ONLY, 0)
# Now merge the files.
call dp_pfmerge (infd, outfd, in_text, out_text, first_file)
# Close the photometry file.
if (in_text == YES)
call close (infd)
else
call tbtclo (infd)
first_file = NO
}
# Close the output file.
if (out_text == YES)
call close (outfd)
else
call tbtclo (outfd)
# Close the input file list.
call fntclsb (plist)
call sfree(sp)
end
|
