ctest/wham_check.sh

   1 #!/bin/bash
   2
   3 if [ "$1" == "1L2Y_wham" ]; then
   4  file=1L2Y_wham.thermal
   5 else
   6  exit 1
   7 fi
   8
   9 # Check if file exist
  10 if [ ! -f $file ]; then
  11     echo "CRITICAL: out $file do not exist"
  12     exit 2
  13 fi
  14
  15 max=`awk '{print $1,$7*1}' $file | sort -n -k 2 | awk 'END{print $1}'`
  16 echo 'T of max Cv(T) ' $max
  17 rms=`awk '{if ($1<260) {a=a+$5;n++}}END{print a/n}' $file`
  18 echo 'average rms for T<260 ' $rms
  19
  20  if [ `echo "a=$max-316;if(0>a)a*=-1;a>10.0"|bc -l` != "0" ]; then
  21   echo 'ERROR wrong T max of Cv(T) by more than 10 K'
  22   exit 1
  23  fi
  24
  25  if [ `echo "a=$rms-3.0;if(0>a)a*=-1;a>0.2"|bc -l` != "0" ]; then
  26   echo 'ERROR wrong average rms for T<260 by more than 0.2'
  27   exit 1
  28  fi