Source
xxxxxxxxxx
casalog.post("NB: Returning dictionary containing fractional Stokes results.")
from __future__ import absolute_import
import os
import numpy as np
from casatasks.private.casa_transition import is_CASA6
if is_CASA6:
from casatools import calibrater
from casatasks import casalog
from .callibrary import *
else:
from callibrary import *
from taskinit import *
calibrater = cbtool
def polcal(vis=None,caltable=None,
field=None,spw=None,intent=None,
selectdata=None,timerange=None,uvrange=None,antenna=None,scan=None,
observation=None, msselect=None,
solint=None,combine=None,preavg=None,refant=None,minblperant=None,minsnr=None,
poltype=None,smodel=None,append=None,
docallib=None,callib=None,
gaintable=None,gainfield=None,interp=None,spwmap=None):
#Python script
casalog.origin('polcal')
try:
mycb=calibrater()
if ((type(vis)==str) & (os.path.exists(vis))):
mycb.open(filename=vis,compress=False,addcorr=False,addmodel=False)
else:
raise ValueError('Visibility data set not found - please verify the name')
# Do data selection according to selectdata
casalog.post("NB: polcal automatically excludes auto-correlations.")
if (selectdata):
# insist no ACs
if len(msselect)>0:
msselect='('+msselect+') && ANTENNA1!=ANTENNA2'
else:
msselect='ANTENNA1!=ANTENNA2'
# pass all data selection parameters in as specified
mycb.selectvis(time=timerange,spw=spw,scan=scan,field=field,
intent=intent, observation=str(observation),
baseline=antenna,uvrange=uvrange,chanmode='none',
msselect=msselect);
else:
# selectdata=F, so time,scan,baseline,uvrange,msselect=''
# using spw and field specifications only
# also insist no ACs
mycb.selectvis(time='',spw=spw,scan='',field=field,intent=intent,
baseline='',uvrange='',chanmode='none',
observation='', msselect='ANTENNA1!=ANTENNA2')
# set the model, if specified
if (len(smodel)>0):
mycb.setptmodel(smodel);
# Arrange applies....