function [FI,TETA,PSI]=lemi_angles_rot_select(current_date)
% Выбор принЯтых на дату MJD углов вращениЯ длЯ LEMI-008
%
%       [FI,TETA,PSI]=lemi_angles_rot_select(current_date);
%
% current_date = mjd или Tser (сериальная дата)
%                mjd - модифицир. юлиантскаЯ дата, длЯ котрой требуютсЯ углы Эйлера
% FI,TETA,PSI - углы на дату MJD (Tser) (в градусах)
%
% !!! Углы выбираютсЯ из файла "LEMI_angles_rot.dat"
%     Дата в файле соответствует дате, с которой (вкл.) принЯты эти углы

% 13.03.2005    Khomutov S.
% 03.06.2008 - дополнено вариантом входной даты как сериальной (кроме mjd)

path_all='D:\Projects\Novosibirsk_new\'; % общий каталог для обсерватории

if ~exist([path_all 'Ini_parameters\LEMI_angles_rot.dat'],'file')
    FI=NaN; TETA=NaN; PSI=NaN;
    disp('!!! Не могу найти файл "LEMI_angles_rot.dat"')
    return
end

load([path_all 'Ini_parameters\LEMI_angles_rot.dat']);
year=LEMI_angles_rot(:,1);
mon=LEMI_angles_rot(:,2);
day=LEMI_angles_rot(:,3);
mjd_rot=cdtmjd(year,mon,day);
T_rot=datenum(year,mon,day);

if current_date>600000  % сериальная дата
    k=max(find(current_date>=T_rot));
else
    k=max(find(current_date>=mjd_rot));
end

FI  =LEMI_angles_rot(k,4);
TETA=LEMI_angles_rot(k,5);
PSI =LEMI_angles_rot(k,6);
