casa5

# alma-24ant-calibration-regression.py

# based on an ALMA QA2 Calibration Script

​

# CASA parallelization test version

# DP, 17 Oct 2012

​

import sys

import casadef

​

#print " *** init.py: QA2-relevant modules will be imported ***"

#sys.path.append('AIV/science/qa2')

#sys.path.append('AIV/science/analysis_scripts')

​

#from QA2_Packaging_module import *

#import analysisUtils as aU

#es = aU.stuffForScienceDataReduction()

​

#print aU.version()

​

​

mynumsubms = 18

​

# Calibration

​

thesteps = []

step_title = {0: 'Partitioning and listobs',

              1: 'A priori flagging',

              2: 'Generation and time averaging of the WVR cal table',

              3: 'Generation of the Tsys cal table',

              4: 'Generation of the antenna position cal table',

              5: 'Application of the WVR, Tsys and antpos cal tables',

              6: 'Split out science SPWs and time average',

              7: 'Listobs, clear pointing table, and save original flags',

              8: 'Initial flagging',

              9: 'Putting a model for the flux calibrator(s)',

              10: 'Save flags before bandpass cal',

              11: 'Bandpass calibration',

              12: 'Save flags before gain cal',

              13: 'Gain calibration',

              14: 'Save flags before applycal cal',

              15: 'Application of the bandpass and gain cal tables'}

​

try:

  print 'List of steps to be executed ...', mysteps

  thesteps = mysteps

except:

  print 'global variable mysteps not set.'

if (thesteps==[]):

  thesteps = range(0,len(step_title))

  print 'Executing all steps: ', thesteps

​

# The Python variable 'mysteps' will control which steps

# are executed when you start the script using

#   execfile('scriptForCalibration.py')

# e.g. setting

#   mysteps = [2,3,4]# before starting the script will make the script execute

# only steps 2, 3, and 4

# Setting mysteps = [] will make it execute all steps.

​

import re

​

#es = aU.stuffForScienceDataReduction() 

​

​

if cu.compare_version('>=',[5]):

 sys.exit('ERROR: PLEASE USE CASA 4.x')

​

​

​

totaltime = 0

inittime = time.time()

ttime = inittime

steptime = []

​

def timing():

    global totaltime

    global inittime

    global ttime

    global steptime

    global step_title

    global mystep

    global thesteps

    thetime = time.time()

    dtime = thetime - ttime

    steptime.append(dtime)

    totaltime += dtime

    ttime = thetime

    casalog.origin('TIMING')

    casalog.post( 'Step '+str(mystep)+': '+step_title[mystep], 'WARN')

    casalog.post( 'Time now: '+str(ttime), 'WARN')

    casalog.post( 'Time used this step: '+str(dtime), 'WARN')