function day=eventtoday(event);
ev=event;
index=isnan(ev);
ev(index)=1000000*ones(size(ev(index)));
mo=floor(ev/1.0e6);
da=floor(ev/1.0e4-mo*100);
hr=floor(ev/100-mo*1.0e4-da*100);
mm=ev-mo*1e6-da*1e4-hr*100;
firstday=[0,31,60,91,121,152,182,213,244,274,305,335]; %leapyear
day=firstday(mo)+da+hr/24+mm/24/60;
day(index)=event(index);
return