diff options
Diffstat (limited to 'triso_study_driver.sh')
| -rwxr-xr-x | triso_study_driver.sh | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/triso_study_driver.sh b/triso_study_driver.sh index 4f8e5d8..b7404c8 100755 --- a/triso_study_driver.sh +++ b/triso_study_driver.sh @@ -1,7 +1,7 @@ #!/bin/bash if [[ $1 == "lcl" ]]; then - prog=openmc-inst01.py hpmr.py + prog=(openmc-inst01.py hpmr.py -l) else prog=./hpmr.py fi @@ -13,14 +13,14 @@ echo "\#" $(date) > results/rrpt_table # 1. Basic calculation without homogenization echo "Running explicit TRISO calculation..." -$prog -t none > results/explicit.out +"${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/rpt/explicit.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2) echo Found $keff +/- $pm for explicit! # 2. VWH calculation echo "Running VWH homogenization calculation..." -$prog -t vwh > results/vwh.out +"${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/rpt/vwh.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2) echo Found $keff +/- $pm for VWH! @@ -29,7 +29,7 @@ echo Found $keff +/- $pm for VWH! # 0.025 step size -> 21 iterations for rad in $(seq 0.65 0.025 1.15); do echo Running RPT with $rad cm... - $prog -t rpt -r $rad > results/rpt/$rad.out + "${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") @@ -41,7 +41,7 @@ done # 0.025 step size -> 35 iterations for rad in $(seq 0.1 0.025 0.95); do echo Running RRPT with $rad cm... - $prog -t rrpt -r $rad > results/rrpt/$rad.out + "${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") |