source: BAORadio/AmasArchive/Amas 21cm/transfert/transfert/Abell85/calib.pic @ 672

defscript getspec
set numcycle $1
set mode     $2
set namedir  $3
set f0 $4
set bw $5
set typefile $6
set source $7
set t0date $8
set t0ttag $9

echo ' -------- processing cycle ' $numcycle


newnt datav0 x0 time
newnt datav1 x1 time

shell find ${source}/${mode}/${namedir}${numcycle}/ -name "${typefile}*" -print | xargs -i basename {} .fits > tmpfile.txt

forinfile fname tmpfile.txt
  del spec
  set var "<${fname}>"
  if ( ${var} == "<>" ) then
     break
  endif
  echo "open file <${fname}>"
  openfits ${source}/${mode}/${namedir}${numcycle}/${fname}.fits
  rename ${fname} spec

#on tag la fenetre par la date du 1er fichier auquel on ajoute l'increment du mean time tag.
#le tout en millisecond. (JEC)
# time tag = t0(date) + mean-tt - mean-tt0

  timetag = ${t0date}
  curMeanTT = 1000. * ${spec.info.MEANTT}
  timetag = ${timetag} + ${curMeanTT}
  timetag = ${timetag} - ${t0ttag}
  timetag = ${timetag}/1000.


  del specv0
  objaoper  spec  row 0 specv0
  del inarr
  cp specv0 inarr
  call func $f0 $bw
  val = ${outarr.sum}/${outarr.size}
  del outarr
#  echo "cycle $numcycle,  file $i x0= $val , timetag= $timetag"
  line2nt  datav0 $val $timetag


  del specv1
  objaoper  spec  row 1 specv1
  del inarr
  cp specv1 inarr
  call func $f0 $bw
  val = ${outarr.sum}/${outarr.size}
  del outarr
  line2nt  datav1 $val $timetag


#clear
  del spec
  del specv0
  del specv1

end
endscript

#c++compile extract
#c++link extract.so func