clc
clear all
close all

load eksperiment.mat
load olga3.mat
load olga2.mat
load matlab_ss_data.mat
load matlab_bif_data.mat
%%
figure

subplot(2,1,1)

plot(100*zs,y1,'k','LineWidth',3)
%plot(OLGA3(:,1),OLGA3(:,2),'r','LineWidth',2)
title('Pressure at inlet of pipeline (\itP_{in}\rm)','FontSize',12,'FontName','Times');
ylabel('\itP\rm_1 [bar]','FontSize',12,'FontName','Times')
xlabel('Valve opening [%]','FontSize',12,'FontName','Times')
hold on

%plot(OLGA2(:,7),OLGA2(:,1),'r','LineWidth',2)
plot(eks(:,1),eks(:,2),'r','LineWidth',2)
legend('Simple Model','Experiment')%,'Experiment')
plot(100*zs,y1min,'k',100*zs,y1max,'k','LineWidth',2)
hold on
%plot(OLGA3(:,1),OLGA3(:,3),'r',OLGA3(:,1),OLGA3(:,4),'r','LineWidth',2)
%plot(OLGA2(:,1),OLGA2(:,3),'r',OLGA2(:,1),OLGA2(:,4),'r','LineWidth',2)

plot(eks(:,1),eks(:,3),'r',eks(:,1),eks(:,4),'r','LineWidth',2)

subplot(2,1,2)
plot(100*zs,y2,'k','LineWidth',3)
%plot(OLGA3(:,1),OLGA3(:,5),'r','LineWidth',2)
title('Pressure at top of riser (\itP_{rt}\rm)','FontSize',12,'FontName','Times');
ylabel('\itP\rm_2 [bar]','FontSize',12,'FontName','Times')
xlabel('Valve opening [%]','FontSize',12,'FontName','Times')
hold on

plot(eks(:,1),eks(:,5),'r','LineWidth',2)

legend('Simple Model','Experiment')
plot(100*zs,y2,'-k')
plot(100*zs,y2min,'-k',100*zs,y2max,'-k','LineWidth',2)
hold on
%plot(OLGA3(:,1),OLGA3(:,6),'r',OLGA3(:,1),OLGA3(:,7),'r','LineWidth',2)

plot(eks(:,1),eks(:,6),'r',eks(:,1),eks(:,7),'r','LineWidth',2)

%C:\Users\Fie\Documents\Fie sine dokumenter\5.klasse\Prosjekt\Bilder

%print -depsc -C:\Users\Fie\Documents\Fie sine dokumenter\5.klasse\Prosjekt\Bilder