-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathmakePointingCont.sh
52 lines (46 loc) · 2.3 KB
/
makePointingCont.sh
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
#!/bin/bash
set -eu
file_freq_set='catalog/B4R_frequency_settings.txt'
file_Pointing_cont='catalog/B4R_PointingCont_list.txt'
original_scr='templates/b4r_PointingCont_3C84_cont02_v01.scr'
original_source_name=`echo ${original_scr} | awk 'BEGIN { FS="_"; OFS="_" } { print $3 }'`
original_freq_name=`echo ${original_scr} | awk 'BEGIN { FS="_"; OFS="_" } { print $4 }'`
datetime=`date -u`
make_dir=products/`basename $original_scr .scr`_all
echo Make all source and freq setting for ${original_scr} in ${make_dir}.
mkdir -p ${make_dir}
cat $file_Pointing_cont | grep -v '^#.*' | while read source ra dec comments; do
if [ $ra != "Planet" ]; then
planet="false"
else
echo $dec
ra="00:00:00.0"
dec="00:00:00.0"
planet="true"
fi
cat $file_freq_set | grep -v '^#.*' | while read freqID line_freq IF_freq comments; do
freqID2=`printf "%02d" $freqID`
printf "making for %s: " $source
printf "freqID %s: " $freqID2
printf "line %.1f GHz: " $line_freq
printf "IF %.1f GHz\n" $IF_freq
cont_or_spec=`echo ${original_scr} | awk 'BEGIN { FS="_"; OFS="_" } { print $4 }' | cut -c 1-4`
cont_or_spec_freqID2=$cont_or_spec$freqID2
tmp_scr=`echo ${original_scr} | awk 'BEGIN { FS="_"; OFS="_" } { $4="'${cont_or_spec_freqID2}'"; $3="'${source}'"; print $0 }'`
tmp_scr=`basename $tmp_scr`
new_scr=./${make_dir}/${tmp_scr}
cat ${original_scr} | awk '
{ for(i=1;i<=NF;i++) { if($i=="-If2Freq[0]") {$(i+1)="'${IF_freq}'"} } }
{ for(i=1;i<=NF;i++) { if($i=="-If2Freq[1]") {$(i+1)="'${IF_freq}'"} } }
{ for(i=1;i<=NF;i++) { if($i=="-LineFreq[0]") {$(i+1)="'${line_freq}'"} } }
{ gsub("'${original_source_name}'", "'${source}'") }
{ gsub("'${original_freq_name}'", "'${cont_or_spec_freqID2}'") }
{ for(i=1;i<=NF;i++) { if($i=="-Dec[0]") {$(i+1)="'${dec}'"} } }
{ for(i=1;i<=NF;i++) { if($i=="-Dec[1]") {$(i+1)="'${dec}'"} } }
{ for(i=1;i<=NF;i++) { if($i=="-Ra[0]") {$(i+1)="'${ra}'"} } }
{ for(i=1;i<=NF;i++) { if($i=="-Ra[1]") {$(i+1)="'${ra}'"} } }
{ if("'${planet}'"=="true") for(i=1;i<=NF;i++) { if($i=="-Planet") {$(i+1)="'${source}'"} } }
{ print $0 }
' > $new_scr
done
done