function uvect = alpdelt2uvect( alpha, delta )
% function uvect = alpdelt2uvect( alpha, delta )
% convert alpha delta to unit vector
%        input  alpha : longitude [deg]
%        input  delta : latitude [deg]
%        output  uvect[3] : derection unit vector 
%        functions used   sin(), cos()                      
%
%     Copyright Tatsuaki Hashimoto  2008.12.5
%
%

PI = 3.141592653589;
RAD = PI/180;

uvect(1) = cos(alpha*RAD) * cos(delta*RAD);
uvect(2) = sin(alpha*RAD) * cos(delta*RAD);
uvect(3) = sin(delta*RAD);

check = uvect(1)^2+uvect(2)^2+uvect(3)^2 ;
if abs(check-1.0) > 1e-6
   fprintf('unit vector error %e \n', check)
end

