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

Removed reversal of energy bins in multigroup mode to agree with energy bin search in source.

Browse source
This commit is contained in:
Brody Bassett 2017-12-09 18:04:02 -05:00
parent 25ce61622d
commit 53434692fa
2 changed files with 1 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/input_xml.F90
View file

@ -4172,7 +4172,6 @@ contains
integer :: n_libraries
logical :: file_exists ! does mgxs.h5 exist?
integer(HID_T) :: file_id
real(8), allocatable :: rev_energy_bins(:)
character(len=MAX_WORD_LEN), allocatable :: names(:)
! Check if MGXS Library exists
@ -4202,7 +4201,6 @@ contains
num_delayed_groups = 0
end if
allocate(rev_energy_bins(num_energy_groups + 1))
allocate(energy_bins(num_energy_groups + 1))
if (attribute_exists(file_id, "group structure")) then
@ -4212,9 +4210,6 @@ contains
call fatal_error("'group structure' attribute must exist!")
end if
! First reverse the order of energy_groups
energy_bins = energy_bins(num_energy_groups + 1:1:-1)
allocate(energy_bin_avg(num_energy_groups))
do i = 1, num_energy_groups
energy_bin_avg(i) = HALF * (energy_bins(i) + energy_bins(i + 1))

2
src/tallies/tally_filter_energy.F90
View file

@ -75,7 +75,7 @@ contains
! ordering of the library and tallying systems).
if (.not. run_CE) then
if (n == num_energy_groups + 1) then
if (all(this % bins == energy_bins(num_energy_groups + 1:1:-1))) &
if (all(this % bins == energy_bins)) &
then
this % matches_transport_groups = .true.
end if