-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot_CD.m
421 lines (383 loc) · 35 KB
/
plot_CD.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
function [] = plot_CD(wingtype, sel_inflation, sel_inflation_double, exp_value, exp_value_hard, exp_value_double, sel_speed, div, chord, kin_viscosity, drag_dir, plot_type, format)
plot_variable = 'CD';
plot_variable_printed_name = 'C_{D}';
double_vel = 0.25;
if wingtype == "hard"
for j = 1:length(sel_speed) % looping over flow speed to create fixed-speed plots
Re = sel_speed(j) * chord / kin_viscosity; % Reynolds number
fprintf('plotting for Re = %.2d (vel = %.2d m/s) \n', Re, sel_speed(j))
figure('Position', [200, 200, 1000, 1000])
if plot_type == "title"
title([plot_variable_printed_name, ' plot # ', num2str(j), '; Flow Speed: ', num2str(sel_speed(j))],'fontweight','bold','fontsize', 24)
end
hold on
grid on
xlabel('AoA [ ˚ ]','fontweight','bold','fontsize', 30);
ylabel(strcat(plot_variable_printed_name, ' [ ]'),'fontweight','bold','fontsize', 30);
ax = gca;
ax.XAxis.LineWidth = 2;
ax.YAxis.LineWidth = 2;
ax.XAxis.TickValues = [-5 0 5 7.5 10 12.5 15];
xlim([-7.5 17.5])
ylim([-1.3 2])
clear k1
for k = 1:length(exp_value.f_avg)
if (exp_value.vel(k) == sel_speed(j))
if exist('k1','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', 'neutral', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#B02013', MarkerEdgeColor = '#B02013')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#B02013', MarkerEdgeColor = '#B02013')
x_vec = [exp_value.aoa(k1), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k1, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--b', 'HandleVisibility','off')
end
k1 = k;
end
end
str_annotation = sprintf('Re = %.2e', Re);
annotation('textbox', [0.696 0.77 0.1 0.1], 'String', str_annotation, 'BackgroundColor','white','LineStyle','-','Fontsize', 16, 'Interpreter','latex' ) % printing Re on plots
hold off
save_plot(gcf, wingtype, plot_type, plot_variable, j, sel_speed(j), format)
disp(sel_speed(j))
end
elseif wingtype == "soft"
for j = 1:length(sel_speed) % looping over flow speed to create fixed-speed plots
Re = sel_speed(j) * chord / kin_viscosity; % Reynolds number
fprintf('plotting for Re = %.2d (vel = %.2d m/s) \n', Re, sel_speed(j))
clear k1 k2 k3 k4 k5
figure('Position', [200, 200, 1000, 1000])
if plot_type == "title"
title([plot_variable_printed_name, ' plot # ', num2str(j), '; Flow Speed: ', num2str(sel_speed(j))],'fontweight','bold','fontsize', 24)
end
legend('Location','north','Orientation','horizontal','fontsize', 22, 'LineWidth', 2)
hold on
grid on
xlabel('AoA [deg]','fontweight','bold','fontsize', 30);
ylabel(strcat(plot_variable_printed_name, ' [ ]'),'fontweight','bold','fontsize', 30);
ax = gca;
ax.XAxis.TickValues = [-5 0 5 7.5 10 12.5 15];
xlim([-7.5 17.5])
ylim([-.1 1.1])
ax.XAxis.LineWidth = 2;
ax.YAxis.LineWidth = 2;
for k = 1:length(exp_value.f_avg)
if (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(1))
if exist('k1','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', 'neutral', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
x_vec = [exp_value.aoa(k1), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k1, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--r', 'HandleVisibility','off')
end
k1 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(2))
if exist('k2','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', '30 mL', 'CapSize', 18, 'MarkerFaceColor', 'k', 'LineWidth', 1, MarkerEdgeColor = 'black')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'MarkerFaceColor', 'k', 'LineWidth', 1, MarkerEdgeColor = 'black')
x_vec = [exp_value.aoa(k2), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k2, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--k', 'HandleVisibility','off')
end
k2 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(3))
if exist('k3','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', '60 mL', 'CapSize', 18, 'MarkerFaceColor', 'm', 'LineWidth', 1, MarkerEdgeColor = 'magenta')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'MarkerFaceColor', 'm', 'LineWidth', 1, MarkerEdgeColor = 'magenta')
x_vec = [exp_value.aoa(k3), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k3, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--m', 'HandleVisibility','off')
end
k3 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(4))
if exist('k4','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', '90 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
x_vec = [exp_value.aoa(k4), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k4, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--b', 'HandleVisibility','off')
end
k4 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(5))
if exist('k5','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'DisplayName', '120 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#A47E3B', MarkerEdgeColor = '#A47E3B')
x_vec = [exp_value.aoa(k5), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k5, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '--g', 'HandleVisibility','off')
end
k5 = k;
end
end
hold off
%saveas(gcf, ['./pic_notitle_paper/CD_plot/','/CD_plot_#', num2str(j), '_flow_speed_0_', num2str(100 * sel_speed(j))], 'png'); % saving plots in desired folder
%saveas(gcf, ['./pic/CD_plot/','/CD_plot_#', num2str(j), '_flow_speed_0_', num2str(100 * sel_speed(j))], 'png'); % saving plots in desired folder
save_plot(gcf, wingtype, plot_type, plot_variable, j, sel_speed(j), format)
end
elseif wingtype == "soft_hard"
for j = 1:length(sel_speed) % looping over flow speed to create fixed-speed plots
Re = sel_speed(j) * chord / kin_viscosity; % Reynolds number
fprintf('plotting for Re = %.2d (vel = %.2d m/s) \n', Re, sel_speed(j))
clear k1 k2 k3 k4 k5 hard1
figure('Position', [200, 200, 1000, 1000])
set(gcf, 'InnerPosition', [1, 1, 1440, 800])
set(gcf, 'OuterPosition', [1, 1, 1440, 880])
if plot_type == "title"
title([plot_variable_printed_name, ' plot # ', num2str(j), '; Flow Speed: ', num2str(sel_speed(j))],'fontweight','bold','fontsize', 24)
end
legend('Location','north','Orientation','horizontal','fontsize', 40, 'LineWidth', 2)
% [h,icons] = legend('Location','north','Orientation','horizontal','fontsize', 35, 'LineWidth', 2)
% set(icons,'MarkerSize',20)
hold on
grid on
xlabel('AoA [ ˚ ]','fontweight','bold','fontsize', 40);
ylabel(strcat(plot_variable_printed_name, ' [ ]'),'fontweight','bold','fontsize', 40);
ax = gca;
ax.FontSize = 40;
ax.FontName = 'Times New Roman';
ax.XAxis.LineWidth = 2;
ax.YAxis.LineWidth = 2;
ax.XAxis.TickValues = [-5 0 5 7.5 10 12.5 15];
xlim([-6 16])
ylim([0 1.3])
ax.InnerPosition = ([0.1014,0.137,0.88,0.83]);
for k = 1:length(exp_value.f_avg)
if k <= length(exp_value_hard.vel)
if (exp_value_hard.vel(k) == sel_speed(j))
if exist('hard1','var') == 0
errorbar(exp_value_hard.aoa(k), (-exp_value_hard.f_avg(k, drag_dir) / div(1, j)), exp_value_hard.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'DisplayName', '–– rigid', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor='#383B3E')
else
errorbar(exp_value_hard.aoa(k), (-exp_value_hard.f_avg(k, drag_dir) / div(1, j)), exp_value_hard.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor='#383B3E')
x_vec = [exp_value_hard.aoa(hard1), exp_value_hard.aoa(k)];
y_vec = [(-exp_value_hard.f_avg(hard1, drag_dir) / div(1, j)), (-exp_value_hard.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, 'Color', '#383B3E', 'HandleVisibility', 'off')
end
hard1 = k;
end
end
% if sel_speed(j) == double_vel
% if k > length(exp_value_double.vel)
% continue
% end
% if exp_value_double.inflation(k) == sel_inflation_double(1)
% if exp_value_double.dir(k) == 1
% if exist('d1','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'DisplayName', '– – neutral', 'CapSize', 18, 'LineWidth', 1, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 1, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
% x_vec = [exp_value_double.aoa(d1), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d1, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#383B3E', 'HandleVisibility','off')
% end
% d1 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d1n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
% x_vec = [exp_value_double.aoa(d1n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d1n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#383B3E', 'HandleVisibility','off')
% end
% d1n = k;
% end
% elseif exp_value_double.inflation(k) == sel_inflation_double(2)
% if exp_value_double.dir(k) == 1
% if exist('d2','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#C9459A','DisplayName', '15 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#C9459A', MarkerEdgeColor = '#C9459A')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#C9459A','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#C9459A', MarkerEdgeColor = '#C9459A')
% x_vec = [exp_value_double.aoa(d2), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d2, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#C9459A', 'HandleVisibility','off')
% end
% d2 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d2n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#C9459A','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#C9459A', MarkerEdgeColor = '#C9459A')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#C9459A', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#C9459A', MarkerEdgeColor = '#C9459A')
% x_vec = [exp_value_double.aoa(d2n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d2n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#C9459A', 'HandleVisibility','off')
% end
% d2n = k;
% end
%
% elseif exp_value_double.inflation(k) == sel_inflation_double(3)
% if exp_value_double.dir(k) == 1
% if exist('d3','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','DisplayName', '30 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
% x_vec = [exp_value_double.aoa(d3), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d3, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#BC3F02', 'HandleVisibility','off')
% end
% d3 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d3n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
% x_vec = [exp_value_double.aoa(d3n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d3n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#BC3F02', 'HandleVisibility','off')
% end
% d3n = k;
% end
%
% elseif exp_value_double.inflation(k) == sel_inflation_double(4)
% if exp_value_double.dir(k) == 1
% if exist('d4','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ED9B00','DisplayName', '45 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ED9B00', MarkerEdgeColor = '#ED9B00')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ED9B00','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ED9B00', MarkerEdgeColor = '#ED9B00')
% x_vec = [exp_value_double.aoa(d4), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d4, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#ED9B00', 'HandleVisibility','off')
% end
% d4 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d4n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ED9B00','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ED9B00', MarkerEdgeColor = '#ED9B00')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ED9B00','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ED9B00', MarkerEdgeColor = '#ED9B00')
% x_vec = [exp_value_double.aoa(d4n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d4n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#ED9B00', 'HandleVisibility','off')
% end
% d4n = k;
% end
% elseif exp_value_double.inflation(k) == sel_inflation_double(5)
% if exp_value_double.dir(k) == 1
% if exist('d5','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','DisplayName', '60 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
% x_vec = [exp_value_double.aoa(d5), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d5, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#ebb81e', 'HandleVisibility','off')
% end
% d5 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d5n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
% x_vec = [exp_value_double.aoa(d5n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d5n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#ebb81e', 'HandleVisibility','off')
% end
% d5n = k;
% end
% elseif exp_value_double.inflation(k) == sel_inflation_double(6)
% if exp_value_double.dir(k) == 1
% if exist('d6','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','DisplayName', '90 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
% x_vec = [exp_value_double.aoa(d6), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d6, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#0A9C6B', 'HandleVisibility','off')
% end
% d6 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d6n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
% x_vec = [exp_value_double.aoa(d6n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d6n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#0A9C6B', 'HandleVisibility','off')
% end
% d6n = k;
% end
% elseif exp_value_double.inflation(k) == sel_inflation_double(7)
% if exp_value_double.dir(k) == 1
% if exist('d7','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','DisplayName', '120 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
% x_vec = [exp_value_double.aoa(d7), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d7, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#1a91c2', 'HandleVisibility','off')
% end
% d7 = k;
% elseif exp_value_double.dir(k) == 0
% if exist('d7n','var') == 0
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
% else
% errorbar(exp_value_double.aoa(k), (-exp_value_double.f_avg(k, drag_dir) / div(1, j)), exp_value_double.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
% x_vec = [exp_value_double.aoa(d7n), exp_value_double.aoa(k)];
% y_vec = [(-exp_value_double.f_avg(d7n, drag_dir) / div(1, j)), (-exp_value_double.f_avg(k, drag_dir) / div(1, j))];
% plot(x_vec, y_vec, '-.', 'Color', '#1a91c2', 'HandleVisibility','off')
% end
% d7n = k;
% end
%
% end
%
% else
if (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(1))
if exist('k1','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'DisplayName', '– – neutral', 'CapSize', 18, 'LineWidth', 2, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#383B3E', 'HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerFaceColor = '#383B3E', MarkerEdgeColor = '#383B3E')
x_vec = [exp_value.aoa(k1), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k1, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '-.', 'Color', '#383B3E', 'HandleVisibility','off')
end
k1 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(2))
if exist('k2','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','DisplayName', '30 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#BC3F02','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#BC3F02', MarkerEdgeColor = '#BC3F02')
x_vec = [exp_value.aoa(k2), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k2, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '-.', 'Color', '#BC3F02', 'HandleVisibility','off')
end
k2 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(3))
if exist('k3','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','DisplayName', '60 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#ebb81e','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#ebb81e', MarkerEdgeColor = '#ebb81e')
x_vec = [exp_value.aoa(k3), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k3, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '-.', 'Color', '#ebb81e', 'HandleVisibility','off')
end
k3 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(4))
if exist('k4','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','DisplayName', '90 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#0A9C6B','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#0A9C6B', MarkerEdgeColor = '#0A9C6B')
x_vec = [exp_value.aoa(k4), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k4, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '-.', 'Color', '#0A9C6B', 'HandleVisibility','off')
end
k4 = k;
elseif (exp_value.vel(k) == sel_speed(j)) && (exp_value.inflation(k) == sel_inflation(5))
if exist('k5','var') == 0
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','DisplayName', '120 mL', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
else
errorbar(exp_value.aoa(k), (-exp_value.f_avg(k, drag_dir) / div(1, j)), exp_value.ci(k, drag_dir) / div(1, j), 'o', 'Color', '#1a91c2','HandleVisibility','off', 'CapSize', 18, 'LineWidth', 2, MarkerSize = 9, MarkerFaceColor = '#1a91c2', MarkerEdgeColor = '#1a91c2')
x_vec = [exp_value.aoa(k5), exp_value.aoa(k)];
y_vec = [(-exp_value.f_avg(k5, drag_dir) / div(1, j)), (-exp_value.f_avg(k, drag_dir) / div(1, j))];
plot(x_vec, y_vec, '-.', 'Color', '#1a91c2', 'HandleVisibility','off')
end
k5 = k;
end
end
hold off
save_plot(gcf, wingtype, plot_type, plot_variable, j, sel_speed(j), format)
end
end
% close all
end