clear all

slutt=8000;
load kaskadep2010_ny;
tid=t(1:slutt)/60;
p=y_m1(1:slutt);
q=y_m2(1:slutt)*60;
u=u(1:slutt);

figure(1)
subplot(3,1,1)
plot(tid,p,tid(2000:8000),0.010,'r')
ylabel('P2 [barg]')
xlabel('time [min]')
legend('P2','Setpoint')
AXIS ([0 30 0 0.15])

subplot(3,1,2)
plot(tid,q)
ylabel('Mass flow [kg/min]')
xlabel('time [min]')
AXIS ([0 30 0 7])

subplot(3,1,3)
plot(tid,u,tid,0.20)
ylabel('Valve position [-]')
xlabel('time [min]')
legend('z','bifurcation point')
AXIS ([0 30 0 1])

break
clear all
slutt=7550;
load kaskadep2015_ny;
tid=t(1:slutt)/60;
p=y_m1(1:slutt);
q=y_m2(1:slutt)*60;
u=u(1:slutt);

figure(2)
subplot(3,1,1)
plot(tid,p,tid,0.015,'r')
ylabel('P2 [barg]')
xlabel('time [min]')
legend('P2','Setpoint')
AXIS ([0 30 0 0.18])

subplot(3,1,2)
plot(tid,q)
ylabel('Mass flow [kg/min]')
xlabel('time [min]')
AXIS ([0 30 0 7])

subplot(3,1,3)
plot(tid,u,tid,0.20,'r')
ylabel('Valve position [-]')
xlabel('time [min]')
legend('z','bifurcation point')
AXIS ([0 30 0 1])
