% pc.m

% Inn: TVN liq, TN vap, TVN ng, Tvegg. Ut: deriverte til ODE-løser via s_pc
% put  liq = u(1)
% put  vap = u(2)
% wng fra forrige = u(3)
% Hng fra forrige = u(4)
% pinn ng  = u(5)
% winn liq = u(6)   [kg/s]
% winn vap = u(7)
% Hinn co2 = u(8)   [J/s] Totalt for l og v

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

y_co2 = y(1:5);   %TVNTN
y_ng  = y(6:8);   %TVN
T_vegg= y(9);     %T

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

dy_co2 = pc_co2(y_co2,[u(1:2);u(6:8);Q_co2]);
% dy_co2 = zeros(5,1);

dy_ng  = pc_ng(y_ng,[u(3:5);Q_ng]);
% dy_ng = zeros(3,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 = 0;

dy = [dy_co2; dy_ng; dT_vegg];