ctest UNRES_M_ss_static_min
[unres.git] / 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