def X(H,S2,VLA,blue,red): #Herschel maps with contours Fig = aplpy.FITSFigure(H) Fig.show_colorscale(vmax=10000,vmin=0, stretch='linear') Fig.add_colorbar() #markers #YSO.w40(Fig) markers = np.loadtxt('/data/damian/run26/SMM/S2-FF-CO/SMM-cores2.tab') ra = markers[:,1] dec = markers[:,2] #Fig.show_markers(ra,dec, edgecolor='k', facecolor='k', marker='*',s=200, alpha=0.25) Fig.show_markers(ra,dec, edgecolor='k', facecolor='k', marker='*',s=200, alpha=1) OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] YSO.w40(Fig) Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1) #contours cset3 = Fig.show_contour(S2, levels=(0.0025,0.0075,0.025), linewidth=2, colors=('black','black','black')) cset4 = Fig.show_contour(VLA, levels=(0.05,0.25), linewidth=2, colors=('white')) #cset1 = Fig.show_contour(red, levels=(5,25,75), linewidth=2, colors=('red')) #cset2 = Fig.show_contour(blue, levels=(20,40,60,80), linewidth=2, colors=('blue')) #figure properties Fig.ticks.show() Fig.ticks.set_color('black') Fig.colorbar.set_axis_label_text('S$_{70\mu m}$ flux (MJy/Sr)') Fig.recenter(277.8615,-2.0700, height=0.30 ,width=0.375) return
def K(before,after,freefree): big = plt.figure(figsize=(18, 12)) #deffine image in APLPY Before = aplpy.FITSFigure(before,figure=big,subplot=[0.1,0.1,0.4,0.8]) After = aplpy.FITSFigure(after,figure=big,subplot=[0.5,0.1,0.4,0.8]) YSO.w40_radio(Before) YSO.w40_radio(After) YSO.w40(Before) YSO.w40(After) #Figure properties Before.show_grayscale(vmax=0.032,vmin=0.0,invert=True,stretch='sqrt') #After.show_grayscale(vmax=0.032,vmin=0.0,invert=True,stretch='sqrt') After.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE') #Deffine contours cset = Before.show_contour(freefree, levels=(0.001,0.005,0.01), linewidth=(6,4,2), colors=('r')) cset2 = After.show_contour(freefree, levels=(0.001,0.005,0.01), linewidth=(6,4,2), colors=('r')) #Before.recenter(277.8925,-2.08200, height=0.1 ,width=0.15) #After.recenter(277.8925,-2.08200, height=0.1 ,width=0.15) Before.recenter(277.8525,-2.10200, height=0.03 ,width=0.035) After.recenter(277.8525,-2.10200, height=0.03 ,width=0.035) Before.tick_labels.set_font(size='medium') After.tick_labels.set_font(size='medium') Before.add_label(0.1, 0.95, r'SCUBA-2 450um',size='large', relative=True) After.add_label(0.1, 0.95, r'VLA 3.6cm',size='large', relative=True) Before.ticks.show() Before.ticks.set_color('black') After.axis_labels.hide_y() After.tick_labels.hide_y() After.ticks.show() After.ticks.set_color('black') #After.add_colorbar() #After.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)') big.canvas.draw() return
def E(fig): #FELLWALKER CLUMPS #deffine image in APLPY Fig = aplpy.FITSFigure(fig) YSO.w40(Fig) #Figure properties Fig.show_colorscale(vmax=82,vmin=0.0, stretch='sqrt') Fig.tick_labels.set_font(size='medium') Fig.ticks.show() Fig.ticks.set_color('black') return
def N(fig,SCUBA2): Fig = aplpy.FITSFigure(fig) YSO.w40(Fig) YSO.w40_radio(Fig) #3.6cm Fig.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE') #21cm #Fig.show_colorscale(vmax=0.025,vmin=0.0015, stretch='linear') cset3 = Fig.show_contour(SCUBA2, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta')) Fig.add_colorbar() Fig.colorbar.set_axis_label_text('3.6 cm (Jy/Beam)') Fig.ticks.show() return