function [pt2,Tt2,m2,Pit,Awdt]=turbina(pt1,Tt1,m1,Pt,et,M_WDT,mch,Tch)

%[Tt2,pt2,Pit]=turbina(Tt1,pt1,Pt,et,m)
%Funkcja oblicza remperaturê i ciœnienie w parmetrach spiêtrzenia za
%turbin¹. Danymi wejœciowymi s¹ temperatura i ciœnienie spiêtrzenia na
%wlocie do turbiny, moc, jak¹ ma wytworzyæ turbina, sprawnoœæ turbiny, 
%strumieñ masowy spalin na wlocie do turbiny. WA¯NE Sprawnoœæ mo¿na
%zdefiniowaæ jako izentropowa i politropowa. Dlatego sprawnoœæ turbiny
%et=[sprawnoœæ, rodzaj sprawnoœci], gdzie sprawnoœæ powinna siê zawieraæ w
%przedziale od 0 do 1, a rodzaj sprawnoœci powinien byæ definiowany jako
%0 dla sprawnoœci izentropowej i 1 dla sprawnooœci politropowej. Dla
%sprawnoœci izentropowej turbiny, wystarczaj¹ce jest aby zdefinowaæ sam¹
%wartoœæ sprawnoœci, a program domyœlnie przyjmie ¿e jest to sprawnoœæ
%izentropowa.
%M_WDT - prêdkoœæ w przekroju minimalnym WDT, Dla przep³ywu krytycznego 
%powinna byæ M_WDT=1 


if exist('mch','var')==0
    mch=0;Tch=0;
end
if isnan(mch)
    mch=0;
    Tch=0;
end

cp=1180;
R=287;
k=1/(1-R/cp);
k=1.33;
cp1=1100;
% Zmienna pomocnicza do sumaowania mass czynnika ch³odz¹cego
z1=ones(length(mch),1);
m2=m1+mch*z1;
Tt2=(m1*cp*Tt1-Pt+cp1*mch*Tch')./(m1*cp+cp1*mch*z1);
if length(et)==1
    et(2)=0;
end
if et(2)==0
     Pit=(et(1)./(et(1)+Tt2./Tt1-1)).^(k./(k-1));
elseif et(2)==1
    Pit=(Tt1./Tt2).^(k/(k-1)/et(1));
else
    Pit=NaN;
    Tt2=NaN;
    disp('¯le zdefinowano rodzaj spreanoœci turbiny')
end
pt2=pt1./Pit;

%% OBLIZCENIA PRZKROJU MINIMALNEGO WDT
TtWDT=(m1*cp*Tt1+mch(1)*Tch(1)*cp1)/(m1+mch(1))/cp;
pkr=pt1*0.985/(1+(k-1)/2*M_WDT^2)^(k/(k-1));
Tkr=TtWDT/(1+(k-1)/2*M_WDT^2);
rokr=pkr/R/Tkr;
ckr=sqrt(k*R*Tkr)*M_WDT;
Awdt=(m1+mch(1))/rokr/ckr;       
%Awdt=m./(pt1*0.98)*sqrt(R./k*Tt1).*((k+1)./2).^((k+1)./2/(k-1));
end

