From 9bbe6a1ead1446930e20b0b7b245f9393ddd48a1 Mon Sep 17 00:00:00 2001 From: David Minton Date: Fri, 21 May 2021 23:42:03 -0400 Subject: [PATCH] Adjusted tolerance value --- src/symba/symba_frag_pos.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index a79542939..b7bd75334 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -113,7 +113,7 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi logical, intent(out) :: lmerge ! Answers the question: Should this have been a merger instead? real(DP), intent(inout) :: Qloss ! Internals - real(DP), parameter :: f_spin = 0.20_DP !! Fraction of pre-impact orbital angular momentum that is converted to fragment spin + real(DP), parameter :: f_spin = 0.03_DP !! Fraction of pre-impact orbital angular momentum that is converted to fragment spin real(DP) :: mscale = 1.0_DP, rscale = 1.0_DP, vscale = 1.0_DP, tscale = 1.0_DP, Lscale = 1.0_DP, Escale = 1.0_DP! Scale factors that reduce quantities to O(~1) in the collisional system real(DP) :: mtot real(DP), dimension(NDIM) :: xcom, vcom @@ -660,7 +660,7 @@ subroutine set_fragment_radial_velocities(lmerge) ! Arguments logical, intent(out) :: lmerge ! Internals - real(DP), parameter :: TOL = 1e-6_DP + real(DP), parameter :: TOL = 1e-10_DP real(DP), dimension(:), allocatable :: vflat logical :: lerr integer(I4B) :: i