crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Cleaned up index_inelastic_scatter a bit to address comments in #963"

Browse source
This commit is contained in:
Adam Nelson 2018-02-13 06:35:50 -05:00
parent 563b165891
commit 718474fb49
1 changed files with 5 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/nuclide_header.F90
View file

@ -303,7 +303,7 @@ contains
real(8) :: temp_actual
type(VectorInt) :: MTs
type(VectorInt) :: temps_to_read
type(VectorInt) :: index_inelastic_scatter_vector
type(VectorInt) :: index_inelastic_scatter
! Get name of nuclide from group
name_len = len(this % name)
@ -478,7 +478,7 @@ contains
MTs % data(i) /= N_3NF .and. MTs % data(i) < 200 .and. &
MTs % data(i) /= N_LEVEL .and. MTs % data(i) /= ELASTIC) then
call index_inelastic_scatter_vector % push_back(i)
call index_inelastic_scatter % push_back(i)
end if
call close_group(rx_group)
@ -486,10 +486,9 @@ contains
call close_group(rxs_group)
! Recast to a regular array to save space
allocate(this % index_inelastic_scatter( &
index_inelastic_scatter_vector % size()))
this % index_inelastic_scatter = index_inelastic_scatter_vector % data(:)
call index_inelastic_scatter_vector % clear()
allocate(this % index_inelastic_scatter(index_inelastic_scatter % size()))
this % index_inelastic_scatter = &
index_inelastic_scatter % data(1: index_inelastic_scatter % size())
! Read unresolved resonance probability tables if present
if (object_exists(group_id, 'urr')) then