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
