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

Prevent use of outgoing energy filters for certain macro tallies.

Browse source
This commit is contained in:
Paul Romano 2011-10-18 11:38:32 -04:00
parent a8efc294cd
commit c6d85b082d
1 changed files with 25 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

25
src/input_xml.f90
View file

@ -745,8 +745,17 @@ contains
select case (trim(word))
case ('flux')
t % macro_bins(j) % scalar = MACRO_FLUX
if (t % n_bins(T_ENERGYOUT) > 0) then
msg = "Cannot tally flux with an outgoing energy filter."
call fatal_error(msg)
end if
case ('total')
t % macro_bins(j) % scalar = MACRO_TOTAL
if (t % n_bins(T_ENERGYOUT) > 0) then
msg = "Cannot tally total reaction rate with an outgoing " &
// "energy filter."
call fatal_error(msg)
end if
case ('scatter')
t % macro_bins(j) % scalar = MACRO_SCATTER
case ('nu-scatter')
@ -767,10 +776,26 @@ contains
t % macro_bins(j) % scalar = MACRO_N_4N
case ('absorption')
t % macro_bins(j) % scalar = MACRO_ABSORPTION
if (t % n_bins(T_ENERGYOUT) > 0) then
msg = "Cannot tally absorption rate with an outgoing " &
// "energy filter."
call fatal_error(msg)
end if
case ('fission')
t % macro_bins(j) % scalar = MACRO_FISSION
if (t % n_bins(T_ENERGYOUT) > 0) then
msg = "Cannot tally fission rate with an outgoing " &
// "energy filter."
call fatal_error(msg)
end if
case ('nu-fission')
t % macro_bins(j) % scalar = MACRO_NU_FISSION
! TODO: Add fission energy transfer and remove this
if (t % n_bins(T_ENERGYOUT) > 0) then
msg = "Cannot tally nu-fission rate with an outgoing " &
// "energy filter."
call fatal_error(msg)
end if
case default
msg = "Unknown macro reaction: " // trim(words(j))
call fatal_error(msg)