|
8 | 8 | % entrainment3D ---does things--- |
9 | 9 | % |
10 | 10 | % Functions called: loadTimestep3D; pulsetitle |
11 | | -% Last edit: Taryn Black, 11 April 2016 |
| 11 | +% Last edit: Taryn Black, 15 April 2016 |
12 | 12 |
|
13 | 13 | % Clear directory of appending files from previous processing attempts |
14 | 14 | cd(savepath) |
|
437 | 437 | 'auto','color','w'); |
438 | 438 | axVol1 = subplot(2,1,1); |
439 | 439 | plot(time,plumevolume) |
440 | | - title(axVol1,{sprintf('%s: Total plume volume',str)},... |
| 440 | + title(axVol1,{sprintf('Total plume volume\n%s',str)},... |
441 | 441 | 'FontWeight','bold') |
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 | 445 | plot(time(2:length(plumevolume)),diff(plumevolume)) |
446 | | - title(axVol2,{sprintf('%s: Change in plume volume',str)},... |
| 446 | + title(axVol2,{sprintf('Change in plume volume\n%s',str)},... |
447 | 447 | 'FontWeight','bold') |
448 | 448 | xlabel(axVol2,{'Time (s)'},'FontWeight','bold') |
449 | 449 | ylabel(axVol2,{'\DeltaVolume (m^3)'},'FontWeight','bold') |
|
474 | 474 | 'DisplayName','Entrainment'); |
475 | 475 | he3 = errorbar(time(2:end),avg_coeff(2:end),std_coeff(2:end),'k',... |
476 | 476 | 'LineWidth',3,'Marker','+','MarkerSize',10,... |
477 | | - 'DisplayName','Total Coefficient'); |
478 | | - tlCoeff1 = 'Plume-averaged entrainment/expansion efficiencies'; |
479 | | - tlCoeff2 = sprintf('%s, overall average entrainment = %.4f',str,mean(avg_entr(2:end))); |
480 | | - title(axCoeff,{tlCoeff1;tlCoeff2},'FontWeight','bold') |
| 477 | + 'DisplayName','Combined'); |
| 478 | + tlCoeff1 = 'EEE averaged over full plume surface'; |
| 479 | + tlCoeff2 = sprintf('%s',str); |
| 480 | + tlCoeff3 = sprintf('Spatiotemporally averaged entrainment = %.4f',mean(avg_entr(2:end))); |
| 481 | + title(axCoeff,{tlCoeff1;tlCoeff2;tlCoeff3},'FontWeight','bold') |
481 | 482 | xlabel(axCoeff,'\bfTime (s)') |
482 | | - ylabel(axCoeff,'\bfEfficiency Coefficient') |
| 483 | + ylabel(axCoeff,'\bfEEE') |
483 | 484 | xlim(axCoeff,[0 time(end)+DT]) |
484 | 485 | ylim(axCoeff,[-1 1]) |
485 | 486 | line(xlim,[0 0],'color',[0.1 0.1 0.1],'LineWidth',0.5); |
|
516 | 517 | 'DisplayName','Entrainment'); |
517 | 518 | hje3 = errorbar(time(2:end),avg_jcoeff(2:end),std_jcoeff(2:end),'k',... |
518 | 519 | 'LineWidth',3,'Marker','+','MarkerSize',10,... |
519 | | - 'DisplayName','Total Coefficient'); |
520 | | - tlJCoeff1 = 'Jet region entrainment/expansion efficiencies'; |
521 | | - tlJCoeff2 = sprintf('%s, overall average entrainment = %.4f',str,mean(avg_jentr(2:end))); |
522 | | - title(axJCoeff,{tlJCoeff1;tlJCoeff2},'FontWeight','bold') |
| 520 | + 'DisplayName','Combined'); |
| 521 | + tlJCoeff1 = 'EEE averaged over jet region'; |
| 522 | + tlJCoeff2 = sprintf('%s',str); |
| 523 | + tlJCoeff3 = sprintf('Spatiotemporally averaged entrainment = %.4f',mean(avg_jentr(2:end))); |
| 524 | + title(axJCoeff,{tlJCoeff1;tlJCoeff2;tlJCoeff3},'FontWeight','bold') |
523 | 525 | xlabel(axJCoeff,'\bfTime (s)') |
524 | | - ylabel(axJCoeff,'\bfEntrainment/Expansion Coefficient') |
| 526 | + ylabel(axJCoeff,'\bfEEE') |
525 | 527 | xlim(axJCoeff,[0 time(end)+DT]) |
526 | 528 | ylim(axJCoeff,[-1 1]) |
527 | 529 | line(xlim,[0 0],'color',[0.1 0.1 0.1],'LineWidth',0.5); |
|
541 | 543 | 'PaperPositionMode','auto','color','w'); |
542 | 544 | axMor = axes('Parent',figMorton,'box','on','FontSize',12); |
543 | 545 | plot(axMor,time,e_Mconic) |
544 | | - title(axMor,sprintf('%s: Morton conic entrainment coefficient',str),... |
| 546 | + title(axMor,sprintf('Morton entrainment coefficient\n%s',str),... |
545 | 547 | 'FontWeight','bold') |
546 | 548 | xlabel(axMor,'\bfTime (s)') |
547 | 549 | ylabel(axMor,'\bfCoefficient') |
|
0 commit comments