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
|
include "sensfunc.h"
# SF_WEIGHTS -- Change weights for point, star, or wavelength.
# The original input weight is permanently lost.
procedure sf_weights (stds, nstds, key, istd, ipt, weight)
pointer stds[nstds] # Standard star data
int nstds # Number of standard stars
int key # Delete point, star, or wavelength
int istd # Index of standard star
int ipt # Index of point
real weight # New weight
int i, j, n
real wave
pointer z, w, iw
begin
switch (key) {
case 'p':
Memr[STD_WTS(stds[istd])+ipt-1] = weight
Memr[STD_IWTS(stds[istd])+ipt-1] = weight
case 's':
n = STD_NWAVES(stds[istd])
w = STD_WTS(stds[istd])
iw = STD_IWTS(stds[istd])
call amovkr (weight, Memr[w], n)
call amovkr (weight, Memr[iw], n)
case 'w':
wave = Memr[STD_WAVES(stds[istd])+ipt-1]
do i = 1, nstds {
if (STD_FLAG(stds[i]) != SF_INCLUDE)
next
n = STD_NWAVES(stds[i])
z = STD_WAVES(stds[i])
w = STD_WTS(stds[i])
iw = STD_IWTS(stds[i])
do j = 1, n {
if (Memr[z] == wave) {
Memr[w] = weight
Memr[iw] = weight
}
w = w + 1
iw = iw + 1
z = z + 1
}
}
}
end
|
