function out=zontide(mjd,par,n_harm);
%Вычисление редукций за влиЯние зональных приливов в UT и LOD
%
%Вызов:  out=zontide(mjd,par,n_harm);
%
% mjd    - массив модифицированных юлианских дат (N);
% par    - параметр, указывающий длЯ чего вычислЯть:
%              1 - длЯ UT, 2 - длЯ LOD;
% n_harm - массив (2x1), задающий номера первой (>=1) и последней (<=62)
%              гармоник, длЯ которых вычислЯютсЯ редукции;
% out    - массив (Nx1), содержащий редукции:
%              dUTR - в сек, dLODR - в мс;
%
%ЗначениЯ коэффициентов и волновых чисел длЯ аргументов содержатсЯ
%в файле ztide.mat
% 1994 IERS Annual Report, Tabl.I-5a, pp. I16 - I17.
% 1000 точек - около 25 с

% Хомутов С.Ю., 19.01.99

load ztide    % файл с волновыми числами аргумента и коэффициентами
              % 10E-4 s for UT (SIN), 10E-5 s for LOD (COS)

amp=ztide(n_harm(1):n_harm(2),7+par);   % коэф-ты длЯ заданного параметра
arg0=ztide(n_harm(1):n_harm(2),2:6);

for k=1:length(mjd),
    if rem(k,500)==0,
       disp(k)
    end
    T=mjd(k)-51544.5;
    EPH(1) = d2r( 134.96 + 13.064993*T );
    EPH(2) = d2r( 357.53 +  0.985600*T );
    EPH(3) = d2r(  93.27 + 13.229350*T );
    EPH(4) = d2r( 297.85 + 12.190749*T );
    EPH(5) = d2r( 125.04 -  0.052954*T );

    if par==1,
       out(k)=sum( amp .* sin( arg0(:,1)*EPH(1)+arg0(:,2)*EPH(2)+...
                               arg0(:,3)*EPH(3)+arg0(:,4)*EPH(4)+...
                               arg0(:,5)*EPH(5) )) / 10000;
    else
       out(k)=sum( amp .* cos( arg0(:,1)*EPH(1)+arg0(:,2)*EPH(2)+...
                               arg0(:,3)*EPH(3)+arg0(:,4)*EPH(4)+...
                               arg0(:,5)*EPH(5) )) / 100;
    end
end

out=out(:);