function [azimuth, elevation] = unit_vector_to_azimuth_elevation( uhat )
    

azimuth   = atan2( uhat(2), uhat(1) );
% elevation = asin( uhat(3) );
elevation = asin( -uhat(3) / sqrt( uhat(3).^2 + uhat(2).^2 + uhat(1).^2) );
