ax.append(pyplot.subplot(gs[0, 0:3])) ax.append(pyplot.subplot(gs[0, 3:6])) ax.append(pyplot.subplot(gs[1, 0:2])) ax.append(pyplot.subplot(gs[1, 2:4])) ax.append(pyplot.subplot(gs[1, 4:6])) ax.append(pyplot.subplot(gs[2, 0:2])) ax.append(pyplot.subplot(gs[2, 2:4])) ax.append(pyplot.subplot(gs[2, 4:6])) plot_indiv_comp(ax[0], ["1e0", "1e1"], ["090"], "Total", "dirty_nz", r"$n_{z}$", fontsize=fontsize) leg = ax[0].legend(fontsize=0.95 * fontsize, loc="lower left") leg.get_frame().set_linewidth(2) ax[0].set_title("Global SED") ax[0].set_xlim(2, 100) plot_converge_slice(ax[1], ["1e0", "1e1"], ["035.11", "151.99"], "090", "dirty_nz", r"$n_{z}$", fontsize=fontsize) ax[1].yaxis.tick_right() ax[1].yaxis.set_label_position("right") leg = ax[1].legend(fontsize=0.95 * fontsize, loc="lower left") leg.get_frame().set_linewidth(2) ax[1].set_title(r"Y Image Slice, $\theta = 090^\circ$") ax[1].set_xlim(2, 100) modname = "dirty_nbinz" cvals = ["2", "5", "10", "20", "50", "100"] custcmap = cubehelix.cmap(reverse=False, rot=1, start=0, startHue=320, sat=2) minmax_y_vals = np.array([1e-4, 20]) for i, cval in enumerate(cvals): tau = "1e0"
r'$\theta = 90^\circ$', r'$\theta = 180^\circ$'], fontsize=fontsize, loc='lower center') leg2.get_frame().set_linewidth(2) #leg = ax[0].legend(fontsize=0.85*fontsize, loc='lower left',ncol=2) #leg.get_frame().set_linewidth(2) ax[0].set_title('Global SED') ax[0].set_ylim([1e-4,1e2]) ax[0].set_xlabel('') ax[0].xaxis.set_ticklabels([]) plot_converge_slice(ax[1], ['1e0','1e1'], ['000.15','000.53','035.11','151.99'], '000', 'dirty_nphot', r'$N$', fontsize=fontsize) #leg = ax[1].legend(fontsize=0.85*fontsize, loc='lower left', ncol=2) #leg.get_frame().set_linewidth(2) ax[1].set_title(r'Y Image Slice, $\theta = 000^\circ$') plot_converge_slice(ax[2], ['1e0','1e1'], ['000.15','000.53','035.11','151.99'], '090', 'dirty_nphot', r'$N$', fontsize=fontsize) #leg = ax[2].legend(fontsize=0.85*fontsize, loc='lower left', ncol=2) #leg.get_frame().set_linewidth(2) ax[2].yaxis.set_ticklabels([]) ax[2].set_ylabel('') ax[2].set_title(r'Y Image Slice, $\theta = 090^\circ$')
ax.append(pyplot.subplot(gs[0,3:6])) ax.append(pyplot.subplot(gs[1,0])) ax.append(pyplot.subplot(gs[1,1])) ax.append(pyplot.subplot(gs[1,2])) ax.append(pyplot.subplot(gs[1,3])) ax.append(pyplot.subplot(gs[1,4])) ax.append(pyplot.subplot(gs[1,5])) plot_indiv_comp(ax[0], ['1e0','1e1'], ['000','180'], 'Total', 'dirty_nxy', r'$n_{xy}$', fontsize=fontsize) leg = ax[0].legend(fontsize=0.95*fontsize, loc='lower left') leg.get_frame().set_linewidth(2) ax[0].set_title('Global SED') ax[0].set_xlim(2,100) plot_converge_slice(ax[1], ['1e0','1e1'], ['035.11','151.99'], '000', 'dirty_nxy', r'$n_{xy}$', fontsize=fontsize) ax[1].yaxis.tick_right() ax[1].yaxis.set_label_position("right") leg = ax[1].legend(fontsize=0.95*fontsize, loc='lower left') leg.get_frame().set_linewidth(2) ax[1].set_title(r'Y Image Slice, $\theta = 000^\circ$') ax[1].set_xlim(2,100) modname = 'dirty_nbinxy' cvals = ['2','5','10','20','50','100'] custcmap = cubehelix.cmap(reverse = False, rot=1, start=0, startHue=320, sat = 2) minmax_y_vals = np.array([0.05,1]) for i, cval in enumerate(cvals):