
function sys=massebal(u,x,li,L,r)
%funksjonen massebal finner de deriverte av vannivået og totalnivået og trykket i separatoren

global M z T;

qinn_v=u(1); % volumstrømmen av vann inn i separatoren
qinn_o=u(2); % volumstrømmen av olje inn i separatoren
winn_g=u(3); % massestrømmen av gass inn i separatoren
qut_v=u(4);  % volumstrømmen av vann ut av separatoren
qut_o=u(5);  % volumstrømmen av olje ut i separatoren
wut_g=u(6);  % massestrømmen av gass ut i separatoren

hv=x(1); % vannivået i separatoren
ht=x(2); % totalt væskenivå i separatoren
Vv=x(3); % volumet av vann i separatoren
Vo=x(4); % volumet av olje i separatoren
p=x(5);  % trykket i separatoren


dhv=((1/(2*li*sqrt(r^2-(r-hv)^2)))*(qinn_v-qut_v));  

dht=((1/(2*L*sqrt(r^2-(r-ht)^2)))*(qinn_v+qinn_o-qut_v-qut_o)); 

dVv=(qinn_v-qut_v); 

dVo=(qinn_o-qut_o); 

Vt=pi*r^2*L; % totalvolumet i separatoren. Antar en sylindrisk tank

Vg=Vt-Vv-Vo; % gassvolumet i separatoren

R=8314;   % [J/kmol*K]
g=9.81;	 % [m/s2]

dP=((((z*R*T)/(Vg*M))*(winn_g-wut_g))-((p/Vg)*(qut_v-qinn_v+qut_o-qinn_o)));   

sys=[dhv;dht;dVv;dVo;dP];

