function [Pt,Tt,ph,Th,V0]=otoczenie(H,dT,VM,Vid)
 
%[Th,ph]=otoczenie_1(H,dT)
%
%funkcja oblicza parametry termodynamiczne ciœnienia i temperatury 
%w zale¿noœci od wysokoœci korzystaj¹c z AMW
%
%Na wejœciu:
%   H - wysokoœæ w [m]
%   dT - poprawka temperatury zale¿na od pory roku (pory dnia)
%        gdy nie jest podana traktowane jest jakby by³a równa 0
%
% zale¿noœci z Muszyñski, Orkisz - Modelowanie Turbinowych Silników Odrzutowych
 
if  exist('dT','var')==0
   dT=0;
end

k=1.4;
R=287;

if H<=11000
   th=288-0.00651*H;
   ph=101325*(th/288)^5.2533;
elseif H>11000
   th=216.5;
   ph=23000*exp((11000-H)/6138);
else
    th=NaN;
    ph=NaN;
end

Th=th+dT;

if Vid=="M"
    Tt=Th*(1+(k-1)/2*VM^2);
    Pt=ph*(1+(k-1)/2*VM^2)^(k/(k-1));
    V0=VM*sqrt(k*R*Th);
else
    M=VM/sqrt(k*R*Th);
    Tt=Th*(1+(k-1)/2*M^2);
    Pt=ph*(1+(k-1)/2*M^2)^(k/(k-1));
    V0=VM;
end
   
end

