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

Moved k_batch back to tally module.

Browse source
This commit is contained in:
Paul Romano 2012-12-04 19:46:04 -05:00
parent 31ba238f89
commit dd87f7ff68
2 changed files with 10 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/eigenvalue.F90
View file

@ -161,10 +161,6 @@ contains
integer :: i ! loop index for state point batches
! Before accumulating results for global_tallies, we need to get the
! current batch estimate of k_analog for displaying to output
k_batch(current_batch) = global_tallies(K_ANALOG) % value
! Collect tallies
if (tallies_on) then
call timer_start(time_tallies)
@ -549,6 +545,8 @@ contains
! =========================================================================
! SINGLE-BATCH ESTIMATE OF K-EFFECTIVE
if (.not. active_batches) k_batch(current_batch) = global_tallies(K_ANALOG) % value
#ifdef MPI
if ((.not. active_batches) .or. (.not. reduce_tallies)) then
! Reduce value of k_batch if running in parallel

8
src/tally.F90
View file

@ -1506,6 +1506,14 @@ contains
curr_ptr => curr_ptr % next
end do
if (run_mode == MODE_EIGENVALUE) then
! Get the current batch estimate of k_analog for displaying to output
! --- this has to be performed after reduce_tally_values and before
! accumulate_result
k_batch(current_batch) = global_tallies(K_ANALOG) % value
end if
! Accumulate results for global tallies
if (active_batches) call accumulate_result(global_tallies)
end if