
% 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=40;
NF=19;
NC=5;
X0=[0.00*ones(NT,1);0.00206825*ones(NT,1);0.22003723*ones(NT,1);0.47840547*ones(NT,1);18; ones(NF-1,1)*1.48;ones(NT-NF-1,1)*1.67;62];
[t,x]=ode15s('multicomAB_lv',[0 20000],X0);

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

% Nominal inputs
%LT=13.5;
LT=8.6; % Reflux kmol/min
%LT=6.0;
%VB=9.9;                        % Boilup kmol/min
VB=11.85;
D=8.1;                               % Distillate
B=3.5;                               % Bottoms
%B=4.11;
F=11.6;                % Feedrate kmol/min
%F=12.21;                % Feedrate kmol/min
qF=0.579;                               % Feed liquid fraction

zF=[0.00 0.00206825 0.22003723 0.47840547 0.29948906];% 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 multicomB_init  % saves in cola_init.m


