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

Fix bug in reimplementation of neighbor_lists.

Browse source
This commit is contained in:
Paul Romano 2015-11-18 07:57:06 -06:00
parent 7661c8902f
commit 0927f25300
1 changed files with 8 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/geometry.F90
View file

@ -903,13 +903,17 @@ contains
do i = 1, n_surfaces
! Copy positive neighbors to Surface instance
j = neighbor_pos(i)%size()
allocate(surfaces(i)%obj%neighbor_pos(j))
surfaces(i)%obj%neighbor_pos(:) = neighbor_pos(i)%data(1:j)
if (j > 0) then
allocate(surfaces(i)%obj%neighbor_pos(j))
surfaces(i)%obj%neighbor_pos(:) = neighbor_pos(i)%data(1:j)
end if
! Copy negative neighbors to Surface instance
j = neighbor_neg(i)%size()
allocate(surfaces(i)%obj%neighbor_neg(j))
surfaces(i)%obj%neighbor_neg(:) = neighbor_neg(i)%data(1:j)
if (j > 0) then
allocate(surfaces(i)%obj%neighbor_neg(j))
surfaces(i)%obj%neighbor_neg(:) = neighbor_neg(i)%data(1:j)
end if
end do
end subroutine neighbor_lists