% Fil init.m
% Eksempel paa oppbygging av en initfil for bruk i 
% programmene nyco7.m, nycobatch.m og idbatch.m 
% Komponentene i matrisene : 1 = metanol 2 = vann 3 = isopropanol

global Xglob Pglob NCglob Rglob Wglob Vglob Aglob Tkglob

Rglob=1.98721; % [cal/K mol]

NCglob=3; % Antall komponenter

Pglob=1;  % Trykk [atm]

% Alle prametere er hentet fra DECHEMA

% Antoine koeffisienter, Psat_i [mmHg] og T [K]

Aglob=[8.08097  1582.271  239.726     
   8.07131  1730.630  233.426  
   8.87829  2010.330  252.636];

% Wilson koeffisienter ved 1 atm., W_ij [cal/mol] 

Wglob=[  0         107.3832 140.1672
  469.5509     0        1238.9966 
   -132.7052   437.9789      0     ];

% Molart volum av rene vaesker [cm^3/mol] 

Vglob=[40.73 18.07 76.92 ];

% Initsiering av start vektor for vaeskesammnesetningen x

NT = 32;
zF = [0.7;0.1;0.2] % Komponenter i vektoren er metanol, vann, isopropanol
X=[];
for i=1:NT
  X=[X;zF];
end

% Initsiering av start temperatur

Tkglob=337*ones(1,NT-1)

% Holdup i koker MB og i kondensator MD [kmol] ved destillasjonsstart

MB = 100;
MD = 1;

% Startvektor for bruk i nyco7.m, nycobatch.m og idbatch.m

%A=[X' MD MB]      %P-regulator, lukket batch

A=[X' MD MB 0]     %PI-regulator, lukket batch

%A=[X' MD MB 0]    %P-regulator, vanlig batch

%A=[X' MD MB 0 0]  %PI-regulator, vanlig batch







