From a72bdc548c22ca4f3cb287f5090e1751226e6fdc Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 6 Aug 2021 18:30:00 -0400 Subject: [PATCH] Moved G out from sqrt --- src/fragmentation/fragmentation.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/fragmentation/fragmentation.f90 b/src/fragmentation/fragmentation.f90 index e1b69dd44..270f48a7a 100644 --- a/src/fragmentation/fragmentation.f90 +++ b/src/fragmentation/fragmentation.f90 @@ -167,7 +167,7 @@ subroutine set_scale_factors() ! Set scale factors Escale = 0.5_DP * (mass(1) * dot_product(v(:,1), v(:,1)) + mass(2) * dot_product(v(:,2), v(:,2))) rscale = sum(radius(:)) - mscale = sqrt(Escale * rscale / param%GU) + mscale = sqrt(Escale * rscale) / param%GU vscale = sqrt(Escale / mscale) tscale = rscale / vscale Lscale = mscale * rscale * vscale