load('dync_tauQ.mat') load('dync_main.mat') load('dync_lcoh.mat') alpha = 1.3 tc_factor = 0.1879; dtc_factor = 0.0259; t0 = 1/(5e-3*22*131.4*2*pi) * 1e3 % in ms that = sqrt(t0*tauQ) % in ms that2 = tauQ.^0.57 * t0^0.43; hbar = 1.054571596e-34; % actual unit wx = 2 * pi * 13, % Hz w_perp = 2 * pi * 131.4 lambda_r = w_perp / wx; mass = 38.17e-27; % kg aho = sqrt(hbar/wx/mass); % meter muf = 22 * w_perp / wx, % dimensionless gamma = 5e-3 % non-dimension xi_0 = (2 * muf)^(-0.5)* aho*1e6; xi_k = (tauQ*1e-3*wx).^0.25 * (4 * muf / gamma)^(-0.25) * aho*1e6; % in micro-meter xi_k2 = xi_0 * (tauQ / t0).^0.29; % Color Codes Color_code{1} = 'r'; Color_code{2} = 'b'; Color_code{3} = [1 1 0] * 3 / 4; Color_code{4} = [1 0 1] * 3 / 4; Color_code{5} = [0.3 0.75 0.93]; Color_code{6} = [0 0.5 0]; figure(8) for jj = 1 : length(tauQ) subplot(3,1,1) color_code = Color_code{jj}; t_compare = abs(time_tc_den{jj}-that(jj)); [t_compare,i_compare] = sort(t_compare); i_ramp = find(time_tc_den{jj} == (1-0.1879)*tauQ(jj)); if isempty(i_ramp) [nearest_tR,i_nearest] = min(abs(time_tc_den{jj}-(1-tc_factor)*tauQ(jj))); i_ramp = i_nearest; end shadedErrorBarG((time_tc_den{jj}-alpha*that(jj))/tauQ(jj),den_cen{jj},[dden_cen{jj};dden_cen{jj}],{'-','color',color_code},1); hold on plot([1 1]*((1-tc_factor)*tauQ(jj)-alpha*that(jj))/tauQ(jj),den_cen{jj}(i_ramp),'o','linewidth',2,'color',color_code,'MarkerFaceColor','w','handlevisibility','off') xlabel('$(t-t_c-1.3\hat{t})/\tau_Q$','interpreter','latex'); ylabel('$N_{cen}$','interpreter','latex') ylim([0 2e4]) xlim([-0.1 1.5]); set(gca,'linewidth',2,'fontsize',14); end for jj = 1: length(tauQ) color_code = Color_code{jj}; subplot(3,1,2) shadedErrorBarG((time_tc{jj}-alpha*that(jj))/tauQ(jj),lcoh{jj},[dlcoh{jj}';dlcoh{jj}']/xi_k(jj),{'-o','color',color_code,'MarkerFaceColor',color_code},1); hold on plot([1 1]*((1-tc_factor)*tauQ(jj)-alpha*that(jj))/tauQ(jj),[0 900],'--','linewidth',1,'color',color_code,'MarkerFaceColor','w','handlevisibility','off') xlabel('$(t-t_c-1.3\hat{t})/\tau_Q$','interpreter','latex'); ylabel('$l_{coh}/\hat{\xi}$','interpreter','latex') ylim([0 500]) xlim([-0.1 1.5]); set(gca,'linewidth',2,'fontsize',14); end for jj = 1 : length(tauQ) color_code = Color_code{jj}; t_compare = abs(t_den_x_front{jj}-tauQ(jj)*(1-tc_factor)); [t_compare,i_compare] = sort(t_compare); y = den_x_front{jj}'*aho*1e6/114; yramp = sum((y(i_compare(1:2)) .* t_compare(1:2))) / sum(t_compare(1:2)); subplot(3,1,3) % longitudinal plot((t_den_x_front{jj}-alpha*that(jj))/tauQ(jj),y,'-','linewidth',2,'color',color_code) hold on plot(( (1 - tc_factor)*tauQ(jj)-alpha*that(jj))/tauQ(jj),yramp,'o','color',color_code,'linewidth',2,'Markerfacecolor','w','handlevisibility','off') % transverse plot((t_den_rho_front{jj}-alpha*that(jj))/tauQ(jj),den_rho_front{jj}*aho*1e6/114*lambda_r,'--','linewidth',2,'color',color_code,'handlevisibility','off') xlim([-0.1 1.5]); xlabel('$(t-t_c-1.3\hat{t})/\tau_Q$','interpreter','latex'); ylabel('scaled density front') set(gca,'linewidth',2,'fontsize',14); end