clc
clear all
folder='EPL3_2009\EPL3\Th0';
M_dir=dir(folder);
M_dir=M_dir(3:end,1);
%Spectra=struct;
%s_i=0;
wavedir=[];
date=[];
for m_i=1:length(M_dir)
    filename=[folder '/' M_dir(m_i,1).name]%;
    [wavedir_day,date_day]=Readwavedir(filename);
    wavedir=[wavedir wavedir_day]; 
    date=[date date_day];
    %pause
end

%%2010
folder='EPL3_2010\EPL3\Th0';
M_dir=dir(folder);
M_dir=M_dir(3:end,1);
%Spectra=struct;
%s_i=0;

for m_i=1:length(M_dir)
    filename=[folder '/' M_dir(m_i,1).name]%;
    [wavedir_day,date_day]=Readwavedir(filename);
    wavedir=[wavedir wavedir_day]; 
    date=[date date_day];
    %pause
end

%%2011
folder='EPL3_2011\EPL3\Th0';
M_dir=dir(folder);
M_dir=M_dir(3:end,1);
%Spectra=struct;
%s_i=0;

for m_i=1:length(M_dir)
    filename=[folder '/' M_dir(m_i,1).name]%;
    [wavedir_day,date_day]=Readwavedir(filename);
    wavedir=[wavedir wavedir_day]; 
    date=[date date_day];
    %pause
end

%%2012
folder='EPL3_2012\EPL3\Th0';
M_dir=dir(folder);
M_dir=M_dir(3:end,1);
%Spectra=struct;
%s_i=0;

for m_i=1:length(M_dir)
    filename=[folder '/' M_dir(m_i,1).name]%;
    [wavedir_day,date_day]=Readwavedir(filename);
    wavedir=[wavedir wavedir_day]; 
    date=[date date_day];
    %pause
end


%%
wavedir(wavedir==999)=NaN;

%%
if 1
    plot(date,wavedir)
    datetick('x')
end

save ('wavedir.mat','date','wavedir')