Changeset 2754 for trunk/python/asapmath.py

Timestamp:

01/31/13 13:14:27 (11 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CSV-2532 (may be related to CSV-1908 and CSV-2161)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed: In tool level, added parameter 'freq_tolsr' to

scantable constructor and function sd.splitant.

Test Programs: test_sdsave, test_importasdm_sd

Put in Release Notes: Yes

Module(s): Module Names change impacts.

Description: Describe your changes here...

In importing MS to Scantable, frequency frame information is
imported as is by default, i.e., base frame in Scantable is
TOPO for ALMA data, which is forcibly converted to LSRK with
wrong time and direction reference.

Some functions have a boolean parameter 'freq_tolsr' that controls
the above behavior. If freq_tolsr is False (default), frequency
is imported as is, while frequency is converted to LSRK (wrongly)
when it is True.

File:

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

trunk/python/asapmath.py

@@ -949 +949 @@
 
 @asaplog_post_dec
-def splitant(filename, outprefix='',overwrite=False):
+def splitant(filename, outprefix='',overwrite=False, getpt=True, freq_tolsr=False):
     """
     Split Measurement set by antenna name, save data as a scantables,
@@ -963 +963 @@
                     The default False is to return with warning
                     without writing the output. USE WITH CARE.
+       getpt        Whether to import direction from MS/POINTING
+                    table or not. Default is True (import direction).
+       freq_tolsr   Whether to convert frequency frame information
+                    to LSRK or not. Default is False (import as is).
 
     """
@@ -1008 +1012 @@
         tb.open(tablename=filename,nomodify=True)
         tbsel=tb.query('ANTENNA1 == %s && ANTENNA2 == %s'%(antid,antid),tmpname)
-        scan=scantable(tmpname,average=False,getpt=True,antenna=int(antid))
+        scan=scantable(tmpname,average=False,antenna=int(antid),getpt=getpt,freq_tolsr=freq_tolsr,)
         outname=prefix+antnames[antid]+'.asap'
         scan.save(outname,format='ASAP',overwrite=overwrite)