CVode fix; force calculations fixed