col2 = FITSFigure(parcubefile,convention='calabretta',slices=[9],figure=pl.figure(2),subplot=topbounds)
col2.show_colorscale(vmin=11,vmax=13.5,cmap=cm_rgb) #cmhot)
col2.recenter(**zoomargs)
col2.colorbar._colorbar_axes.set_ylabel('log$_{10}$(N(H$_2$) cm$^{-2}$)')
col2.hide_xaxis_label()
col2.hide_xtick_labels()
savefig('W51_H2CO_2parfittry11_v2_columnvelocity%s.png' % extrastr,bbox_inches='tight')
pl.figure(3)
pl.clf()
pl.rc('font',size=24)
dens1 = FITSFigure(parcubefile,convention='calabretta',slices=[0],figure=pl.figure(3))
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.set_tick_labels_format('dd.d','d.dd')
#dens1.hide_xaxis_label()
#dens1.hide_xtick_labels()
dens1.refresh()
savefig('W51_H2CO_2parfittry11_v1_justdensity%s.png' % extrastr,bbox_inches='tight')
dens1.show_markers(np.array(ysos['_Glon'][ysos['Cl1']=='I']),
np.array(ysos['_Glat'][ysos['Cl1']=='I']),
edgecolor='k',marker='x')
dens1.show_markers(np.array(masers['GLON']),
np.array(masers['GLAT']),
marker='+',edgecolor='k')
dens1.show_markers(glon,glat,marker='+',edgecolor='k')
dens1.refresh()
savefig('W51_H2CO_2parfittry11_v1_justdensity%s_withYSOs.png' % extrastr,bbox_inches='tight')
column1 = 10**parcubefile[0].data[1] * parcubefile[0].data[5]
