fixed some indices; added computeNewVelocity(..., MultiIndexType& noLD) (TB, AH)