%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% PRESA D'ARIA %%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Tt_2, Pt_2] = presa_daria(Ta, Pa)

global gamma_aria eta_d M

Tt_2 = Ta*(1+(gamma_aria-1)*(M^2)/2)
Pt_2 = Pa*(1+eta_d*(gamma_aria-1)*(M^2)/2)^(gamma_aria/(gamma_aria-1))


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% COMPRESSORE %%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Tt_3, Pt_3] = compressore(Tt_2, Pt_2, pi_c)

global gamma_aria eta_c

Tt_3 = Tt_2+(1/eta_c)*(Tt_2*(pi_c^((gamma_aria-1)/gamma_aria)-1))
Pt_3 = Pt_2*pi_c


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%% CAMERA COMBUSTIONE %%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Tt_4, Pt_4] = camera_combustione(Tt_3, Pt_3)

global eta_cc, eta_m, delta_p, cpf, cpc, Tt_5, Tt_2

pi_b = 1-delta_p/Pt_3