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

Changed isotropic in lab scattering tag to iso-in-lab

Browse source
This commit is contained in:
Will Boyd 2015-11-08 16:06:34 -05:00
parent 0a0c4c4d5f
commit a11ccce8a3
3 changed files with 9 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

12
docs/source/usersguide/input.rst
View file

@ -1138,12 +1138,12 @@ Each ``material`` element can have the following attributes or sub-elements:
An optional attribute/sub-element for each nuclide is ``scattering``. This
attribute may be set to "ace" to use the scattering laws specified in the
ACE files (default). Alternatively, when set to "lab", the ACE scattering
laws are used to sample the outgoing energy but an isotropic-in-lab
distribution is used to sample the outgoing angle at each scattering
interaction. The ``scattering`` attribute may be most useful when using
OpenMC to compute multi-group cross-sections for deterministic transport
codes and to quantify the effects of anisotropic scattering.
ACE files (default). Alternatively, when set to "iso-in-lab", the ACE
scattering laws are used to sample the outgoing energy but an
isotropic-in-lab distribution is used to sample the outgoing angle at each
scattering interaction. The ``scattering`` attribute may be most useful
when using OpenMC to compute multi-group cross-sections for deterministic
transport codes and to quantify the effects of anisotropic scattering.
*Default*: None

4
openmc/nuclide.py
View file

@ -102,9 +102,9 @@ class Nuclide(object):
@scattering.setter
def scattering(self, scattering):
if not scattering in ['ace', 'lab']:
if not scattering in ['ace', 'iso-in-lab']:
msg = 'Unable to set scattering for Nuclide to {0} ' \
'which is not "ace" or "lab"'.format(scattering)
'which is not "ace" or "iso-in-lab"'.format(scattering)
raise ValueError(msg)
self._scattering = scattering

2
src/input_xml.F90
View file

@ -1938,7 +1938,7 @@ contains
! Check enforced isotropic lab scattering
if (check_for_node(node_nuc, "scattering")) then
call get_node_value(node_nuc, "scattering", temp_str)
if (trim(adjustl(to_lower(temp_str))) == "lab") then
if (trim(adjustl(to_lower(temp_str))) == "iso-in-lab") then
call list_iso_lab % append(1)
else if (trim(adjustl(to_lower(temp_str))) == "ace") then
call list_iso_lab % append(0)