---

idx = find(imag(Wjw(1:end-1)).*imag(Wjw(2:end)) < 0, 1);
w180 = w(idx)
P180 = real(Wjw(idx))
k1kr = k1/abs(P180)