crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Moved reading <sab> tag outside do loop for nuclides.

Browse source
This commit is contained in:
Paul Romano 2012-07-01 10:43:31 -04:00
parent 3e46473104
commit 09ede9cb90
1 changed files with 18 additions and 18 deletions
Download patch file Download diff file Expand all files Collapse all files

36
src/input_xml.F90
View file

@ -775,29 +775,29 @@ contains
else
mat % atom_percent(j) = -nuc % wo
end if
end do
! Read S(a,b) table information
if (size(material_(i) % sab) == 1) then
! Get pointer to S(a,b) table
sab => material_(i) % sab(1)
! Read S(a,b) table information
if (size(material_(i) % sab) == 1) then
! Get pointer to S(a,b) table
sab => material_(i) % sab(1)
! Determine name of S(a,b) table
name = trim(sab % name) // "." // trim(sab % xs)
mat % sab_name = name
! Determine name of S(a,b) table
name = trim(sab % name) // "." // trim(sab % xs)
mat % sab_name = name
! Check that this nuclide is listed in the cross_sections.xml file
if (.not. dict_has_key(xs_listing_dict, name)) then
message = "Could not find S(a,b) table " // trim(name) // &
" in cross_sections.xml file!"
call fatal_error()
end if
mat % has_sab_table = .true.
elseif (size(material_(i) % sab) > 1) then
message = "Cannot have multiple S(a,b) tables on a single material."
! Check that this nuclide is listed in the cross_sections.xml file
if (.not. dict_has_key(xs_listing_dict, name)) then
message = "Could not find S(a,b) table " // trim(name) // &
" in cross_sections.xml file!"
call fatal_error()
end if
end do
mat % has_sab_table = .true.
elseif (size(material_(i) % sab) > 1) then
message = "Cannot have multiple S(a,b) tables on a single material."
call fatal_error()
end if
! Check to make sure either all atom percents or all weight percents are
! given