% tank.m
% Gir entalpien inn på PC for CO2 og NG, og tetthet for CO2

function tank(b,d)
% CO2:
T_co2 = d.T_co2tank;
p_co2 = b.p_co2tank;
V_co2 = d.V_co2tank;
N0_co2 = d.rho_co2liq*V_co2/b.MW_pc(1);


% Væske:
N_co2 = N0_co2;
dN_co2= N_co2;
while abs(dN_co2)>1e-6*N_co2
    A_co2 = co2_sw([T_co2,V_co2,N_co2]);
    dN_co2 = (p_co2 + A_co2.g(2))/(-A_co2.H(2,3));
    N_co2 = N_co2 + dN_co2;
end
rho_co2 = N_co2/V_co2*b.MW_pc(1); % [kg/m3]
H_co2 = (-A_co2.g(1)*T_co2 + A_co2.g(3)*N_co2)/N_co2;   % [J/mol]

% Naturgass:
T_ng = d.T_pc_inn(2);
V_ng = d.V_ngtank;
N0_ng= b.p_pc_inn(2)*V_ng/(b.R*T_ng);

N_ng = N0_ng;
dN_ng= N_ng;
while abs(dN_ng)>1e-6*N_ng
    A_ng = ng_pr([T_ng, V_ng, N_ng, N_ng*1e-6, N_ng*1e-6]');
    dN_ng = (b.p_ngtank + A_ng.g(2))/(-A_ng.H(2,3));
    N_ng = N_ng + dN_ng;
end
rho_ng = N_ng/V_ng*b.MW_pc(2); % [kg/m3]
H_ng = (-A_ng.g(1)*T_ng + A_ng.g(3)*N_ng)/N_ng; %[J/mol]

%Resultat:
co2_tank = [rho_co2, H_co2];
ng_tank  = [rho_ng, H_ng];
save tankverdier co2_tank ng_tank