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
>>>>>>> master #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