From aa5d093a592f641720671b9e832f799fd78a4262 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Fri, 22 Dec 2017 13:40:24 +0700 Subject: [PATCH] Remove unused components of the Material derived type --- src/cross_section.F90 | 22 ---------------------- src/material_header.F90 | 7 ------- 2 files changed, 29 deletions(-) diff --git a/src/cross_section.F90 b/src/cross_section.F90 index 0bdadd87e..46cbb93ab 100644 --- a/src/cross_section.F90 +++ b/src/cross_section.F90 @@ -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 diff --git a/src/material_header.F90 b/src/material_header.F90 index d77220415..a3ce32db3 100644 --- a/src/material_header.F90 +++ b/src/material_header.F90 @@ -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