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

Fixed reading of elastic cross sections on S(a,b) tables and added attributes for thresholds.

Browse source
This commit is contained in:
Paul Romano 2011-10-27 12:09:21 -04:00
parent cab8cc629d
commit 7c3c8ee0f9
2 changed files with 18 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/cross_section.f90
View file

@ -1207,6 +1207,9 @@ contains
table % inelastic_e_in = get_real(NE_in)
table % inelastic_sigma = get_real(NE_in)
! set threshold value
table % threshold_inelastic = table % inelastic_e_in(NE_in)
! allocate space for outgoing energy/angle for inelastic
! scattering
NE_out = NXS(4)
@ -1239,9 +1242,18 @@ contains
allocate(table % elastic_e_in(NE_in))
allocate(table % elastic_P(NE_in))
! read elastic energies and P
XSS_index = JXS4 + 1
table % elastic_e_in = get_real(NE_in)
table % elastic_P = get_real(NE_in)
! set threshold
table % threshold_elastic = table % elastic_e_in(NE_in)
! determine whether sigma=P or sigma = P/E
table % n_elastic_type = NXS(5)
table % elastic_type = NXS(5)
else
table % threshold_elastic = ZERO
table % n_elastic_e_in = 0
end if

6
src/cross_section_header.f90
View file

@ -123,6 +123,10 @@ module cross_section_header
real(8) :: awr
real(8) :: temp
! threshold for S(a,b) treatment (usually ~4 eV)
real(8) :: threshold_inelastic
real(8) :: threshold_elastic = 0.0
! Inelastic scattering data
integer :: n_inelastic_e_in
integer :: n_inelastic_e_out
@ -133,8 +137,8 @@ module cross_section_header
real(8), allocatable :: inelastic_mu(:,:,:)
! Elastic scattering data
integer :: elastic_type
integer :: n_elastic_e_in
integer :: n_elastic_type
integer :: n_elastic_mu
real(8), allocatable :: elastic_e_in(:)
real(8), allocatable :: elastic_P(:)