def waitForCASA():
    os.system('sleep 2')
    clearstat()
ones = [1,1,1,1,1,1,1,1,1,1,
        1,1,1,1,1,1,1,1,1,1,
        1,1,1,1,1,1,1,1,1,1,
        1,1,1,1,1,1,1,1,1,1,
        1,1,1,1,1,1,1,1,1]
#execfile('smaImportFix.py')
#waitForCASA()
#os.system('rm -r -f Upper*')
#os.system('rm -r -f Lower*')
vis ='MyDataUpper'
ms.open(vis)
ms.split(outputms=vis+'.tile',tileshape=[1,256,54])
ms.done()
vis ='MyDataUpper.tile'
#flagdata(vis=vis,timerange='18:00:00~23:59:59',antenna='SMA8',mode='manual',spw='')
#waitForCASA()
#flagdata(vis=vis,timerange='05:00:00~06:00:00',antenna='',mode='manual',spw='',field='3c454.3')
setjy(vis=vis,field='3c279',standard='manual',fluxdensity=[13.2,0,0,0],scalebychan=False,usescratch=True)
waitForCASA()
setjy(vis=vis,field='1911-201',standard='manual',fluxdensity=[2.3,0,0,0],scalebychan=False,usescratch=True)
waitForCASA()
setjy(vis=vis,field='1924-292',standard='manual',fluxdensity=[4.4,0,0,0],scalebychan=False,usescratch=True)
waitForCASA()
gaincal(vis=vis,caltable=vis+'.bpself.gcal',field='3c279',spw='1~24',refant='SMA8',calmode='p',solint='int',minblperant=3,minsnr=2.0,combine='spw')

waitForCASA()
bandpass(vis=vis,caltable=vis+'.bandpass.bcal',field='3c279',refant='SMA8',spwmap=[ones],spw='1~24',solint=['inf','3.25MHz'],combine='scan',solnorm=False,gaintable=[vis+'.bpself.gcal'],minsnr=2.0)
waitForCASA()
applycal(vis=vis,gaintable=vis+'.bandpass.bcal')
waitForCASA()
gaincal(vis=vis,spw='1~24',solint='inf',combine='spw',refant='SMA8',gaintype='G',calmode='ap',caltable=vis+'.gaincal',field='1911-201,1924-292')
waitForCASA()
smoothcal(vis=vis,tablein=vis+'.gaincal',caltable=vis+'.gaincal.smoothed',smoothtype='mean',smoothtime=7200.0)
waitForCASA()
applycal(vis=vis,gaintable=[vis+'.bandpass.bcal', vis+'.gaincal.smoothed'],spwmap=[[],ones])
waitForCASA()
split(vis=vis,outputvis='gomezUSB',datacolumn= 'corrected',field='gomez')
vis='gomezUSB'