summaryrefslogtreecommitdiff
path: root/hp_height_study_driver.sh
blob: 234bbaff9200b5bbe5c88acd7b2ac84e9ead29f1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
#!/bin/bash

mkdir -p results/outputs
touch results/keff_table
echo "\#" $(date) > results/keff_table


if [[ $1 == "lcl" ]]; then
    prog=openmc-inst01.py krusty.py
else
    prog=./krusty.py
fi

green="\e[0;32m"
reset="\e[0m"

for pool_height in $(seq 1 45); do
    echo -e "$green Starting runs for $pool_height cm pool height... $reset"
    for sat_height in $(seq $pool_height 45); do
        output_file=results/outputs/p$pool_height\_s$sat_height.out
        $prog -p $pool_height -s $sat_height > $output_file
        keff=$(cat $output_file | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
        pm=$(cat $output_file | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
        md5=$(cat $output_file | grep "MD5:" | cut -d ":" -f 2 | sed "s/ //g")
        echo K-eff $keff +/- $pm found for p=$pool_height and s=$sat_height \($md5\)
        echo $pool_height $sat_height $keff $pm $md5 \($(date)\) >> results/keff_table
    done
done