Python scripts to automate SCALE actions

[Arielle Miller]

Does anyone have any python scripts they have been using to modify the input files to change small things like the radius of a sphere over multiple files?  Or a script that pulls a data string out of multiple output files and dumps it into a text file (like keff and signma)?  I'm not familiar with how to do this but i know it can be done.  Any help would be really appreciated.

[Alex]

Arielle,

Here is a script that I use to create multiple input files with different temperatures based on a single base input file. Don't know whether you'll find it a useful example or not.

Alex

from shutil import copyfile
import fileinput
import re
from math import exp

# These functions for density are based on volumetric expansion coefficients
# given on page 40 of MSRE Design and Operations Report, Part iii, Nuclear
# Analysis

def fdens(temperature):
    alpha_fuel_faeh = 1.18e-4
    alpha_fuel_kel = 1.8 * alpha_fuel_faeh
    rho0_fuel = 2.146
    T0 = 922
    return rho0_fuel * exp(-alpha_fuel_kel * (temperature - T0))

def mdens(temperature):
    alpha_graph_faeh = 1.0e-5
    alpha_graph_kel = 1.8 * alpha_graph_faeh
    rho0_graph = 1.86
    T0 = 922
    return rho0_graph * exp(-alpha_graph_kel * (temperature - T0))

start = 922
stop = 1022
step = 50
root = "calc_buckling"
for temp in range(start, stop + step, step):
    file_name = root + "_" + str(temp) + ".inp"
    copyfile(root + ".inp", file_name)
    for line in fileinput.input(file_name, inplace=True):
        if re.search(r'2\.146.*922', line):
            print(line.replace(line, re.sub(r'2\.146(.*)922', '{:.3f}'.format(fdens(temp)) + r'\g<1>' + str(temp), line)), end="")
        elif re.search(r'1\.86.*922', line):
            print(line.replace(line, re.sub(r'1\.86(.*)922', '{:.3f}'.format(mdens(temp)) + r'\g<1>' + str(temp), line)), end="")
        else:
            print(line, end="")