From f658bf85ebdc379d157757b5807fa529cbe49750 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Sat, 8 Apr 2017 14:30:28 -0500 Subject: [PATCH] Show error message if 0K scattering not present --- src/input_xml.F90 | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/src/input_xml.F90 b/src/input_xml.F90 index bf981e9f5f..a58cd3af84 100644 --- a/src/input_xml.F90 +++ b/src/input_xml.F90 @@ -5303,6 +5303,13 @@ contains do i = 1, size(res_scat_nuclides) if (nuc % name == res_scat_nuclides(i)) then + ! Make sure nuclide has 0K data + if (.not. allocated(nuc % energy_0K)) then + call fatal_error("Cannot treat " // trim(nuc % name) // " as a & + &resonant scatterer because 0 K elastic scattering data is not & + &present.") + end if + ! Set nuclide to be resonant nuc % resonant = .true.