Skip to content

Commit ec46e3b

Browse files
committed
Improved makie plotting
1 parent 66ac77d commit ec46e3b

File tree

1 file changed

+29
-1
lines changed

1 file changed

+29
-1
lines changed

ext/VortexStepMethodMakieExt.jl

Lines changed: 29 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1037,7 +1037,7 @@ function VortexStepMethod.plot_combined_analysis(
10371037
set_axes_equal_makie!(ax_geo, panels; zoom=0.5)
10381038
axislegend(ax_geo; position=:lt)
10391039

1040-
# [1,2] Polar Data Surfaces
1040+
# [1,2] Polar Data Surfaces or Curves
10411041
if body_aero.panels[1].aero_model == POLAR_MATRICES
10421042
alphas = collect(deg2rad.(-5:0.3:25))
10431043
delta_tes = collect(deg2rad.(-5:0.3:25))
@@ -1062,6 +1062,30 @@ function VortexStepMethod.plot_combined_analysis(
10621062
wireframe!(ax, delta_tes, alphas, interp_matrix;
10631063
color=:blue, linewidth=0.5, transparency=true)
10641064
end
1065+
elseif body_aero.panels[1].aero_model == POLAR_VECTORS
1066+
alphas_deg = collect(-5:0.5:25)
1067+
alphas = deg2rad.(alphas_deg)
1068+
1069+
ax_cl_curve = Axis(fig[1, 2][1, 1];
1070+
title="Cl vs α",
1071+
xlabel="α [°]",
1072+
ylabel="Cl")
1073+
ax_cd_curve = Axis(fig[1, 2][1, 2];
1074+
title="Cd vs α",
1075+
xlabel="α [°]",
1076+
ylabel="Cd")
1077+
ax_cm_curve = Axis(fig[1, 2][1, 3];
1078+
title="Cm vs α",
1079+
xlabel="α [°]",
1080+
ylabel="Cm")
1081+
1082+
cl_vals = [body_aero.panels[1].cl_interp(a) for a in alphas]
1083+
cd_vals = [body_aero.panels[1].cd_interp(a) for a in alphas]
1084+
cm_vals = [body_aero.panels[1].cm_interp(a) for a in alphas]
1085+
1086+
lines!(ax_cl_curve, alphas_deg, cl_vals; color=:blue, linewidth=2)
1087+
lines!(ax_cd_curve, alphas_deg, cd_vals; color=:red, linewidth=2)
1088+
lines!(ax_cm_curve, alphas_deg, cm_vals; color=:green, linewidth=2)
10651089
end
10661090

10671091
# [2,1] Spanwise Distributions (3×3 grid)
@@ -1158,6 +1182,10 @@ function VortexStepMethod.plot_combined_analysis(
11581182
label=label_with_re, marker=:star5, markersize=12)
11591183
axislegend(ax_polar, position=:lt)
11601184

1185+
# Set column widths: left column wider for 3x3 grid
1186+
colsize!(fig.layout, 1, Relative(0.6))
1187+
colsize!(fig.layout, 2, Relative(0.4))
1188+
11611189
if is_show
11621190
display(fig)
11631191
end

0 commit comments

Comments
 (0)