From e2fcd68d1f4dffaa739fe3f6bda98d6913146e2b Mon Sep 17 00:00:00 2001 From: Sterling Harper Date: Mon, 10 Jul 2017 14:43:18 -0400 Subject: [PATCH] Fatal error for nuclides in multiple S(a,b) tables --- src/input_xml.F90 | 13 +++++++++++++ 1 file changed, 13 insertions(+) diff --git a/src/input_xml.F90 b/src/input_xml.F90 index 2e7219a31..ef742ed45 100644 --- a/src/input_xml.F90 +++ b/src/input_xml.F90 @@ -5193,6 +5193,19 @@ contains end if end do ASSIGN_SAB + ! Make sure each nuclide only appears in one table. + do j = 1, i_sab_nuclides % size() + do k = j+1, i_sab_nuclides % size() + if (i_sab_nuclides % data(j) == i_sab_nuclides % data(k)) then + call fatal_error(trim( & + nuclides(mat % nuclide(i_sab_nuclides % data(j))) % name) & + // " in material " // trim(to_str(mat % id)) // " was found & + &in multiple S(a,b) tables. Each nuclide can only appear in & + &one S(a,b) table per material.") + end if + end do + end do + ! Update i_sab_tables and i_sab_nuclides deallocate(mat % i_sab_tables) deallocate(mat % sab_fracs)