• Print
• .

DoProcessCalOld

  #!/bin/csh
  # converts dtacq file into wera cal file then calls plotlsq and launch standard CAL processing
  # do_process_CAL.sh  20123661630_acq164_091.dta

  # stupid GKS hack
  unsetenv DISPLAY

  date

  set sitename='kok'

  # Make the output filenames

  set infile=$1
  set timetag=`echo $infile | sed -e s/_acq164_....dta//`
  set dtaname=`echo $infile | sed -e s/201........_acq/acq/ -e s/.dta//`
  set hdrfile=`echo $infile | sed -e s/.dta/.hdr/`
  set filename="$timetag"_$sitename
  set outfile=$filename.DAT
  set calfile=$filename.CAL

  cd /home/wera/data/$sitename

  (echo filein="'$infile'"; echo fileout="'$outfile'";cat /home/wera/wera_import/dtacq2wera.m) | /usr/local/bin/matlab -nodisplay
  #/bin/gzip $infile &
  #/bin/rm $infile

  # prepend the standard header
  (echo -n '2048 SAMPLES   ';cat $hdrfile /home/wera/wera_import/$sitename.cal.hdr $outfile) > $calfile
  /bin/rm $outfile
  echo $calfile made...

  exit
  # make individual channels spectrum

  (echo g; echo 1; echo y; echo g; echo 2; echo y; echo g; echo 3; echo y;   echo g; echo 4; echo y; echo g; echo 5; echo y; echo g; echo 6; echo y;   echo g; echo 7; echo y; echo g; echo 8; echo y; echo g; echo 9; echo y;   echo g; echo 10; echo y; echo g; echo 11; echo y; echo g; echo 12; echo y;   echo g; echo 13; echo y; echo g; echo 14; echo y; echo g; echo 15; echo y;   echo g; echo 16; echo n) | plotlsq $calfile
  gifsicle --delay 100 $filename.CAL_ant0{1,2,3,4,5,6,7,8}.gif --loop  >  $filename.CAL_ant.gif
  /bin/rm $filename.CAL_ant??.gif
  echo $filename.CAL_ant.gif made...

  # cleanup