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

Added plot_xs to mgxs-part-iv notebook

Browse source
This commit is contained in:
Adam Nelson 2016-11-20 20:13:40 -05:00
parent f76ee0867c
commit 823d6fe1be
2 changed files with 146 additions and 65 deletions
Download patch file Download diff file Expand all files Collapse all files

192
docs/source/pythonapi/examples/mgxs-part-iv.ipynb
View file

File diff suppressed because one or more lines are too long

19
openmc/plotter.py
View file

@ -197,6 +197,10 @@ def plot_xs(this, types, divisor_types=None, temperature=294., axis=None,
plot_func(E, data[i, :], label=types[i])
ax.set_xlabel('Energy [eV]')
if plot_CE:
ax.set_xlim(1.E-5, 20.E6)
else:
ax.set_xlim(E[-1], E[0])
if divisor_types:
if data_type == 'nuclide':
ylabel = 'Nuclidic Microscopic Data'
@ -214,7 +218,10 @@ def plot_xs(this, types, divisor_types=None, temperature=294., axis=None,
ax.set_ylabel(ylabel)
ax.legend(loc='best')
if this.name is not None:
ax.set_title('Cross Section for ' + this.name)
if len(types) > 1:
ax.set_title('Cross Sections for ' + this.name)
else:
ax.set_title('Cross Section for ' + this.name)
return fig
@ -649,11 +656,19 @@ def calculate_mgxs(this, types, orders=None, temperature=294.,
# Convert the data to the format needed
data = np.zeros((len(types), 2 * library.energy_groups.num_groups))
energy_grid = np.zeros(2 * library.energy_groups.num_groups)
i = 0
for g in range(library.energy_groups.num_groups):
energy_grid[i: i + 2] = library.energy_groups.group_edges[g: g + 2]
i += 2
# Ensure the energy will show on a log-axis by replacing 0s with a
# sufficiently small number
if energy_grid[0] <= 0.:
energy_grid[0] = 1.E-5
for line in range(len(types)):
i = 0
for g in range(library.energy_groups.num_groups):
data[line, i: i + 2] = mgxs[line, g]
energy_grid[i: i + 2] = library.energy_groups.group_edges[g: g + 2]
i += 2
return np.flipud(energy_grid), data