
% Generate steady-state data for column T300
% Saves data in file colT300_init.mat

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

clear all


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

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

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

% Nominal inputs

feed=78900;		      %kg/h
refluks=202811;		%kg/h 
oppkok=233540+600;	%kg/h 
h=60;			         %min/h

LT=(refluks/MW)/h;		 % Reflux
VB=(oppkok/MW)/h;                % Boilup

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

Uinit300 = [ 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 colT300_init  % saves in colT300_init.m


