diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index cda9d3f4b..3e3bca5d4 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -626,8 +626,7 @@ function ke_objective_function(v_r_mag, v_r_unit, v_t_mag, v_t_unit, x_frag, m_f v_shift(:, i) = v_shift(:, i) + v_t_mag(i) * v_t_unit(:, i) + vcom(:) fval = fval + 0.5_DP * m_frag(i) * dot_product(v_shift(:, i), v_shift(:, i)) end do - !write(*,*) 'fval: ',fval - fval = fval**2 + fval = fval**4 return