diff options
Diffstat (limited to 'triso_study_driver.sh')
| -rwxr-xr-x | triso_study_driver.sh | 18 |
1 files changed, 12 insertions, 6 deletions
diff --git a/triso_study_driver.sh b/triso_study_driver.sh index 835d819..4c99c1e 100755 --- a/triso_study_driver.sh +++ b/triso_study_driver.sh @@ -1,20 +1,26 @@ #!/bin/bash +if [[ $1 == "lcl" ]]; then + prog=openmc-inst01.py +else + prog=./hpmr.py +fi + mkdir -p results/ results/rpt results/rrpt touch results/rpt_table results/rrpt_table -echo "\#" `date` > results/rpt_table -echo "\#" `date` > results/rrpt_table +echo "\#" $(date) > results/rpt_table +echo "\#" $(date) > results/rrpt_table # 1. Basic calculation without homogenization echo "Running explicit TRISO calculation..." -./hpmr.py > 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..." -./hpmr.py > 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! @@ -23,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... - ./hpmr.py -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) echo Found $keff +/- $pm for RPT with r=$rad cm... @@ -34,7 +40,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... - ./hpmr.py -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) echo Found $keff +/- $pm for RRPT with r=$rad cm... |