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()