% routine to label the yearday axis with mth/day
xtl=get(gca,'xticklabels');
[nr,mn]=size(xtl);
nl=floor(nr/3);

% intitialize new xticklabels
for k=1:nr;
  nxtl(k,:)='     ';
end;
ln=0;

for k=1:nl:nr;
  ln=ln+1;
  [date]=calday(95,str2num(xtl(k,:)));

% take care of the four combinations of string lengths
  if date(2)<10;
    if date(3)<10;
      nxtl(k,:)=[' ',num2str(date(2)),'/0',num2str(date(3))];
    else;
      nxtl(k,:)=[' ',num2str(date(2)),'/',num2str(date(3))];
    end;
  elseif date(2)>=10;
    if date(3)<10;
      nxtl(k,:)=[num2str(date(2)),'/0',num2str(date(3))];
    else;
      nxtl(k,:)=[num2str(date(2)),'/',num2str(date(3))];
    end;
  end;
end;
set(gca,'xticklabels',nxtl);