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] column2 = 10**parcubefile[0].data[9] * parcubefile[0].data[13] totalcolumn = np.log10(column1 + column2) totalcolumn[totalcolumn == 0] = np.nan