casa5

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

# $Id:$

# Test Name:                                                                #

# alma-titan-ephemeris-regression.py                                        #

#                                                                           # 

# Rationale for Inclusion:                                                  #

#    Need test of analysis of MS with ephemeris attached                    #

#                                                                           # 

# Input data:                                                               #

#     two MSs                                                               #

#                                                                           #

# Author: DP, Sept. 2013                                                    #

#                                                                           #

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

​

​

step_title = { 0 : 'fixplanets with ephemeris',

           1 : 'concat',

               2 : 'selfcal',

           3 : 'cvel to outframe SOURCE',

               4 : 'imaging',

               5 : 'verification of the regression results'

               }

​

# global defs

basename=['X19d','X346']

makeplots=False

​

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

​

# Some infrastructure to make repeating individual parts

#   of this workflow more convenient.

​

thesteps = []

​

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 'mysteps empty. Executing all steps: ', thesteps

​

# The Python variable 'mysteps' will control which steps

# are executed when you start the script using

#   execfile('alma-titan-ephemeris-regression.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.

​

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')

    casalog.post( 'Total time used so far: ' + str(totaltime), 'WARN')

    casalog.post( 'Step  Time used (s)     Fraction of total time (percent) [description]', 'WARN')

    for i in range(0, len(steptime)):

        casalog.post( '  '+str(thesteps[i])+'   '+str(steptime[i])+'  '+str(steptime[i]/totaltime*100.)

                      +' ['+step_title[thesteps[i]]+']', 'WARN')

​

# default ASDM dataset name

myasdm_dataset_name = "X19d.ms.split.cal.titan.spw0"

myasdm_dataset2_name = "X346.ms.split.cal.titan.spw0"

myephemeris = 'Titan_55197-59214dUTC_J2000.tab'

​

# get the dataset name from the wrapper if possible

mydict = locals()

if mydict.has_key("asdm_dataset_name"):

    if(myasdm_dataset_name != mydict["asdm_dataset_name"]):