#!/bin/bash

#PRMTOP="../1SVT_chainA.prmtop"
PRMTOP="../1XCK_chainA_MGATP_sol.prmtop"
PREFIX="126_1XCK_chainA_MGATP_prod"
SUFFIX_ORIG=".mdcrd"
SUFFIX_REIMG="_reimg.nc"

C=250
END=310
#OFFSET=' 1 500 1 '
#OFFSET=' 1 10000 1 '
#type="radgyr"
#type="distance"
#type=project
#type="snapshots"
#type="rmsd"
type="short"
#type="cluster"
#type="mmtsb_clustering"
#type=qha
#type=correl
#type=proj
ptraj_file=${type}.ptraj

if [ -f $ptraj_file ]; then
	rm $ptraj_file
fi

if [ ! $WHOLE == "" ]; then
	echo 'trajin '${WHOLE} $OFFSET >> $ptraj_file
else 
	let END=${END}+1
	while [ $C -lt $END ]; do
		traj_orig=${PREFIX}${C}${SUFFIX_ORIG}
		traj_reimg=${PREFIX}${C}${SUFFIX_REIMG}

                if [ -f $traj ]; then
                        echo 'trajin ../mdcrd/'${traj_orig} $OFFSET >> $ptraj_file
                fi
	
		let C=${C}+1
	done
fi

echo `cat $type >> $ptraj_file`

ptraj $PRMTOP $ptraj_file