#!/bin/bash

INFILE="./input/params.in"
OUT="./output"
SIM="./simulation"

for rate in `seq 100 150 400`
do
for loc in `seq 60 10 60`
do
for sites in `seq 1 2 5`
do  
  sed -i 's/^GEN_ENZYME_RATE=.*$/GEN_ENZYME_RATE='$rate'/g' $INFILE
  sed -i 's/^GEN_LOC_X=.*$/GEN_LOC_X=50/g' $INFILE
  sed -i 's/^GEN_LOC_Y=.*$/GEN_LOC_Y='$loc'/g' $INFILE
  sed -i 's/^GEN_LOC_Z=.*$/GEN_LOC_Z='$loc'/g' $INFILE
  sed -i 's/^CUT_SITES_TERM=.*$/CUT_SITES_TERM='$sites'/g' $INFILE
  rep=102
  cores=6 #Cores must divide rep
  echo "Simulation: $rate-50-$loc-$loc-$sites for "$rep" repetitions"
  for i in `seq $((rep/cores))`
  do
    for j in `seq $cores`
    do
       if [ 1 -lt $i ] 
       then
         wait ${pids[$j]}
       fi
       nice $SIM >> $OUT/$rate-50-$loc-$loc-$sites &
       pids[$j]=$!
       echo -ne $(((i-1)*cores+j)) \\r
    done
  done
done
done
done  
