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
|
include "../lib/display.h"
include "../lib/center.h"
include "../lib/fitsky.h"
include "../lib/phot.h"
include "../lib/polyphot.h"
# AP_CEPARS -- Procedure to write out the current CENTERPARS parameters
# to the current CENTERPARS parameter file.
procedure ap_cepars (ap)
pointer ap # pointer to apphot structure
pointer sp, str, cp
bool itob()
int apstati()
pointer clopset()
real apstatr()
begin
# Open the parameter set.
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
cp = clopset ("centerpars")
# Write the centering parameters.
call apstats (ap, CSTRING, Memc[str], SZ_FNAME)
call clppset (cp, "calgorithm", Memc[str])
call clppsetr (cp, "cbox", 2.0 * apstatr (ap, CAPERT))
call clppsetr (cp, "cthreshold", apstatr (ap, CTHRESHOLD))
call clppsetr (cp, "minsnratio", apstatr (ap, MINSNRATIO))
call clppseti (cp, "cmaxiter", apstati (ap, CMAXITER))
call clppsetr (cp, "maxshift", apstatr (ap, MAXSHIFT))
call clppsetb (cp, "clean", itob (apstati (ap, CLEAN)))
call clppsetr (cp, "rclean", apstatr (ap, RCLEAN))
call clppsetr (cp, "rclip", apstatr (ap, RCLIP))
call clppsetr (cp, "kclean", apstatr (ap, SIGMACLEAN))
call clppsetb (cp, "mkcenter", itob (apstati (ap, MKCENTER)))
# Close the pset file.
call clcpset (cp)
call sfree (sp)
end
# AP_SAPARS -- Procedure to write out the current FITSKYPARS parameters
# to the FITSKYPARS file.
procedure ap_sapars (ap)
pointer ap # pointer to apphot structure
pointer mp, str, sp
bool itob()
int apstati()
pointer clopset()
real apstatr()
begin
# Open the parameter sets.
call smark (mp)
call salloc (str, SZ_LINE, TY_CHAR)
sp = clopset ("fitskypars")
# Set the sky fitting parameters.
call apstats (ap, SSTRING, Memc[str], SZ_FNAME)
call clppset (sp, "salgorithm", Memc[str])
call clppsetr (sp, "annulus", apstatr (ap, ANNULUS))
call clppsetr (sp, "dannulus", apstatr (ap, DANNULUS))
call clppsetr (sp, "skyvalue", apstatr (ap, SKY_BACKGROUND))
call clppsetr (sp, "khist", apstatr (ap, K1))
call clppsetr (sp, "binsize", apstatr (ap, BINSIZE))
call clppsetb (sp, "smooth", itob (apstati (ap, SMOOTH)))
call clppsetr (sp, "sloclip", apstatr (ap, SLOCLIP))
call clppsetr (sp, "shiclip", apstatr (ap, SHICLIP))
call clppseti (sp, "smaxiter", apstati (ap, SMAXITER))
call clppseti (sp, "snreject", apstati (ap, SNREJECT))
call clppsetr (sp, "sloreject", apstatr (ap, SLOREJECT))
call clppsetr (sp, "shireject", apstatr (ap, SHIREJECT))
call clppsetr (sp, "rgrow", apstatr (ap, RGROW))
call clppsetb (sp, "mksky", itob (apstati (ap, MKSKY)))
# Close up the pset files.
call clcpset (sp)
call sfree (mp)
end
# AP_PHPARS -- Procedure to write out the PHOTPARS parameters to the
# PHOTPARS task.
procedure ap_phpars (ap)
pointer ap # pointer to apphot structure
pointer mp, str, pp
bool itob()
int apstati()
pointer clopset()
real apstatr()
begin
# Open the parameter set.
call smark (mp)
call salloc (str, SZ_LINE, TY_CHAR)
pp = clopset ("photpars")
# Set the photometry parameters.
call apstats (ap, APERTS, Memc[str], SZ_LINE)
call clppset (pp, "apertures", Memc[str])
call clppsetr (pp, "zmag", apstatr (ap, ZMAG))
call apstats (ap, PWSTRING, Memc[str], SZ_FNAME)
call clppset (pp, "weighting", Memc[str])
call clppsetb (pp, "mkapert", itob (apstati (ap, MKAPERT)))
# Close the pset file.
call clcpset (pp)
call sfree (mp)
end
# AP_POPARS -- Procedure to write the current POLYPARS parameters to the
# current POLYPARS parameter file.
procedure ap_popars (ap)
pointer ap # pointer to apphot structure
pointer pp
bool itob()
int apstati()
pointer clopset()
real apstatr()
begin
# Open the psets.
pp = clopset ("polypars")
# Set the photometry parameters.
call clppsetr (pp, "zmag", apstatr (ap, PYZMAG))
call clppsetb (pp, "mkpolygon", itob (apstati (ap, MKPOLYGON)))
# Close the pset files.
call clcpset (pp)
end
|
