#!/bin/bash
#$ -cwd
#$ -V
#$ -P bsg.prjb -q long.qb
#$ -e ErrFiles
#$ -o OutFiles
#$ -N process_data.sh
#$ -t 1-32
#$ -j y

replicate=1000
case="Constantpopsize"

#joblist=("ms" \
"scrmwindow1000000" "scrmwindow700000" "scrmwindow500000" "scrmwindow300000" \
#"scrmwindow100000" "scrmwindow70000" "scrmwindow50000" "scrmwindow30000" \
#"scrmwindow10000" "scrmwindow7000" "scrmwindow5000" "scrmwindow3000" \
#"scrmwindow1000" "scrmwindow500" "scrmwindow0" \
#"macsretain1000000" "macsretain700000" "macsretain500000" "macsretain300000" \
#"macsretain100000" "macsretain70000" "macsretain50000" "macsretain30000" \
#"macsretain10000" "macsretain7000" "macsretain5000" "macsretain3000" \
#"macsretain1000" "macsretain500" "macsretain0" \
#"fastsimcoal")

joblist=("scrmwindow1000000" "scrmwindow700000" "scrmwindow500000" "scrmwindow300000" \
"scrmwindow100000" "scrmwindow70000" "scrmwindow50000" "scrmwindow30000" \
"scrmwindow10000" "scrmwindow7000" "scrmwindow5000" "scrmwindow3000" \
"scrmwindow1000" "scrmwindow500" "scrmwindow0" )

#joblist=("scrmwindow1000000" "scrmwindow700000" "scrmwindow500000" "scrmwindow300000" \
#"macsretain1000000" "macsretain700000" "macsretain500000" "macsretain300000" )

rep=$(expr $SGE_TASK_ID - 1)
#rep=0
#for rep in $(seq 0 1 16)
    #do
    Job=${case}${joblist[${rep}]}
    JobParamFile=${Job}param
echo -e "case: ${case}\n\
nsam: 20\n\
replicate: ${replicate}\n\
seqlen: 10000001\n\
rho: 4000\n\
job: ${Job}_" > ${JobParamFile}
    python process_data.py ${JobParamFile}
    #done

