From 8afe3e8c41266e1ece6135e627dbb5cfbf2bfe14 Mon Sep 17 00:00:00 2001 From: Adam Nelson Date: Sat, 14 Nov 2015 09:00:36 -0500 Subject: [PATCH] Updated manual changed inverse velocities units --- docs/source/usersguide/index.rst | 1 + docs/source/usersguide/input.rst | 2 +- src/input_xml.F90 | 13 ++++++------- 3 files changed, 8 insertions(+), 8 deletions(-) diff --git a/docs/source/usersguide/index.rst b/docs/source/usersguide/index.rst index 5a7e7addf..0338c5cef 100644 --- a/docs/source/usersguide/index.rst +++ b/docs/source/usersguide/index.rst @@ -14,6 +14,7 @@ essential aspects of using OpenMC to perform simulations. beginners install input + mgxs_library output/index processing troubleshoot diff --git a/docs/source/usersguide/input.rst b/docs/source/usersguide/input.rst index 3da2382a8..d792851a0 100644 --- a/docs/source/usersguide/input.rst +++ b/docs/source/usersguide/input.rst @@ -1977,7 +1977,7 @@ attributes/sub-elements: automatically assumes a one energy group calculation over the entire energy range. - .. note:: When running in the multi-group :ref:`energy_mode`, these + .. note:: When running in the multi-group :ref:`energy_mode`, the energy bins must match the data library's group boundaries. :albedo: diff --git a/src/input_xml.F90 b/src/input_xml.F90 index b95f9f483..2cc829454 100644 --- a/src/input_xml.F90 +++ b/src/input_xml.F90 @@ -4453,7 +4453,8 @@ contains ! assumed to be the midpoint do i = 1, energy_groups inverse_velocities(i) = & - (sqrt(TWO * energy_bin_avg(i) / (MASS_NEUTRON_MEV)) * C_LIGHT) + (sqrt(TWO * energy_bin_avg(i) / (MASS_NEUTRON_MEV)) * & + C_LIGHT * 100.0_8) end do end if @@ -4485,14 +4486,12 @@ contains if (check_for_node(node_xsdata, "zaid")) then call get_node_value(node_xsdata, "zaid", listing % zaid) else - listing % zaid = 100 + listing % zaid = -1 end if - if (check_for_node(node_xsdata, "kT")) & - call get_node_value(node_xsdata, "kT", listing % kT) - if (check_for_node(node_xsdata, "awr")) then - call get_node_value(node_xsdata, "awr", listing % awr) + if (check_for_node(node_xsdata, "kT")) then + call get_node_value(node_xsdata, "kT", listing % kT) else - listing % awr = ONE + listing % kT = 2.53E-8_8 end if ! determine type of cross section