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)
col1.recenter(**zoomargs)
col1.colorbar._colorbar_axes.set_ylabel('log$_{10}$(N(H$_2$) cm$^{-2}$)')
