casa6

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

#                                             #

# Regression/Benchmarking Script for L02D     #

#             26 JAN               (U Cam)    #

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

import time

import os

# Data must be pre-filled

# l.ms, l2.ms

​

pathname=os.environ.get('CASAPATH').split()[0]

datapath=pathname+'/data/regression/ATST1/L02D/l*.ms'

​

startTime=time.time()

startProc=time.clock()

​

print '--Copy MS--'

# Fill - really copy over pristine MS

os.system('rm -rf l.* l2.* all.* l1.* *.ms* *.image *.model *.residual')

copystring='cp -r '+datapath+' .'

os.system(copystring)

clearcal(vis='l.ms')

clearcal(vis='l2.ms')

copytime=time.time()

​

print '--Flag data--'

default('flagdata')

#flagdata(vis='l.ms',

#    timerange=['2002/01/27/05:45:47.0', '2002/01/27/07:00:00.0'],

#    fieldid=-1)

flagdata(vis='l.ms',

     timerange='2002/01/27/05:45:47.0~2002/01/27/07:00:00.0')

flagtime=time.time()

​

#

# l1 calibration 

#

print '--Calibration phase/bandpass (3mm)--'

#setjy

default('setjy')

setjy(vis='l.ms',field='0',spw='6',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='14',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='18',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='7',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='15',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='19',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='3',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='11',scalebychan=False,standard='manual',fluxdensity=[10.87,0.,0.,0.])

setjy(vis='l.ms',field='0',spw='26',scalebychan=False,standard='manual',fluxdensity=[6.0,0,0,0])

setjy(vis='l.ms',field='0',spw='30',scalebychan=False,standard='manual',fluxdensity=[6.0,0,0,0])

setjy(vis='l.ms',field='0',spw='27',scalebychan=False,standard='manual',fluxdensity=[6.0,0,0,0])

setjy(vis='l.ms',field='0',spw='31',scalebychan=False,standard='manual',fluxdensity=[6.0,0,0,0])

​

# 3mm USB

#preliminary time-dependent phase solutions to improve coherent average for bandpass solution

default('gaincal')

gaincal(vis='l.ms',caltable='l.3mmUSB.gcal0',

    field='0',spw='7,15,19:9~118',

    gaintype='GSPLINE',calmode='p',splinetime=10000.,

    refant='1',phasewrap=260,

    preavg=0)

​

print '--bandpass (3mm)--'

#derive bandpass calibration for 3mm (C34S J=2-1) using USB continuum phase solns

default('bandpass')

bandpass(vis='l.ms',caltable='l.3mmC34S.bpoly',

     field='0',spw='3',

     bandtype='BPOLY',degamp=2,degphase=2,solnorm=False,

     maskcenter=2,maskedge=0,refant='1',

     gaintable='l.3mmUSB.gcal0')

​

# Band pass for CH3OH,3mm continuum USB/LSB

default('bandpass')

bandpass(vis='l.ms',caltable='l.3mmch3oh.bpoly',

     field='0',spw='11',

     bandtype='BPOLY',degamp=2,degphase=4,solnorm=False,

     maskcenter=2,maskedge=0,refant='1',

     gaintable='l.3mmUSB.gcal0')

default('bandpass')

bandpass(vis='l.ms',caltable='l.3mmcont.bpoly',

     field='0',spw='7,15,19',

     bandtype='BPOLY',degamp=10,degphase=25,solnorm=False,

     maskcenter=2,maskedge=0,refant='1',

     gaintable='l.3mmUSB.gcal0')

default('bandpass')

bandpass(vis='l.ms',caltable='l.3mmcont.bpoly',

     field='0',spw='6,14,18',

     bandtype='BPOLY',degamp=10,degphase=25,solnorm=False,

     maskcenter=2,maskedge=0,refant='1',

     gaintable='l.3mmUSB.gcal0',append=True)

​