#!/bin/bash

# Driver script for parametric pin size analysis of the 
# 52-element Apollonian fuel bundle. AnnCon 2026 paper.
# Connor Moore, 2026. <connor.moore@ontariotechu.net>

green="\033[32m"
reset="\033[0m"

mkdir -p 52-results/
touch 52-results/52_keff_table

for fmr in $(seq 0.1 0.05 0.95); do
    echo -e ${green}Running 52-element Apollonian bundle with fmr = $fmr... ${reset}
    filename=52-results/$fmr.log
    python3 apollonian-52.py -fmr $fmr > $filename
    keff=$(cat $filename | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
    pm=$(cat $filename | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
    flow_area=$(cat $filename | grep "Flow area" | cut -d "=" -f 2)
    clad_circ=$(cat $filename | grep "Cladding circumference" | cut -d = -f 2)
    fuel_mass=$(cat $filename | grep "Fuel mass" | cut -d "=" -f 2)
    clad_mass=$(cat $filename | grep "Cladding mass" | cut -d "=" -f 2)
    dtu=$(cat $filename | grep "DTU ratio" | cut -d "=" -f 2)

    echo $keff +/- $pm, $flow_area cm² flow, $fuel_mass g/cm fuel, $clad_mass g/cm clad, $dtu DTU ratio | cowsay
    echo $keff $pm $flow_area $clad_circ $fuel_mass $clad_mass $dtu >> 52-results/keff_table
done