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

Make sure source is in geometry before fiss. check

Browse source
This commit is contained in:
Sterling Harper 2017-04-19 19:13:07 -04:00
parent bc4683be1c
commit c6b67e6443
1 changed files with 12 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

22
src/source.F90
View file

@ -130,7 +130,7 @@ contains
! Repeat sampling source location until a good site has been found
found = .false.
do while (.not.found)
do while (.not. found)
! Set particle defaults
call p % initialize()
@ -145,16 +145,18 @@ contains
call find_cell(p, found)
! Check if spatial site is in fissionable material
select type (space => external_source(i) % space)
type is (SpatialBox)
if (space % only_fissionable) then
if (p % material == MATERIAL_VOID) then
found = .false.
elseif (.not. materials(p % material) % fissionable) then
found = .false.
if (found) then
select type (space => external_source(i) % space)
type is (SpatialBox)
if (space % only_fissionable) then
if (p % material == MATERIAL_VOID) then
found = .false.
elseif (.not. materials(p % material) % fissionable) then
found = .false.
end if
end if
end if
end select
end select
end if
! Check for rejection
if (.not. found) then