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