From 9b1e72e7bc4176d0898db90a0887ca138356449b Mon Sep 17 00:00:00 2001
From: Martin Diehl <m.diehl@mpie.de>
Date: Thu, 15 Feb 2018 11:26:19 +0100
Subject: [PATCH] more sensible cutoff for scatter

---
 src/math.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/math.f90 b/src/math.f90
index e88a5ed03..3ede28beb 100644
--- a/src/math.f90
+++ b/src/math.f90
@@ -1783,7 +1783,7 @@ function math_sampleGaussOri(center,noise)
  real(pReal), dimension(3), parameter :: ORIGIN = 0.0_pReal
  real(pReal), dimension(5) :: rnd
 
- noScatter: if (abs(noise) < tol_math_check) then
+ noScatter: if (noise < 0.5_pReal*INRAD) then
    math_sampleGaussOri = center
  else noScatter
   ! Helming uses different distribution with Bessel functions