function ind = ang2ind(a,hz)
% what index < wavelength is angle a (<2pi)?
global N Fs;
wl = hztowlind(hz);
ind = (a/(2*pi))*wl;
end