% Plotter simuleringsdata for TE-problemet

subplot(321)
plot(tout,simout(:,7))
hold on,
plot([0 tout(length(tout))],[Y(7) Y(7)],'kd-')
xlabel('time(hours)')
ylabel('Reactor pressure(kPa gauge)')
axis([0 100 2700 2900])

subplot(3,2,2)
plot(tout,simout(:,17))
hold on,
plot(tout,-Y17setpkt+Y(17),'k')
xlabel('time(hours)')
ylabel('Production rate(m3/h)')
axis([0 100 20 28 ])


subplot(3,2,3)
plot(tout,simout(:,40))
hold on,
plot(tout,-Y40setpkt+Y(40),'k')
xlabel('time(hours)')
ylabel('mole% G in product')
axis([0 100 40 60])


subplot(3,2,4)
plot(tout,simout(:,31))
hold on,
plot([0 tout(length(tout))],[Y(31) Y(31)],'kd-')
xlabel('time(hours)')
ylabel('mole% C in purge')
axis([0 100 0 30])


subplot(3,2,5)
plot(tout,podrag(:,3))
hold on,
plot([0 tout(length(tout))],[U(3) U(3)],'kd-')
xlabel('time(hours)')
ylabel('A feedflow(% of max flow)')
axis([0 100 0 100])


subplot(3,2,6)
plot(tout,podrag(:,4))
hold on,
plot([0 tout(length(tout))],[U(4) U(4)],'kd-')
xlabel('time(hours)')
ylabel('C feedflow(% of max flow)')
axis([0 100 0 100])




