Changeset 1781 for branches/alma


Ignore:
Timestamp:
07/30/10 11:32:50 (14 years ago)
Author:
Malte Marquarding
Message:

Fix for handling multiple restfrequencies. list of list didn't work and was using fixed units

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/alma/python/scantable.py

    r1757 r1781  
    7575            tmpstr = ''
    7676            for i in range( len(antenna) ):
    77                 if type(antenna[i]) == int: 
    78                     tmpstr = tmpstr + ('%s,'%(antenna[i])) 
     77                if type(antenna[i]) == int:
     78                    tmpstr = tmpstr + ('%s,'%(antenna[i]))
    7979                elif type(antenna[i]) == str:
    8080                    tmpstr=tmpstr+antenna[i]+','
     
    8282                    asaplog.push('Bad antenna selection.')
    8383                    print_log('ERROR')
    84                     return 
     84                    return
    8585            antenna = tmpstr.rstrip(',')
    8686        parallactify = parallactify or rcParams['scantable.parallactify']
     
    485485        if stat.lower().endswith('_abc'): rtnabc = True
    486486        getchan = False
    487         if stat.lower().startswith('min') or stat.lower().startswith('max'): 
     487        if stat.lower().startswith('min') or stat.lower().startswith('max'):
    488488            chan = self._math._minmaxchan(self, mask, stat)
    489489            getchan = True
     
    522522            if self.npol(-1) > 1: out +=  ' Pol[%d] ' % (self.getpol(i))
    523523            #outvec.append(callback(i))
    524             #out += ('= %'+form) % (outvec[i]) +'   '+refstr+'\n' 
    525             out += ('= %'+form) % (statvals[i]) +'   '+refstr+'\n' 
     524            #out += ('= %'+form) % (outvec[i]) +'   '+refstr+'\n'
     525            out += ('= %'+form) % (statvals[i]) +'   '+refstr+'\n'
    526526            out +=  sep+"\n"
    527527
     
    927927            else: raise
    928928        self._add_history("flag_row", varlist)
    929        
     929
    930930    def clip(self, uthres=None, dthres=None, clipoutside=True, unflag=False):
    931931        """
     
    12211221            # list values are scalars
    12221222            if isinstance(freqs[-1], int) or isinstance(freqs[-1], float):
    1223                 self._setrestfreqs(freqs, [""],unit)
     1223                self._setrestfreqs(freqs, [""], unit)
    12241224            # list values are tuples, (value, name)
    12251225            elif isinstance(freqs[-1], dict):
     
    12331233                #self._setselection(savesel)
    12341234                self._setrestfreqs(freqs["value"],
    1235                                    freqs["name"], "MHz")
     1235                                   freqs["name"], unit)
    12361236            elif isinstance(freqs[-1], list) or isinstance(freqs[-1], tuple):
    12371237                sel = selector()
     
    12431243                    sel.set_ifs(iflist[i])
    12441244                    self._setselection(sel)
    1245                     self._setrestfreqs(freqs[i]["value"],
    1246                                        freqs[i]["name"], "MHz")
     1245                    self._setrestfreqs(freqs[i], [""], unit)
    12471246                self._setselection(savesel)
    12481247        # freqs are to be taken from a linecatalog
     
    15401539        """
    15411540        Return a scan where all spectra have been binned up.
    1542        
     1541
    15431542        Parameters:
    15441543            width:       The bin width (default=5) in pixels
     
    17341733        if mask is None:
    17351734            mask = [True for i in xrange(self.nchan(-1))]
    1736        
     1735
    17371736        from asap.asapfitter import fitter
    17381737        try:
     
    17461745            if len(rows) > 0:
    17471746                self.blpars = []
    1748            
     1747
    17491748            for r in rows:
    17501749                # take into account flagtra info (CAS-1434)
     
    17791778            print_log()
    17801779            if insitu: self._assign(workscan)
    1781             else: return workscan 
     1780            else: return workscan
    17821781        except RuntimeError:
    17831782            msg = "The fit failed, possibly because it didn't converge."
     
    19131912            f.data = None
    19141913            f.fit()
    1915            
     1914
    19161915            # Show mask list
    19171916            masklist=workscan.get_masklist(fl.get_mask(),row=r)
     
    24162415        print_log()
    24172416        return s
    2418 
    2419        
Note: See TracChangeset for help on using the changeset viewer.