def get_data(CL):
make_avl_file()
output_avl = lift_distribution(CL)
x_pos, cl, cdi = get_correct_data(output_avl)
# print(x_pos, cl, cdi)
PolyFitCurveCl = interp1d(x_pos, cl, kind="cubic", fill_value="extrapolate")
PolyFitCurveidrag = interp1d(x_pos, cdi, kind='cubic', fill_value='extrapolate')
return PolyFitCurveCl, PolyFitCurveidrag
def calc_polars(CL):
output_avl = lift_distribution(CL)
x_pos, cl, cdi = get_correct_data(output_avl)
# Lift Code
PolyFitCurveCl = interpolate.interp1d(x_pos,
cl,
kind="cubic",
fill_value="extrapolate")
# Drag Code
PolyFitCurveidrag = interpolate.interp1d(x_pos,
cdi,
kind='cubic',
fill_value='extrapolate')
return PolyFitCurveCl, PolyFitCurveidrag
#print(V)
n_ult= 3.75
n = n_ult/1.5#4.7/1.5#4.21/1.5#4.4/1.5#4.4/1.5
W = 1536496 #- weights["W_F"]#1801946.31#1510125.47#1549762.26#1806203.58
def input_CL(S,V,rho,W):
input_CL = W/(0.5*rho*V**2*S)
return input_CL
print(input_CL(S, V, rho, W), "cl")
#print(input_CL(n, S, V, rho, W))
make_avl_file()
## Import File List:
output_avl = lift_distribution(input_CL(S, V, rho, W))
x_pos = get_correct_data(output_avl)[0]
#print(x_pos)
##Lift Code:
cl = get_correct_data(output_avl)[1]
PolyFitCurveCl = sp.interpolate.interp1d(x_pos, cl, kind="cubic", fill_value="extrapolate")
##Drag Code:
cdi = get_correct_data(output_avl)[2]
PolyFitCurveidrag = sp.interpolate.interp1d(x_pos, cdi, kind='cubic', fill_value='extrapolate')
#print(cdi)
>>>>>>> master
### Define your functions at the beginning of the program
def c(z):
c = Cr - ((Cr - Ct) / (b / 2)) * z
return c