%% loading the data load('fig2.mat'), %% Figure Frame Setup figure(2), set(gcf,'position',[4 100 1904 672]); font_size = 13, colormap([hsv(128);jet(128)]); %% Subplot Setup A_1 = subplot('position',[0.056 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_2 = subplot('position',[0.056 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_3 = subplot('position',[0.056 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_4 = subplot('position',[0.056 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_5 = subplot('position',[0.056 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_6 = subplot('position',[0.056 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_7 = subplot('position',[0.056 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); A_8 = subplot('position',[0.056 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_1 = subplot('position',[0.156 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_2 = subplot('position',[0.156 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_3 = subplot('position',[0.156 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_4 = subplot('position',[0.156 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_5 = subplot('position',[0.156 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_6 = subplot('position',[0.156 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_7 = subplot('position',[0.156 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); B_8 = subplot('position',[0.156 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_1 = subplot('position',[0.256 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_2 = subplot('position',[0.256 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_3 = subplot('position',[0.256 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_4 = subplot('position',[0.256 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_5 = subplot('position',[0.256 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_6 = subplot('position',[0.256 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_7 = subplot('position',[0.256 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); C_8 = subplot('position',[0.256 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_1 = subplot('position',[0.356 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_2 = subplot('position',[0.356 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_3 = subplot('position',[0.356 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_4 = subplot('position',[0.356 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_5 = subplot('position',[0.356 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_6 = subplot('position',[0.356 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_7 = subplot('position',[0.356 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); D_8 = subplot('position',[0.356 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_1 = subplot('position',[0.456 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_2 = subplot('position',[0.456 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_3 = subplot('position',[0.456 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_4 = subplot('position',[0.456 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_5 = subplot('position',[0.456 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_6 = subplot('position',[0.456 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_7 = subplot('position',[0.456 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); E_8 = subplot('position',[0.456 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_1 = subplot('position',[0.556 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_2 = subplot('position',[0.556 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_3 = subplot('position',[0.556 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_4 = subplot('position',[0.556 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_5 = subplot('position',[0.556 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_6 = subplot('position',[0.556 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_7 = subplot('position',[0.556 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); F_8 = subplot('position',[0.556 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_1 = subplot('position',[0.656 0.834 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_2 = subplot('position',[0.656 0.726 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_3 = subplot('position',[0.656 0.640 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_4 = subplot('position',[0.656 0.499 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_5 = subplot('position',[0.656 0.381 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_6 = subplot('position',[0.656 0.273 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_7 = subplot('position',[0.656 0.122 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); G_8 = subplot('position',[0.656 0.041 0.09 0.05]); set(gca,'fontsize',font_size,'TickDir','out'); d = subplot('position',[0.844 0.781 0.15 0.204]); set(gca,'fontsize',font_size,'TickDir','out'); e = subplot('position',[0.844 0.557 0.15 0.204]); set(gca,'fontsize',font_size,'TickDir','out'); f = subplot('position',[0.833 0.078 0.161 0.384]); set(gca,'fontsize',font_size,'TickDir','out'); f1 = subplot('position',[1 0.232 0.097 0.22]); set(gca,'fontsize',font_size,'TickDir','out'); %% dx = diff(fig2(1).x); dx = dx(1); %% Plotting Column 1 k=1 axes(A_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); ylabel({['Rb'];['x (\mum)']}); title('t=0.06s','fontsize',font_size) set(gca,'position',[0.056 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(A_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); ylabel({['Cs'];['x (\mum)']}); set(gca,'position',[0.056 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(A_3); plot(fig2(1).z,fig2(8).uc1(3,:),'r','linewidth',2); hold on; plot(fig2(1).z,fig2(8).uc1(4,:),'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.056 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'TickDir','out','clim',[0 62]); ylabel({['n_{int}'];['(\mum^{-1})']},'fontsize',font_size+2) axes(A_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); ylabel({['Rb'];['x (\mum)']}); set(gca,'position',[0.056 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(A_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); ylabel({['Cs'];['x (\mum)']}); set(gca,'position',[0.056 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(A_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.056 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'TickDir','out','clim',[0 62]); ylabel({['n_{int}'];['(\mum^{-1})']},'fontsize',font_size+2) axes(A_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); ylabel({['Rb'];['x (\mum)']}); set(gca,'position',[0.056 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(A_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); ylabel({['Cs'];['x (\mum)']}); xlabel('z (\mum)'); set(gca,'position',[0.056 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); %% Plotting Column 2 k=2 axes(B_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=0.1s','fontsize',font_size) set(gca,'position',[0.156 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(B_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.156 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(B_3); plot(fig2(1).z,fig2(9).uc1(3,:),'r','linewidth',2); hold on; plot(fig2(1).z,fig2(9).uc1(4,:),'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.156 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(B_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.156 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(B_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.156 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(B_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.156 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(B_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.156 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(B_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.156 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 3 k=3 axes(C_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=0.25s','fontsize',font_size) set(gca,'position',[0.256 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(C_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.256 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(C_3); plot(fig2(1).z,sum(fig2(k).uc1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uc2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.256 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(C_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.256 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(C_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.256 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(C_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.256 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(C_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.256 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(C_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.256 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 4 k=4 axes(D_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=0.5s','fontsize',font_size) set(gca,'position',[0.356 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(D_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.356 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(D_3); plot(fig2(1).z,sum(fig2(k).uc1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uc2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.356 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(D_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.356 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(D_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.356 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(D_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.356 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(D_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.356 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(D_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.356 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 5 k=5 axes(E_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=0.7s','fontsize',font_size) set(gca,'position',[0.456 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(E_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.456 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(E_3); plot(fig2(1).z,sum(fig2(k).uc1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uc2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.456 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(E_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.456 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(E_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.456 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(E_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.456 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(E_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.456 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(E_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.456 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 6 k=6 axes(F_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=1.0s','fontsize',font_size) set(gca,'position',[0.556 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(F_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.556 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(F_3); plot(fig2(1).z,sum(fig2(k).uc1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uc2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.556 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(F_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.556 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(F_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.556 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(F_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.556 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(F_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.556 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(F_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.556 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 7 k=7 axes(G_1); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc1); title('t=2.0s','fontsize',font_size) set(gca,'position',[0.656 0.834 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); colorbar('position',[0.75 0.764 0.01 0.195],'fontsize',font_size,'ylim',[0 62]) axes(G_2); imagesc(fig2(1).z,fig2(1).x,fig2(k).uc2); set(gca,'position',[0.656 0.726 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); ylabel('n_{col} (\mum^{-2})','position',[67.158 -15.03 1.002]) axes(G_3); plot(fig2(1).z,sum(fig2(k).uc1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uc2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.656 0.640 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(G_4); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO1); set(gca,'position',[0.656 0.499 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); ylabel('n_{0,col} (\mum^{-2})','position',[67.158 11.216 1.002]) colorbar('position',[0.75 0.417 0.01 0.2],'fontsize',font_size,'ylim',[0 62]) axes(G_5); imagesc(fig2(1).z,fig2(1).x,fig2(k).uPO2); set(gca,'position',[0.656 0.381 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 1]*62); axes(G_6); plot(fig2(1).z,sum(fig2(k).uPO1,1)*dx,'r','linewidth',2); hold on; plot(fig2(1).z,sum(fig2(k).uPO2,1)*dx,'--b','linewidth',2); xlim([-1 1]*33); ylim([0 450]); set(gca,'position',[0.656 0.298 0.09 0.1],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'YtickLabel',[],'TickDir','out','clim',[0 62]); axes(G_7); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase1); set(gca,'position',[0.656 0.157 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'XtickLabel',[],'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); axes(G_8); imagesc(fig2(1).z,fig2(1).x,fig2(k).phase2); set(gca,'position',[0.656 0.041 0.09 0.156],'Xtick',[-1 -0.5 0 0.5 1]*30,'Ytick',[-1 0 1]*10,'YtickLabel',[],'TickDir','out','clim',[0 62]); axis equal; xlim([-1 1]*33); ylim([-12 12]); caxis([-1 3]*pi); xlabel('z (\mum)'); %% Plotting Column 8 axes(d); plot(fig2(1).NPO(:,1),fig2(1).NPO(:,2),'.-r','linewidt',2,'markersize',20); hold on; plot(fig2(1).NPO(:,1),fig2(1).NPO(:,3),'.--b','linewidt',2,'markersize',20); set(gca,'XtickLabel',[],'Ytick',[0 2.5 5.0 7.5]*1e3,'YtickLabel',[0 2.5 5.0 7.5]); ylabel('N_{0}(10^{3})'); ylim([0 8e3]); legend('Rb','Cs') axes(e); plot(fig2(1).r(:,1),fig2(1).r(:,2),'.-r','linewidt',2,'markersize',20); hold on; plot(fig2(1).r(:,1),fig2(1).r(:,3),'.--b','linewidth',2,'markersize',20); ylim([0 28]); xlabel('t (s)'); ylabel('r'); axes(f); errorbar(fig2(1).Nd(:,1),fig2(1).Nd(:,2),fig2(1).Nd(:,3),fig2(1).Nd(:,4),'linewidth',2,'Color',[0 0.5 0],'Marker','.','MarkerSize',20); ylabel('N_{sol}'); xlabel('t (s)'); axes(f1); errorbar(fig2(1).Ndt(:,1),fig2(1).Ndt(:,2),fig2(1).Ndt(:,3),fig2(1).Ndt(:,4),'LineStyle','none','Color',[0 0.5 0],'Marker','.','MarkerSize',20); set(gca,'position',[0.893 0.232 0.097 0.22],'Xscale','log','Yscale','log','Ytick',[1 5 10],'YtickLabel',[1 5 10],'Xtick',[10 100]); ylim([1 12]); xlim([6 700]); hold on; plot(fig2(2).Ndt(1:2),fig2(2).Ndt(3:4),'--k','linewidth',2); plot(fig2(3).Ndt(1:2),fig2(3).Ndt(3:4),':','linewidth',2,'color',[153 51 0]/255); plot(fig2(4).Ndt(1:2),fig2(4).Ndt(3:4),'-.','linewidth',2,'color',[11 132 199]/255); xlabel('t_{R} (ms)'); ylabel('\langleN_{sol}(t_{R})\rangle'); %% a = annotation('textbox',[0 0.926 0.05 0.068],'string',['(a)'],'FontSize',font_size+2,'LineStyle','none'); b = annotation('textbox',[0 0.572 0.05 0.068],'string',['(b)'],'FontSize',font_size+2,'LineStyle','none'); c = annotation('textbox',[0 0.233 0.05 0.068],'string',['(c)'],'FontSize',font_size+2,'LineStyle','none'); d = annotation('textbox',[0.847 0.908 0.05 0.068],'string',['(d)'],'FontSize',font_size+2,'LineStyle','none'); e = annotation('textbox',[0.847 0.684 0.05 0.068],'string',['(e)'],'FontSize',font_size+2,'LineStyle','none'); f = annotation('textbox',[0.836 0.394 0.05 0.061],'string',['(f)'],'FontSize',font_size+2,'LineStyle','none'); tR_13 = annotation('textbox',[0.936 0.369 0.053 0.08],'string',['t_{R}^{-1/3}'],'FontSize',font_size+2,'LineStyle','none'); tR_23 = annotation('textbox',[0.964 0.314 0.053 0.08],'string',['t_{R}^{-2/3}'],'FontSize',font_size+2,'LineStyle','none','Color',[153 51 0]/255); tR_12 = annotation('textbox',[0.925 0.266 0.053 0.08],'string',['t_{R}^{-1/2}'],'FontSize',font_size+2,'LineStyle','none','Color',[11 132 199]/255);