from __future__ import absolute_import
from casatasks.private.casa_transition import is_CASA6
from casatools import table, measures, quanta, ms
from casatasks import casalog
from taskinit import tbtool as table
from taskinit import metool as measures
from taskinit import qatool as quanta
from taskinit import mstool as ms
def converttopoephem2geo(tablename='', outtablename='', overwrite=True):
Convert the given topo ephemeris table to the geocentric ref frame.
Converted are the RA, DEC, and RadVel columns only
tablename -- name of the TOPO frame ephemeris table (in canonical CASA format)
outtablename -- name of the output GEO frame ephemeris table
if type(tablename)!=str or tablename=='':
casalog.post('Invalid parameter tablename', 'WARN')
if type(outtablename)!=str or outtablename=='':
casalog.post('Invalid parameter outtablename', 'WARN')
radvel = tbt.getcol('RadVel')
tmpkw = tbt.getcolkeywords('RadVel')
radvelunit = tmpkw['UNIT']
elif 'QuantumUnits' in tmpkw:
radvelunit = tmpkw['QuantumUnits'][0]
casalog.post('Cannot determine units of radial velocity column. Assuming km/s.', 'WARN')