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

Remove unused components of the Material derived type

Browse source
This commit is contained in:
Paul Romano 2017-12-22 13:40:24 +07:00
parent b10b7204b2
commit aa5d093a59
2 changed files with 0 additions and 29 deletions
Download patch file Download diff file Expand all files Collapse all files

22
src/cross_section.F90
View file

@ -573,28 +573,6 @@ contains
end subroutine calculate_urr_xs
!===============================================================================
! FIND_ENERGY_INDEX determines the index on the union energy grid at a certain
! energy
!===============================================================================
pure function find_energy_index(mat, E) result(i)
type(Material), intent(in) :: mat ! pointer to current material
real(8), intent(in) :: E ! energy of particle
integer :: i ! energy grid index
! if the energy is outside of energy grid range, set to first or last
! index. Otherwise, do a binary search through the union energy grid.
if (E <= mat % e_grid(1)) then
i = 1
elseif (E > mat % e_grid(mat % n_grid)) then
i = mat % n_grid - 1
else
i = binary_search(mat % e_grid, mat % n_grid, E)
end if
end function find_energy_index
!===============================================================================
! MULTIPOLE_EVAL evaluates the windowed multipole equations for cross
! sections in the resolved resonance regions

7
src/material_header.F90
View file

@ -42,13 +42,6 @@ module material_header
! it is not present in the material, the entry is set to zero.
integer, allocatable :: mat_nuclide_index(:)
! Energy grid information
integer :: n_grid ! # of union material grid points
real(8), allocatable :: e_grid(:) ! union material grid energies
! Unionized energy grid information
integer, allocatable :: nuclide_grid_index(:,:) ! nuclide e_grid pointers
! S(a,b) data
integer :: n_sab = 0 ! number of S(a,b) tables
integer, allocatable :: i_sab_nuclides(:) ! index of corresponding nuclide