
% Generate steady-state data for column T200
% Saves data in file colT200_init.mat

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

clear all

MW=58.124;		%kg/kmol
mt=280;			%kg vaeske pr trinn

start1=0.5*ones(1,94)'; start2=(mt/MW)*ones(1,94)';
start=[start1;start2]; 

[t,x]=ode15s('colT200_lv',[0 20000],start); 
%Xinit = x(sel(size(x),1,1),:)';
lengthx=size(x); Xinit200 = x(lengthx(1),:)';

feed=22795+000;		%kg/h
refluks=88680+0000;		%kg/h
oppkok=95290+0000;		%kg/h
h=60;			%min/h

% Nominal inputs

LT=(refluks/MW)/h;    	         	 % Reflux
VB=(oppkok/MW)/h;   	                 % Boilup
                                                                  
F=(feed/MW)/h;                           % Feedrate
zF=0.29;                             	 % Feed composition
qF=1;                             	 % Feed liquid fraction
B=LT+F-VB;			       	%Bottom
D=F-B;					%Distillate


Uinit200 = [ LT VB D B F zF qF]';

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

save colT200_init  % saves in colT200_init.m


