clear all
close all
clc

load('D:\werf_je\Werk\1.201.819-Proefsuppletie_Galgeplaat_2e_evaluatie\Data_analyse\Condities\zw_data.mat');

for i = 1:length(zw_data.naam)
    clear T Zw dum1 dum2 dum1b dum2b
    T = zw_data.t{i};
    Zw = zw_data.zw{i};
    k1 = 0;
    k2 = 0;    
    for j = 3:length(T)-2
        if Zw(j) >= Zw(j-1) & Zw(j) >= Zw(j-2) & Zw(j) > Zw(j+1) & Zw(j) > Zw(j+2) & Zw(j) > 0
            k1 = k1 +1;
            dum1(k1) = Zw(j);            
            dum1b(k1) = T(j);                        
        end
        if Zw(j) <= Zw(j-1) & Zw(j) <= Zw(j-2) & Zw(j) < Zw(j+1) & Zw(j) < Zw(j+2) & Zw(j) < 0.5
            k2 = k2 +1;
            dum2(k2) = Zw(j);            
            dum2b(k2) = T(j); 
        end        
    end
    vloed(i) = {dum1};
    Tvloed(i) = {dum1b};
    eb(i) = {dum2};
    Teb(i) = {dum2b};
end

% figure
% plot(zw_data.t{1},zw_data.zw{1},'k-')
% hold on
% plot(Tvloed{1},vloed{1},'r.')
% plot(Teb{1},eb{1},'r.')
% 
% figure
% plot(zw_data.t{2},zw_data.zw{2},'k-')
% hold on
% plot(Tvloed{2},vloed{2},'r.')
% plot(Teb{2},eb{2},'r.')
% 
% figure
% plot(zw_data.t{3},zw_data.zw{3},'k-')
% hold on
% plot(Tvloed{3},vloed{3},'r.')
% plot(Teb{3},eb{3},'r.')
% 
% figure
% plot(zw_data.t{4},zw_data.zw{4},'k-')
% hold on
% plot(Tvloed{4},vloed{4},'r.')
% plot(Teb{4},eb{4},'r.')

length(Tvloed{3}) + length(Tvloed{4})

%Vb plot waterstand (oktober 2009)
figure
subplot(2,1,1)
plot(zw_data.t{4},zw_data.zw{4},'b-','LineWidth',1)
hold on
xlabel('dag')
ylabel('z_w (m NAP)')
datum_begin = datenum('01-06-2009','dd-mm-yyyy');
datum_einde = datenum('01-07-2009','dd-mm-yyyy');
xData = datum_begin:2:datum_einde;        
set(gca,'XTick',xData);
datetick('x','dd','keepticks')
ylim([-3 3])
grid on
title('Waterstanden, juni 2009')
set(gca,'ytick',[-3:0.5:3])    
saveas(gcf,'EPS\zw_vb.eps','psc2');
saveas(gcf,'PNG\zw_vb.png');

levels = [-3.5:0.05:3.5]';

%Aangenomen duur getijcyclus in [uren]
T = 12.6;

%"Levels of interest" (loi)
loi = [-3,-2,-1.5:0.25:1.5,2,3];

jj = 0;
%Loop voor meetstations
for i = [1,3]
    jj = jj + 1;
    clear Z
    Z = [zw_data.zw{i},zw_data.zw{i+1}];
    zmax(:,jj) = max(Z);
    zmin(:,jj) = min(Z);
    zlength = sum(~isnan(Z));
    for j = 1:length(levels)
        clear dum1
        dum1 = find(Z < levels(j));
        if ~isempty(dum1)
            onder(j,jj) = length(dum1)/zlength;
        else
            onder(j,jj) = 0;
        end        
    end
end
% 
%Maken figuren
figure
% plot(levels,onder(:,1)*100,'k-','LineWidth',1.0)
% hold on
plot(levels,onder(:,2)*100,'b-','LineWidth',1.0)
hold on
grid on
xlabel('z_w (m NAP)')
ylabel('onderschrijding (%)')
%legend('Marollegat','Stavenisse',2)
xlim([-3.5 3.5])
saveas(gcf,'EPS\onderschrijding2.eps','psc2');
saveas(gcf,'PNG\onderschrijding2.png');
% 
% figure
% plot(levels,onder(:,1)*T,'k-','LineWidth',1.0)
% hold on
% plot(levels,onder(:,2)*T,'r-','LineWidth',1.0)
% grid on
% xlabel('z_w (m NAP)')
% ylabel('duur onderschrijding (uren)')
% legend('Marollegat','Stavenisse',2)
% xlim([-3.5 3.5])
% saveas(gcf,'EPS\droogvalduur.eps','psc2');
% saveas(gcf,'PNG\droogvalduur.png');
% 
% %Maken data struct en wegschrijven data
% dvd_data.droogvalduur = {T*onder(:,1) T*onder(:,2)};
% dvd_data.getijperiode = {T};
% dvd_data.levels = {levels};
% dvd_data.naam = {'Marollegat' 'Stavenisse'};
% dvd_data.onder = {onder(:,1) onder(:,2)};
% dvd_data.zwmax = {zmax};
% dvd_data.zwmin = {zmin};
% save dvd_data dvd_data
% 
% %Wegschrijven droogvalduur in txt-file voor bepaalde waterstanden
% for i = 1:length(loi)
%     M(i,1) = loi(i);
%     dum = find(levels == loi(i));
%     M(i,2) = T*onder(dum,1);
%     M(i,3) = T*onder(dum,2);    
% end
% dlmwrite('dvd.txt', M,'delimiter', '\t','precision','%0.2f');