summaryrefslogtreecommitdiff
path: root/triso_study_driver.sh
diff options
context:
space:
mode:
Diffstat (limited to 'triso_study_driver.sh')
-rwxr-xr-xtriso_study_driver.sh84
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