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

Another fix for mesh tracklength tallies.

Browse source
This commit is contained in:
Paul Romano 2012-03-26 19:51:37 -04:00
parent 1466eca5fc
commit 54feb3fe20
2 changed files with 12 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/physics.F90
View file

@ -81,9 +81,6 @@ contains
! Select smaller of the two distances
distance = min(d_boundary, d_collision)
! Save coordinates for tallying purposes
p % last_xyz = p % coord0 % xyz
! Advance particle
coord => p % coord0
do while (associated(coord))
@ -131,6 +128,9 @@ contains
p % surface = NONE
call collision()
! Save coordinates for tallying purposes
p % last_xyz = p % coord0 % xyz
! Set last material to none since cross sections will need to be
! re-evaluated
p % last_material = NONE

13
src/tally.F90
View file

@ -552,7 +552,7 @@ contains
! since multiple bins can be scored to with a single track
if (t % n_filter_bins(FILTER_MESH) > 0) then
call score_tl_on_mesh(tracklength_tallies(i))
call score_tl_on_mesh(tracklength_tallies(i), distance)
cycle
end if
@ -621,9 +621,10 @@ contains
! these tallies, it is possible to score to multiple mesh cells for each track.
!===============================================================================
subroutine score_tl_on_mesh(index_tally)
subroutine score_tl_on_mesh(index_tally, d_track)
integer, intent(in) :: index_tally
real(8), intent(in) :: d_track
integer :: i ! loop index for filter/score bins
integer :: j ! loop index for direction
@ -657,8 +658,8 @@ contains
! CHECK IF THIS TRACK INTERSECTS THE MESH
! Copy starting and ending location of particle
xyz0 = p % last_xyz
xyz1 = p % coord0 % xyz
xyz0 = p % coord0 % xyz - (d_track - TINY_BIT) * p % coord0 % uvw
xyz1 = p % coord0 % xyz - TINY_BIT * p % coord0 % uvw
! Determine indices for starting and ending location
m => meshes(t % mesh)
@ -671,6 +672,10 @@ contains
if (.not. mesh_intersects(m, xyz0, xyz1)) return
end if
! Reset starting and ending location
xyz0 = p % coord0 % xyz - d_track * p % coord0 % uvw
xyz1 = p % coord0 % xyz
! =========================================================================
! CHECK FOR SCORING COMBINATION FOR FILTERS OTHER THAN MESH