% %Initial point

%This script calculates the initial point for the equations oriented
%simulation.

%Load the steady state results 
filename='SSresults.mat';
load (filename)

%Initialize compressibility 
Zv=zeros(1,15); % Vapour phase
Zl=zeros(1,15); % Liquid phase

% Calculate compressibilities for each stream on each phase
for i=1:15
    [Zl(i),phi,H,V]= srks(z(i,:),Ts(i),P(i),1);
    [Zv(i),phi,H,V]= srks(z(i,:),Ts(i),P(i),2);
end