function ang = ind2ang(i,hz,p)
%what angle at hz hertz are we at index i and phase p?
global N Fs;
wl = hz2wlind(hz);
cyc = i/wl;
ang = (cyc-fix(cyc))*2*pi + p;
end

