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):
