Skip to content

Commit e11a436

Browse files
committed
added CL plot
1 parent 768a5c5 commit e11a436

File tree

1 file changed

+40
-0
lines changed

1 file changed

+40
-0
lines changed

src/plotting.jl

Lines changed: 40 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -556,6 +556,46 @@ function plot_polars(
556556

557557
# Number of computational results (excluding literature)
558558
n_solvers = length(solver_list)
559+
for (i, (polar_data, label)) in enumerate(zip(polar_data_list, label_list))
560+
if i < n_solvers
561+
linestyle = "-"
562+
marker = "*"
563+
markersize = 7
564+
else
565+
linestyle = "-"
566+
marker = "."
567+
markersize = 5
568+
end
569+
if contains(label, "LLT")
570+
label = replace(label, "e5" => raw"\cdot10^5")
571+
label = replace(label, " " => "~")
572+
label = replace(label, "LLT" => raw"\mathrm{LLT}{~}")
573+
label = raw"$" * label * raw"$"
574+
else
575+
label = replace(label, "e5" => raw"\cdot10^5")
576+
label = replace(label, " " => "~")
577+
label = replace(label, "VSM" => raw"\mathrm{VSM}")
578+
label = raw"$" * label * raw"$"
579+
end
580+
axs[1, 1].plot(
581+
polar_data[1],
582+
polar_data[2],
583+
label=label,
584+
linestyle=linestyle,
585+
marker=marker,
586+
markersize=markersize,
587+
)
588+
# Limit y-range if CL > 10
589+
if maximum(polar_data[2]) > 10
590+
axs[1, 1].set_ylim([-0.5, 2])
591+
end
592+
println(label)
593+
title = raw"$C_L" * raw"$"
594+
# axs[1, 1].set_title(L"C_L \textrm{ vs } %$angle_type")
595+
axs[1, 1].set_xlabel("$angle_type [°]")
596+
axs[1, 1].set_ylabel(L"$C_L$")
597+
axs[1, 1].legend()
598+
end
559599

560600
# # Plot CL vs angle
561601
# plot!(res[1])

0 commit comments

Comments
 (0)