ctest/csa_energy_check.sh

   1 #!/bin/bash
   2
   3 if [ "$1" == "1l2y_csa_4P" ]; then
   4  file=1l2y_csa_4P.CSA.history
   5 elif [ "$1" == "1l2y_csa_local_4P" ]; then
   6  file=1l2y_csa_local_4P.CSA.history
   7 elif [ "$1" == "1l2y_csa" ]; then
   8  file=1l2y_csa.CSA.history
   9 elif [ "$1" == "1l2y_csa_local" ]; then
  10  file=1l2y_csa_local.CSA.history
  11 else
  12  exit 1
  13 fi
  14
  15 # Check if file exist
  16 if [ ! -f $file ]; then
  17     echo "CRITICAL: out $file do not exist"
  18     exit 2
  19 fi
  20
  21 emin=`grep jlee $file |awk '{print $10}'`
  22 echo 'Emin= ' $emin
  23
  24 if [ "$1" == "1l2y_csa_4P" ]; then
  25  echo 'difference ' `echo "a=$emin+215.41208;if(0>a)a*=-1;a"|bc -l` ' from reference emin -215.41208'
  26  if [ `echo "a=$emin+215.41208;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
  27   echo 'ERROR difference from reference emin -215.41208 greater than 5.0'
  28   exit 1
  29  fi
  30 elif [ "$1" == "1l2y_csa_local_4P" ]; then
  31  echo 'difference ' `echo "a=$emin+200.13486;if(0>a)a*=-1;a"|bc -l` ' from reference emin -200.13486'
  32  grep rms 1l2y_csa_local_4P*.pdb
  33  if [ `echo "a=$emin+200.13486;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
  34   echo 'ERROR difference from reference emin 200.13486 greater than 5.0'
  35   exit 1
  36  fi
  37 elif [ "$1" == "1l2y_csa" ]; then
  38  echo 'difference ' `echo "a=$emin-29.92824;if(0>a)a*=-1;a"|bc -l` ' from reference emin 29.92824'
  39  if [ `echo "a=$emin-29.92824;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
  40   echo 'ERROR difference from reference emin 29.92824 greater than 5.0'
  41   exit 1
  42  fi
  43 elif [ "$1" == "1l2y_csa_local" ]; then
  44  echo 'difference ' `echo "a=$emin-31.52677;if(0>a)a*=-1;a"|bc -l` ' from reference emin 31.52677'
  45  grep rms 1l2y_csa_local_4P*.pdb
  46  if [ `echo "a=$emin-31.52677;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
  47   echo 'ERROR difference from reference emin 31.52677 greater than 5.0'
  48   exit 1
  49  fi
  50
  51 fi