% loading data (+2 days for growth continuying for another 2 days maximum)
disp('loading data');
dmpsfitdata=[];dmpsdata=[];MNCdata=[];sizes=[];
for i=1:period(end)-period(1)+1
    disp(['Adding data:', datestr(period(i))])
    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']);
        varkoll=tempdata; %min koll för att se hur datum och tid fungerar
        tempdata(:,1)=(tempdata(:,1)-floor(tempdata(:,1)))+period(i);  %to transform the date and time to Matalab time
        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 
    if exist([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.MNC'])>0;
        tempdata=[];
        tempdata=[i load([stationname,datestr(period(i),'YY'),datestr(period(i),'mm'),datestr(period(i),'dd'),'.MNC'])];
        MNCdata=[MNCdata; tempdata];
        disp(['Loading data from : ',datestr(period(i),'yymmdd')]);
    end
end
%Add the next two days outside the chosen time period to be able to plot
%three days (for the last day in the time period)
if exist([stationname,datestr(period(i)+1,'YY'),datestr(period(i)+1,'mm'),datestr(period(i)+1,'dd'),'.fit'])>0;
    tempdata=[];        
    tempdata=load([stationname,datestr(period(i)+1,'YY'),datestr(period(i)+1,'mm'),datestr(period(i)+1,'dd'),'.fit']);
    tempdata(:,1)=((tempdata(:,1)-floor(tempdata(:,1)))+period(i)+1);
    tempdata=[((i+1)*ones(1,length(tempdata(:,1))))' tempdata];
    dmpsfitdata=[dmpsfitdata; tempdata];    
end
if exist([stationname,datestr(period(i)+1,'YY'),datestr(period(i)+1,'mm'),datestr(period(i)+1,'dd'),'.sum'])>0;
    tempdata=[];            
    tempdata=load([stationname,datestr(period(i)+1,'YY'),datestr(period(i)+1,'mm'),datestr(period(i)+1,'dd'),'.sum']);
    tempdata(:,1)=(tempdata(:,1)-floor(tempdata(:,1)))+period(i)+1;
    tempdata=[((i+1)*ones(1,length(tempdata(:,1))-1))' tempdata(2:end,:)];
    dmpsdata=[dmpsdata; tempdata];
end
if exist([stationname,datestr(period(i)+2,'YY'),datestr(period(i)+2,'mm'),datestr(period(i)+2,'dd'),'.fit'])>0;
    tempdata=[];        
    tempdata=load([stationname,datestr(period(i)+2,'YY'),datestr(period(i)+2,'mm'),datestr(period(i)+2,'dd'),'.fit']);
    tempdata(:,1)=((tempdata(:,1)-floor(tempdata(:,1)))+period(i)+2);
    tempdata=[((i+2)*ones(1,length(tempdata(:,1))))' tempdata];
    dmpsfitdata=[dmpsfitdata; tempdata];    
end
if exist([stationname,datestr(period(i)+2,'YY'),datestr(period(i)+2,'mm'),datestr(period(i)+2,'dd'),'.sum'])>0;
    tempdata=[];            
    tempdata=load([stationname,datestr(period(i)+2,'YY'),datestr(period(i)+2,'mm'),datestr(period(i)+2,'dd'),'.sum']);
    tempdata(:,1)=(tempdata(:,1)-floor(tempdata(:,1)))+period(i)+2;
    tempdata=[((i+2)*ones(1,length(tempdata(:,1))-1))' tempdata(2:end,:)];
    dmpsdata=[dmpsdata; tempdata];    
end