Added coordinate transform step before resolving collisions to get al…

…l the coordinates consistenent between barycentric and heliocentric before starting the fragment code
MintonGroup · Aug 16, 2021 · f9af5d7 · f9af5d7
1 parent dc4a86a
commit f9af5d7
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90
@@ -1083,6 +1083,8 @@ module subroutine symba_collision_resolve_plplenc(self, system, param, t)
                if (plplcollision_list%nenc == 0) return ! No collisions to resolve
 
                write(*, *) "Collision between massive bodies detected at time t = ", t
+               call pl%vb2vh(system%cb)
+               call pl%h2b(system%cb)
                if (param%lfragmentation) then
                   call plplcollision_list%resolve_fragmentations(system, param)
                else