correction: Jdot = J*DIV_VEL not just DIV_VEL