#!/bin/bash

#fuels=("uo2" "uco5" "uco10" "uco20" "un")
fuels=("uo2")

for mat in $fuels; do

    echo Beginning runs for $mat fuel

    prog=(./pin.py -e)

    mkdir -p results/ results/rpt results/rrpt
    touch results/rpt_table results/rrpt_table
    echo "\#" $(date) > results/rpt_table
    echo "\#" $(date) > results/rrpt_table

    # 1. Basic calculation without homogenization
    echo -n "Running explicit TRISO calculation... "
    "${prog[@]}" -t none -f $mat > results/explicit.out
    keff=$(cat results/explicit.out| grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
    pm=$(cat results/explicit.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
    echo found $keff +/- $pm

    # 2. VWH calculation
    echo -n "Running VWH homogenization calculation... "
    "${prog[@]}" -t vwh -f $mat > results/vwh.out
    keff=$(cat results/vwh.out| grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
    pm=$(cat results/vwh.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
    echo found $keff +/- $pm

    # 3. RPT
    for rad in $(seq 0.25 0.05 0.75); do
        echo -n "Running RPT with $rad cm... "
        "${prog[@]}" -t rpt -r $rad -f $mat > results/rpt/$rad.out
        keff=$(cat results/rpt/$rad.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
        pm=$(cat results/rpt/$rad.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
        md5=$(cat results/rpt/$rad.out | grep "MD5:" | cut -d ":" -f 2 | sed "s/ //g")
        echo found $keff +/- $pm
        echo $rad $keff $pm $(hostname) $md5 \($(date)\) >> results/rpt_table
    done

    # 4. RRPT
    for rad in $(seq 0.1 0.05 0.95); do
        echo -n "Running RRPT with $rad cm... "
        "${prog[@]}" -t rrpt -r $rad -f $mat > results/rrpt/$rad.out
        keff=$(cat results/rrpt/$rad.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
        pm=$(cat results/rrpt/$rad.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
        md5=$(cat results/rrpt/$rad.out | grep "MD5:" | cut -d ":" -f 2 | sed "s/ //g")
        echo found $keff +/- $pm
        echo $rad $keff $pm $(hostname) $md5 \($(date)\) >> results/rrpt_table
    done
    mv results/ results_$mat/
done