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
|
include "../lib/apphotdef.h"
include "../lib/radprofdef.h"
include "../lib/apphot.h"
include "../lib/radprof.h"
# define the #N, #U and #K radprof strings
define RP_NSTR1 "#N%4tPFWHM%14tINORM%29tTINORM%66tRIER%71tRERROR%80t\\\n"
define RP_USTR1 "#U%4tScale%14tcounts%29tcounts%66t##%71trerrors%80t\\\n"
define RP_FSTR1 "#F%4t%%-13.3f%14t%%-15.7f%29t%%-36.7f%66t%%-5d%71t\
%%-9s%80t \n"
define RP_WSTR1 "%4t%-10.3f%-15.7g%-36.7g%-5d%-9.9s%80t\\\n"
define RP_NSTR2 "#N%4tPRADIUS%14tINTENSITY%29tTINTENSITY%80t\\\n"
define RP_USTR2 "#U%4tscale%14tcounts%29tcounts%80t\\\n"
define RP_FSTR2 "#F%4t%%-13.3f%14t%%-15.7f%29t%%-15.7f%80t \n"
define RP_WSTR2 "%4t%-10.3f%-15.7g%-15.7g%79t%2s\n"
# AP_RHDR -- Print the radprof column header strings.
procedure ap_rhdr (ap, fd)
pointer ap # apphot descriptor (unused)
int fd # output file descriptor
begin
if (fd == NULL)
return
call fprintf (fd, RP_NSTR1)
call fprintf (fd, RP_USTR1)
call fprintf (fd, RP_FSTR1)
call fprintf (fd, "#\n")
call fprintf (fd, RP_NSTR2)
call fprintf (fd, RP_USTR2)
call fprintf (fd, RP_FSTR2)
call fprintf (fd, "#\n")
end
# AP_WRRES -- Write the results of the radprof task to the output file.
procedure ap_wrres (ap, fd, ier)
pointer ap # pointer to apphot structure
int fd # output text file descriptor
int ier # radial profile error
int i, nrpts
pointer sp, str, rprof
real apstatr()
begin
# Initialize.
if (fd == NULL)
return
rprof = AP_RPROF(ap)
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
# Print the radprof parameters.
call ap_srerrors (ier, Memc[str], SZ_LINE)
call fprintf (fd, RP_WSTR1)
call pargr (apstatr (ap, RPFWHM) / apstatr (ap, SCALE))
call pargr (apstatr (ap, INORM))
call pargr (apstatr (ap, TNORM))
call pargi (ier)
call pargstr (Memc[str])
# Print the radial profile.
nrpts = apstatr (ap, RPRADIUS) / apstatr (ap, RPSTEP) + 1
if (nrpts == 0) {
call fprintf (fd, RP_WSTR2)
call pargr (INDEFR)
call pargr (INDEFR)
call pargr (INDEFR)
call pargstr (" ")
} else {
do i = 1, nrpts {
if (nrpts == 1) {
call fprintf (fd, RP_WSTR2)
call pargr (Memr[AP_RPDIST(rprof)+i-1] / AP_SCALE(ap))
call pargr (Memr[AP_INTENSITY(rprof)+i-1])
call pargr (Memr[AP_TINTENSITY(rprof)+i-1])
call pargstr (" ")
} if (i == nrpts) {
call fprintf (fd, RP_WSTR2)
call pargr (Memr[AP_RPDIST(rprof)+i-1] / AP_SCALE(ap))
call pargr (Memr[AP_INTENSITY(rprof)+i-1])
call pargr (Memr[AP_TINTENSITY(rprof)+i-1])
call pargstr ("* ")
} else {
call fprintf (fd, RP_WSTR2)
call pargr (Memr[AP_RPDIST(rprof)+i-1] / AP_SCALE(ap))
call pargr (Memr[AP_INTENSITY(rprof)+i-1])
call pargr (Memr[AP_TINTENSITY(rprof)+i-1])
call pargstr ("*\\")
}
}
}
call sfree (sp)
end
# AP_SRERRORS -- Encode the radial profile error message in a string.
procedure ap_srerrors (ier, str, maxch)
int ier # error code
char str[ARB] # encoded error string
int maxch # maximum number of characters
begin
switch (ier) {
case AP_RP_NOPROFILE:
call sprintf (str, maxch, "%s")
call pargstr ("OffImage")
case AP_RP_OUTOFBOUNDS:
call sprintf (str, maxch, "%s")
call pargstr ("EdgeImage")
case AP_RP_NPTS_TOO_SMALL:
call sprintf (str, maxch, "%s")
call pargstr ("TooFewPts")
case AP_RP_SINGULAR:
call sprintf (str, maxch, "%s")
call pargstr ("Singular")
default:
call sprintf (str, maxch, "%s")
call pargstr ("NoError")
}
end
|
