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

Changed count_bank_sites to use weight of bank sites.

Browse source
This commit is contained in:
Paul Romano 2012-03-31 12:59:44 -04:00
parent a47c7ab91c
commit 1fb88c1360
2 changed files with 10 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/intercycle.F90
View file

@ -4,7 +4,7 @@ module intercycle
use error, only: fatal_error, warning
use global
use mesh, only: count_fission_sites
use mesh, only: count_bank_sites
use mesh_header, only: StructuredMesh
use output, only: write_message
use random_lcg, only: prn, set_particle_seed, prn_skip
@ -320,7 +320,7 @@ contains
end if
! count number of fission sites over mesh
call count_fission_sites(m, fission_bank, entropy_p, total, &
call count_bank_sites(m, fission_bank, entropy_p, total, &
size_bank=n_bank, sites_outside=sites_outside)
! display warning message if there were sites outside entropy box

16
src/mesh.F90
View file

@ -145,13 +145,13 @@ contains
end subroutine bin_to_mesh_indices
!===============================================================================
! COUNT_FISSION_SITES determines the number of fission bank sites in each cell
! of a given mesh as well as an optional energy group structure. This can be
! used for a variety of purposes (Shannon entropy, CMFD, uniform fission source
! COUNT_BANK_SITES determines the number of fission bank sites in each cell of a
! given mesh as well as an optional energy group structure. This can be used for
! a variety of purposes (Shannon entropy, CMFD, uniform fission source
! weighting)
!===============================================================================
subroutine count_fission_sites(m, bank_array, cnt, total, &
subroutine count_bank_sites(m, bank_array, cnt, total, &
energies, size_bank, sites_outside)
type(StructuredMesh), pointer :: m ! mesh to count sites
@ -206,7 +206,7 @@ contains
end if
! add weight
weight = weight + ONE
weight = weight + bank_array(i) % wgt
! determine energy bin
if (present(energies)) then
@ -222,8 +222,8 @@ contains
end if
! add to appropriate mesh box
! TODO: if tracking weight through bank, add weight instead
cnt(e_bin,ijk(1),ijk(2),ijk(3)) = cnt(e_bin,ijk(1),ijk(2),ijk(3)) + 1
cnt(e_bin,ijk(1),ijk(2),ijk(3)) = cnt(e_bin,ijk(1),ijk(2),ijk(3)) + &
bank_array(i) % wgt
end do FISSION_SITES
#ifdef MPI
@ -253,7 +253,7 @@ contains
sites_outside = outside
#endif
end subroutine count_fission_sites
end subroutine count_bank_sites
!===============================================================================
! MESH_INTERSECT determines if a line between xyz0 and xyz1 intersects the outer