From 24230347c763f01e028d0da7854ea81fe2b27841 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Thu, 23 Feb 2023 14:24:24 -0500 Subject: [PATCH] Improved the bracketing method for SFD slope --- src/fraggle/fraggle_util.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/fraggle/fraggle_util.f90 b/src/fraggle/fraggle_util.f90 index 23c0a1a21..3b7ee4fe2 100644 --- a/src/fraggle/fraggle_util.f90 +++ b/src/fraggle/fraggle_util.f90 @@ -182,8 +182,8 @@ module subroutine fraggle_util_set_mass_dist(self, param) y1 = y1 - (i-1)**(-x1/3.0_DP) end do if (y0*y1 < 0.0_DP) exit - x1 = x1 * 1.1_DP - x0 = x0 * 1.1_DP + if (y0 > 0.0_DP) x0 = x0 * 1.6_DP + if (y1 < 0.0_DP) x1 = x1 * 1.6_DP end do ! Find the mass scaling factor with bisection