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

Fixed sample_angles interpolation of mu. The old way was resulting in mus outside the range of [-1,1] and was effectively using the slope of the randomly selected bin but extrapolating to the next lower bin.

Browse source
This commit is contained in:
Adam Nelson 2014-05-19 12:19:14 -04:00
parent e22b452b5f
commit d1ca3519f2
1 changed files with 1 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/physics.F90
View file

@ -1164,7 +1164,7 @@ contains
! calculate cosine
mu0 = rxn % adist % data(lc + k)
mu1 = rxn % adist % data(lc + k+1)
mu = mu0 + (32.0_8 * xi - k) * (mu1 - mu0)
mu = mu0 + (32.0_8 * xi - k + ONE) * (mu1 - mu0)
elseif (type == ANGLE_TABULAR) then
interp = int(rxn % adist % data(lc + 1))