function val = testPhase(rot, zci)
% find the "annoyance of fit" of rotation rot
out = mod(zci+rot, 2*pi);
out1 = out(out<pi);
out2 = 2*pi - out(out>=pi);
clf
plot(out1, 'rx');
hold on
plot(out2, 'b+');