close all;
clear;
clc;
y=1;
s = char ('C');
delete ('Verlaeufe.xls');
for x = (1:.1:1.3);  %Zaehler für Druck


In{1} = 1500;       %Drehzahl
In{2} = 0.0;        %Kraftstoffeinspritzung
In{3} = 10.15;      %Epsilon
In{4} = x;          %Druck
parameterSetzen( ...
    'MTU396.xml', ...
    'Motormodell.xml', ...
    In, ...
    true );

!cDyn

load Ausgabe.mat;

Matrix.Grad(1) = {Ausgabe_KW.phi_gradKW};
Matrix.Druck(1,y) = {Ausgabe_KW.Zylinder1.p};
Matrix.Temperatur(1,y) = {Ausgabe_KW.Zylinder1.T};
Matrix.Waerme(1,y) = {Ausgabe_KW.Zylinder1.Q};


figure(1);

%Druckverlauf
h(1) = subplot(131); 
hold on;
title(h(1), sprintf('Druckverlauf %s'),...
            'Interpreter', 'none' );
yLabel(h(1), 'Druck [bar]');
xLabel(h(1), 'Kurbelwinkel [°]');
%         for m = 1:length(Matrix.Druck)
plot( h(1), Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.p);
%             l_str{m} = num2str(Druck.n_Line(m).n);
%         end
        
%         legend(h(1), l_str, 'FontName','Arial',...
%             'FontSize',8,...
%             'Position',[0.3674 0.5906 0.09453 0.15]);



%Temperaturverlauf
h(2) =subplot(132);   
hold on;
title(h(2), sprintf('Temperaturverlauf %s'),...
            'Interpreter', 'none' );
yLabel(h(2), 'Temperatur [Kelvin]');
xLabel(h(2), 'Kurbelwinkel [°]');
plot( h(2), Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.T);

%Waermeverlauf
h(3) = subplot(133);
hold on;
title(h(3), sprintf('Waermeverlauf %s'),...
            'Interpreter', 'none' );
yLabel(h(3), 'Q [J]');
xLabel(h(3), 'Kurbelwinkel [°]');
plot( h(3), Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.Q);

%Druckverlauf in Excell
k = {'Kurbelwinkel [°]'};
l ={'Druck [bar]'};
xlswrite('Verlaeufe', k, 'Druck', 'B2');
xlswrite ('Verlaeufe', l, 'Druck', 'B1');
xlswrite('Verlaeufe', x, 'Druck', [s '1']);
xlswrite('Verlaeufe',Matrix.Druck{1,y} , 'Druck', [s '3']);
xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Druck', 'B3'  );
% Temperaturverlauf in Excell
k = {'Kurbelwinkel [°]' 'Temperatur [K]'};
l ={'Druck [bar]'};
xlswrite('Verlaeufe', k, 'Temperatur', 'B2');
xlswrite ('Verlaeufe', l, 'Temperatur', 'B1');
xlswrite('Verlaeufe', x , 'Temperatur', [s '1']);
xlswrite('Verlaeufe', Matrix.Temperatur{1,length(Matrix.Temperatur)}, 'Temperatur', [s '3']);
xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Temperatur', 'B3'  );
% Waermeverlauf in Excell
k = {'Kurbelwinkel [°]' 'Waerme [J]'};
l ={'Druck [bar]'};
xlswrite('Verlaeufe', k, 'Waerme', 'B2');
xlswrite ('Verlaeufe', l, 'Waerme', 'B1');
xlswrite('Verlaeufe', x , 'Waerme', [s '1']);
xlswrite('Verlaeufe', Matrix.Waerme{1,y}, 'Waerme', [s '3']);
xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Waerme', 'B3'  );


y=y+1;                      %Zaehlvariable für Spalte in Matrix
s=s+1;                      %Zaehlvariable für Spalte in Excell
end


hold off;
%Druckverlauf in Excell
% k = {'Kurbelwinkel [°]'};
% l ={'Druck [bar]'};
% xlswrite('Verlaeufe', k, 'Druck', 'B2');
% xlswrite ('Verlaeufe', l, 'Druck', 'B1');
% xlswrite('Verlaeufe', x, 'Druck', [s '1']);
% xlswrite('test',Matrix.Druck , 'Druck');
% xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Druck', 'B3'  );
% % Temperaturverlauf in Excell
% k = {'Kurbelwinkel [°]' 'Temperatur [K]'};
% l ={'Druck [bar]'};
% xlswrite('Verlaeufe', k, 'Temperatur', 'B2');
% xlswrite ('Verlaeufe', l, 'Temperatur', 'B1');
% xlswrite('Verlaeufe', x , 'Temperatur', [s '1']);
% xlswrite('Verlaeufe', Matrix.Temperatur{1,length(Matrix.Temperatur)}, 'Temperatur', [s '3']);
% xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Temperatur', 'B3'  );
% % Waermeverlauf in Excell
% k = {'Kurbelwinkel [°]' 'Waerme [J]'};
% l ={'Druck [bar]'};
% xlswrite('Verlaeufe', k, 'Waerme', 'B2');
% xlswrite ('Verlaeufe', l, 'Waerme', 'B1');
% xlswrite('Verlaeufe', x , 'Waerme', [s '1']);
% xlswrite('Verlaeufe', Matrix.Waerme{1,y}, 'Waerme', [s '3']);
% xlswrite('Verlaeufe', Matrix.Grad{1,1}, 'Waerme', 'B3'  );
% subplot(131);
% title('Druckverlauf');
% xlabel('Kurbelwinkel [°]');
% ylabel('Druck [bar]');
% plot(Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.p);
% h(1) = subplot (131);
% title('Temperaturverlauf');
% xlabel('Kurbelwinkel [°]');
% ylabel('Temperatur [Kelvin]');
% plot(Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.T);
% hold all;
% title('Waermeverlauf');
% xlabel('Kurbelwinkel [°]');
% ylabel('Q');
% plot(Ausgabe_KW.phi_gradKW, Ausgabe_KW.Zylinder1.Q);
% hold all;