% test Gthurst2
disp(' ')
disp(' Start here')
disp(' ')
Vtrueknots=[150,150]
Hpress=[5000, 5000]
[TempdegR,press,rho,Hgeopvector]=atmosphere4(Hpress,0);
%TempdegR=
Nblades=2
BHP0=300
RPM0=2700
BAF=113
Din=7*12
C=.12
Z=.4
istractor=1

[BHP,thrust]=Gthrust2(Vtrueknots,Hpress,TempdegR,Nblades,BHP0,RPM0,BAF,Din,Z,istractor,C)

i=1
[BHP,thrust,Mtip]=Gthrust(Vtrueknots(i),Hpress(i),RPM0,TempdegR(i),Nblades)


Nblades=3
BAF=95
[BHP,thrust]=Gthrust2(Vtrueknots,Hpress,TempdegR,Nblades,BHP0,RPM0,BAF,Din,Z,istractor,C)

i=1
[BHP,thrust,Mtip]=Gthrust(Vtrueknots(i),Hpress(i),RPM0,TempdegR(i),Nblades)
