casa6

from math import *

from tasks import *

from parameter_dictionary import *

from task_help import *

import numpy as np

from casa_stack_manip import stack_frame_find

​

from parameter_check import *

​

## allow globals for taskby default

casaglobals=True

​

​

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

def checkgeodetic():

    """

    Verify the contents of the most important Measures tables

    """

    import os

    rval = True

    casadatadir = os.getenv("CASAPATH").split()[0]+'/data'

    if not os.path.isdir(casadatadir):

        casalog.post('Data repository directory does not exist. Expected at '+casadatadir, 'WARN')

        rval = False

    else:

        if not os.path.isdir(casadatadir+'/geodetic'):

            casalog.post('Data repository sub-directory \"geodetic\" does not exist. Expected at '+casadatadir+'/geodetic', 'WARN')

            rval = False

        else:

            casalog.post('\n', 'INFO')

            casalog.post('Checking measures tables in data repository sub-directory '+casadatadir+'/geodetic', 'INFO')

            mytb = tbtool()

            mytables=['IERSeop2000', 'IERSeop97', 'IERSpredict', 'TAI_UTC']

            for mytable in mytables:

                if not os.path.isdir(casadatadir+'/geodetic/'+mytable):

                    casalog.post('Measures table '+mytable+' does not exist. Expected at '+casadatadir+'/geodetic/'+mytable, 'WARN')

                    rval = False

                else:

                    mytb.open(casadatadir+'/geodetic/'+mytable)

                    vsdate = mytb.getkeyword('VS_DATE')

                    mjd = mytb.getcol('MJD')

                    if len(mjd)>0:

                        mydate = qa.time({'value': mjd[-1], 'unit': 'd'}, form='ymd')[0]

                        casalog.post('  '+mytable+' (version date, last date in table (UTC)): '+vsdate+', '+mydate, 'INFO')

                    else:

                        casalog.post(mytable+' contains no entries.', 'WARN')

                        rval = False

                    mytb.close()

    return rval

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

def go(taskname=None):

    """ Execute taskname: """

    myf=stack_frame_find( )

    if taskname==None: taskname=myf['taskname']

    oldtaskname=taskname

    if(myf.has_key('taskname')):

        oldtaskname=myf['taskname']

    #myf['taskname']=taskname

    if type(taskname)!=str:

        taskname=taskname.__name__

        myf['taskname']=taskname

    try:

        parameter_checktype(['taskname'],[taskname],str)

    except TypeError, e:

        print "go -- TypeError: ",e

        return

    fulltaskname=taskname+'()'

    print 'Executing: ',fulltaskname

    exec(fulltaskname)

    myf['taskname']=oldtaskname

​

​

# from RR

def announce_async_task(taskname):

    """Use the logger to herald the beginning of an asynchronous task."""

    casalog.origin(taskname)

    casalog.post('')

    casalog.post('###############################################')

    casalog.post('###  Begin Task: ' + taskname + ' ' * (27 - len(taskname)) + '###')

    casalog.post("")

    casalog.post("Use: ")

    casalog.post("      tm.retrieve(return_value) # to retrieve the status")

    casalog.post("")

def write_task_obit(taskname):

    """Eulogize the task in the logger."""

    casalog.post('###  End Task: ' + taskname + ' ' * (29 - len(taskname)) + '###')

    casalog.post('###############################################')

    casalog.post('')

​