• Print
• .

DoProcessCalNew

 #!/bin/csh
 # converts dtacq file into wera raw file then calls wera16_raw_sort and launch standard SORT processing
 # new version to be used with new version of dtacq2wera.m 2015.03.17
 # do_process_DTA.sh  20123661630_acq164_091.dta

 # stupid GKS hack
 unsetenv DISPLAY

 date

 set sitename='kkh'

 # 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 timefile=`echo $infile | sed -e s/.dta/.hdr/`
 set filename="$timetag"_$sitename
 set outfile=$filename.mat
 set rawfile=$filename.CAL

 cd /home/wera/data/$sitename

 (echo filein="'$infile'"; echo fileout="'$rawfile'"; echo filemat="'$outfile'"; echo filetime="'$timefile'"; \
 echo filehead="'/home/wera/wera_import/$sitename.cal.hdr'"; \
 cat /home/wera/wera_import/dtacq2wera.m) | /usr/local/bin/matlab -nodisplay
 #/bin/gzip $infile &
 /bin/rm $infile
 echo $rawfile made...
 echo $outfile made...

 (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
 gifsicle --delay 100 $filename.CAL_ant{01,02,03,04,05,06,07,08,09,10,11,12,13,14,15,16}.gif --loop  >  $filename.CAL_ant.gif
 /bin/rm $filename.CAL_ant??.gif
 echo $filename.CAL_ant.gif made...