load Figure7.mat
 
  
figure('position',[0,0,1200,500])
h1=subplot(121);hold on
h2=subplot(122);hold on
 
set(gcf,'CurrentAxes',h1)
semilogx( k_fst,   KE , 'linewidth', 1); %kinetic energy flux, generated by"fig7and*_fst_snapshot_and_vs_k"
set(gcf,'CurrentAxes',h2)
semilogx( k_fst,   PE , 'linewidth', 1); %kinetic energy flux
 
set(gcf,'CurrentAxes',h1 )
%     set(gca,'ylim',[10^4,10^8],'fontsize',15,'xlim',[2, 300],'fontsize',12,'xScale','log','yScale','log')
xlabel('$k$', 'fontsize',25,'interpreter','latex');
ylabel('$\bar{\Pi}_u(k)$', 'fontsize',25,'interpreter','latex');
set(gca,'fontsize',18,'xScale','log','xlim',[3 256],'ylim',[-1.8e-3 2.e-4]);grid on;box on
text(1,1e-4,'$(a)$','fontsize',30,'interpreter','latex')


set(gcf,'CurrentAxes',h2)
xlabel('$k$', 'fontsize',25,'interpreter','latex');
ylabel('$\bar{\Pi}_b(k)$', 'fontsize',25,'interpreter','latex');
 set(gca,'fontsize',18,'xScale','log','xlim',[3 256]);grid on;box on;
legend(lg,'location','northwest','box','off','interpreter','latex','fontsize',18)
text(1,9.5e-4,'$(b)$','fontsize',30,'interpreter','latex')

