diff --git a/examples/Fragmentation/Fragmentation_Movie.py b/examples/Fragmentation/Fragmentation_Movie.py index cf99bf67d..0917abef4 100755 --- a/examples/Fragmentation/Fragmentation_Movie.py +++ b/examples/Fragmentation/Fragmentation_Movie.py @@ -59,9 +59,9 @@ } vel_vectors = {"disruption_headon" : [np.array([ 0.00, 6.280005, 0.0]), - np.array([ 0.00, 4.28, 0.0])], + np.array([ 0.00, 3.28, 0.0])], "disruption_off_axis" : [np.array([ 0.00, 6.280005, 0.0]), - np.array([ 0.05, 4.28, 0.0])], + np.array([ 0.05, 3.28, 0.0])], "supercatastrophic_headon": [np.array([ 0.00, 6.28, 0.0]), np.array([ 0.00, 4.28, 0.0])], "supercatastrophic_off_axis": [np.array([ 0.00, 6.28, 0.0]), diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90 index bc5332138..6c63d446b 100644 --- a/src/fraggle/fraggle_generate.f90 +++ b/src/fraggle/fraggle_generate.f90 @@ -717,7 +717,7 @@ module subroutine fraggle_generate_vel_vec(collider, nbody_system, param, lfailu ! Remove a constant amount of velocity from the bodies so we don't shift the center of mass and screw up the momentum ke_avail = 0.0_DP do i = fragments%nbody, 1, -1 - ke_avail = ke_avail + 0.5_DP * fragments%mass(i) * max(fragments%vmag(i) - vesc,0.0_DP)**2 + ke_avail = ke_avail + 0.5_DP * fragments%mass(i) * max(fragments%vmag(i) - vesc/try,0.0_DP)**2 end do ke_remove = min(E_residual, ke_avail)