#!/bin/bash

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

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

mkdir -p 37-results/
touch 37-results/37_keff_table

for r in $(seq 0.1 0.05 0.8); do
    echo -e ${green}Running rf = $r cm... ${reset}
    filename=37-results/$r.log
    python3 candu-37.py -rf $r > $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)
    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
done