From 25081276a3642e192b243ed061f588c18657cfcd Mon Sep 17 00:00:00 2001 From: David Minton Date: Thu, 29 Aug 2019 09:33:19 -0400 Subject: [PATCH] Updated ray model to be the narrow ray model from Minton et al. (2019) --- src/ejecta/ejecta_ray_pattern.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_ray_pattern.f90 b/src/ejecta/ejecta_ray_pattern.f90 index 2e15dc73..32982add 100644 --- a/src/ejecta/ejecta_ray_pattern.f90 +++ b/src/ejecta/ejecta_ray_pattern.f90 @@ -186,7 +186,7 @@ pure function pattern(theta,r,rmin,rmax,thetari,ej) result(ans) ans = 0.0_DP end if else - rw0 = rmin * pi / Nraymax / 1 + rw0 = rmin * pi / Nraymax / 2 rw1 = 2 * pi / Nraymax rw = rw0 * (1._DP - (1.0_DP - rw1 / rw0) * exp(1._DP - (r / rmin)**2)) n = max(min(floor((Nraymax**rayp - (Nraymax**rayp - 1) * log(r/minray) / log(rray/minray))**(1._DP/rayp)),Nraymax),1) ! Exponential decay of ray number with distance