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

Ability to sample multiple secondary energy distributions for single reaction. Closes gh-5.

Browse source
This commit is contained in:
Paul Romano 2011-11-11 13:10:41 -05:00
parent d311721e5d
commit 58a5f871c4
3 changed files with 23 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/endf_header.f90
View file

@ -7,7 +7,7 @@ module endf_header
!===============================================================================
type Tab1
integer :: n_regions ! # of interpolation regions
integer :: n_regions = 0 ! # of interpolation regions
integer, allocatable :: nbt(:) ! values separating interpolation regions
integer, allocatable :: int(:) ! interpolation scheme
integer :: n_pairs ! # of pairs of (x,y) values

2
src/interpolation.f90
View file

@ -138,7 +138,7 @@ contains
! determine number of interpolation regions and pairs
n_regions = obj % n_regions
n_pairs = obj % n_pairs
n_pairs = obj % n_pairs
! find which bin the abscissa is in -- if the abscissa is outside the
! tabulated range, the first or last point is chosen, i.e. no interpolation

29
src/physics.f90
View file

@ -1425,7 +1425,7 @@ contains
! SAMPLE_ENERGY
!===============================================================================
subroutine sample_energy(edist, E_in, E_out, mu_out, A, Q)
recursive subroutine sample_energy(edist, E_in, E_out, mu_out, A, Q)
type(DistEnergy), pointer :: edist
real(8), intent(in) :: E_in
@ -1447,6 +1447,8 @@ contains
integer :: ND ! number of discrete lines
integer :: NP ! number of points in distribution
real(8) :: p_valid ! probability of law validity
real(8) :: E_i_1, E_i_K ! endpoints on outgoing grid i
real(8) :: E_i1_1, E_i1_K ! endpoints on outgoing grid i+1
real(8) :: E_1, E_K ! endpoints interpolated between i and i+1
@ -1479,15 +1481,26 @@ contains
real(8) :: x, y, v ! intermediate variables for n-body dist
real(8) :: r1, r2, r3, r4, r5, r6
! TODO: If there are multiple scattering laws, sample scattering law
! ==========================================================================
! SAMPLE ENERGY DISTRIBUTION IF THERE ARE MULTIPLE
! Check for multiple interpolation regions
if (edist % p_valid % n_regions > 0) then
message = "Multiple interpolation regions not supported while &
&attempting to sample secondary energy distribution."
call fatal_error()
if (associated(edist % next)) then
if (edist % p_valid % n_regions > 0) then
p_valid = interpolate_tab1(edist % p_valid, E_in)
if (rang() > p_valid) then
if (edist % law == 44 .or. edist % law == 61) then
call sample_energy(edist%next, E_in, E_out, mu_out)
elseif (edist % law == 66) then
call sample_energy(edist%next, E_in, E_out, A=A, Q=Q)
else
call sample_energy(edist%next, E_in, E_out)
end if
return
end if
end if
end if
! Determine which secondary energy distribution law to use
select case (edist % law)
case (1)