print "Split the calibrated data into single source measurementSet -"
 os.system('rm -rf '+BPname+'.ms')
 print "split: "+BPname
 split(vis=datain+'_'+str(bw)+'.ms',outputvis=BPname+'.ms',spw=allspw,
      field=str(BP), correlation='', width=str(width0))
 clearcal(vis=BPname+'.ms')
 listobs(vis=BPname+'.ms')
 os.system('rm -rf '+FLname+'.ms')
 print "split: "+FLname
 split(vis=datain+'_'+str(bw)+'.ms',outputvis=FLname+'.ms',spw=allspw,
      field=str(FL),  correlation='', width=str(width0))
 clearcal(vis=FLname+'.ms')
 listobs(vis=FLname+'.ms')
 os.system('rm -rf '+T2name+'.ms')
 print "split: "+T2name
 split(vis=datain+'_'+str(bw)+'.ms',outputvis=T2name+'.ms',spw=allspw,
      field=str(T2),  correlation='', width=str(width0))
 clearcal(vis=T2name+'.ms')
 listobs(vis=T2name+'.ms')
 os.system('rm -rf '+T1name+'.ms')
 print "split: "+T1name
 split(vis=datain+'_'+str(bw)+'.ms',outputvis=T1name+'.ms',spw=allspw,
      field=str(T1),  correlation='', width=str(width7))
 clearcal(vis=T1name+'.ms')
 listobs(vis=T1name+'.ms')
 os.system('rm -rf '+CGname+'.ms')
 print "split: "+CGname
 split(vis=datain+'_'+str(bw)+'.ms',outputvis=CGname+'.ms',spw=allspw,
      field=str(CG), correlation='', width=str(width0))
 clearcal(vis=CGname+'.ms')
 listobs(vis=CGname+'.ms')
 print "Done with split!"