function rotMat   = calRotMat(pos);
%---------------------------------------------------------
% calculate rotation matrix for given position
% Ju-Yong Do,  Feb 2, 2009
%---------------------------------------------------------
% input  : pos
% output : rotMat
%=========================================================
%=========================================================
%rotation matrix for coordinate change
lla             = wgsxyz2lla(pos);
rotMat          = rot(90-lla(1), 1)*rot(90+lla(2), 3);

return;
%=========================================================
%=========================================================