% Visualisering av resultater
%--------------------------------------------------------------------------
% Funksjonene mottar løsningen på det dynamiske problemet og plotter
% typiske størrelser mot tidsaksen. I tillegg plottes pådrag og strømmer og
% trykk som blir tatt inn som globale størrelser
%--------------------------------------------------------------------------
function grafer = Plotter(t,T,d)
global u P w Tventil y

% Gir lokale navn
%--------------------------------------------------------------------------
T_NG1A = [T(:,1),T(:,d.n+2:2*d.n+1),T(:,2*d.n+2:3*d.n+1),T(:,3*d.n+2:4*d.n+1),T(:,4*d.n+2:5*d.n+1)];
T_NG1B = [T(:,5*d.n+13),T(:,6*d.n+14:7*d.n+13),T(:,7*d.n+14:8*d.n+13),T(:,8*d.n+14:9*d.n+13),T(:,9*d.n+14:10*d.n+13)];
Tw_NG1A = T(:,2:d.n+1);
Tw_NG1B = T(:,5*d.n+14:6*d.n+13);
T_SW1 = [T(:,5*d.n+2),T(:,5*d.n+3)];
T_SW2 = [T(:,5*d.n+4),T(:,5*d.n+5)];
T_tank = [T(:,5*d.n+10),T(:,5*d.n+11),T(:,5*d.n+12)];
mg_tank = T(:,10*d.n+17:10*d.n+19);
ml_tank3 = T(:,10*d.n+20);
mg_SW1 = T(:,10*d.n+16);
ml_1a = T(:,5*d.n+6); 
mg_1a = T(:,5*d.n+7); 
ml_SW2 = T(:,5*d.n+8);
mg_SW2 = T(:,5*d.n+9);
ml_1b = T(:,10*d.n+14);
mg_1b = T(:,10*d.n+15);
%--------------------------------------------------------------------------

z = [0:1/19:1]';
%--------------------------------------------------------------------------
% Plotter resultater
%--------------------------------------------------------------------------

% Temperaturer inn og ut av NG1A og NG1B
%--------------------------------------------------------------------------
subplot(2,2,1)
plot(t,T_NG1A(:,1),'g',t,Tw_NG1A(:,1),'*',t,T_NG1A(:,2),'r',t,T_NG1A(:,d.n+2),'k',t,T_NG1A(:,2*d.n+2),'m',t,T_NG1A(:,3*d.n+2),'b')
title('Temperaturer i NG1A ved inngang')

subplot(2,2,2)
plot(t,Tventil(:,1),'g',t,Tw_NG1A(:,d.n),'*',t,T_NG1A(:,d.n+1),'r',t,T_NG1A(:,2*d.n+1),'k',t,T_NG1A(:,3*d.n+1),'m',t,T_NG1A(:,4*d.n+1),'b')
title('Temperaturer ut av NG1A')

subplot(2,2,3)
plot(t,T_NG1B(:,1),'g',t,Tw_NG1B(:,1),'*',t,T_NG1B(:,2),'r',t,T_NG1B(:,d.n+2),'k',t,T_NG1B(:,2*d.n+2),'m',t,T_NG1B(:,3*d.n+2),'b')
title('Temperaturer i NG1B ved inngang')

subplot(2,2,4)
plot(t,Tventil(:,2),'g',t,Tw_NG1B(:,d.n),'*',t,T_NG1B(:,d.n+1),'r',t,T_NG1B(:,2*d.n+1),'k',t,T_NG1B(:,3*d.n+1),'m',t,T_NG1B(:,4*d.n+1),'b')
title('Temperaturer ut av NG1B')
legend('1','Vegg','2','3','4','5')

% Temperaturer i SW1 og SW2
%--------------------------------------------------------------------------
figure
plot(t,T_SW1(:,1),'y',t,T_SW1(:,2),'b',t,T_SW2(:,1),'r',t,T_SW2(:,2),'k')
title(' Temperaturer ut av SW1 og SW2')
legend('SW1(1)','SW1(2)','SW2(1)','SW2(2)')

% Temperaturer i Tank 1-3
%--------------------------------------------------------------------------
figure
plot(t,T_tank(:,1),'r',t,T_tank(:,2),'k',t,T_tank(:,3),'b')
title('Tank temperaturer')
legend('1','2','3')

% Trykk
%--------------------------------------------------------------------------
figure
plot(t,P(:,7),'r',t,P(:,8),'b',t,P(:,9),'g')
title(' Trykk i tankene')
legend('1','2','3')

% Gjennomstrømning
%--------------------------------------------------------------------------
figure
plot(t,w(:,1),'r',t,w(:,6),'g',t,w(:,11),'k')
title('Gjennomstrømning')
legend('NG1A','NG1B','Totalt')

% Pådrag
%--------------------------------------------------------------------------
figure
plot(t,u.n(1:length(t)),t,u.z(1:length(t),1),t,u.z(1:length(t),2))
title('Pådrag')
legend('Normalisert turtall på kompressorene','Ventilåpning for ventil 1','Ventilåpning for ventil 2')
%-------------------------------------------------------------------------- 
