From 4db72c90af2dd62d52a199f4e86ef26676abab53 Mon Sep 17 00:00:00 2001 From: Will Boyd Date: Sun, 8 Nov 2015 20:53:18 -0500 Subject: [PATCH] Fixed issues in iso-in-lab scattering in physics.F90 per comments by @walshjon --- src/physics.F90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/physics.F90 b/src/physics.F90 index 00db08950c..9f8fb535b1 100644 --- a/src/physics.F90 +++ b/src/physics.F90 @@ -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