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
|
include "../lib/apphotdef.h"
include "../lib/polyphot.h"
include "../lib/fitpsf.h"
include "../lib/radprof.h"
# AP_WPROFS -- Procedure to print out the profile fitting parameters.
procedure ap_wprofs (ap, out)
pointer ap # apphot structure pointer
int out # output file descriptor
pointer sp, str
int apstati()
real apstatr()
begin
if (out == NULL)
return
call smark (sp)
call salloc (str, SZ_FNAME, TY_CHAR)
if (AP_RPROF(ap) != NULL) {
call ap_rparam (out, KY_RPRADIUS, apstatr (ap, RPRADIUS),
UN_RSCALEUNIT, "fitting radius")
call ap_rparam (out, KY_RPSTEP, apstatr (ap, RPSTEP),
UN_RSCALEUNIT, "step size in radius")
call ap_iparam (out, KY_RPORDER, apstati (ap, RPORDER), UN_RNUMBER,
"number of splines pieces")
call ap_rparam (out, KY_RPKSIGMA, apstatr (ap, RPKSIGMA),
UN_RSIGMA, "k-sigma rejection criterion")
call ap_iparam (out, KY_RPNREJECT, apstati (ap, RPNREJECT),
UN_RNUMBER, "maximum number of rejection cycles")
call fprintf (out, "#\n")
}
call sfree (sp)
end
# AP_WPOLY -- Procedure to write out the polyphot parameters.
procedure ap_wpoly (ap, out)
pointer ap # apphot structure pointer
int out # output file descriptor
pointer sp, str
real apstatr()
begin
if (out == NULL)
return
call smark (sp)
call salloc (str, SZ_FNAME, TY_CHAR)
if (AP_POLY(ap) != NULL) {
call ap_sparam (out, "WEIGHTING", "constant", "model", "")
call ap_rparam (out, KY_PYZMAG, apstatr (ap, PYZMAG), UN_PYZMAG,
"zero point of magnitdue scale")
call fprintf (out, "#\n")
}
call sfree (sp)
end
# AP_WPSF -- Procedure to write the psf fitting parameters .
procedure ap_wpsf (ap, out)
pointer ap # apphot strucuture pointer
int out # output file descriptor
pointer sp, str
int apstati()
real apstatr()
begin
if (out == NULL)
return
call smark (sp)
call salloc (str, SZ_FNAME, TY_CHAR)
if (AP_PPSF(ap) != NULL) {
call apstats (ap, PSFSTRING, Memc[str], SZ_FNAME)
call ap_sparam (out, KY_PSFUNCTION, Memc[str], UN_PSFMODEL,
"fitting function")
call ap_rparam (out, KY_PSFAPERT, 2.0 * apstatr (ap, PSFAPERT),
UN_PSFSCALEUNIT, "width of the fitting box")
call ap_iparam (out, KY_PMAXITER, apstati (ap, PMAXITER),
UN_PSFSCALEUNIT, "maximum number of iterations")
call ap_rparam (out, KY_PK2, apstatr (ap, PK2), UN_PSFSIGMA,
"k-sigma rejection limit for the fit")
call ap_iparam (out, KY_PNREJECT, apstati (ap, PNREJECT),
UN_PSFNUMBER, "maximum number of rejection cycles")
call fprintf (out, "#\n")
}
call sfree (sp)
end
|
