clear all
close all
clc

load('D:\werf_je\Werk\1.201.819-Proefsuppletie_Galgeplaat_2e_evaluatie\Data_analyse\Condities\zw_data.mat');
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,'r-','LineWidth',1.0)
grid on
xlabel('z_w (m NAP)')
ylabel('onderschrijding (%)')
legend('Marollegat','Stavenisse',2)
xlim([-3.5 3.5])
saveas(gcf,'EPS\onderschrijding.eps','psc2');
saveas(gcf,'PNG\onderschrijding.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');