if (nd==2 && i>=itcross(1) && i<=itcross(3) && rem(i-itcross(1),itcross(2))==1)

[xy,nCross] = fixEdgeCrossings(xy,adjAll,ecthresh);

if (plt(3)==2) axis auto; axis equal; a2 = axis; end

end

[maxVel,imax] = max(abs(velocity));

[maxVel,jmax] = max(maxVel);

imax = imax(jmax);