% This file contain the data for the case
global coldata dyndata
coldata.Nt=8; % Number of trays in each section
coldata.alfa=[9 3 1]; % Relative volality
coldata.M0=1; % Holdup reboil and condense
coldata.yset=[0.97;0.97;0.97];
coldata.d0=[1;1/3;1/3;1;coldata.yset]; % Nominal disturbance
%coldata.Tboil=[0;50;100];
coldata.Tboil=[300;342.85;400];
coldata.P=1;%Pressure 1 atm
% Data for the dynamic simulations
load nyx
load nyu
xnomopt=X;
unomopt=u;
xnom=xnomopt;
unom=unomopt; %u=[L,V,S,Rl,Rv,F,z(1),z(2),q]
clear xnomopt unomopt X u
dyndata.L0=unom(1);
dyndata.V0=unom(2);
dyndata.S0=unom(3);
dyndata.Rl0=unom(4);
dyndata.Rv0=unom(5);
dyndata.F0=unom(6);
dyndata.zF=unom(7:8);
dyndata.q=unom(9);
dyndata.yset=coldata.yset
dyndata.yselect=1; % Compositions and temperatures

%Get states, and calculate Antoine coeff.
[xp]=genpet([],xnom,unom,3,coldata.Nt,coldata.alfa,coldata.M0,3,[]); %ynom=[xA,xB,xC,T1..T48]
[Temp,coeff]=temp_antoine(xp);%Based on Antoines eq.
coldata.A=coeff(1);
coldata.B=coeff(2);
coldata.C=coeff(3);
clear unom xnom coldata.Tboil
