function [azim,z_dist]=equ2hor(t_angle,delta,lat)
%Перевод экваториальных координат (t_angle,delta) 
%    в горизонтальные (a,z)
%
%Call:	[azim,z_dist]=equ2hor(t_angle,delta,lat);
%
%Вход: 	t_angle	- часовой угол (в часах);
%	  delta - склонение (в градусах);
%	    lat - широта места наблюдени?, в град. (положительна? к северу);
%
%Выход:  azim - геодезический азимут в град. (от N через E) ;
%	    z_dist - зенитное рассто?ние в град. (от 0 до 180);
%
%Внимание !!!  Все углы д.б. в град. иди часах

% Sergey Khomutov, 08.05.97,  Novosibirsk
%     06.06.2002 - перевод входных величин в градусы/часы

% промежуточные параметры - тригонометрические функции
st=sin(d2r(t_angle*15)); ct=cos(d2r(t_angle*15));
sd=sin(d2r(delta)); cd=cos(d2r(delta));
sl=sin(d2r(lat)); cl=cos(d2r(lat));

% Вычисление азимута
azim = atan2(-cd.*st, sd.*cl-cd.*sl.*ct);
k=azim<0; azim(k)=azim(k)+2*pi;		% transform to (0...2*pi);

% Вычисление зенитного рассто?ни?
z_dist = atan2(-cd.*st, (sd.*sl+cd.*cl.*ct).*sin(azim));
k=z_dist<0; z_dist(k)=z_dist(k)+pi;	% transform to (0...pi);

azim=r2d(azim);
z_dist=r2d(z_dist);
