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

added coremap routine for no reflector acceleration

Browse source
This commit is contained in:
Bryan Herman 2012-01-03 14:27:53 -08:00
parent a43db54057
commit 3bc1f734fc
1 changed files with 49 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

49
src/cmfd_utils.F90
View file

@ -414,4 +414,53 @@ contains
end subroutine neutron_balance
!===============================================================================
! SET_COREMAP is a routine that sets the core mapping information
!===============================================================================
subroutine set_coremap
integer :: kount=1 ! counter for unique fuel assemblies
integer :: nx ! number of mesh cells in x direction
integer :: ny ! number of mesh cells in y direction
integer :: nz ! number of mesh cells in z direction
integer :: ng ! number of energy groups
integer :: i ! iteration counter for x
integer :: j ! iteration counter for y
integer :: k ! iteration counter for z
! extract spatial indices from object
nx = cmfd % indices(1)
ny = cmfd % indices(2)
nz = cmfd % indices(3)
! begin loops over spatial indices
ZLOOP: do i = 1,nz
YLOOP: do j = 1,ny
XLOOP: do k = 1,nx
! check for reflector
if (cmfd % coremap(i,j,k) == 1) then
! reset value to 99999
cmfd % coremap(i,j,k) = 99999
else
! must be a fuel --> give unique id number
cmfd % coremap(i,j,k) = kount
kount = kount + 1
end if
end do XLOOP
end do YLOOP
end do ZLOOP
end subroutine set_coremap
end module cmfd_utils