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

print a warning if the maximum neutron energy is < 20 MeV, close #1138

Browse source
This commit is contained in:
liangjg 2018-12-19 11:40:15 -05:00
parent 493341dbe1
commit 8ae7c783b5
2 changed files with 8 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/input_xml.F90
View file

@ -2258,6 +2258,8 @@ contains
call write_message("Maximum neutron transport energy: " // &
trim(to_str(energy_max(NEUTRON))) // " eV for " // &
trim(adjustl(nuclides(i) % name)), 7)
if (master .and. energy_max(NEUTRON) < 20.E6) call warning("Maximum &
&neutron energy is below 20 MeV. This may bias the results.")
exit
end if
end if
@ -2272,9 +2274,10 @@ contains
exit
end if
end do
if (.not. mp_found) call warning("Windowed multipole functionality is &
&turned on, but no multipole libraries were found. Make sure that &
&windowed multipole data is present in your cross_sections.xml file.")
if (master .and. .not. mp_found) call warning("Windowed multipole &
&functionality is turned on, but no multipole libraries were found. &
&Make sure that windowed multipole data is present in your &
&cross_sections.xml file.")
end if
call already_read % clear()

4
src/nuclide_header.F90
View file

@ -685,8 +685,8 @@ contains
subroutine nuclide_init_grid(this, E_min, E_max, M)
class(Nuclide), intent(inout) :: this
real(8), intent(in) :: E_min ! Minimum energy in MeV
real(8), intent(in) :: E_max ! Maximum energy in MeV
real(8), intent(in) :: E_min ! Minimum energy in eV
real(8), intent(in) :: E_max ! Maximum energy in eV
integer, intent(in) :: M ! Number of equally log-spaced bins
integer :: i, j, k ! Loop indices