
clear all
close all
clc

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

load('dvd_data.mat');
load('dva_data.mat')

levels = dvd_data.levels{:};

% figure
% plot(levels,10^-4*dva_data.areaal{1},'k-','LineWidth',1)
% hold on
% plot(levels,10^-4*dva_data.areaal{2},'r-','LineWidth',1)
% plot(levels,10^-4*dva_data.areaal{5},'b-','LineWidth',1)
% plot(levels,10^-4*dva_data.areaal{8},'g-','LineWidth',1)
% plot(levels,10^-4*dva_data.areaal{10},'m-','LineWidth',1)
% plot(levels,10^-4*dva_data.areaal{13},'c-','LineWidth',1)
% grid on
% legend(datestr(dva_data.datum{1}),datestr(dva_data.datum{2}),datestr(dva_data.datum{5}),...
%        datestr(dva_data.datum{8}),datestr(dva_data.datum{10}),datestr(dva_data.datum{13}));
%  xlabel('waterstand (m NAP)')
%  ylabel('droogvallend areaal (ha)')
%  title('Gebied: cirkel bestek Galgeplaat')
%  xlim([-1 1])
% saveas(gcf,'EPS\dva.eps','psc2');
% saveas(gcf,'PNG\dva.png');

figure
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{1},'k-','LineWidth',1)
hold on
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{2},'r-','LineWidth',1)
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{5},'b-','LineWidth',1)
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{8},'g-','LineWidth',1)
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{10},'m-','LineWidth',1)
plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{13},'c-','LineWidth',1)
grid on
 legend(datestr(dva_data.datum{1}),datestr(dva_data.datum{2}),datestr(dva_data.datum{5}),...
        datestr(dva_data.datum{8}),datestr(dva_data.datum{10}),datestr(dva_data.datum{13}));
xlabel('droogvalduur (uren/getij)')
ylabel('droogvallend areaal (ha)')
title('Gebied: cirkel bestek Galgeplaat')
xlim([0 13])
set(gca,'xtick',[0:13])
saveas(gcf,'EPS\dva2.eps','psc2');
saveas(gcf,'PNG\dva2.png');



% figure
% plot(10^-4*dva_data.areaal{1},dvd_data.onder{2}*12.4,'k-','LineWidth',1)
% hold on
% plot(10^-4*dva_data.areaal{2},dvd_data.onder{2}*12.4,'r-','LineWidth',1)
% % plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{5},'b-','LineWidth',1)
% % plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{8},'g-','LineWidth',1)
% % plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{10},'m-','LineWidth',1)
% % plot(dvd_data.onder{2}*12.4,10^-4*dva_data.areaal{13},'c-','LineWidth',1)
% grid on
% %  legend(datestr(dva_data.datum{1}),datestr(dva_data.datum{2}),datestr(dva_data.datum{5}),...
% %         datestr(dva_data.datum{8}),datestr(dva_data.datum{10}),datestr(dva_data.datum{13}));
% % xlabel('droogvalduur (uren/getij)')
% % ylabel('droogvallend areaal (ha)')
% % title('Gebied: cirkel bestek Galgeplaat')
% ylim([0 13])
% set(gca,'ytick',[0:13])
% % saveas(gcf,'EPS\dva2.eps','psc2');
% % saveas(gcf,'PNG\dva2.png');

% % figure
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{1},'k-','LineWidth',1)
% % hold on
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{2},'r-','LineWidth',1)
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{5},'b-','LineWidth',1)
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{8},'g-','LineWidth',1)
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{10},'m-','LineWidth',1)
% % plot(dvd_data.droogvalduur{2},10^-4*dva_data.areaal{13},'c-','LineWidth',1)
% % grid on
% % legend(datestr(dva_data.datum{1}),datestr(dva_data.datum{2}),datestr(dva_data.datum{5}),...
% %        datestr(dva_data.datum{8}),datestr(dva_data.datum{10}),datestr(dva_data.datum{13}));
% % xlabel('droogvalduur (uren)')
% % ylabel('droogvallend areaal (ha)')
% % title('Gebied: cirkel bestek Galgeplaat')
% % saveas(gcf,'EPS\dva2.eps','psc2');
% % saveas(gcf,'PNG\dva2.png');
% 
for i = 1:length(loi)
    M(i,1) = loi(i);
    dum = find(levels == loi(i));
    for j = 1:1:length(dva_data.areaal)
        clear dum2
        dum2 = 10^-4*dva_data.areaal{j};
        M(i,j+1) = dum2(dum);        
    end    
end
dlmwrite('dva.txt', M,'delimiter', '\t','precision','%0.2f');