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

Fixed bug in tallies that depend on post-collision parameters.

Browse source
This commit is contained in:
Paul Romano 2011-10-12 21:59:39 -04:00
parent 1fb0e4fe99
commit 4edbee8f71
3 changed files with 14 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

1
src/mpi_routines.f90
View file

@ -356,6 +356,7 @@ contains
p % xyz_local = temp_bank(i) % xyz
p % uvw = temp_bank(i) % uvw
p % E = temp_bank(i) % E
p % last_E = p % E
! set defaults
call initialize_particle(p)

1
src/source.f90
View file

@ -89,6 +89,7 @@ contains
! set particle energy
p % E = E
p % last_E = E
end do
end if
end do

16
src/tally.f90
View file

@ -217,6 +217,14 @@ contains
end do
do i = 1, n_cells
do j = 1, size(tally_maps(T_CELL) % items(i) % elements)
print *, cells(i) % uid, &
tally_maps(T_CELL) % items(i) % elements(j) % index_tally, &
tally_maps(T_CELL) % items(i) % elements(j) % index_bin
end do
end do
end subroutine create_tally_map
!===============================================================================
@ -259,10 +267,10 @@ contains
! SCORE_TALLY contains the main logic for scoring user-specified tallies
!===============================================================================
subroutine score_tally(p, scatter)
subroutine score_tally(p, scattered)
type(Particle), pointer :: p
logical, intent(in) :: scatter
logical, intent(in) :: scattered
integer :: i
integer :: j
@ -290,7 +298,7 @@ contains
do i = 1, n_tallies
t => tallies(i)
! =======================================================================
! DETERMINE SCORING BIN COMBINATION
@ -395,7 +403,7 @@ contains
! If this tally has an outgoing energy filter, the only supported
! reaction is scattering. For all other reactions, about
if (.not. scatter) return
if (.not. scattered) cycle
! Make sure bin is scattering -- since scattering has already
! occured, we do not need to multiply by the scattering cross