ডিপথট:Flagging tiles

ডিপথট থেকে
import os, elementflag
o=os.system
ef = elementflag.station()

def visibilise(n):
	global o, ef
	
	# create MS
	pref = '3C196.tilex.'+str(n)
	ms = 'ms/'+pref+'.MS'
	o('rm -rf '+ms)
	print '... creating '+ms
	o('cp -r ms/3C196.root.MS '+ms)
	
	# predict
	print '... making tile configuration default'
	ef.ElementFlagDefault('ms/3C196.root.MS', ms)
	print '... predicting in DATA of '+ms
	o('calibrate-stand-alone -v -f '+ms+' pset/sim.parset sky/3C196.sky &> log/'+pref+'.sim')
	print '... turning '+str(n)+' tiles of all the stations off in '+ms
	ef.WriteElementFlag(ms, 'all', '', n)
	print '... predicting with turned off tiles in MODEL_DATA of '+ms
	o('calibrate-stand-alone -v -f '+ms+' pset/sim.MDATA.parset sky/3C196.sky &> log/'+pref+'.sim')
	
	# calibrate
	print '... calibrating and saving in CORRECTED_DATA of '+ms
	o('calibrate-stand-alone -v -f '+ms+' pset/cal.parset sky/3C196.sky &> log/'+pref+'.cal')
	
	# imaging with lwimager
	comm = "lwimager ms="+ms+" image=dimg/"+pref+".DATA.img fits=fits/"+pref+".DATA.fits npix=900 cellsize=40arcsec data=DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
	print "... executing %s"%comm
	o(comm)
	comm = "lwimager ms="+ms+" image=dimg/"+pref+".MODEL_DATA.img fits=fits/"+pref+".MOEL_DATA.fits npix=900 cellsize=40arcsec data=MODEL_DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
	print "... executing %s"%comm
	o(comm)
	comm = "lwimager ms="+ms+" image=dimg/"+pref+".CORRECTED_DATA.img fits=fits/"+pref+".CORRECTED_DATA.fits npix=900 cellsize=40arcsec data=CORRECTED_DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
	print "... executing %s"%comm
	o(comm)

if __name__=='__main__':
	ns = [1,3,8,12,16,20]
	for i in ns:
		print 'predicting, turning '+str(i)+' tiles off, again predicting and finally calibrating'
		visibilise(i)