#!/bin/sh

# This is for running on EXPANSE

source $HOME/.bashrc

module purge
module load cpu/0.15.4
module load gcc/9.2.0
module load mafft/7.453
module load slurm

MAFFTBINARY=mafft
/expanse/projects/ngbt/opt/expanse/slurm-aff $MAFFTBINARY &

case $1 in
                -qinsi)
                        shift; MAFFTBINARY=mafft-qinsi;;
                        #check that seed is a positive int
                -xinsi)
                        shift; MAFFTBINARY=mafft-xinsi;;
                        #check that num_chains is a positive int.
                *)
                        ;;
esac

echo "CIPRES_THREADSPP=$CIPRES_THREADSPP" 1>&2
echo "CIPRES_NP=$CIPRES_NP" 1>&2
echo "running:" 1>&2
echo "${MAFFTBINARY} --progress ${PWD}/progress.mafft $*" 1>&2

${MAFFTBINARY} --progress ${PWD}/progress.mafft $*

treeoutfile=*tree
if [ -f $treeoutfile ]
then chmod g+r *tree
fi