
% 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=92;
NF=45;
NC=5;
X0=[0.00*ones(NT,1);0.00300052*ones(NT,1);0.31695238*ones(NT,1);0.67604645*ones(NT,1);66.4; ones(NF-1,1)*3.5793777;ones(NT-NF-1,1)*5.4250431;88.1];
[t,x]=ode15s('multicomAS_lv',[0 20000],X0);

%Xinit = x(sel(size(x),1,1),:)';
lengthx=size(x); Xinit = x(lengthx(1),:)';
  
% Nominal inputs
%LT=25.901199;                              % Reflux kmol/min
LT=32;
VB=34.5;                               % Boilup kmol/min
D=2.5;                               % Distillate
B=5.6;                               % Bottoms
F=8.1;                               % Feedrate kmol/min
qF=1;                               % Feed liquid fraction

zF=[0.00 0.00300052 0.31695238  0.67604645 0.00400065];% 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 multicomS_init  % saves in cola_init.m


