|
425 | 425 |
|
426 | 426 |
|
427 | 427 | % ------------------- ENTRAINMENT TIME SERIES PLOTS ------------------- % |
428 | | - % Total plume volume and change in plume volume |
429 | 428 | if strcmp(PULSE,'T') == 1 |
430 | 429 | str = sprintf('%s: Unsteady flow %g Hz',titlerun,FREQ); |
431 | 430 | elseif strcmp(PULSE,'F') == 1 |
432 | 431 | str = sprintf('%s: Steady flow',titlerun); |
433 | 432 | end |
434 | 433 |
|
| 434 | + % Total plume volume and change in plume volume |
435 | 435 | figPlumeVol = figure('Name','Plume Volume','visible',vis,'units',... |
436 | 436 | 'centimeters','outerposition',[0 0 33.33 18.75],'PaperPositionMode',... |
437 | 437 | 'auto','color','w'); |
|
442 | 442 | xlabel(axVol1,{'Time (s)'},'FontWeight','bold') |
443 | 443 | ylabel(axVol1,{'Volume (m^3)'},'FontWeight','bold') |
444 | 444 | axVol2 = subplot(2,1,2); |
445 | | - plot(time(2:length(plumevolume)),diff(plumevolume)) |
| 445 | + plot(time,[0 diff(plumevolume)]) |
446 | 446 | title(axVol2,{sprintf('Change in plume volume\n%s',str)},... |
447 | 447 | 'FontWeight','bold') |
448 | 448 | xlabel(axVol2,{'Time (s)'},'FontWeight','bold') |
|
462 | 462 | for t = 2:length(time) |
463 | 463 | coeff_all = load(sprintf('ecoeff_all_t%03d.txt',time(t))); |
464 | 464 | hs = scatter(time(t)*ones(1,length(coeff_all)),coeff_all,... |
465 | | - 'MarkerEdgeColor',[0.65 0.65 0.65]); |
| 465 | + 'MarkerEdgeColor',[0.65 0.65 0.65],'SizeData',16); |
466 | 466 | hold on |
467 | 467 | end |
468 | 468 | hs.DisplayName = 'Full plume'; |
469 | 469 | he1 = errorbar(time(2:end),avg_expn(2:end),std_expn(2:end),'r',... |
470 | | - 'LineWidth',6,'LineStyle','none','Marker','+','MarkerSize',10,... |
| 470 | + 'LineWidth',4,'LineStyle','none','Marker','+','MarkerSize',10,... |
471 | 471 | 'DisplayName','Expansion'); |
472 | 472 | he2 = errorbar(time(2:end),avg_entr(2:end),std_entr(2:end),'c',... |
473 | | - 'LineWidth',6,'LineStyle','none','Marker','+','MarkerSize',10,... |
| 473 | + 'LineWidth',4,'LineStyle','none','Marker','+','MarkerSize',10,... |
474 | 474 | 'DisplayName','Entrainment'); |
475 | 475 | he3 = errorbar(time(2:end),avg_coeff(2:end),std_coeff(2:end),'k',... |
476 | | - 'LineWidth',3,'Marker','+','MarkerSize',10,... |
| 476 | + 'LineWidth',1.5,'Marker','+','MarkerSize',10,... |
477 | 477 | 'DisplayName','Combined'); |
478 | 478 | tlCoeff1 = 'EEE averaged over full plume surface'; |
479 | 479 | tlCoeff2 = sprintf('%s',str); |
|
505 | 505 | for t = 2:length(time) |
506 | 506 | jcoeff_all = load(sprintf('jcoeff_all_t%03d.txt',time(t))); |
507 | 507 | hjs = scatter(time(t)*ones(1,length(jcoeff_all)),jcoeff_all,... |
508 | | - 'MarkerEdgeColor',[0.65 0.65 0.65]); |
| 508 | + 'MarkerEdgeColor',[0.65 0.65 0.65],'SizeData',16); |
509 | 509 | hold on |
510 | 510 | end |
511 | 511 | hjs.DisplayName = 'Jet region'; |
512 | 512 | hje1 = errorbar(time(2:end),avg_jexpn(2:end),std_jexpn(2:end),'r',... |
513 | | - 'LineWidth',6,'LineStyle','none','Marker','+','MarkerSize',10,... |
| 513 | + 'LineWidth',4,'LineStyle','none','Marker','+','MarkerSize',10,... |
514 | 514 | 'DisplayName','Expansion'); |
515 | 515 | hje2 = errorbar(time(2:end),avg_jentr(2:end),std_jentr(2:end),'c',... |
516 | | - 'LineWidth',6,'LineStyle','none','Marker','+','MarkerSize',10,... |
| 516 | + 'LineWidth',4,'LineStyle','none','Marker','+','MarkerSize',10,... |
517 | 517 | 'DisplayName','Entrainment'); |
518 | 518 | hje3 = errorbar(time(2:end),avg_jcoeff(2:end),std_jcoeff(2:end),'k',... |
519 | | - 'LineWidth',3,'Marker','+','MarkerSize',10,... |
| 519 | + 'LineWidth',1.5,'Marker','+','MarkerSize',10,... |
520 | 520 | 'DisplayName','Combined'); |
521 | 521 | tlJCoeff1 = 'EEE averaged over jet region'; |
522 | 522 | tlJCoeff2 = sprintf('%s',str); |
|
0 commit comments