aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/apphot/aplib/apset2.x
blob: 85200dafbf3edbeb64b3295bc97537f6d53ddba9 (plain) (blame)
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
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
include "../lib/apphotdef.h"
include "../lib/apphot.h"
include "../lib/displaydef.h"
include "../lib/display.h"
include "../lib/noisedef.h"
include "../lib/noise.h"
include "../lib/polyphotdef.h"
include "../lib/polyphot.h"
include "../lib/radprofdef.h"
include "../lib/radprof.h"
include "../lib/finddef.h"
include "../lib/find.h"

# AP2SETS -- Procedure to set an apphot string parameter.

procedure ap2sets (ap, param, str)

pointer	ap		# pointer to apphot structure
int	param		# parameter
char	str[ARB]	# string parameter

pointer	nse, ply

begin
	nse = AP_NOISE(ap)
	ply = AP_POLY(ap)

	switch (param) {
	case PYNAME:
	    call strcpy (str, AP_PYNAME(ply), SZ_FNAME)
	case PYROOT:
	    call strcpy (str, AP_PYROOT(ply), SZ_FNAME)
	case GAIN:
	    call strcpy (str, AP_GAIN(nse), SZ_FNAME)
	case NSTRING:
	    call strcpy (str, AP_NSTRING(nse), SZ_FNAME)
	case CCDREAD:
	    call strcpy (str, AP_CCDREAD(nse), SZ_FNAME)
	default:
	    call error (0, "APSETS: Unknown apphot string parameter")
	}
end


# AP2SETI -- Procedure to set an integer apphot parameter.

procedure ap2seti (ap, param, ival)

pointer	ap		# pointer to apphot structure
int	param		# parameter
int	ival		# integer value

pointer	dsp, nse, ply, rprof

begin
	nse = AP_NOISE(ap)
	dsp = AP_PDISPLAY(ap)
	ply = AP_POLY(ap)
	rprof = AP_RPROF(ap)

	switch (param) {
	case RPORDER:
	    AP_RPORDER(rprof) = ival
	case RPNREJECT:
	    AP_RPNREJECT(rprof) = ival
	case PYNVER:
	    AP_PYNVER(ply) = ival
	case PYBADPIX:
	    AP_PYBADPIX(ply) = ival
	case MKSKY:
	    AP_MKSKY(dsp) = ival
	case MKCENTER:
	    AP_MKCENTER(dsp) = ival
	case MKAPERT:
	    AP_MKAPERT(dsp) = ival
	case MKPOLYGON:
	    AP_MKPOLYGON(dsp) = ival
	case MKDETECTIONS:
	    AP_MKDETECTIONS(dsp) = ival
	case NOISEFUNCTION:
	    AP_NOISEFUNCTION(nse) = ival
	case MKPSFBOX:
	    AP_MKPSFBOX(dsp) = ival
	case RADPLOTS:
	    AP_RADPLOTS(dsp) = ival
	case RPNPTS:
	    AP_RPNPTS(rprof) = ival
	case RPNDATA:
	    AP_RPNDATA(rprof) = ival
	case RPNDATAREJ:
	    AP_RPNDATAREJ(rprof) = ival
	default:
	    call error (0, "APSETI: Unknown apphot integer parameter")
	}
end


# AP2SETR -- Procedure to set a real apphot parameter.

procedure ap2setr (ap, param, rval)

pointer	ap		# pointer to apphot structure
int	param		# parameter
real	rval		# real value

pointer	nse, ply, rprof, fnd

begin
	nse = AP_NOISE(ap)
	ply = AP_POLY(ap)
	rprof = AP_RPROF(ap)
	fnd = AP_PFIND(ap)

	switch (param) {
	case RPFWHM:
	    AP_RPFWHM(rprof) = rval
	case INORM:
	    AP_INORM(rprof) = rval
	case TNORM:
	    AP_TINORM(rprof) = rval
	case DNORM:
	    AP_DNORM(rprof) = rval
	case RPXCUR:
	    AP_RPXCUR(rprof) = rval
	case RPYCUR:
	    AP_RPYCUR(rprof) = rval
	case ORPXCUR:
	    AP_ORPXCUR(rprof) = rval
	case ORPYCUR:
	    AP_ORPYCUR(rprof) = rval
	case RPRADIUS:
	    AP_RPRADIUS(rprof) = rval
	    AP_RPNPTS(rprof) = int (AP_RPRADIUS(rprof) / AP_RPSTEP(rprof)) + 1
	    call realloc (AP_RPDIST(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_INTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_DINTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_TINTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	case RPSTEP:
	    AP_RPSTEP(rprof) = rval
	    AP_RPNPTS(rprof) = int (AP_RPRADIUS(rprof) / AP_RPSTEP(rprof)) + 1
	    call realloc (AP_RPDIST(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_INTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_DINTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	    call realloc (AP_TINTENSITY(rprof), AP_RPNPTS(rprof), TY_REAL)
	case RPKSIGMA:
	    AP_RPKSIGMA(rprof) = rval
	case PYZMAG:
	    AP_PYZMAG(ply) = rval
	case PYMAG:
	    AP_PYMAG(ply) = rval
	case PYMAGERR:
	    AP_PYMAGERR(ply) = rval
	case PYX:
	    AP_PYX(ply) = rval
	case PYY:
	    AP_PYY(ply) = rval
	case PYMINRAD:
	    AP_PYMINRAD(ply) = rval
	case PYCX:
	    AP_PYCX(ply) = rval
	case PYCY:
	    AP_PYCY(ply) = rval
	case OPYCX:
	    AP_OPYCX(ply) = rval
	case OPYCY:
	    AP_OPYCY(ply) = rval
	case PYXMEAN:
	    AP_PYXMEAN(ply) = rval
	case PYYMEAN:
	    AP_PYYMEAN(ply) = rval
	case OPYXMEAN:
	    AP_OPYXMEAN(ply) = rval
	case OPYYMEAN:
	    AP_OPYYMEAN(ply) = rval
	case SKYSIGMA:
	    AP_SKYSIGMA(nse) = rval
	case EPADU:
	    AP_EPADU(nse) = rval
	case READNOISE:
	    AP_READNOISE(nse) = rval
	case THRESHOLD:
	    AP_THRESHOLD(fnd) = rval
	case RATIO:
	    AP_RATIO(fnd) = rval
	case THETA:
	    AP_THETA(fnd) = rval
	case NSIGMA:
	    AP_NSIGMA(fnd) = rval
	case SHARPLO:
	    AP_SHARPLO(fnd) = rval
	case SHARPHI:
	    AP_SHARPHI(fnd) = rval
	case ROUNDLO:
	    AP_ROUNDLO(fnd) = rval
	case ROUNDHI:
	    AP_ROUNDHI(fnd) = rval
	default:
	    call error (0, "APSETR: Unknown apphot real parameter")
	}
end


# AP2SETD -- Procedure to set a double apphot parameter.

procedure ap2setd (ap, param, dval)

pointer	ap		# pointer to apphot structure
int	param		# parameter
double	dval		# double value

pointer	nse, ply, rprof, fnd

begin
	nse = AP_NOISE(ap)
	ply = AP_POLY(ap)
	rprof = AP_RPROF(ap)
	fnd = AP_PFIND(ap)

	switch (param) {
	case PYNPIX:
	    AP_PYNPIX(ply) = dval
	case PYFLUX:
	    AP_PYFLUX(ply) = dval
	default:
	    call error (0, "APSETD: Unknown apphot double parameter")
	}
end