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

Add energy checks suggested by @walshjon

Browse source
This commit is contained in:
Paul Romano 2017-03-25 16:58:37 -05:00
parent f49e0f9cf1
commit 60e4d6c539
2 changed files with 13 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/input_xml.F90
View file

@ -885,8 +885,9 @@ contains
if (check_for_node(node_res_scat, "energy_max")) then
call get_node_value(node_res_scat, "energy_max", res_scat_energy_max)
end if
if (res_scat_energy_max < ZERO) then
call fatal_error("Upper resonance scattering energy bound is negative")
if (res_scat_energy_max < res_scat_energy_min) then
call fatal_error("Upper resonance scattering energy bound is below the &
&lower resonance scattering energy bound.")
end if
! Get nuclides that resonance scattering should be applied to
@ -5301,6 +5302,15 @@ contains
end do
end associate
! Check to make sure resonance scattering upper energy is below URR
! minimum, if present
if (nuc % urr_present) then
if (res_scat_energy_max > nuc % urr_data(1) % energy(1)) then
call fatal_error("Resonance scattering maximum energy is above the &
&bottom of the unresolved resonance region.")
end if
end if
exit
end if
end do

2
tests/test_resonance_scattering/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
1.439342E+00 1.224486E-02
1.441041E+00 8.832991E-03