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

Fixed issues in iso-in-lab scattering in physics.F90 per comments by @walshjon

Browse source
This commit is contained in:
Will Boyd 2015-11-08 20:53:18 -05:00
parent 4b70c4abc1
commit 4db72c90af
1 changed files with 2 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/physics.F90
View file

@ -316,7 +316,6 @@ contains
type(Reaction), pointer :: rxn
real(8) :: uvw_new(3) ! outgoing uvw for iso-in-lab scattering
real(8) :: uvw_old(3) ! incoming uvw for iso-in-lab scattering
real(8) :: mu_lab ! polar angle cosine for iso-in-lab scattering
real(8) :: phi ! azimuthal angle for iso-in-lab scattering
! copy incoming direction
@ -407,10 +406,10 @@ contains
phi = TWO * PI * prn()
uvw_new(2) = cos(phi) * sqrt(ONE - uvw_new(1)*uvw_new(1))
uvw_new(3) = sin(phi) * sqrt(ONE - uvw_new(1)*uvw_new(1))
mu_lab = dot_product(uvw_old, uvw_new)
p % mu = dot_product(uvw_old, uvw_new)
! change direction of particle
p % coord(1) % uvw = rotate_angle(uvw_new, mu_lab)
p % coord(1) % uvw = uvw_new
end if
end subroutine scatter