function t_month=f_month(mjd_day,month,year)
%Вычисление даты в долях месяца по модифицированной юлианской 
%или календарной дате
%
%Вызов:     t_month=f_month(mjd_day,month,year)
%
% mjd_day - модифицированная юлианская дата (другие аргументы
%           должны отсутствовать) или день месяца (обязательно
%           должны быть аргументы month и year)
%   month - месяц
%    year - год
% t_month - дата в долях месяца.
%
% 1-е число каждого месяца имеет дату 0, т.е. 1 марта -> 3.0

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


if nargin==1   % вычисления по MJD
   [year,month,day]=mjdcdt(mjd_day);
else           % вычисления по календарной дате
   if nargin~=3
      error ('!!! Мало аргументов !!!')
   end
end

t_month=zeros(length(mjd_day),1);
for km=1:12,
   k=month==km;
   if sum(k)>0,   % есть данные для этого месяца
      t_month(k)=month(k)+(day(k)-1)/max(day(k));
   end
end
