Waterfall plots

Subscript to plot the fields as funtion of pump power using waterfall plots.

plotlimit = floor(Param.t_steps/2-Param.t_steps/18):floor(Param.t_steps/2+Param.t_steps/18);
powerplot = 1:powersteps;
FontSize  = 20;
if( Param.t_FWHM*1e9 / sqrt(2) > 1)
    plotName  = [int2str(Param.t_FWHM*1e9 / sqrt(2)) 'ns'];
elseif( Param.t_FWHM*1e12 / sqrt(2) > 1)
    plotName  = [int2str(Param.t_FWHM*1e12 / sqrt(2)) 'ps'];
else
    plotName  = [num2str(Param.t_FWHM / sqrt(2),'%10.3e\n') 's'];
end

    figure(1)
    clearvars plot1data;
    plot1data(1:size(powerplot,2),1:Param.z_steps) = ...
        squeeze( Signal.I.depleted(powerplot,:,end/2) );
    w1 = waterfall(axis.z.*1e3,power(powerplot)./1000,plot1data.*1e-9);
    title('$I_{signal}^{+dep.}\,\,(\tau=0,z)$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$z~[mm]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w1,'edgecolor','k');
    %zlim([0 2e11])
    view([-40 50])
    export_fig(['CenterPropagation-' plotName '-Signal-Depleted'], '-pdf', '-transparent');%%

    figure(2)
    clearvars plot2data;
    plot2data(1:size(powerplot,2),1:size(plotlimit,2)) = ...
        squeeze( Signal.I.depleted(powerplot,end,plotlimit) );
    w2 = waterfall(axis.t(plotlimit).*1e12,power(powerplot)./1000,plot2data.*1e-9);
    title('$I_{signal}^{+dep.}\,\,(\tau,z=z_0)\,$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$\tau~[ps]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w2,'edgecolor','k');
    %zlim([0 2e11])
    view([-140 50])
    export_fig(['CrystalEnd-' plotName '-Signal-Depleted'], '-pdf', '-transparent');

    figure(3)
    clearvars plot3data;
    plot3data(1:size(powerplot,2),1:size(plotlimit,2)) = ...
        squeeze( Pump.I.depleted(powerplot,end,plotlimit) );
    w3 = waterfall(axis.t(plotlimit).*1e12,power(powerplot)./1000,plot3data.*1e-9);
    title('$I_{pump}^{+dep.}\,\,(\tau,z=z_0)\,$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$\tau~[ps]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w3,'edgecolor','k');
    %zlim([0 2e11])
    view([-140 50])
    export_fig(['CrystalEnd-' plotName '-Pump-Depleted'], '-pdf', '-transparent');

    figure(4)
    clearvars plot4data;
    plot4data(1:size(powerplot,2),1:Param.z_steps) = ...
        squeeze( Signal.I.undepleted(powerplot,:,end/2) );
    w4 = waterfall(axis.z.*1e3,power(powerplot)./1000,plot4data.*1e-9);
    title('$I_{signal}^{-dep.}\,\,(\tau=0,z)$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$z~[mm]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w4,'edgecolor','k');
    %zlim([0 2e11])
    view([-40 50])
    export_fig(['CenterPropagation-' plotName '-Signal-Undepleted'], '-pdf', '-transparent');

    figure(5)
    clearvars plot5data;
    plot5data(1:size(powerplot,2),1:size(plotlimit,2)) = ...
        squeeze( Signal.I.undepleted(powerplot,end,plotlimit) );
    w5 = waterfall(axis.t(plotlimit).*1e12,power(powerplot)./1000,plot5data.*1e-9);
    title('$I_{signal}^{-dep.}\,\,(\tau,z=z_0)\,$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$\tau~[ps]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w5,'edgecolor','k');
    %zlim([0 2e11])
    view([-140 50])
    export_fig(['CrystalEnd-' plotName '-Signal-Undepleted'], '-pdf', '-transparent');

    figure(6)
    clearvars plot6data;
    plot6data(1:size(powerplot,2),1:size(plotlimit,2)) = ...
        squeeze( Pump.I.undepleted(powerplot,end,plotlimit) );
    w6 = waterfall(axis.t(plotlimit).*1e12,power(powerplot)./1000,plot6data.*1e-9);
    title('$I_{pump}^{-dep.}\,\,(\tau,z=z_0)\,$','Interpreter','latex','FontSize',FontSize+4);
    xlabel('$\tau~[ps]$','Interpreter','latex','FontSize',FontSize);
    ylabel('$P_p~[kW]$','Interpreter','latex','FontSize',FontSize);
    zlabel('$I~[GW/m^2]$','Interpreter','latex','FontSize',FontSize);
    set(w6,'edgecolor','k');
    %zlim([0 2e11])
    view([-140 50])
    export_fig(['CrystalEnd-' plotName '-Pump-Undepleted'], '-pdf', '-transparent');

Published with MATLAB® R2013a