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

Added check for particle on upper_right of lattice

Browse source
This commit is contained in:
Nicholas Horelik 2013-04-18 22:01:17 -03:00
parent 40543168ae
commit c9d0bb94f7
1 changed files with 12 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/geometry.F90
View file

@ -78,6 +78,7 @@ contains
integer :: n ! number of cells to search
integer :: index_cell ! index in cells array
real(8) :: xyz(3) ! temporary location
real(8) :: upper_right(3) ! lattice upper_right
logical :: use_search_cells ! use cells provided as argument
logical :: outside_lattice ! if particle is not inside lattice bounds
logical :: lattice_edge ! if particle is on a lattice edge
@ -187,12 +188,21 @@ contains
i_z < 1 .or. i_z > n_z) then
! Check for when particle is on lattice edge
upper_right(1) = lat % lower_left(1) + &
lat % width(1) * dble(lat % dimension(1))
upper_right(2) = lat % lower_left(2) + &
lat % width(2) * dble(lat % dimension(2))
if ( abs(xyz(1) - lat % lower_left(1)) < FP_COINCIDENT .or. &
abs(xyz(2) - lat % lower_left(2)) < FP_COINCIDENT) then
abs(xyz(2) - lat % lower_left(2)) < FP_COINCIDENT .or. &
abs(upper_right(1) - xyz(1)) < FP_COINCIDENT .or. &
abs(upper_right(2) - xyz(2)) < FP_COINCIDENT) then
lattice_edge = .true.
end if
if (lat % n_dimension == 3) then
if (abs(xyz(3) - lat % lower_left(3)) < FP_COINCIDENT) then
upper_right(3) = lat % lower_left(3) + &
lat % width(3) * dble(lat % dimension(3))
if (abs(xyz(3) - lat % lower_left(3)) < FP_COINCIDENT .or. &
abs(upper_right(3) - xyz(3)) < FP_COINCIDENT) then
lattice_edge = .true.
end if
end if