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

Fixed bug in write_tallies(): filter_matches index j --> t % filter(j)

Browse source
This commit is contained in:
amandalund 2017-05-09 12:09:55 -05:00
parent cf6b71f719
commit cd58aaefa7
1 changed files with 6 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/output.F90
View file

@ -816,18 +816,18 @@ contains
if (size(t % filter) == 0) exit find_bin
! Increment bin combination
filter_matches(j) % bins % data(1) = &
filter_matches(j) % bins % data(1) + 1
filter_matches(t % filter(j)) % bins % data(1) = &
filter_matches(t % filter(j)) % bins % data(1) + 1
! =================================================================
! REACHED END OF BINS FOR THIS FILTER, MOVE TO NEXT FILTER
if (filter_matches(j) % bins % data(1) > &
if (filter_matches(t % filter(j)) % bins % data(1) > &
filters(t % filter(j)) % obj % n_bins) then
! If this is the first filter, then exit
if (j == 1) exit print_bin
filter_matches(j) % bins % data(1) = 0
filter_matches(t % filter(j)) % bins % data(1) = 0
j = j - 1
indent = indent - 2
@ -841,7 +841,7 @@ contains
! Print current filter information
write(UNIT=unit_tally, FMT='(1X,2A)') repeat(" ", indent), &
trim(filters(t % filter(j)) % obj % &
text_label(filter_matches(j) % bins % data(1)))
text_label(filter_matches(t % filter(j)) % bins % data(1)))
indent = indent + 2
j = j + 1
end if
@ -852,7 +852,7 @@ contains
if (size(t % filter) > 0) then
write(UNIT=unit_tally, FMT='(1X,2A)') repeat(" ", indent), &
trim(filters(t % filter(j)) % obj % &
text_label(filter_matches(j) % bins % data(1)))
text_label(filter_matches(t % filter(j)) % bins % data(1)))
end if
! Determine scoring index for this bin combination -- note that unlike