Updated driver file and added gnuplot script

1 files changed, 46 insertions, 38 deletions

diff --git a/triso_study_driver.sh b/triso_study_driver.sh
index 7b7e7a1..e7b72e7 100755
--- a/triso_study_driver.sh
+++ b/triso_study_driver.sh
@@ -1,44 +1,52 @@
 #!/bin/bash
 
-prog=(./pin.py -e -f un)
-
-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 > 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 > 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.30 0.05 0.75); do
-    echo -n "Running RPT with $rad cm... "
-    "${prog[@]}" -t rpt -r $rad > 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")
+fuels=("uo2" "uco10" "uco20" "un")
+
+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
-    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 > 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")
+    # 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
-    echo $rad $keff $pm $(hostname) $md5 \($(date)\) >> results/rrpt_table
+
+    # 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