| 1 | defscript scan |
|---|
| 2 | |
|---|
| 3 | set numcycle $1 |
|---|
| 4 | set namecycle $2 |
|---|
| 5 | set firstf $3 |
|---|
| 6 | set lastf $4 |
|---|
| 7 | set f0 $5 |
|---|
| 8 | set bw $6 |
|---|
| 9 | set spectype $7 |
|---|
| 10 | |
|---|
| 11 | set dirname "./On" |
|---|
| 12 | |
|---|
| 13 | size = $lastf - $firstf + 1 |
|---|
| 14 | |
|---|
| 15 | newvec dataCh0 $size |
|---|
| 16 | newvec dataCh1 $size |
|---|
| 17 | |
|---|
| 18 | set index 0 |
|---|
| 19 | for f $firstf:$lastf |
|---|
| 20 | del spec${numcycle}f${f} |
|---|
| 21 | openfits $dirname/${namecycle}${numcycle}/${spectype}${f}.fits |
|---|
| 22 | rename ${spectype}${f} spec${numcycle}f${f} |
|---|
| 23 | del spec${numcycle}f${f}ch0 spec${numcycle}f${f}ch1 |
|---|
| 24 | objaoper spec${numcycle}f${f} row 0 spec${numcycle}f${f}ch0 |
|---|
| 25 | objaoper spec${numcycle}f${f} row 1 spec${numcycle}f${f}ch1 |
|---|
| 26 | del inarr |
|---|
| 27 | cp spec${numcycle}f${f}ch0 inarr |
|---|
| 28 | call func $f0 $bw 0 |
|---|
| 29 | x = ${outarr.sum}/${outarr.size} |
|---|
| 30 | ####### x = ${medianFilterFreq.sum}/${medianFilterFreq.size} |
|---|
| 31 | ####### echo " Cha 0: $x" |
|---|
| 32 | dataCh0[$index] = $x |
|---|
| 33 | ####### del medianFilterFreq |
|---|
| 34 | del outarr |
|---|
| 35 | |
|---|
| 36 | del inarr |
|---|
| 37 | cp spec${numcycle}f${f}ch1 inarr |
|---|
| 38 | call func $f0 $bw 1 |
|---|
| 39 | x = ${outarr.sum}/${outarr.size} |
|---|
| 40 | ####### x = ${medianFilterFreq.sum}/${medianFilterFreq.size} |
|---|
| 41 | ####### echo " Cha 1: $x" |
|---|
| 42 | dataCh1[$index] = $x |
|---|
| 43 | ####### del medianFilterFreq |
|---|
| 44 | del outarr |
|---|
| 45 | |
|---|
| 46 | index = $index + 1 |
|---|
| 47 | #clean |
|---|
| 48 | ####### disp spec${numcycle}f${f}ch0 "blue" |
|---|
| 49 | ####### disp spec${numcycle}f${f}ch1 "red same" |
|---|
| 50 | |
|---|
| 51 | # del spec${numcycle}f${f} |
|---|
| 52 | # del spec${numcycle}f${f}ch0 |
|---|
| 53 | # del spec${numcycle}f${f}ch1 |
|---|
| 54 | end |
|---|
| 55 | line2vec dataCh0$numcycle $dataCh0 |
|---|
| 56 | line2vec dataCh1$numcycle $dataCh1 |
|---|
| 57 | ##### arrtoascii dataCh0$numcycle dataCh0${numcycle}.dat |
|---|
| 58 | ##### arrtoascii dataCh1$numcycle dataCh1${numcycle}.dat |
|---|
| 59 | endscript |
|---|
| 60 | |
|---|
| 61 | c++compile fcalibextrac |
|---|
| 62 | c++link fcalibextrac.so func |
|---|
