Initial push with working study

1 files changed, 44 insertions, 0 deletions

diff --git a/triso_study_driver.sh b/triso_study_driver.sh
new file mode 100755
index 0000000..7b7e7a1
--- /dev/null
+++ b/triso_study_driver.sh
@@ -0,0 +1,44 @@
+#!/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")
+    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")
+    echo found $keff +/- $pm
+    echo $rad $keff $pm $(hostname) $md5 \($(date)\) >> results/rrpt_table
+done