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('droogvalduur (uren)')
% ylabel('droogvallend areaal (ha)')
% title('Gebied: cirkel bestek Galgeplaat')
% saveas(gcf,'EPS\dva.eps','psc2');
% saveas(gcf,'PNG\dva.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');