Changeset 1880 for trunk/python

Timestamp:

08/20/10 12:53:05 (14 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

An implementation of asapmath.splitant() is changed since
the new filler doesn't support antenna parameter at the
moment.

The default behavior of MS filler also changed.

File:

• trunk/python/asapmath.py (modified) (2 diffs)

trunk/python/asapmath.py

@@ -945 +945 @@
     """
     # Import the table toolkit from CASA
-
     import casac
     tbtool = casac.homefinder.find_home_by_name('tableHome')
     tb = tbtool.create()
-    tb2 = tbtool.create()
     # Check the input filename
     if isinstance(filename, str):
@@ -973 +971 @@
         prefix=filename.rstrip('/')
     # Now do the actual splitting.
+    # 2010/08/20 TN
+    # The antenna parameter for scantable constructor is ineffective at the moemnt.
+    # Thus, we should go back to original implementation of splitant.
+##     outfiles=[]
+##     tb.open(tablename=filename+'/ANTENNA',nomodify=True)
+##     nant=tb.nrows()
+##     antnames=tb.getcol('NAME',0,nant,1)
+##     antpos=tb.getcol('POSITION',0,nant,1).transpose()
+##     tb.close()
+##     tb.open(tablename=filename,nomodify=True)
+##     ant1=tb.getcol('ANTENNA1',0,-1,1)
+##     tb.close()
+##     for antid in set(ant1):
+##         scan=scantable(filename,average=False,getpt=True,antenna=int(antid))
+##         outname=prefix+antnames[antid]+'.asap'
+##         scan.save(outname,format='ASAP',overwrite=overwrite)
+##         del scan
+##         outfiles.append(outname)
     outfiles=[]
     tb.open(tablename=filename+'/ANTENNA',nomodify=True)
     nant=tb.nrows()
     antnames=tb.getcol('NAME',0,nant,1)
-    antpos=tb.getcol('POSITION',0,nant,1).transpose()
     tb.close()
     tb.open(tablename=filename,nomodify=True)
     ant1=tb.getcol('ANTENNA1',0,-1,1)
-    tb.close()
+    tmpname='asapmath.splitant.tmp'
     for antid in set(ant1):
-        scan=scantable(filename,average=False,getpt=True,antenna=int(antid))
+        tbsel=tb.query('ANTENNA1 == %s && ANTENNA2 == %s'%(antid,antid))
+        tbsel.copy(tmpname,deep=True)
+        #tbsel=tb.query('ANTENNA1 == %s && ANTENNA2 == %s'%(antid,antid),tmpname)
+        tbsel.close()
+        del tbsel
+        scan=scantable(tmpname,average=False,getpt=True,antenna=int(antid))
         outname=prefix+antnames[antid]+'.asap'
         scan.save(outname,format='ASAP',overwrite=overwrite)
         del scan
         outfiles.append(outname)
-    del tb, tb2
+        os.system('rm -rf '+tmpname)
+    tb.close()
+    del tb
     return outfiles