From b9d815ee484d635eb1584489857083e61953408f Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 16 Aug 2021 19:14:40 -0400 Subject: [PATCH] Fixed energy term in the merger case --- src/symba/symba_collision.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90 index 2e4d50b58..05ac54fd4 100644 --- a/src/symba/symba_collision.f90 +++ b/src/symba/symba_collision.f90 @@ -277,7 +277,7 @@ module function symba_collision_casemerge(system, param, family, x, v, mass, rad pe = 0.0_DP do j = 1, nfamily do i = j + 1, nfamily - pe = pe - pl%mass(i) * pl%mass(j) / norm2(pl%xb(:, i) - pl%xb(:, j)) + pe = pe - pl%Gmass(i) * pl%mass(j) / norm2(pl%xb(:, i) - pl%xb(:, j)) end do end do system%Ecollisions = system%Ecollisions + pe