% Script to generate rate of climb data for Warrior 3
% Use data on page 5-19 of the Pilot's Information Manual
VtrueKTS=[117.02013,116.75445,116.48888,116.18462,115.88029,115.49863,115.07826,114.50304,113.73432,112.69470,111.03595]; %Data fro PIM page 5-20b
% True maximum level flight airspeed knots
hp=0:1000:10000; % pressure altitude, feet
W=2300*ones(size(hp)); 

[temp,press,rho,Hgeopvector]=atmosphere4(hp,0)
TdegF=temp-459.67


% Write data to command window
n=length(hp);
str2=['% Hp(ft)   OAT (degF)    Wt(lbf) Vtrue(knots)  for Warrior 3'];
disp(' '); disp(str2)
for i=1:n
    strx=[num2str(hp(i)),'       ',num2str(TdegF(i)),'       ',num2str(W(i)),'       ',num2str(VtrueKTS(i))];
    disp(strx)
end

disp(' '); disp('Copy the above data into a text file called Warrior3SPEEDdata.txt.')
disp('This file is the input file for Step_3_Speed_Warrior3.m.')
format loose