Updated ray model to be the narrow ray model from Minton et al. (2019)

MintonGroup · Aug 29, 2019 · 2508127 · 2508127
1 parent d58fb91
commit 2508127
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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