close
clear

%Auswählen der auszuwertenden *.xls Datei
[FileNames,PathName] = uigetfile('*.xls','Open xls files (*.xls)','MultiSelect','on');

%Einlesen der *.xls Datei
[nums, txt] = xlsread(FileNames,'Regen1');


[typ, desc] = xlsfinfo(FileNames)

numSheets = length(desc);  % Anzahl Arbeitsblätter 

ab= nums(:,1) %=H gemessen
b= nums(:,2)  %=Q gemessen
c= nums(:,3)  %=Q berechnet
d= nums(:,4)  %=H berechnet
Regensumme = nums(1,10)
Flaeche = nums(1,11)
Trockenwetterabfluss = nums(1,12)
Drift=nums(1,8) %Drifthöhe 
DriftZeile=nums(1,9)  %Driftzeile
  
H=nums(1,7) %Stationshöhe
filenamesmod=txt(:,1) %Bezeichnung X-Achse
z=(length(ab))

a=ab+H %H Stationshöhe

xa=(1:z)
x=xa'
Nullhgem = find(ab==0) %Ist ein gemessener H-Wert =0?
NullQgem = find(b==0)  %Ist ein gemessener Q-Wert =0?

w = a;

 if Drift~=0
    w(DriftZeile:end) = a(DriftZeile:end) + Drift; 
 else
 end

Qmaxber = max(b)
Qmaxgem = max(c)
Vgem=sum(c)*300
Vber=sum(b)*300
Prozentabweichung=Vgem/Vber

f = -Qmaxber/10
g = f/2


figure(1)
set(gca,'position',[0.2 0.2 0.8 0.8])
%set(gcf, 'paperposition',[0,8, 20, 20])
set(gcf, 'papersize',[21,29])
hold on
plot(x,c,'r')
plot(x,b,'b')
title('Durchfluss Q in l/s')
grid
legend('Qber','Qgem')
xlabel('Zeit'); ylabel('Q in l/s')

set(gca,'xtick',x,'xticklabel',filenamesmod)
hold off

text(0.2,f,['Qmax berechnet = ',num2str(Qmaxber) ' l/s'],'FontSize',10)
text(0.2,f+g,['Qmax gemessen = ',num2str(Qmaxgem) ' l/s'],'FontSize',10)
text(0.2,f+2*g,['Abflussvolumen berechnet = ',num2str(Vber) ' m³'],'FontSize',10)
text(0.2,f+3*g,['Abflussvolumen gemessen = ',num2str(Vgem) ' m³'],'FontSize',10)
saveas(1,'Ergebnisse_Regen1.pdf')


HMaxber = max(a)
HMaxgem = max(d)


figure(2)
hold on
plot(x,d,'r')
plot(x,w,'b')
title('Wassterstand')
grid
legend('Hber','Hgem')
xlabel('Zeit'); ylabel('Wasserstand in müNN')
set(gca,'xtick',x,'xticklabel',filenamesmod)
hold off