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

Changed tallies.out to use mesh indices rather than bins. Changed get_uid subroutine to get_label.

Browse source
This commit is contained in:
Paul Romano 2011-09-29 22:37:10 -04:00
parent 56a425ca9a
commit 7c8196b6a0
1 changed files with 25 additions and 19 deletions
Download patch file Download diff file Expand all files Collapse all files

44
src/tally.f90
View file

@ -4,7 +4,7 @@ module tally
use cross_section, only: get_macro_xs
use error, only: fatal_error
use global
use mesh, only: get_mesh_bin
use mesh, only: get_mesh_bin, bin_to_mesh_indices
use output, only: message
use search, only: binary_search
use string, only: int_to_str, real_to_str
@ -610,7 +610,7 @@ contains
if (has_filter(j)) then
! Print current filter information
write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') space(1:indent), &
trim(filter_name(j)), trim(get_uid(t, j, bins(j)))
trim(filter_name(j)), trim(get_label(t, j, bins(j)))
indent = indent + 2
end if
j = j + 1
@ -621,7 +621,7 @@ contains
! Print filter information
write(UNIT=UNIT_TALLY, FMT='(1X,2A,1X,A)') space(1:indent), &
trim(filter_name(j)), trim(get_uid(t, j, bins(j)))
trim(filter_name(j)), trim(get_label(t, j, bins(j)))
! Determine scoring index for this bin combination -- note that unlike
! in the score_tally subroutine, we have to use max(bins,1) since all
@ -648,50 +648,56 @@ contains
end subroutine write_tallies
!===============================================================================
! GET_UID returns the uid for a cell/surface/etc given a tally, filter type, and
! corresponding bin
! GET_LABEL returns a label for a cell/surface/etc given a tally, filter type,
! and corresponding bin
!===============================================================================
function get_uid(t, filter_type, bin) result(uid)
function get_label(t, filter_type, bin) result(label)
type(TallyObject), pointer :: t ! tally object
integer, intent(in) :: filter_type ! type of filter
integer, intent(in) :: bin ! bin in filter array
character(30) :: uid ! user-specified identifier
character(30) :: label ! user-specified identifier
integer :: index ! index in cells/surfaces/etc array
real(8) :: E0 ! lower bound for energy bin
real(8) :: E1 ! upper bound for energy bin
integer :: index ! index in cells/surfaces/etc array
integer, allocatable :: ijk(:) ! indices in mesh
real(8) :: E0 ! lower bound for energy bin
real(8) :: E1 ! upper bound for energy bin
type(StructuredMesh), pointer :: m
select case(filter_type)
case (T_UNIVERSE)
index = t % universe_bins(bin) % scalar
uid = int_to_str(universes(index) % uid)
label = int_to_str(universes(index) % uid)
case (T_MATERIAL)
index = t % material_bins(bin) % scalar
uid = int_to_str(materials(index) % uid)
label = int_to_str(materials(index) % uid)
case (T_CELL)
index = t % cell_bins(bin) % scalar
uid = int_to_str(cells(index) % uid)
label = int_to_str(cells(index) % uid)
case (T_CELLBORN)
index = t % cellborn_bins(bin) % scalar
uid = int_to_str(cells(index) % uid)
label = int_to_str(cells(index) % uid)
case (T_SURFACE)
index = t % surface_bins(bin) % scalar
uid = int_to_str(surfaces(index) % uid)
label = int_to_str(surfaces(index) % uid)
case (T_MESH)
uid = "Bin " // int_to_str(bin)
m => meshes(t % mesh)
allocate(ijk(m % n_dimension))
call bin_to_mesh_indices(m, bin, ijk)
label = "Index (" // trim(int_to_str(ijk(1))) // ", " // &
trim(int_to_str(ijk(2))) // ", " // trim(int_to_str(ijk(3))) // ")"
case (T_ENERGYIN)
E0 = t % energy_in(bin)
E1 = t % energy_in(bin + 1)
uid = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
label = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
case (T_ENERGYOUT)
E0 = t % energy_out(bin)
E1 = t % energy_out(bin + 1)
uid = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
label = "[" // trim(real_to_str(E0)) // ", " // trim(real_to_str(E1)) // ")"
end select
end function get_uid
end function get_label
!===============================================================================
! TALLY_STATISTICS computes the mean and standard deviation of the mean of each