small changes to the collision resolver for mergers

MintonGroup · Dec 22, 2022 · d62cec5 · d62cec5
1 parent 0896ffb
commit d62cec5
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/collision/collision_regime.f90 b/src/collision/collision_regime.f90
@@ -36,9 +36,9 @@ module subroutine collision_regime_impactors(self, nbody_system, param)
          case("MERGE")
             impactors%regime = COLLRESOLVE_REGIME_MERGE
             mtot = sum(impactors%mass(:))
+            if (allocated(impactors%mass_dist)) deallocate(impactors%mass_dist)
+            allocate(impactors%mass_dist(1))
             impactors%mass_dist(1) = mtot
-            impactors%mass_dist(2) = 0.0_DP
-            impactors%mass_dist(3) = 0.0_DP
             impactors%rbcom(:) = (impactors%mass(1) * impactors%rb(:,1) + impactors%mass(2) * impactors%rb(:,2)) / mtot
             impactors%vbcom(:) = (impactors%mass(1) * impactors%vb(:,1) + impactors%mass(2) * impactors%vb(:,2)) / mtot
          case default