add CD plot

Author: ufechner7

src/plotting.jl

@@ -567,9 +567,9 @@ function plot_polars(
             markersize = 5
         end
         if contains(label, "LLT")
-            label = replace(label, "e5" => raw"\cdot10^5")
-            label = replace(label, " " => "~")
-            label = replace(label, "LLT" => raw"\mathrm{LLT}{~}")
+            label = replace(label, "e5"  => raw"\cdot10^5")
+            label = replace(label, " "   => raw"~")
+            label = replace(label, "LLT" => raw"\mathrm{LLT}{~\,}")
             label = raw"$" * label * raw"$"
         else
             label = replace(label, "e5" => raw"\cdot10^5")
@@ -590,44 +590,54 @@ function plot_polars(
             axs[1, 1].set_ylim([-0.5, 2])
         end
         println(label)
-        title = raw"$C_L" * raw"$"
-        # axs[1, 1].set_title(L"C_L \textrm{ vs } %$angle_type")
+        title = raw"$C_L" * raw"$" * " vs $angle_type [°]"
+        axs[1, 1].set_title(title)
         axs[1, 1].set_xlabel("$angle_type [°]")
         axs[1, 1].set_ylabel(L"$C_L$")
         axs[1, 1].legend()
     end
-    
-#     # Plot CL vs angle
-#     plot!(res[1])
-#     for (i, (polar_data, label)) in enumerate(zip(polar_data_list, label_list))
-#         style = i ≤ n_solvers ? (:solid, :star, 7) : (:solid, :circle, 5)
-#         plot!(res[1], polar_data[1], polar_data[2],
-#               label=label, linestyle=style[1], marker=style[2], markersize=style[3])
-        
-#         # Limit y-range if CL > 10
-#         if maximum(polar_data[2]) > 10
-#             ylims!(res[1], (-0.5, 2.0))
-#         end
-#     end
-#     title!(res[1], L"C_L \textrm{ vs } %$angle_type")
-#     xlabel!(res[1], "$angle_type [deg]")
-#     ylabel!(res[1], L"C_L")
 
-#     # Plot CD vs angle
-#     plot!(res[2])
-#     for (i, (polar_data, label)) in enumerate(zip(polar_data_list, label_list))
-#         style = i ≤ n_solvers ? (:solid, :star, 7) : (:solid, :circle, 5)
-#         plot!(res[2], polar_data[1], polar_data[3],
-#               label=label, linestyle=style[1], marker=style[2], markersize=style[3])
-        
-#         # Limit y-range if CD > 10
-#         if maximum(polar_data[3]) > 10
-#             ylims!(res[2], (-0.2, 0.5))
-#         end
-#     end
-#     title!(res[2], L"C_D \textrm{ vs } %$angle_type")
-#     xlabel!(res[2], "$angle_type [deg]")
-#     ylabel!(res[2], L"C_D")
+    for (i, (polar_data, label)) in enumerate(zip(polar_data_list, label_list))
+        if i < n_solvers
+            linestyle = "-"
+            marker = "*"
+            markersize = 7
+        else
+            linestyle = "-"
+            marker = "."
+            markersize = 5
+        end
+        if contains(label, "LLT")
+            label = replace(label, "e5"  => raw"\cdot10^5")
+            label = replace(label, " "   => raw"~")
+            label = replace(label, "LLT" => raw"\mathrm{LLT}{~\,}")
+            label = raw"$" * label * raw"$"
+        else
+            label = replace(label, "e5" => raw"\cdot10^5")
+            label = replace(label, " " => "~")
+            label = replace(label, "VSM" => raw"\mathrm{VSM}")
+            label = raw"$" * label * raw"$"
+        end
+        axs[1, 2].plot(
+            polar_data[1],
+            polar_data[3],
+            label=label,
+            linestyle=linestyle,
+            marker=marker,
+            markersize=markersize,
+        )
+        # Limit y-range if CL > 10
+        if maximum(polar_data[2]) > 10
+            axs[1, 1].set_ylim([-0.5, 2])
+        end
+        println(label)
+        title = raw"$C_D" * raw"$" * " vs $angle_type [°]"
+        axs[1, 2].set_title(title)
+        axs[1, 2].set_xlabel("$angle_type [°]")
+        axs[1, 2].set_ylabel(L"$C_D$")
+        axs[1, 2].legend()
+    end
+    
 
 #     # Plot CS vs angle (if available)
 #     plot!(res[3])