From c18a6eafa4dadf8ac998b04b7ad3e3502cadef23 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Sun, 1 Sep 2013 22:26:11 -0400 Subject: [PATCH] Check for valid secondary mode on S(a,b) tables. --- src/ace.F90 | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/src/ace.F90 b/src/ace.F90 index 2617bcbed4..f83921b858 100644 --- a/src/ace.F90 +++ b/src/ace.F90 @@ -1178,8 +1178,14 @@ contains integer :: NMU ! number of outgoing angles integer :: JXS4 ! location of elastic energy table - ! read secondary energy mode for inelastic scattering + ! read secondary energy mode for inelastic scattering and check table % secondary_mode = NXS(7) + if (table % secondary_mode /= SAB_SECONDARY_EQUAL .and. & + table % secondary_mode /= SAB_SECONDARY_SKEWED) then + message = "Unsupported secondary mode on S(a,b) table " // & + trim(adjustl(table % name)) // ": " // to_str(table % secondary_mode) + call fatal_error() + end if ! read number of inelastic energies and allocate arrays NE_in = int(XSS(JXS(1)))