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
