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

Fix bug in openmc_tally_set_scores (for depletion reactions)

Browse source
This commit is contained in:
Paul Romano 2018-02-09 14:51:36 -06:00
parent 37f552a5dc
commit 0a772252cc
1 changed files with 9 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/tallies/tally_header.F90
View file

@ -732,8 +732,10 @@ contains
integer :: MT
character(C_CHAR), pointer :: string(:)
character(len=:, kind=C_CHAR), allocatable :: score_
logical :: depletion_rx
err = E_UNASSIGNED
depletion_rx = .false.
if (index >= 1 .and. index <= size(tallies)) then
associate (t => tallies(index) % obj)
if (allocated(t % score_bins)) deallocate(t % score_bins)
@ -757,10 +759,13 @@ contains
t % score_bins(i) = SCORE_NU_SCATTER
case ('(n,2n)')
t % score_bins(i) = N_2N
depletion_rx = .true.
case ('(n,3n)')
t % score_bins(i) = N_3N
depletion_rx = .true.
case ('(n,4n)')
t % score_bins(i) = N_4N
depletion_rx = .true.
case ('absorption')
t % score_bins(i) = SCORE_ABSORPTION
case ('fission', '18')
@ -829,8 +834,10 @@ contains
t % score_bins(i) = N_NC
case ('(n,gamma)')
t % score_bins(i) = N_GAMMA
depletion_rx = .true.
case ('(n,p)')
t % score_bins(i) = N_P
depletion_rx = .true.
case ('(n,d)')
t % score_bins(i) = N_D
case ('(n,t)')
@ -839,6 +846,7 @@ contains
t % score_bins(i) = N_3HE
case ('(n,a)')
t % score_bins(i) = N_A
depletion_rx = .true.
case ('(n,2a)')
t % score_bins(i) = N_2A
case ('(n,3a)')
@ -879,6 +887,7 @@ contains
end do
err = 0
t % depletion_rx = depletion_rx
end associate
else
err = E_OUT_OF_BOUNDS