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
    cmhot.set_under('white')
    cmhot.set_over('black')

    cmjet = mpl.cm.jet
    cmjet = mpl.cm.RdYlBu_r
    cmjet.set_bad('white')
    cmjet.set_under('white')
    cmjet.set_over('black')

    topbounds = [0.1,0.50,0.8,0.45]
    bottombounds = [0.1,0.05,0.8,0.464]
    # centerx, centery, radius, width, height

    pl.figure(1,figsize=(12,12))
    pl.clf()
    dens1 = FITSFigure(parcubefile,convention='calabretta',slices=[0],figure=pl.figure(1),subplot=topbounds)
    dens1.show_colorscale(vmin=2,vmax=6,cmap=cmhot)
    dens1.recenter(**zoomargs)
    dens1.colorbar._colorbar_axes.set_ylabel('log$_{10}$(n(H$_2$) cm$^{-3}$)')
    dens1.hide_xaxis_label()
    dens1.hide_xtick_labels()
    velo1 = FITSFigure(parcubefile,convention='calabretta',slices=[4],figure=pl.figure(1),subplot=bottombounds)
    velo1.show_colorscale(vmin=52,vmax=66,cmap=cmjet)
    velo1.recenter(**zoomargs)
    velo1.colorbar._colorbar_axes.set_ylabel('Velocity ($V_{LSR}$ km s$^{-1}$)')
    savefig('W51_H2CO_2parfittry11_v1_densityvelocity%s.png' % extrastr,bbox_inches='tight')

    dens1.remove_colorbar()
    dens1.hide_colorscale()
    col1 = FITSFigure(parcubefile,convention='calabretta',slices=[1],figure=pl.figure(1),subplot=topbounds)
    col1.show_colorscale(vmin=11,vmax=13.5,cmap=cmhot)
        tbl.add_column(Column(data=tbl['{0} Flux Over Threshold'.format(stat)] / total_co_slab_reg.value,
                              dtype='float', name='{0} F_total(slab)'.format(stat)))

    pkfrac = cube13ss.with_mask(cubemask).spectral_slab(vrange1[0], vrange2[1])[slices].sum().value / total_co_slab_reg.value
    tbl.meta['Peak Fraction'] = pkfrac

    tables[region] = tbl


if __name__ == "__main__":
    from aplpy_figure_maker import FITSFigure
    import aplpy
    import pylab as pl
    fig = pl.figure(1, figsize=(12,12))
    fig.clf()
    F1 = FITSFigure(cube13ss_slab3_masked_mom0.hdu, subplot=[0.05,0.5,0.4,0.4], figure=fig)
    F1.show_grayscale()
    F1._ax1.set_title("$^{13}$CO masked with H$_2$CO")
    F2 = FITSFigure(cube13_slab3_masked_mom0.hdu, subplot=[0.50,0.5,0.4,0.4], figure=fig)
    F2.show_grayscale()
    F2._ax1.set_title("$^{13}$CO masked by S/N")
    F3 = FITSFigure(high5e4dens_co_slab3_mom0.hdu, subplot=[0.05,0.05,0.4,0.4], figure=fig)
    F3.show_grayscale()
    F3._ax1.set_title("$^{13}$CO masked by $n>5\\times10^4$")
    F4 = FITSFigure(high1e4dens_co_slab3_mom0.hdu, subplot=[0.50,0.05,0.4,0.4], figure=fig)
    F4.show_grayscale()
    F4._ax1.set_title("$^{13}$CO masked by $n>1\\times10^4$")
    F13 = FITSFigure(high1e5dens_co_slab3_mom0.hdu, subplot=[0.05,0.05,0.4,0.4], figure=fig)
    F13.show_grayscale()

    pl.matplotlib.rc_file('pubfiguresrc')