bond_analyze_script

   1 #!/bin/bash
   2
   3 if [ ! -d $1 ]; then
   4         echo "$1 is not a directory"
   5         exit 1
   6 fi
   7
   8 echo "processing $1 ..."
   9
  10 rm -f $1/bond_analyze.txt
  11
  12 for i in $1/*.save; do
  13         bai="`./bond_analyze $i | grep '\[bond' | grep quality`"
  14         time=`echo $bai | awk '{ print $4 }'`
  15         quality4=`echo $bai | awk '{ print $8 }'`
  16         quality3=`echo $bai | awk '{ print $10 }'`
  17         count=`echo $bai | awk '{ print $12 }'`
  18         echo "$time $quality4 $quality3 $count" >> $1/bond_analyze.txt
  19 done
  20
  21 echo "done"