function   dVdHpersec=dVdH_Warrior3(ViVector,flapVector,HpVector)
% function   dVdHpersec=dVdH_Warrior3(ViVector,flapVector,HpVector)
% Compute ratio of the increment of true speed 
% divided by increment of altitude (1/sec)
% assuming constant indicated airspeed climb
% Input units: knots, nondimensionsl, geopotential feet
%                    flap=0   for no flap
%                    flap=0.5 for half flap
%                    flap=1   for full flap
% Output units: per second
% Input units are knots indicated, ND, ft

dH=10; % Increment of altitude (ft)
VcVector=AS1Warrior3(ViVector,flapVector);
HcVector=HpVector; % a vector same length as Hp2
[tVector,pVector,rho,Hgeopvector]=atmosphere4(HcVector,0); % altitude is geopotential, p is a vector
[VeVector,qcOverp0,qcOverp]=AS2(VcVector,pVector);
[VtrueKnots,M,aknots,rho,sigma]=AS3(VeVector,tVector,HcVector);
VtFtperSec=VtrueKnots*1.687808; % Convert knots to ft/sec



HcVector2=HpVector+dH; % a vector same length as Hp2
[tVector2,pVector2,rho2,Hgeopvector2]=atmosphere4(HcVector2,0); % altitude is geopotential, p is a vector
[VeVector2,qcOverp02,qcOverp2]=AS2(VcVector,pVector2);
[VtrueKnots2,M2,aknots2,rho2,sigma2]=AS3(VeVector2,tVector2,HcVector2);
VtFtperSec2=VtrueKnots2*1.687808; % Convert knots to ft/sec
dVdHpersec=(VtFtperSec2-VtFtperSec)/dH; % Increment of true speed divided by increment of altitude (1/sec)
