crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Check for neighboring energy points that are the same on nuclide energy grid.

Browse source
This commit is contained in:
Paul Romano 2011-12-18 17:59:29 -05:00
parent f31c7e5f0a
commit ffdcead9e4
1 changed files with 6 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

7
src/physics.F90
View file

@ -235,8 +235,13 @@ contains
! TODO: If not using unionized energy grid, we need to find the index on the
! nuclide energy grid using lethargy mapping or whatever other technique
! search nuclide energy grid
! get index on nuclide energy grid
IE = nuc % grid_index(p % IE)
! check for rare case where two energy points are the same
if (nuc % energy(IE) == nuc % energy(IE+1)) IE = IE + 1
! calculate interpolation factor
f = (p%E - nuc%energy(IE))/(nuc%energy(IE+1) - nuc%energy(IE))
micro_xs(index_nuclide) % index_grid = IE