casa6

###############################################

#                                             #

# Regression/Benchmarking Script for G192     #

#                                             #

###############################################

​

import time

import os

​

os.system('rm -rf g192*.ms g192_a* g192.*.im')

​

startTime = time.time()

startProc = time.clock()

​

print '--Import--'

datapath=os.environ.get('CASAPATH').split()[0] +'/data/regression/ATST1/G192/'

​

default('importvla')

importvla(archivefiles=[datapath + 'AS758_C030425.xp1',datapath+'AS758_C030425.xp2',datapath+'AS758_C030425.xp3',datapath+'AS758_C030426.xp4',datapath+'AS758_C030426.xp5'],

      vis='g192_a.ms',bandname='K',frequencytol=10000000.0)

importtime = time.time() 

print '--Observation summary--'

default('listobs')

listobs(vis='g192_a.ms')

#listtime = time.time()

print '--Flag auto-correlations--'

default('flagdata')

flagdata(vis='g192_a.ms',mode='manual',autocorr=True)

flagtime = time.time()

print '--Setjy--'

default('setjy')

setjy(vis='g192_a.ms',field='4',scalebychan=False,standard='Perley-Taylor 99') #set flux density for 1331+305 (3C286)

setjytime = time.time()

​

print '--Gencal(opac)--'

# make opacity caltable

default('gencal')

gencal(vis='g192_a.ms',caltable='g192_a.opac',caltype='opac',parameter=[0.062])

gencaltime = time.time()

​

print '--Gaincal--'

#Select data for gain calibrators and drop outer channesl that may bias the solution

default('gaincal')

gaincal(vis='g192_a.ms',caltable='g192_a.gcal',

    field='0,2,3,4',spw='0:3~117', gaintype='G',

    solint='inf',combine='',refant='VA05',

    gaintable=['g192_a.opac'])

gaintime = time.time()

print '--Bandpass--'

#Select bandpass calibrator. Arrange to solve for a single solution over the entire observation 

default('bandpass')

bandpass(vis='g192_a.ms',caltable='g192_a.bcal',

     field='3',  

     gaintable=['g192_a.opac','g192_a.gcal'],gainfield=['','3'],interp=['','nearest'],

     solint='inf',combine='scan',

     refant='VA05')

bptime = time.time()

print '--Fluxscale--'

#Transfer the flux density scale from the flux calibrator to the gain calibrators

#Solutions are written to a table (g192_a.fluxcal) on disk.

default('fluxscale')

fluxscale(vis='g192_a.ms',caltable='g192_a.gcal',fluxtable='g192_a.fluxcal',

      reference=['1331+305'],transfer=['0530+135','05309+13319'])

fstime = time.time()

​

print '--Correct--'

#Apply calibration solutions to the data

default('applycal')

applycal(vis='g192_a.ms',

     field='0,1,2', 

     gaintable=['g192_a.opac','g192_a.fluxcal','g192_a.bcal'],gainfield=['','0,2']);

correcttime = time.time()

​

print '--Split (Cal/src data)--'

#split out the data of interest into a new visibility data set

default('split')

split(vis='g192_a.ms',outputvis='g192_cal.split.ms',

 #     field=4,spw=0,nchan=100,start=9,step=1,datacolumn='CORRECTED_DATA')

    field='4',spw='0:9~108',datacolumn='corrected')

splitcaltime = time.time()

default('split')

split(vis='g192_a.ms',outputvis='g192_src.split.ms',

#      field=1,spw=0,nchan=100,start=9,step=1,datacolumn='CORRECTED_DATA')

    field='1',spw='0:9~108',datacolumn='corrected')

splitsrctime = time.time()

​

print '--Flag bad time range--'

#flag data in the specified time range for the source and spw

flagdata(vis="g192_src.split.ms", field="0", spw="0", 

      timerange="2003/04/26/02:45:00.0~2003/04/26/02:46:30.0")

flagsrctime=time.time()