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
|
#!/usr/bin/env tcsh -f
#******************************************************************************
#* Johns Hopkins University
#* Center For Astrophysical Sciences
#* FUSE
#******************************************************************************
#*
#* Synopsis: cf_make_900_obs.csh association_file
#*
#* Description: This routine is derived from cf_make_all_obs.csh, but modified
#* to consider only 900+ (airglow) files. Output is a single
#* quick-look image of the airglow spactrum.
#*
#* Extracted spectra are combined, not IDF files.
#*
#* When an expected data set is missing, the script stops, cleans
#* the directory, and returns 1.
#*
#* History: 08/08/08 1.00 wvd Create separate quick-look image
#* for 900+ (airglow) exposures.
#* Always combine extracted spectra. If
#* cf_xcorr fails, combine with no shift.
#* 08/08/08 1.01 wvd Rename airglow quick-look file to
#* M112580100000airgttagf.gif
#*
#*****************************************************************************/
# Delete files after processing? (Default is no.)
#set DELETE_IDF # Delete intermediate data files
set DELETE_BPM # Delete bad-pixel map files
# Set program path
set rm = "rm -f"
set cf_xcorr = cf_xcorr
set cf_combine = cf_combine
set cf_pack = cf_pack
set cf_nvo = cf_nvo
set idl_obsplot = idl_obsplot.pl
set modhead = modhead
# Init var list
set detector = (1a 2b)
set channel = (lif sic)
set resolution = (2 3 4)
set obsmod = (hist ttag)
# Determine the root name and the program ID
set asnf = $1
set rn = ${asnf:s/000asnf.fit//}
set pid = `echo $rn | awk '{print substr($1, 1, 4)}'`
# Determine the object class
set tmp_file = `ls ${rn}*fcal.fit | awk '{if (NR == 1) print}'`
set tmp_buf = `$modhead $tmp_file OBJCLASS`
set objclass = $tmp_buf[2]
# Clean tmp files that the script will create (safe)
$rm tmp_xcorr.res tmp_combine.lis
$rm tmp_good_exp.lis tmp_exp.lis tmp_seg_dn.lis
$rm DN_${rn}*.fit
foreach om ($obsmod)
foreach res ($resolution)
foreach chan ($channel)
foreach det ($detector)
set ignore_exp_stat = ''
if ($om == hist) set ignore_exp_stat = -a
if ($pid == S100) set ignore_exp_stat = -a
if ($pid == M106) set ignore_exp_stat = -a
# Find exposures that match the current segment
set seg = $det$chan$res$om
set readfiles = 0
ls ${rn}9[0-9][0-9]${seg}fcal.fit |& grep -v 000${seg} > tmp_exp.lis # Keep only "9xx" exposures
if ($? == 0) then
set readfiles = 2
set ignore_exp_stat = -a
endif
if ($readfiles >= 1) then # There are one or more exposures
echo " "
echo "*** Processing: $seg ***"
# [1a][lif] -> [1b][lif], [2b][lif] -> [2a][lif] etc...
if ($det == 1a) set det2 = 1b
if ($det == 2b) set det2 = 2a
set seg2 = $det2$chan$res$om
echo "----- cf_xcorr input -----"
cat tmp_exp.lis
echo "----------------------------"
$cf_xcorr tmp_exp.lis tmp_xcorr.res # Compute shift and sigma_shift
echo "----- cf_xcorr results -----"
cat tmp_xcorr.res
echo "----------------------------"
awk '{if ($3 >= 0) {print $6, $2} else {print $6}}' tmp_xcorr.res > tmp_good_exp.lis
#
# --- Path 1: Optimize resolution ---
#
echo "Optimize resolution..."
set n_good = `cat tmp_good_exp.lis | wc -l`
# Extract [Day + Night] spectra
set s = DN_${rn}900${seg}fcal.fit
awk '{print "'$rn'"$1"'$seg'fcal.fit",$2}' tmp_good_exp.lis > tmp_combine.lis # Combine $seg [dn]
# echo "----- Combining Files ----- "
# cat tmp_combine.lis
# echo "Output: " $s
$cf_combine -k $ignore_exp_stat tmp_combine.lis $s
# $modhead "${s}[1]" NUM_EXP $n_good
$modhead "${s}[1]" COMBMETH XCORR
echo $s >> tmp_seg_dn.lis
set s = DN_${rn}900${seg2}fcal.fit
awk '{print "'$rn'"$1"'$seg2'fcal.fit",$2}' tmp_good_exp.lis > tmp_combine.lis # Combine $seg2 [dn]
# echo "----- Combining Files ----- "
# cat tmp_combine.lis
# echo "Output: " $s
$cf_combine -k $ignore_exp_stat tmp_combine.lis $s
# $modhead "${s}[1]" NUM_EXP $n_good
$modhead "${s}[1]" COMBMETH XCORR
echo $s >> tmp_seg_dn.lis
endif
end
$rm tmp_xcorr.res tmp_combine.lis
$rm tmp_good_exp.lis tmp_exp.lis
end
# Pack the 8 [detector][channel] pairs together
if (-e tmp_seg_dn.lis) then
set fcal_all = ${rn}00900all$res${om}fcal.fit # Final output name
$rm $fcal_all # Clean (safe)
set n_segs = `cat tmp_seg_dn.lis | wc -l`
if (!($n_segs == 8)) then
@ mseg = 8 - $n_segs
echo "ERROR: $mseg (day + night) segments are missing"
$rm tmp_seg_dn.lis
goto crash
else
$cf_pack tmp_seg_dn.lis $fcal_all
$rm tmp_seg_dn.lis
# Plot figures, delete unwanted files
$idl_obsplot {$rn} airglow
mv ${rn}00900spec${om}f.gif ${rn}00000airg${om}f.gif
$rm ${rn}00900lif*.gif ${rn}00900sic*.gif
$rm $fcal_all
endif
endif
end
end
# Clean [dn] files.
$rm DN_${rn}*.fit
# Delete IDF files
if $?DELETE_IDF then
echo "NOTE: Deleting intermediate data files."
$rm ${rn}*idf.fit
endif
# Delete bad-pixel-map (bpm) files
if $?DELETE_BPM then
echo "NOTE: Deleting bad pixel map (bpm) files."
$rm ${rn}*bpm.fit
endif
exit(0)
crash: # Procedure when script crashes
# Clean directory
$rm tmp_xcorr.res tmp_combine.lis
$rm tmp_good_exp.lis tmp_exp.lis tmp_seg_dn.lis
$rm DN_*.fit
# Return 1
exit(1)
|
