function mjd_axis_code(mjd_step,code);
% Разметка временной оси по NJD с шагом dmjd суток
%
% Вызов: mjd_axis_code(dmjd,code)
%
% mjd_step - шаг маркировки длЯ ticks
% code     - тип маркировки оси:
%       1 - шаг в сутках от начала оси
%       2 - шаг в месЯцах (напр., "1" - 1-е число месЯца, 
%           "2" - 1-е число каждого второго месЯца)

% 06.01.2004 - Khomutov

mon_str=['Jan';'Feb';'Mar';'Apr';'May';'Jun';
         'Jul';'Aug';'Sep';'Oct';'Nov';'Dec'];

% пределы оси
tmp=axis; mjd1=tmp(1); mjd2=tmp(2);

if code==1, % разметка оси через dmjd суток
    mjd_tick=(mjd1:mjd_step:mjd2)';
    [year,mon,day]=mjdcdt(mjd_tick);
    if year(1)<1900,
        error('Error in MJD (year<1900)')
    end
    mjd_label=[];
    for k=1:length(mjd_tick),
        mjd_label=strvcat(mjd_label,sprintf('%02d/%02d',mon(k),day(k)));
    end
    set(gca,'Xlim',[mjd1 mjd2],'Xtick',mjd_tick,'Xticklabel',mjd_label)
elseif code==2,
    tmp=(mjd1:mjd2)';
    [year,mon,day]=mjdcdt(tmp);
    mjd_tmp=tmp(day==1);
    mjd_tick=mjd_tmp(1:mjd_step:length(mjd_tmp));
    [year,mon,day]=mjdcdt(mjd_tick);
    if year(1)<1900,
        error('Error in MJD (year<1900)')
    end
    mjd_label=[];
    for k=1:length(mjd_tick),
        mjd_label=strvcat(mjd_label,sprintf('%s',mon_str(mon(k),:)));
    end
    set(gca,'Xlim',[mjd1 mjd2],'Xtick',mjd_tick,'Xticklabel',mjd_label)
end
