#! /bin/sh

# Accessible surface area and volume calculation.

set -e

surface xyzin1 $CEXAM/toxd/toxd.pdb xyzout $CCP4_SCR/toxd.surface <<+ 
NFILES 1
ZSTEP 0.25
RUN
ALLATM
NEXT
STOP
+

# Nb XtalView also has a program called volume
# Specifically use the version in the CCP4 bin directory
${CBIN}/volume xyzin $CCP4_SCR/toxd.surface xyzout $CCP4_SCR/calc_shell.vol \
       shellfile $CCP4_SCR/toxd.shell <<+ 
SHELL
NOCALC
RUN
+

${CBIN}/volume xyzin $CCP4_SCR/toxd.surface xyzout $CCP4_SCR/with_shell.vol \
       shellfile $CCP4_SCR/toxd.shell <<+ 
NOSHELL
CALC
WATER
RUN
+

# find awk executable
for myawk in nawk gawk awk
do
  if ( `type $myawk >/dev/null 2>&1` ); then
    break
  fi
done

$myawk <$CCP4_SCR/with_shell.vol  '
BEGIN {first = 1;}
/^ BEGIN/ {
    started = 1;
    print "\nResidue, number, volume (>1000 probably indicates ignored atom)\n";
    next;
}
/^ TOTAL/ {exit;}
{
    if (!started) { next; }
    seqnum = substr ($0, 29, 5);
    volume = substr ($0, 68, 8);
    res3 = substr ($0, 23, 6);
    if (old_seqnum == seqnum) {
      volume_sum += volume ;
    } else if (first) {
      volume_sum = volume;
      old_seqnum = seqnum ;
      old_res3 = res3;
      first = 0;
    } else {
      printf "%6s%5d %8.3f\n", old_res3,old_seqnum,volume_sum;
      old_res3 = res3;
      volume_sum = volume;
      old_seqnum = seqnum ;
    }      
}
'
echo ' Normal Termination'
#