From 018ed9f58a0c98525f981a18bc04a873ea965988 Mon Sep 17 00:00:00 2001 From: yardasol Date: Thu, 21 Jul 2022 17:47:50 -0500 Subject: [PATCH] add flux parameter to SII integrator --- openmc/deplete/abc.py | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/openmc/deplete/abc.py b/openmc/deplete/abc.py index c9e0d566b9..24859151f7 100644 --- a/openmc/deplete/abc.py +++ b/openmc/deplete/abc.py @@ -526,9 +526,9 @@ class Integrator(ABC): initial heavy metal inventory to get total power if ``power`` is not speficied. flux : float or iterable of float, optional - Neutron flux in [neur/s-cm^2] for each interval in :attr: `timesteps` + neutron flux in [neut/s-cm^2] for each interval in :attr:`timesteps` source_rates : float or iterable of float, optional - Source rate in [neutron/sec] for each interval in :attr:`timesteps` + source rate in [neutron/sec] for each interval in :attr:`timesteps` .. versionadded:: 0.12.1 timestep_units : {'s', 'min', 'h', 'd', 'MWd/kg'} @@ -866,6 +866,8 @@ class SIIntegrator(Integrator): Power density of the reactor in [W/gHM]. It is multiplied by initial heavy metal inventory to get total power if ``power`` is not speficied. + flux : float or iterable of float, optional + neutron flux in [neut/s-cm^2] for each interval in :attr:`timesteps` source_rates : float or iterable of float, optional Source rate in [neutron/sec] for each interval in :attr:`timesteps` @@ -922,12 +924,12 @@ class SIIntegrator(Integrator): """ def __init__(self, operator, timesteps, power=None, power_density=None, - source_rates=None, timestep_units='s', n_steps=10, + flux=None, source_rates=None, timestep_units='s', n_steps=10, solver="cram48"): check_type("n_steps", n_steps, Integral) check_greater_than("n_steps", n_steps, 0) super().__init__( - operator, timesteps, power, power_density, source_rates, + operator, timesteps, power, power_density, flux, source_rates, timestep_units=timestep_units, solver=solver) self.n_steps = n_steps