print " D : ", AsUnit(Prop.D, 'in')
print " Pitch : ", AsUnit(Prop.Pitch, 'in')
Vmax = 50
h = 0 * FT
N = npy.linspace(1000, 6800, 6) * RPM
Alpha = npy.linspace(-25, 25, 41) * ARCDEG
V = npy.linspace(0, Vmax, 30) * FT / SEC
Prop.CoefPlot(Alpha, fig=1)
Prop.PTPlot(N, V, h, 'V', fig=2)
#
# N = npy.linspace(0, 13000,31)*RPM
# V = npy.linspace(0,Vmax,5)*FT/SEC
#
# Prop.PTPlot(N,V,h,'N', fig = 3)
Prop.PlotTestData(fig=4)
N = 2414 * RPM
print
print "Static Thrust : ", AsUnit(Prop.T(N, 0 * FT / SEC, h), 'lbf')
print "Measured Thrust : ", AsUnit(max(npy.array(Prop.ThrustData)[:, 1]),
'lbf')
print
print "Static Torque : ", AsUnit(
Prop.P(N, 0 * FT / SEC, h) / N, 'in*ozf')
# print "Measured Torque : ", AsUnit( max(npy.array(Prop.TorqueData)[:,1]), 'in*ozf' )
pyl.show()
print " Pitch : ", AsUnit(Prop.Pitch, 'in')
Vmax = 50
h = 0 * FT
N = npy.linspace(1000, 6800, 5) * RPM
Alpha = npy.linspace(-25, 25, 41) * ARCDEG
V = npy.linspace(0, Vmax, 30) * FT / SEC
Prop.CoefPlot(Alpha, fig=1)
Prop.PTPlot(N, V, h, 'V', fig=2)
# N = npy.linspace(0, 13000,31)*RPM
# V = npy.linspace(0,Vmax,5)*FT/SEC
# Prop.PTPlot(N,V,h,'N', fig = 3)
Prop.PlotTestData(fig=4)
N = 6024 * RPM
print 'Static Thrust : ', AsUnit(Prop.T(N, 0 * FT / SEC, h), 'lbf')
print 'Measured Thrust : ', AsUnit(max(npy.array(Prop.ThrustData)[:, 1]),
'lbf')
N = 6410 * RPM
print
print "Static Torque : ", AsUnit(
Prop.P(N, 0 * FT / SEC, h) / N, 'in*ozf')
print "Measured Torque : ", AsUnit(max(npy.array(Prop.TorqueData)[:, 1]),
'in*ozf')
pyl.show()
