#!/bin/csh -f

#Tests chamber to ensure it gives the correct output and creates the expected files

#1) Test Chamber
/bin/rm -f *.prmtop *.inpcrd *.out *.restrt

# Test for the presence of the 'chamber' program

if ( -x $AMBERHOME/bin/chamber ) then
   #CHARMM CRD FILE
   $AMBERHOME/bin/chamber \
             -top ./top_cbh1_cellulose.rtf \
             -param ./par_cbh1_cellulose.prm \
             -psf gluc.psf -nocmap\
             -crd glucmin.pdb \
             -p prmtop \
             -inpcrd inpcrd >& chamber.out || goto error

   ../../dacdif chamber.out.save chamber.out
   ../../dacdif -r 10E-15 prmtop.save prmtop
   ../../dacdif inpcrd.save inpcrd

   exit(0)

else
   echo "Could not find the 'chamber' program. "
   echo "It may happen if you do not compile the serial"
   echo "version of amber."
   echo "=============================================================="
   exit(0)
endif   
error:
echo "  ${0}:  Program error"
exit (1)