%% Put Compare Heatmaps Together
for probii=8:4:16
    for std=[3 6 12]
        str1=sprintf('Zoom_compare_state0_probii%1.0f_std%1.0f_REP.csv', probii, std);
        str2=sprintf('Zoom_compare_state1_probii%1.0f_std%1.0f_REP.csv', probii, std);
        
        heatmap1=csvread(str1);
        heatmap2=csvread(str2);
        
        heatmap1_adj=heatmap1(heatmap1>.3);
        heatmap2_adj=heatmap2(heatmap2>.3);
        total1=sum(heatmap1_adj);
        total2=sum(heatmap2_adj);
        
        figure('units','normalized','position',[0 0 1 1])
        colormap jet
                
        subaxis(1,2,1, 'Spacing', 0.01, 'Padding', 0.05, 'Margin', 0.05, 'PaddingBottom', 0.12, 'PaddingTop', 0.075, 'PaddingLeft', 0.07, 'PaddingRight', 0.01);
        pcolor(heatmap1)
        caxis([-.1 .9])
        colorbar
        
        xlabels=0.25:.5:3;
        ylabels=150:20:400;
        
        set(gca,'TickDir','out')
        set(gca,'box','off')
        set(gca, 'XTick', 1.5:2:12.5);
        set(gca, 'XTickLabel', xlabels);
        set(gca, 'YTick', 1.5:4:51.5);
        set(gca, 'YTickLabel', ylabels);
        xtickangle(45)
        ytickangle(45)
        set(gca, 'FontSize', 20);
        % ylabel({'Inter-connectivity Density'; 'Compared to E-E Density'}, 'FontSize', 26)
        % xlabel('Inhibitory Synaptic Weight (nS)', 'FontSize', 26)
        ylabel('External Applied Current (pA)', 'FontSize', 26)
        strtit=sprintf('Sync: Post-Pre, State=0, Total=%1.2f', total1);
        title(strtit, 'FontSize', 30);
        
        subaxis(1,2,2, 'Spacing', 0.01, 'Padding', 0.05, 'Margin', 0.05, 'PaddingBottom', 0.12, 'PaddingTop', 0.075, 'PaddingRight', 0.01, 'PaddingLeft', 0.05);
        pcolor(heatmap2)
        caxis([-.1 .9])
        colorbar
        
        set(gca,'TickDir','out')
        set(gca,'box','off')
        set(gca, 'XTick', 1.5:2:20.5);
        set(gca, 'XTickLabel', xlabels);
        set(gca, 'YTick', []);
        set(gca, 'YTickLabel', ylabels);
        xtickangle(45)
        ytickangle(45)
        set(gca, 'FontSize', 20);
        % ylabel({'Inter-connectivity Density'; 'Compared to E-E Density'}, 'FontSize', 26)
        % xlabel('Inhibitory Synaptic Weight (nS)', 'FontSize', 26)
        % ylabel('External Applied Current (pA)', 'FontSize', 26)
        strtit=sprintf('Sync: Post-Pre, State=0, Total=%1.2f', total2);
        title(strtit, 'FontSize', 30);
        
        supAxis=[.0750 .10 .85 .85];
        str1= sprintf('Inhibitory Synaptic Weight (nS)');
        str3= sprintf('Connectivity Density: %1.2f; STD: %1.2f; Post-Pre', probii/100, std);
        [ax,h]=suplabel(str1, 'x', supAxis);
        set(h, 'FontSize', 26);
        [ax3,h3]=suplabel(str3, 't', supAxis);
        set(h3, 'FontSize', 20);
        
        
        
        str1=sprintf('Zoom_compare_bistabilitymeasure_probii%1.0f_std%1.0f_REP.png', probii, std);
        %     saveas(gcf, str1)
        set(gcf,'PaperPositionMode','auto')
        print(str1, '-dpng', '-r0');
        
        
    end
end