def step1():
config = ConfigParser.ConfigParser()
config.read(config_filename)
structured = config.getboolean("general", "structured")
sub_voxel = config.getboolean("general", "sub_voxel")
meshtal = config.get("step1", "meshtal")
tally_num = config.getint("step1", "tally_num")
flux_tag = config.get("step1", "flux_tag")
decay_times = config.get("step2", "decay_times").split(",")
geom = config.get("step1", "geom")
reverse = config.getboolean("step1", "reverse")
num_rays = config.getint("step1", "num_rays")
grid = config.getboolean("step1", "grid")
load(geom)
# get meshtal info from meshtal file
flux_mesh = resolve_mesh(meshtal, tally_num, flux_tag)
# create the cell_fracs array before irradiation_steup
if flux_mesh.structured:
cell_fracs = discretize_geom(flux_mesh, num_rays=num_rays, grid=grid)
# tag cell fracs for both default and subvoxel r2s modes
flux_mesh.tag_cell_fracs(cell_fracs)
else:
cell_fracs = discretize_geom(flux_mesh)
cell_mats = cell_materials(geom)
irradiation_setup(
flux_mesh,
cell_mats,
cell_fracs,
alara_params_filename,
tally_num,
num_rays=num_rays,
grid=grid,
reverse=reverse,
flux_tag=flux_tag,
decay_times=decay_times,
sub_voxel=sub_voxel,
)
# create a blank mesh for step 2:
ves = list(flux_mesh.iter_ve())
tags_keep = (
"cell_number",
"cell_fracs",
"cell_largest_frac_number",
"cell_largest_frac",
)
for tag in flux_mesh.get_all_tags():
if tag.name not in tags_keep and isinstance(tag, NativeMeshTag):
tag.delete()
flux_mesh.write_hdf5("blank_mesh.h5m")
print("The file blank_mesh.h5m has been saved to disk.")
print("Do not delete this file; it is needed by r2s.py step2.\n")
print("R2S step1 complete, run ALARA with the command:")
print(">> alara alara_inp > output.txt")
def step1():
config = ConfigParser.ConfigParser()
config.read(config_filename)
structured = config.getboolean("general", "structured")
meshtal = config.get("step1", "meshtal")
tally_num = config.getint("step1", "tally_num")
flux_tag = config.get("step1", "flux_tag")
decay_times = config.get("step2", "decay_times").split(",")
geom = config.get("step1", "geom")
reverse = config.getboolean("step1", "reverse")
num_rays = config.getint("step1", "num_rays")
grid = config.getboolean("step1", "grid")
responses = config.get("general", "responses").split(",")
wdr_file = config.get("general", "wdr_file")
load(geom)
# get meshtal info from meshtal file
flux_mesh = resolve_mesh(meshtal, tally_num, flux_tag)
# create the cell_fracs array before irradiation_steup
if flux_mesh.structured:
cell_fracs = discretize_geom(flux_mesh, num_rays=num_rays, grid=grid)
else:
cell_fracs = discretize_geom(flux_mesh)
cell_mats = cell_materials(geom)
irradiation_setup(
flux_mesh,
cell_mats,
cell_fracs,
alara_params_filename,
tally_num,
num_rays=num_rays,
grid=grid,
reverse=reverse,
output_mesh="activation_responses_step1.h5m",
flux_tag=flux_tag,
decay_times=decay_times,
sub_voxel=False,
responses=responses,
wdr_file=wdr_file,
)
# create a blank mesh for step 2:
ves = list(flux_mesh.iter_ve())
for tag in flux_mesh.get_all_tags():
tag.delete()
flux_mesh.write_hdf5("blank_mesh.h5m")
print("The file blank_mesh.h5m has been saved to disk.")
print(
"Do not delete this file; it is needed by activation_responses.py step2.\n"
)
print("Decay_heat step1 complete, run ALARA with the command:")
print(">> alara alara_inp > output.txt")
def step1():
config = ConfigParser.ConfigParser()
config.read(config_filename)
structured = config.getboolean('general', 'structured')
meshtal = config.get('step1', 'meshtal')
tally_num = config.getint('step1', 'tally_num')
flux_tag = config.get('step1', 'flux_tag')
decay_times = config.get('step2', 'decay_times').split(',')
geom = config.get('step1', 'geom')
reverse = config.getboolean('step1', 'reverse')
num_rays = config.getint('step1', 'num_rays')
grid = config.getboolean('step1', 'grid')
responses = config.get('general', 'responses').split(',')
wdr_file = config.get('general', 'wdr_file')
load(geom)
# get meshtal info from meshtal file
flux_mesh = resolve_mesh(meshtal, tally_num, flux_tag)
# create the cell_fracs array before irradiation_steup
if flux_mesh.structured:
cell_fracs = discretize_geom(flux_mesh, num_rays=num_rays, grid=grid)
else:
cell_fracs = discretize_geom(flux_mesh)
cell_mats = cell_materials(geom)
irradiation_setup(flux_mesh,
cell_mats,
cell_fracs,
alara_params_filename,
tally_num,
num_rays=num_rays,
grid=grid,
reverse=reverse,
output_mesh="activation_responses_step1.h5m",
flux_tag=flux_tag,
decay_times=decay_times,
sub_voxel=False,
responses=responses,
wdr_file=wdr_file)
# create a blank mesh for step 2:
ves = list(flux_mesh.iter_ve())
for tag in flux_mesh.get_all_tags():
tag.delete()
flux_mesh.write_hdf5('blank_mesh.h5m')
print('The file blank_mesh.h5m has been saved to disk.')
print(
'Do not delete this file; it is needed by activation_responses.py step2.\n'
)
print('Decay_heat step1 complete, run ALARA with the command:')
print('>> alara alara_inp > output.txt')
def step1():
config = ConfigParser.ConfigParser()
config.read(config_filename)
structured = config.getboolean('general', 'structured')
sub_voxel = config.getboolean('general', 'sub_voxel')
meshtal = config.get('step1', 'meshtal')
tally_num = config.getint('step1', 'tally_num')
flux_tag = config.get('step1', 'flux_tag')
decay_times = config.get('step2', 'decay_times').split(',')
geom = config.get('step1', 'geom')
reverse = config.getboolean('step1', 'reverse')
num_rays = config.getint('step1', 'num_rays')
grid = config.getboolean('step1', 'grid')
load(geom)
# get meshtal info from meshtal file
flux_mesh = resolve_mesh(meshtal, tally_num, flux_tag)
# create the cell_fracs array before irradiation_steup
if flux_mesh.structured:
cell_fracs = discretize_geom(flux_mesh, num_rays=num_rays, grid=grid)
# tag cell fracs for both default and subvoxel r2s modes
flux_mesh.tag_cell_fracs(cell_fracs)
else:
cell_fracs = discretize_geom(flux_mesh)
cell_mats = cell_materials(geom)
irradiation_setup(flux_mesh,
cell_mats,
cell_fracs,
alara_params_filename,
tally_num,
num_rays=num_rays,
grid=grid,
reverse=reverse,
flux_tag=flux_tag,
decay_times=decay_times,
sub_voxel=sub_voxel)
# create a blank mesh for step 2:
ves = list(flux_mesh.iter_ve())
tags_keep = ("cell_number", "cell_fracs", "cell_largest_frac_number",
"cell_largest_frac")
for tag in flux_mesh.get_all_tags():
if tag.name not in tags_keep and isinstance(tag, NativeMeshTag):
tag.delete()
flux_mesh.write_hdf5('blank_mesh.h5m')
print('The file blank_mesh.h5m has been saved to disk.')
print('Do not delete this file; it is needed by r2s.py step2.\n')
print('R2S step1 complete, run ALARA with the command:')
print('>> alara alara_inp > output.txt')