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