crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Updating comments

Browse source
This commit is contained in:
Patrick Shriwise 2022-03-10 15:20:40 -06:00
parent 7abb11f7c4
commit 7b8bc4ae90
1 changed files with 11 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

23
openmc/weight_windows.py
View file

@ -570,21 +570,20 @@ def wwinp_to_wws(path):
e_bounds *= 1E6
# create an array for weight window lower bounds
# # MCNP ordering for weight windows matches that of OpenMC
# # ('xyz' with x changing fastest)
ww_lb = np.zeros((*mesh.dimension, ne))
for ijk in mesh.indices:
idx = tuple([v - 1 for v in ijk] + [slice(None)])
ww_lb[idx] = [float(next(wwinp)) for _ in range(ne)]
# for e in range(ne):
# # MCNP ordering for weight windows matches that of OpenMC
# # ('xyz' with x changing fastest)
# ww_vals = [float(next(wwinp)) for _ in range(n_elements)]
# ww_lb[e, :] = np.asarray(ww_vals).reshape(mesh.dimension)
settings = WeightWindows(id=None,
mesh=mesh,
lower_ww_bounds=ww_lb.flatten(),
upper_bound_ratio=5.0,
energy_bins=e_bounds,
particle_type=particle)
wws.append(settings)
# create a WeightWindows object and add it to the output list
ww = WeightWindows(id=None,
mesh=mesh,
lower_ww_bounds=ww_lb.flatten(),
upper_bound_ratio=5.0,
energy_bins=e_bounds,
particle_type=particle)
wws.append(ww)
return wws