Fix example stall_model

.gitignore

@@ -1,2 +1,3 @@
 Manifest.toml
 .vscode/settings.json
+results/TUDELFT_V3_LEI_KITE/polars/$C_L$ vs $C_D$.pdf

examples/testing_stall_model.jl → examples/stall_model.jl

@@ -76,47 +76,48 @@ results = solve(VSM, wing_aero_CAD_19ribs)
 @time results_with_stall = solve(VSM_with_stall_correction, wing_aero_CAD_19ribs)
 @time results_with_stall = solve(VSM_with_stall_correction, wing_aero_CAD_19ribs)
 
-# CAD_y_coordinates = [panel.aerodynamic_center[2] for panel in wing_aero_CAD_19ribs.panels]
+CAD_y_coordinates = [panel.aerodynamic_center[2] for panel in wing_aero_CAD_19ribs.panels]
 
-# plot_distribution(
-#     [CAD_y_coordinates, CAD_y_coordinates],
-#     [results, results_with_stall],
-#     ["VSM", "VSM with stall correction"];
-#     title="CAD_spanwise_distributions_alpha_$(round(aoa, digits=1))_beta_$(round(side_slip, digits=1))_yaw_$(round(yaw_rate, digits=1))_Umag_$(round(Umag, digits=1))",
-#     data_type=".pdf",
-#     save_path=joinpath(save_folder, "spanwise_distributions"),
-#     is_save=false,
-#     is_show=true
-# )
+plot_distribution(
+    [CAD_y_coordinates, CAD_y_coordinates],
+    [results, results_with_stall],
+    ["VSM", "VSM with stall correction"];
+    title="CAD_spanwise_distributions_alpha_$(round(aoa, digits=1))_beta_$(round(side_slip, digits=1))_yaw_$(round(yaw_rate, digits=1))_Umag_$(round(Umag, digits=1))",
+    data_type=".pdf",
+    save_path=joinpath(save_folder, "spanwise_distributions"),
+    is_save=false,
+    is_show=true
+)
 
-# # Plotting polar
-# save_path = joinpath(root_dir, "results", "TUD_V3_LEI_KITE")
-# path_cfd_lebesque = joinpath(
-#     root_dir,
-#     "data",
-#     "TUDELFT_V3_LEI_KITE",
-#     "literature_results",
-#     "V3_CL_CD_RANS_Lebesque_2024_Rey_300e4.csv"
-# )
+# Plotting polar
+save_path = joinpath(root_dir, "results", "TUD_V3_LEI_KITE")
+path_cfd_lebesque = joinpath(
+    root_dir,
+    "data",
+    "TUDELFT_V3_LEI_KITE",
+    "literature_results",
+    "V3_CL_CD_RANS_Lebesque_2024_Rey_300e4.csv"
+)
 
-# plot_polars(
-#     [VSM, VSM_with_stall_correction],
-#     [wing_aero_CAD_19ribs, wing_aero_CAD_19ribs],
-#     [
-#         "VSM CAD 19ribs",
-#         "VSM CAD 19ribs , with stall correction",
-#         "CFD_Lebesque Rey 30e5"
-#     ];
-#     literature_path_list=[path_cfd_lebesque],
-#     angle_range=range(0, 25, length=25),
-#     angle_type="angle_of_attack",
-#     angle_of_attack=0,
-#     side_slip=0,
-#     yaw_rate=0,
-#     Umag=10,
-#     title="tutorial_testing_stall_model_n_panels_$(n_panels)_distribution_$(spanwise_panel_distribution)",
-#     data_type=".pdf",
-#     save_path=joinpath(save_folder, "polars"),
-#     is_save=true,
-#     is_show=true
-# )
+plot_polars(
+    [VSM, VSM_with_stall_correction],
+    [wing_aero_CAD_19ribs, wing_aero_CAD_19ribs],
+    [
+        "VSM CAD 19ribs",
+        "VSM CAD 19ribs , with stall correction",
+        "CFD_Lebesque Rey 30e5"
+    ];
+    literature_path_list=[path_cfd_lebesque],
+    angle_range=range(0, 25, length=25),
+    angle_type="angle_of_attack",
+    angle_of_attack=0,
+    side_slip=0,
+    yaw_rate=0,
+    Umag=10,
+    title="tutorial_testing_stall_model_n_panels_$(n_panels)_distribution_$(spanwise_panel_distribution)",
+    data_type=".pdf",
+    save_path=joinpath(save_folder, "polars"),
+    is_save=true,
+    is_show=true
+)
+nothing