#!/bin/csh -f
#TEST-PROGRAM sander
#TEST-DESCRIP TO_BE_DEtermined
#TEST-PURPOSE regression, basic
#TEST-STATE   undocumented

set sander = "../../bin/sander"
if( $?TESTsander ) then
   set sander = $TESTsander
endif

if( ! $?DO_PARALLEL ) then
        setenv DO_PARALLEL " "
else
  set numprocs=`echo $DO_PARALLEL | awk -f ../numprocs.awk `
  if ( $numprocs > 19 ) then
      echo " DO_PARALLEL set to $DO_PARALLEL"
      goto toomany
  else if ( $?MP_PROCS)then
      if ($MP_PROCS > 19)then
#               --- For IBM AIX ---
          echo "    MP_PROCS set to "$MP_PROCS
          goto toomany
      endif
  endif
endif

cat <<EOF > mdin
  test of rdc and pseudo-csa restraints for an RNA stem, allowing
  the molecule to move
 &cntrl                                                                        
    pencut=0.05, nmropt=0,                                     
    imin=1, maxcyc=2000, igb=1, rgbmax=15.,
    ntpr=50,
    ipnlty=1, cut=25.0, ntb=0, iscale=0, scalm=300.,                            
    ntc=1, ntb=0,
 /
 &wt type='END'  /                                                          
LISTOUT=POUT
DIPOLE=RST.dip3
CSA=RST.csa
END                                                                            
EOF

$DO_PARALLEL $sander -O -i mdin -o stem.csa2.o < /dev/null || goto error

../dacdif -s -w stem.csa2.o.save stem.csa2.o
#/bin/rm -f mdin restrt mdinfo
exit(0)

error:
echo "  ${0}:  Program error"
exit(1)

toomany:
    echo " too many processors for this test, exiting (Max = 19)"
    echo "============================================================"
    exit(0)