Moved G out from sqrt

MintonGroup · Aug 6, 2021 · a72bdc5 · a72bdc5
1 parent 6706e16
commit a72bdc5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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