Skip to content

Commit 68d6a03

Browse files
committed
succesful plotting of stall model
1 parent 653ddfc commit 68d6a03

File tree

3 files changed

+23
-22
lines changed

3 files changed

+23
-22
lines changed

examples/testing_stall_model.jl

Lines changed: 21 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -99,24 +99,24 @@ path_cfd_lebesque = joinpath(
9999
"V3_CL_CD_RANS_Lebesque_2024_Rey_300e4.csv"
100100
)
101101

102-
# plot_polars(
103-
# [VSM, VSM_with_stall_correction],
104-
# [wing_aero_CAD_19ribs, wing_aero_CAD_19ribs],
105-
# [
106-
# "VSM CAD 19ribs",
107-
# "VSM CAD 19ribs , with stall correction",
108-
# "CFD_Lebesque Rey 30e5"
109-
# ];
110-
# literature_path_list=[path_cfd_lebesque],
111-
# angle_range=range(0, 25, length=25),
112-
# angle_type="angle_of_attack",
113-
# angle_of_attack=0,
114-
# side_slip=0,
115-
# yaw_rate=0,
116-
# Umag=10,
117-
# title="tutorial_testing_stall_model_n_panels_$(n_panels)_distribution_$(spanwise_panel_distribution)",
118-
# data_type=".pdf",
119-
# save_path=joinpath(save_folder, "polars"),
120-
# is_save=true,
121-
# is_show=true
122-
# )
102+
plot_polars(
103+
[VSM, VSM_with_stall_correction],
104+
[wing_aero_CAD_19ribs, wing_aero_CAD_19ribs],
105+
[
106+
"VSM CAD 19ribs",
107+
"VSM CAD 19ribs , with stall correction",
108+
"CFD_Lebesque Rey 30e5"
109+
];
110+
literature_path_list=[path_cfd_lebesque],
111+
angle_range=range(0, 25, length=25),
112+
angle_type="angle_of_attack",
113+
angle_of_attack=0,
114+
side_slip=0,
115+
yaw_rate=0,
116+
Umag=10,
117+
title="tutorial_testing_stall_model_n_panels_$(n_panels)_distribution_$(spanwise_panel_distribution)",
118+
data_type=".pdf",
119+
save_path=joinpath(save_folder, "polars"),
120+
is_save=true,
121+
is_show=true
122+
)
572 KB
Binary file not shown.

src/plotting.jl

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -557,8 +557,9 @@ function plot_polars(
557557
# Read all data first
558558
data = readdlm(path, ',')
559559
# Skip the header row by taking data from row 2 onwards
560+
@show data[1, :]
560561
data = data[2:end, :]
561-
push!(polar_data_list, [data[:,3], data[:,1], data[:,2]])
562+
push!(polar_data_list, [data[:,3], data[:,1], data[:,2], zeros(length(data[:,1]))])
562563
end
563564
end
564565

0 commit comments

Comments
 (0)