def S(fig,H70,HII1,HII2):#,mask):
#W40 temp map with Herschel Contours
#deffine image in APLPY
Fig = aplpy.FITSFigure(fig)
#Deffine Markers
#YSO.GBS(Fig)
#YSO.w40(Fig)
YSO.w40_radio(Fig)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
#Deffine contours
cset1 = Fig.show_contour(HII1, levels=(0.05,0.075,0.1,0.2), linewidth=4, colors=('blue')) #Units: Jy/Beam, Resolution: 45''
#cset2 = Fig.show_contour(HII2, levels=(0.001,0.005,0.0075,0.01), linewidth=2, colors=('magenta')) #Units: Jy/Beam, Resolution: 45''
#cset3 = Fig.show_contour(COmask, levels=1, linewidth=2, colors=('black')) #Mask of VanDerWeil2012 CO
#cset4 = Fig.show_contour(H70,levels=(300,1200,4800,12000),colors=('k')) #Herschel 70um data.
#cset3 = Fig.show_contour(flux, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta'))
#cset5 = Fig.show_contour(YSO,levels=(15,45,75),colors=('magenta'))
#Figure properties
Fig.show_colorscale(vmax=30.0,vmin=10, stretch='sqrt')
Fig.tick_labels.set_font(size='medium')
Fig.add_colorbar()
Fig.colorbar.set_axis_label_text('$T_{d}\,(K)$')
Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42)
Fig.ticks.show()
Fig.ticks.set_color('black')
return
def T(fig,H70,s21):
#poster image - 850 + H70/21cm contours + CO
#deffine image in APLPY
Fig = aplpy.FITSFigure(fig)
YSO.w40_radio(Fig)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
#cset4 = Fig.show_contour(H70,levels=(300,1200,4800,12000),colors=('blue')) #Herschel 70um data.
cset1 = Fig.show_contour(s21, levels=(0.01,0.03,0.05,0.1), linewidth=4, colors=('red')) #Units: Jy/Beam, Resolution: 45''
#Figure properties
Fig.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE')
Fig.tick_labels.set_font(size='medium')
Fig.add_colorbar()
#Fig.colorbar.set_axis_label_text('$850\,(K)$')
Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42)
Fig.ticks.show()
Fig.ticks.set_color('black')
return
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 R(fig,surfaceA,surfaceB,H,mask):
#produce surface map of YSO distributions
Fig = aplpy.FITSFigure(fig)
levels = (15,45,75,100,150)
cset_H = Fig.show_contour(H, levels=(200,10000),colors='blue')
#cset1 = Fig.show_contour(surfaceA, levels=levels,colors='green')
#cset2 = Fig.show_contour(surfaceB, levels=levels,colors='k')
cset_mask = Fig.show_contour(mask, levels=(0.5,10),colors='k')
#YSO.w40(Fig)
YSO.c2dGBS(Fig)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=250, alpha=1)
ra = 277.299583333
dec = -2.06388888889
Fig.show_markers(ra,dec, layer='W40a', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1)
Fig.show_grayscale(vmax=0.027,vmin=0.0015, stretch='sqrt',invert='TRUE')
#Fig.recenter(277.8995,-2.09200, height=0.6 ,width=0.6)
Fig.recenter(277.4406666667,-1.90194444444, height=0.85 ,width=0.5)
Fig.add_colorbar()
#Fig.colorbar.set_location('top')
Fig.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)')
def I(fig,red,blue):
#deffine image in APLPY
Fig = aplpy.FITSFigure(fig)
#Deffine contours
cset_red = Fig.show_contour(red, levels=(5,10,15,20,25,30,35), linewidth=3, colors=('r'))
cset_red2 = Fig.show_contour(red, levels=(30,35,100), linewidth=3, filled=('TRUE'), colors=('r'))
cset_blue = Fig.show_contour(blue, levels=(5,10,15,20,25,30,35), linewidth=3, colors=('b'))
cset_blue2 = Fig.show_contour(blue, levels=(30,35,100), linewidth=3, filled=('TRUE'), colors=('b'))
cset = Fig.show_contour(fig, levels=(0.0027,0.0081,0.027,0.081), linewidth=3, colors=('black'))
#Figure properties
Fig.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE')
#Fig.add_label(0.07, 0.95, r'$850\mu m$',size='large', relative=True)
#Deffine Markers
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=25, alpha=1)
YSO.w40_radio(Fig)
Fig.tick_labels.set_font(size='medium')
Fig.add_colorbar()
#Fig.colorbar.set_location('top')
Fig.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)')
Fig.recenter(277.86083,-2.0885,height=0.1 ,width=0.1)
Fig.ticks.show()
Fig.ticks.set_color('white')
return
def A(s450,s850,temp,sig450,sig850,surface):
#tri panel map - should tile all maps at once. 1
big = plt.figure(figsize=(12, 17))
#deffine image in APLPY
Fig_850 = aplpy.FITSFigure(s850,figure=big,subplot=[0.1,0.64,0.8,0.26])
Fig_450 = aplpy.FITSFigure(s450,figure=big,subplot=[0.1,0.37,0.8,0.26])
Fig_temp = aplpy.FITSFigure(temp,figure=big,subplot=[0.1,0.1,0.8,0.26])
#Deffine contours
cset1 = Fig_850.show_contour(s850, levels=(sig850,300*sig850), linewidth=3, colors=('black','black'))
cset2 = Fig_450.show_contour(s450, levels=(sig450,300*sig450), linewidth=3, colors=('black','black'))
#cset3 = Fig_temp.show_contour(temp, levels=(12,20,30), linewidth=3, colors=('black'))
#Figure properties
Fig_450.show_grayscale(vmax=10*sig450,vmin=0.0,invert=True, stretch='sqrt')
Fig_450.add_label(0.07, 0.95, r'b) $450\mu m$',size='large', relative=True)
Fig_850.show_grayscale(vmax=10*sig850,vmin=0.0,invert=True, stretch='sqrt')
Fig_850.add_label(0.07, 0.95, r'a) $850\mu m$',size='large', relative=True)
Fig_temp.show_colorscale(vmax=35.0,vmin=10.0, stretch='sqrt',cmap='gist_heat')
Fig_temp.add_label(0.05, 0.95, r'c) $T_{d}$',size='large', relative=True)
#PLOT ysos
YSO.c2dGBS(Fig_850)
#PLOT contours
cset850 = Fig_temp.show_contour(s850, levels=(sig850,100), linewidth=3, colors=('black'))
csetYSO = Fig_temp.show_contour(surface, levels=(20,60,100), linewidth=3, colors=('red'))
Fig_450.tick_labels.set_font(size='medium')
Fig_450.add_colorbar()
#Fig_450.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_850.tick_labels.set_font(size='medium')
Fig_850.add_colorbar()
#Fig_850.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_temp.tick_labels.set_font(size='medium')
Fig_temp.add_colorbar()
#Fig_temp.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_450.hide_xaxis_label()
Fig_850.hide_xaxis_label()
Fig_450.hide_xtick_labels()
Fig_850.hide_xtick_labels()
Fig_850.hide_yaxis_label()
Fig_temp.hide_yaxis_label()
Fig_450.ticks.show()
Fig_450.ticks.set_color('black')
Fig_850.ticks.show()
Fig_850.ticks.set_color('black')
Fig_temp.ticks.show()
Fig_temp.ticks.set_color('black')
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
def U(small,large):
#deffine image in APLPY
Fig = aplpy.FITSFigure(large)
cset1 = Fig.show_contour(small, levels=(0.001,0.005,0.0075,0.01), linewidth=4, colors=('red')) #Units: Jy/Beam, Resolution: 45''
cset2 = Fig.show_contour(large, levels=(0.001,0.005,0.0075,0.01), linewidth=4, colors=('lime')) #Units: Jy/Beam, Resolution: 45''
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=30, alpha=1)
YSO.w40_radio(Fig)
Fig.ticks.set_color('black')
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 Q(fig,SCUBA2):
#VLA 3.6cm map
Fig = aplpy.FITSFigure(fig)
YSO.w40_radio(Fig)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=200, alpha=1)
cset850 = Fig.show_contour(SCUBA2, levels=(0.0025,0.0075,0.025),colors='k')
Fig.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE')
Fig.add_colorbar()
Fig.recenter(277.8525,-2.09200, height=0.06 ,width=0.07)
Fig.colorbar.set_axis_label_text('$S_{3.6\,cm}$ (Jy/beam)')
Fig.add_beam(9.97*u.arcsecond,9.97*u.arcsecond,24.42*u.degree,color='black')
Fig.ticks.show()
Fig.ticks.set_color('black')
return
def B(s850,temp,sig850,pos):
Fig_temp = aplpy.FITSFigure(temp)
Fig_temp.show_colorscale(vmax=35.0,vmin=10.0, stretch='linear',cmap='gist_heat')
#PLOT ysos
YSO.c2dGBS(Fig_temp)
data = np.loadtxt('OBstars_list.txt')
RA_OB = data[:,0]
DEC_OB = data[:,1]
Fig_temp.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
#PLOT contours
cset850 = Fig_temp.show_contour(s850, levels=(sig850,1000), linewidth=3, colors=('b'))
Fig_temp.ticks.show()
Fig_temp.ticks.set_color('black')
Fig_temp.tick_labels.set_font(size='medium')
Fig_temp.add_colorbar()
Fig_temp.colorbar.set_axis_label_text('T$_{d}$ (K)')
def B(fig,contourA,contourB,contourC):
#deffine image in APLPY
Fig = aplpy.FITSFigure(fig)
#Deffine contours
cset1 = Fig.show_contour(contourA, levels=(0.00655,0.03328), linewidth=3, colors=('red'))
cset2 = Fig.show_contour(contourB, levels=(92.4,462), linewidth=3, colors=('blue'))
cset3 = Fig.show_contour(contourC, levels=(92.4,462), linewidth=3, colors=('blue'))
#Deffine Markers
#YSO.oth(Fig)
YSO.GBS(Fig)
YSO.mwc297(Fig)
#YSO.mwc297_oth(Fig)
#Figure properties
Fig.show_grayscale(vmax=35.0,vmin=5.0,invert=True, stretch='sqrt')
#Fig.add_label(0.1, 0.95, r'$450\mu m$', relative=True)
Fig.tick_labels.set_font(size='medium')
Fig.add_colorbar()
Fig.ticks.show()
Fig.ticks.set_color('black')
return
def L(fig,HII1):#,mask):
#W40 with CO observstion areas
#deffine image in APLPY
Fig = aplpy.FITSFigure(fig)
#Deffine Markers
#YSO.GBS(Fig)
#YSO.oth(Fig)
YSO.w40_radio(Fig)
#YSO.w40(Fig)
#Deffine contour
cset1 = Fig.show_contour(HII1, levels=(0.01,0.05,0.1), linewidth=4, colors=('blue')) #Units: Jy/Beam, Resolution: 45''
#cset2 = Fig.show_contour(HII2, levels=(0.0015,0.015,0.0175,0.02), linewidth=2, colors=('cyan')) #Units: Jy/Beam, Resolution: 45''
cset3 = Fig.show_contour(fig, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta'))
#Figure properties
Fig.show_grayscale(vmax=0.054,vmin=0.0,invert=True)
Fig.add_label(0.07, 0.95, r'$850\mu m$',size='large', relative=True)
Fig.tick_labels.set_font(size='medium')
Fig.add_colorbar()
Fig.ticks.show()
Fig.ticks.set_color('black')
Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42)
return
def P(s450,s850,s450ff,s850ff):
big = plt.figure(figsize=(16, 12))
#deffine image in APLPY
S450 = aplpy.FITSFigure(s450,figure=big,subplot=[0.1,0.1,0.4,0.8])
S850 = aplpy.FITSFigure(s850,figure=big,subplot=[0.5,0.1,0.4,0.8])
#Figure properties
S450.show_grayscale(vmax=0.185,vmin=0., stretch='linear',invert='TRUE')
S850.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE')
#Deffine contours
cset = S450.show_contour(s450ff, levels=(0.0104,0.0173,0.0519,100), colors=('r','k'),filled='TRUE')
cset2 = S850.show_contour(s850ff, levels=(0.0015,0.0025,0.0075,100), colors=('r','k'),filled='TRUE')
cset450 = S450.show_contour(s450, levels=(0.0104,0.0173,0.0519,0.1110), colors=('k'))
cset850 = S850.show_contour(s850, levels=(0.0015,0.0025,0.0075,0.0161), colors=('k'))
#recenter
S450.recenter(277.8525,-2.10200, height=0.025 ,width=0.05)
S850.recenter(277.8525,-2.10200, height=0.025 ,width=0.05)
YSO.w40_radio(S450)
#YSO.w40_radio(S850)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
S850.show_markers(ra,dec, layer='W40a', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1)
S450.show_markers(ra,dec, layer='W40b', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1)
ra = 277.849937
dec = -2.091536
S450.show_markers(ra,dec, edgecolor='k', facecolor='blue', marker='*',s=100, alpha=1)
S850.show_markers(ra,dec, edgecolor='k', facecolor='blue', marker='*',s=100, alpha=1)
#850 set up
S850.axis_labels.hide_y()
S850.tick_labels.hide_y()
S850.ticks.show()
S850.ticks.set_color('black')
S850.add_colorbar()
S850.colorbar.set_font(size='small')
S850.colorbar.set_location('top')
S850.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)')
S850.add_beam(14.5*u.arcsecond,14.5*u.arcsecond,0,color='black')
S850.tick_labels.set_font(size='small')
#450 set up
S450.add_colorbar()
S450.colorbar.set_location('top')
S450.colorbar.set_axis_label_text('450$\mu m$ (Jy/pixel)')
S450.colorbar.set_font(size='small')
S450.add_beam(9.8*u.arcsecond,9.8*u.arcsecond,0,color='black')
S450.ticks.show()
S450.ticks.set_color('black')
S450.tick_labels.set_font(size='small')
big.canvas.draw()
return
RA_OB = data1[:,0]
DEC_OB = data1[:,1]
data2 = np.loadtxt('OBstars.txt',dtype='string')
RA_OBx = map(float, data2[:,2])
DEC_OBx = map(float, data2[:,3])
Fig1 = aplpy.FITSFigure(temp)
Fig1.show_colorscale(vmax=40.0,vmin=10.0, stretch='linear',cmap='CMRmap')
#PLOT contours
#SCUBA-2
Fig1.show_contour(s850, levels=(sig850,1000), linewidth=3, colors=('k'))
#Plot YSOs
YSO.c2dGBS_v2(Fig1)
#PLOT OB stars - primary stars in cyan
Fig1.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='none', marker='*',s=250, alpha=1)
Fig1.show_markers(RA_OBx,DEC_OBx, edgecolor='k', facecolor='cyan', marker='*',s=250, alpha=1)
#Recenter
Fig1.recenter(pos1[0],pos1[1],pos1[2],pos1[3])
#Label
Fig1.add_label(0.65, 0.49, 'A', relative=True)
Fig1.add_label(0.22, 0.42, 'B', relative=True)
Fig1.add_label(0.45, 0.26, 'C', relative=True)
Fig1.add_label(0.4, 0.15, 'E', relative=True)
Fig1.add_label(0.25, 0.08, 'F', relative=True)
"> /dev/null",
)
os.system(cmd)
ndf2fits.ndf2fits(mask)
Fig1.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))
Fig2.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))
Fig3.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))
cmd = "rm %s %s %s" % (thresh, mask, magic)
os.system(cmd)
# Plot YSOs
YSO.c2dGBS_v2(Fig1)
YSO.c2dGBS_v2(Fig2)
YSO.c2dGBS_v2(Fig3)
# PLOT OB stars - primary stars in cyan
Fig1.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)
Fig2.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)
Fig3.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)
Fig1.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)
Fig2.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)
Fig3.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)
# Recenter
Fig1.recenter(pos1[0], pos1[1], pos1[2], pos1[3])
Fig2.recenter(pos2[0], pos2[1], pos2[2], pos2[3])
big = plt.figure(figsize=(25, 6))
Fig = aplpy.FITSFigure(s850,figure=big,subplot=[0.1,0.1,0.375,0.8])
Fig.show_grayscale(vmax=10*sig850,vmin=(1./5.)*sig850,invert=True, stretch='sqrt')
###############
#PLOT-contours#
###############
#OPEN up the region specific clump list
data = fits.open(clumpcat)
DATA = data[1].data
#PLOT ysos
YSO.c2dGBS(Fig)
#OB star position
RA_OB = 328.3708
DEC_OB = 47.2667
Fig.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
Fig.show_contour(s850, levels=(sig850,300*sig850), linewidth=3, colors=('black','black'))
###########
#PLOT-misc#
###########
Fig.ticks.show()
Fig.ticks.set_color('black')
