def show_con(dgmf, fignum=3, label="$f(n>10^4$ cm$^{-3})$"):
figN = pl.figure(fignum)
figN.clf()
gray = copy.copy(pl.cm.gray)
gray.set_bad('black')
gray.set_under('black')
FF = FITSFigure(cube13ss_slab3_masked_mom0.hdu, figure=figN, colorbar=False,
cmap=gray, transparent_nan=False)
#FF = FITSFigure(dgmf1e4, figure=fig4)
#FF.show_grayscale()
FF.show_contour(dgmf, levels=[0.1,0.3,0.5,0.7,0.9], zorder=1000, smooth=1)
FF.add_colorbar()
cax = FF.colorbar._colorbar_axes
cax.collections[0].set_visible(False)
FF.colorbar = pl.colorbar(FF._layers['contour_set_1'], cax=cax)
FF.colorbar.set_label("Dense Fraction", rotation=270, labelpad=30)
for lines in FF.colorbar.lines:
lines.set_linewidth(76)
FF._ax1.set_title(label)
#FF.scalebar._scalebar.set_zorder(40)
#FF.refresh()
return FF
fig7 = pl.figure(7)
fig7.clf()
gray = copy.copy(pl.cm.gray_r)
gray.set_bad('white')
gray.set_under('white')
FMM = FITSFigure(sfmasshdu1e4, figure=fig7, cmap=gray, stretch='log', vmax=99.95)
FMM.colorbar._colorbar.set_label("Star-Forming ($n>10^4$ cm$^{-3}$) Gas Surface Density\n[$M_{\odot}$ pc$^{-2}$]",
rotation=270, labelpad=50)
FMM.save(paths.fpath('StarFormingMassMap.pdf'))
FMM.show_markers(cl1lt10coords.l, cl1lt10coords.b, marker='x', edgecolor='r', zorder=1100, layer='kangx')
FMM.show_markers(cl1gt10coords.l, cl1gt10coords.b, marker='o', edgecolor='r', zorder=1100, layer='kango')
FMM.save(paths.fpath('StarFormingMassMap_Kang2009ClassIYSOs.pdf'))
FMM.remove_layer('kangx')
FMM.remove_layer('kango')
FMM.show_contour(sfmasshdu1e4, levels=[100,1000], zorder=1000, smooth=1, colors=['#AA1111','#FF4411','#BB6666'], linewidths=[2,3], alpha=0.7)
FMM.show_contour(sfmasshdu1e5, levels=[100,1000], zorder=1000, smooth=1, colors=['#1111BB','#1144FF','#6699BB'], linewidths=[2,3], alpha=0.9)
FMM.save(paths.fpath('StarFormingMassMap_thresholdcontours.pdf'))
fig8 = pl.figure(8)
fig8.clf()
ax8 = fig8.gca()
ax8.loglog(totalmassmap[~badsfmass1e4 & ok2cm], sfr2cmd[~badsfmass1e4 & ok2cm], 'r.', alpha=0.1)
ax8.loglog(sfmassmap1e4[~badsfmass1e4 & ok2cm], sfr2cmd[~badsfmass1e4 & ok2cm], 'k.', alpha=0.1)
ax8.loglog([1,1e4],np.array([1,1e4])*1.2e-8*1e6,'b--', alpha=0.5)
ax8.loglog([700,2900],[2,900],'g:', alpha=0.7) # by-eye fit to some of the data
ax8.set_xlabel("Star Forming Gas Surface Density ($M_{\\odot}$ pc$^{-2}$)")
ax8.set_ylabel("2 cm continuum SFR ($M_{\\odot}$ yr$^{-1}$ kpc$^{-2}$)")
ax8.set_ylim(0.5,1e3)
fig8.savefig(paths.fpath('ksrelation_2cm_densegas.pdf'))
