
% Generate steady-state data for column A
% Saves data in file multicom_init.mat

%% Do this by simulating 20000 min with stabilized LV-model:

clear all
NT=32;
NF=20;
NC=5;
X0=[0.4057*ones(NT,1);0.35375*ones(NT,1);0.05642*ones(NT,1);0.1134*ones(NT,1);121; ones(NF-1,1)*6.32;ones(NT-NF-1,1)*3.97;226];
[t,x]=ode15s('multicomAE_lv',[0 20000],X0);

%Xinit = x(sel(size(x),1,1),:)';
lengthx=size(x); Xinit = x(lengthx(1),:)';

% Nominal inputs
%LT=71.46;                              % Reflux kmol/min
LT=63.5;
VB=88.62;                              % Boilup kmol/min
D=30.32;                               % Distillate
B=45.32;                               % Bottoms
F=75.64;                               % Feedrate kmol/min
qF=0.93;                               % Feed liquid fraction

zF=[0.40570 0.35375 0.05642 0.1134 0.07073];% Feed compositions for component A to component NC (CO2+C1+C2,C3,iC4,nC4,C5+)
Uinit = [ LT VB D B F qF zF]';

clear x; clear t; clear LT; clear VB; clear D; clear B;
clear F; clear zF; clear qF;

save multicomE_init  % saves in cola_init.m


