% loading data
disp('loading data');

dmpsfitdata=[];dmpsdata=[];sizes=[];

for i=1:period(end)-period(1)+1
    if exist([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.fit'])>0;
    tempdata=[];        
    tempdata=load([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.fit']);
    tempdata(:,1)=((tempdata(:,1)-floor(tempdata(:,1)))+period(i));
    tempdata=[(i*ones(1,length(tempdata(:,1))))' tempdata];
    dmpsfitdata=[dmpsfitdata; tempdata];    
    end
    if exist([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.sum'])>0;
    tempdata=[];            
    tempdata=load([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.sum']);
    tempdata(:,1)=(tempdata(:,1)-floor(tempdata(:,1)))+period(i);
    if size(sizes)<1
        % sizes of dmps data
        sizes=tempdata(1,3:end);
    end
    tempdata=[(i*ones(1,length(tempdata(:,1))-1))' tempdata(2:end,:)];
    dmpsdata=[dmpsdata; tempdata];
    end 
end
   