#!/bin/bash

mkdir -p results/outputs
touch results/keff_table
echo "#" $(date) > results/keff_table


if [[ $1 == "lclrs" ]]; then
    prog=(openmc-inst01.py krusty.py -l -r)
    pool_seq=$(seq 1 45)
elif [[ $1 == "merlin" ]]; then
    prog=(../krusty.py)
    pool_seq=$2
else
    prog=(./krusty.py -r)
    pool_seq=$(seq 1 45)
fi

green="\e[0;32m"
reset="\e[0m"

for pool_height in $pool_seq; do
    echo -e "$green Starting runs for $pool_height cm pool height... $reset"
    for sat_height in $(seq $pool_height 45); do
        output_file=results/outputs/p$pool_height\_s$sat_height.out
        "${prog[@]}" -p $pool_height -s $sat_height > $output_file
        keff=$(cat $output_file | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
        pm=$(cat $output_file | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
        md5=$(cat $output_file | grep "MD5:" | cut -d ":" -f 2 | sed "s/ //g")
        echo K-eff $keff +/- $pm found for p=$pool_height and s=$sat_height \($md5\)
        echo $pool_height $sat_height $keff $pm $(hostname) $md5 \($(date)\) >> results/keff_table
    done
done