% pc.m

% Inn: TVN liq, TN vap, TVN ng, Tvegg. Ut: deriverte til ODE-løser via s_pc
% winn liq = u(1)   [kg/s]
% wut  vap = u(2)
% Hinn co2 = u(3)   [J/s] Totalt for l og v (beregnet i tank.m)
% w NG     = u(4)
% Hinn NG  = u(5)

function dy=pc(t,y,u,ni)
global b

y_co2 = y(1:5*b.n);         %TVNTN
y_ng  = y(5*b.n+1:8*b.n);   %TVN
T_vegg= y(8*b.n+1:end);     %T

Q_co2 = b.UA_pc(1)*(T_vegg-y_co2(1:b.n));   %[J/s]
Q_ng  = b.UA_pc(2)*(T_vegg-y_ng(1:b.n));    

%Q_co2 = zeros(1,b.n);

dy_co2 = pc_co2(y_co2, [u(1:3);Q_co2]);
% dy_co2 = zeros(5*b.n,1);

dy_ng  = pc_ng(y_ng, [u(4:5);Q_ng]);
% dy_ng = zeros(3*b.n,1);

% Qco2 > 0. Ut fra vegg--> minus foran
% Qng  < 0. Inn i vegg--> minus foran:
dT_vegg = b.mCp_vegg*(-Q_ng-Q_co2);
% dT_vegg = zeros(b.n,1);


dy = [dy_co2; dy_ng; dT_vegg];