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)
col1.recenter(**zoomargs)
col1.colorbar._colorbar_axes.set_ylabel('log$_{10}$(N(H$_2$) cm$^{-2}$)')
col1.hide_xaxis_label()
col1.hide_xtick_labels()
savefig('W51_H2CO_2parfittry11_v1_columnvelocity%s.png' % extrastr,bbox_inches='tight')
pl.figure(2,figsize=(12,12))
pl.clf()
dens2 = FITSFigure(parcubefile,convention='calabretta',slices=[8],figure=pl.figure(2),subplot=topbounds)
dens2.show_colorscale(vmin=2,vmax=6,cmap=cmhot)
dens2.recenter(**zoomargs)
dens2.colorbar._colorbar_axes.set_ylabel('log$_{10}$(n(H$_2$) cm$^{-3}$)')
dens2.hide_xaxis_label()
