%% Data Plotter for Fig. 3 % This script loads the data package, fig3.mat, which contains the (i) utz % the c-field density at x=y=0 for b1.1.1, b1.1.2, b1.1.6 and b1.1.4 for the % structure label 1, 2, 3 and 4 respectively; (ii) t and z are the vertical % axis and horizontal for pannel (a) to (d); (iii) upoA and upoB are the PO % integrated densiy for pannel (e) to (h) and the corresponding insets, % with structure label 1, 2, 3 and 4 as well. %% loading the data load('fig3.mat') %% Set Fontsize font_size = 14, %% Figure Frame Configuration % Set the position and the size for pannel (a) to (h) without the plots of % insetn, which will be dealt in the bottom. figure(3), a_1 = subplot('position',[0.090 0.732 0.18 0.215]); set(gca,'fontsize',font_size); a_2 = subplot('position',[0.285 0.732 0.18 0.215]); set(gca,'fontsize',font_size); b_1 = subplot('position',[0.520 0.732 0.18 0.215]); set(gca,'fontsize',font_size); b_2 = subplot('position',[0.715 0.732 0.18 0.215]); set(gca,'fontsize',font_size); c_1 = subplot('position',[0.090 0.464 0.18 0.215]); set(gca,'fontsize',font_size); c_2 = subplot('position',[0.285 0.464 0.18 0.215]); set(gca,'fontsize',font_size); d_1 = subplot('position',[0.520 0.464 0.18 0.215]); set(gca,'fontsize',font_size); d_2 = subplot('position',[0.715 0.464 0.18 0.215]); set(gca,'fontsize',font_size); e_1 = subplot('position',[0.101 0.225 0.4 0.142]); set(gca,'fontsize',font_size); f_1 = subplot('position',[0.52 0.228 0.4 0.142]); set(gca,'fontsize',font_size); g_1 = subplot('position',[0.101 0.06 0.4 0.142]); set(gca,'fontsize',font_size); h_1 = subplot('position',[0.52 0.06 0.4 0.142]); set(gca,'fontsize',font_size); %% Plotting the main plots axes(a_1); imagesc(fig3(1).z,fig3(1).t,fig3(1).utz(:,1:length(fig3(1).z))); ... axis xy; xlim([-35 35]); ylim([0 1.5]);... ylabel('t (s)','fontsize',font_size);... title('n_{Rb}^{z} (\mum^{-3})','fontsize',font_size);... set(gca,'TickDir','out','XtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(a_2); imagesc(fig3(1).z,fig3(1).t,fig3(1).utz(:,length(fig3(1).z)+1:end));... axis xy; xlim([-35 35]); ylim([0 1.5]);... title('n_{Cs}^{z} (\mum^{-3})','fontsize',font_size);... set(gca,'TickDir','out','XtickLabel',[],'YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(b_1); imagesc(fig3(2).z,fig3(2).t,fig3(2).utz(:,1:length(fig3(2).z))); ... axis xy; xlim([-35 35]); ylim([0 1.5]);... title('n_{Rb}^{z} (\mum^{-3})','fontsize',font_size);... set(gca,'TickDir','out','XtickLabel',[],'YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(b_2); imagesc(fig3(2).z,fig3(2).t,fig3(2).utz(:,length(fig3(2).z)+1:end));... axis xy; xlim([-35 35]); ylim([0 1.5]);... title('n_{Cs}^{z} (\mum^{-3})','fontsize',font_size);... set(gca,'TickDir','out','XtickLabel',[],'YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(c_1); imagesc(fig3(3).z,fig3(3).t,fig3(3).utz(:,1:length(fig3(3).z))); ... axis xy; xlim([-35 35]); ylim([0 1.5]);... ylabel('t (s)','fontsize',font_size);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(c_2); imagesc(fig3(3).z,fig3(3).t,fig3(3).utz(:,length(fig3(3).z)+1:end));... axis xy; xlim([-35 35]); ylim([0 1.5]);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(d_1); imagesc(fig3(4).z,fig3(4).t,fig3(4).utz(:,1:length(fig3(4).z))); ... axis xy; xlim([-35 35]); ylim([0 1.5]);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); axes(d_2); imagesc(fig3(4).z,fig3(4).t,fig3(4).utz(:,length(fig3(4).z)+1:end));... axis xy; xlim([-35 35]); ylim([0 1.5]);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','YtickLabel',[],'Xtick',[-30 -15 0 15 30],'Ytick',[0 0.5 1 1.5]); caxis([0 16.3]); caxis([0 16.3]); colorbar('position',[0.919 0.466 0.027 0.475],'Ytick',[0:2:16],'ylim',[0 16],'fontsize',font_size); axes(e_1); plot(fig3(1).z,fig3(1).upoA(1:length(fig3(1).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(1).upoA(length(fig3(1).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... ylabel('n_{0,int} (\mum^{-1})','fontsize',font_size);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[]); axes(f_1); plot(fig3(1).z,fig3(2).upoA(1:length(fig3(2).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(2).upoA(length(fig3(2).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[],'YtickLabel',[]); axes(g_1); plot(fig3(1).z,fig3(3).upoA(1:length(fig3(3).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(3).upoA(length(fig3(3).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... ylabel('n_{0,int} (\mum^{-1})','fontsize',font_size);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30]); axes(h_1); plot(fig3(1).z,fig3(4).upoA(1:length(fig3(3).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(4).upoA(length(fig3(3).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... %ylabel('n_{0,int} (\mum^{-1})','fontsize',font_size);... xlabel('z (\mum)','fontsize',font_size);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'YtickLabel',[]); %axes(h_1); plot(fig3(1).z,fig3(4).upoB(1:216),'r'); hold on; plot(fig3(1).z,fig3(4).upoB(217:end),'--b') %% The plotting for the insent e_2 = subplot('position',[1 0.297 0.131 0.065]); set(gca,'fontsize',font_size); axes(e_2); plot(fig3(1).z,fig3(1).upoB(1:length(fig3(2).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(1).upoB(length(fig3(2).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[],'YtickLabel',[],'position',[0.354 0.297 0.131 0.065]); f_2 = subplot('position',[1 0.296 0.131 0.065]); set(gca,'fontsize',font_size); axes(f_2); plot(fig3(1).z,fig3(2).upoB(1:length(fig3(2).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(2).upoB(length(fig3(2).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[],'YtickLabel',[],'position',[0.784 0.297 0.131 0.065]); g_2 = subplot('position',[1 0.132 0.131 0.065]); set(gca,'fontsize',font_size); axes(g_2); plot(fig3(1).z,fig3(3).upoB(1:length(fig3(2).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(3).upoB(length(fig3(2).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[],'YtickLabel',[],'position',[0.354 0.132 0.131 0.065]); h_2 = subplot('position',[1 0.132 0.131 0.065]); set(gca,'fontsize',font_size); axes(h_2); plot(fig3(1).z,fig3(4).upoB(1:length(fig3(2).z)),'r','linewidth',2);... hold on; plot(fig3(1).z,fig3(4).upoB(length(fig3(2).z)+1:end),'--b','linewidth',2);... xlim([-35 35]); ylim([0 450]);... set(gca,'TickDir','out','Xtick',[-30 -15 0 15 30],'XtickLabel',[],'YtickLabel',[],'position',[0.784 0.132 0.131 0.065]); %% Reset the frame of the figure and Add labels set(gcf,'position',[200 50 880 864]); a = annotation('textbox',[0.011 0.953 0.06 0.039],'string',['(a)'],'FontSize',font_size+2,'LineStyle','none'); b = annotation('textbox',[0.479 0.953 0.06 0.039],'string',['(b)'],'FontSize',font_size+2,'LineStyle','none'); c = annotation('textbox',[0.011 0.685 0.06 0.039],'string',['(c)'],'FontSize',font_size+2,'LineStyle','none'); d = annotation('textbox',[0.479 0.685 0.06 0.039],'string',['(d)'],'FontSize',font_size+2,'LineStyle','none'); e = annotation('textbox',[0.106 0.318 0.061 0.039],'string',['(e)'],'FontSize',font_size+2,'LineStyle','none'); f = annotation('textbox',[0.525 0.318 0.061 0.039],'string',['(f)'],'FontSize',font_size+2,'LineStyle','none'); g = annotation('textbox',[0.106 0.158 0.061 0.039],'string',['(g)'],'FontSize',font_size+2,'LineStyle','none'); h = annotation('textbox',[0.525 0.158 0.061 0.039],'string',['(h)'],'FontSize',font_size+2,'LineStyle','none'); %% Input the notation in (e) to (h) eN = annotation('textbox',[0.142 0.27 0.148 0.086],'string',fig3(1).text,'FontSize',font_size,'LineStyle','none'); fN = annotation('textbox',[0.562 0.27 0.148 0.086],'string',fig3(2).text,'FontSize',font_size,'LineStyle','none'); gN = annotation('textbox',[0.142 0.111 0.148 0.086],'string',fig3(3).text,'FontSize',font_size,'LineStyle','none'); hN = annotation('textbox',[0.562 0.111 0.148 0.086],'string',fig3(4).text,'FontSize',font_size,'LineStyle','none');