Adjusted initial distribution of angular momentum value to help KE converge better

MintonGroup · May 17, 2021 · cc18f2d · cc18f2d
1 parent e945044
commit cc18f2d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
@@ -506,7 +506,7 @@ subroutine set_fragment_tangential_velocities()
             call util_crossproduct(v_r_unit(:, i), v_t_unit(:, i), L(:))
             A(4:6, i) = m_frag(i) * rmag(i) * L(:)
          else !For the remining bodies, distribute the angular momentum equally amongs them
-            v_t_mag(i) = 0.25_DP * L_orb_mag / (m_frag(i) * rmag(i) * nfrag)
+            v_t_mag(i) = 0.9_DP * L_orb_mag / (m_frag(i) * rmag(i) * nfrag)
             v_frag(:, i) = v_t_mag(i) * v_t_unit(:, i)
             L_lin_others(:) = L_lin_others(:) + m_frag(i) * v_frag(:, i)
             call util_crossproduct(x_frag(:, i), v_frag(:, i), L(:))