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
|
include <error.h>
include <smw.h>
# SMW_SWATTRS -- Set spectrum attribute parameters
# This routine has the feature that if the coordinate system of a single
# spectrum in an EQUISPEC WCS is changed then the image WCS is changed
# to a MULTISPEC WCS.
procedure smw_swattrs (smw, index1, index2, ap, beam, dtype, w1, dw, nw, z,
aplow, aphigh, coeff)
pointer smw # SMW pointer
int index1 # Spectrum index
int index2 # Spectrum index
int ap # Aperture number
int beam # Beam number
int dtype # Dispersion type
double w1 # Starting coordinate
double dw # Coordinate interval
int nw # Number of valid pixels
double z # Redshift factor
real aplow[2], aphigh[2] # Aperture limits
char coeff[ARB] # Nonlinear coeff string
bool fp_equald()
int i, j, sz_val, strlen()
double a, b
pointer sp, str, val, mw
errchk smw_mw
define start_ 10
begin
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
start_
switch (SMW_FORMAT(smw)) {
case SMW_ND:
if (!IS_INDEFI(SMW_DTYPE(smw)) && (!fp_equald(w1,SMW_W1(smw)) ||
!fp_equald(dw,SMW_DW(smw)) || !fp_equald(z,SMW_Z(smw)))) {
call malloc (val, 15, TY_DOUBLE)
mw = SMW_MW(smw,0)
i = SMW_PDIM(smw)
j = SMW_PAXIS(smw,1)
call mw_gwtermd (mw, Memd[val], Memd[val+i], Memd[val+2*i], i)
Memd[val+j-1] = 1.
switch (dtype) {
case DCNO, DCLINEAR:
a = w1 / (1 + z)
b = dw / (1 + z)
case DCLOG:
a = log10 (w1 / (1 + z))
b = log10 ((w1 + (nw - 1) * dw) / w1) / (nw - 1)
case DCFUNC:
call error (1,
"Nonlinear functions not allowed for NSPEC format")
}
Memd[val+i+j-1] = a
Memd[val+2*i+(i+1)*(j-1)] = b
call mw_swtermd (mw, Memd[val], Memd[val+i], Memd[val+2*i], i)
call mfree (val, TY_DOUBLE)
}
SMW_DTYPE(smw) = dtype
SMW_NW(smw) = nw
SMW_W1(smw) = w1
SMW_DW(smw) = dw
SMW_Z(smw) = z
case SMW_ES:
# Check for any changes to the dispersion system.
if (dtype == DCFUNC) {
call smw_esms(smw)
goto start_
}
if (!IS_INDEFI(SMW_DTYPE(smw)) && (dtype != SMW_DTYPE(smw) ||
nw != SMW_NW(smw) || !fp_equald(w1,SMW_W1(smw)) ||
!fp_equald(dw,SMW_DW(smw)) || !fp_equald(z,SMW_Z(smw)))) {
if (SMW_NSPEC(smw) > 1 && index1 > 0) {
call smw_esms(smw)
goto start_
}
if (!fp_equald(w1,SMW_W1(smw)) || !fp_equald(dw,SMW_DW(smw)) ||
!fp_equald(z,SMW_Z(smw))) {
call malloc (val, 15, TY_DOUBLE)
mw = SMW_MW(smw,0)
i = SMW_PDIM(smw)
j = SMW_PAXIS(smw,1)
call mw_gwtermd (mw, Memd[val], Memd[val+i],
Memd[val+2*i], i)
Memd[val+j-1] = 1.
switch (dtype) {
case DCNO, DCLINEAR:
a = w1 / (1 + z)
b = dw / (1 + z)
case DCLOG:
a = log10 (w1 / (1 + z))
b = log10 ((w1 + (nw - 1) * dw) / w1) / (nw - 1)
}
Memd[val+i+j-1] = a
Memd[val+2*i+(i+1)*(j-1)] = b
call mw_swtermd (mw, Memd[val], Memd[val+i],
Memd[val+2*i], i)
call mfree (val, TY_DOUBLE)
}
}
SMW_DTYPE(smw) = dtype
SMW_NW(smw) = nw
SMW_W1(smw) = w1
SMW_DW(smw) = dw
SMW_Z(smw) = z
if (index1 > 0) {
Memi[SMW_APS(smw)+index1-1] = ap
Memi[SMW_BEAMS(smw)+index1-1] = beam
Memr[SMW_APLOW(smw)+2*index1-2] = aplow[1]
Memr[SMW_APHIGH(smw)+2*index1-2] = aphigh[1]
Memr[SMW_APLOW(smw)+2*index1-1] = aplow[2]
Memr[SMW_APHIGH(smw)+2*index1-1] = aphigh[2]
}
case SMW_MS:
# We can't use SPRINTF for the whole string because it can only
# handle a limited length and trucates long coefficient strings.
# Use STRCAT instead.
call smw_mw (smw, index1, index2, mw, i, j)
sz_val = strlen (coeff) + SZ_LINE
call salloc (val, sz_val, TY_CHAR)
call sprintf (Memc[str], SZ_LINE, "spec%d")
call pargi (i)
call sprintf (Memc[val], sz_val,
"%d %d %d %.14g %.14g %d %.14g %.2f %.2f")
call pargi (ap)
call pargi (beam)
call pargi (dtype)
if (dtype == DCLOG) {
call pargd (log10 (w1))
call pargd (log10 ((w1+(nw-1)*dw)/w1)/(nw-1))
} else {
call pargd (w1)
call pargd (dw)
}
call pargi (nw)
call pargd (z)
call pargr (aplow[1])
call pargr (aphigh[1])
if (coeff[1] != EOS) {
call strcat (" ", Memc[val], sz_val)
call strcat (coeff, Memc[val], sz_val)
}
call mw_swattrs (mw, 2, Memc[str], Memc[val])
if (SMW_APS(smw) != NULL)
Memi[SMW_APS(smw)+index1-1] = ap
}
call sfree (sp)
end
|
