source: branches/alma/plugins/vanvleck_flag.py@ 3089

Last change on this file since 3089 was 1313, checked in by mar637, 18 years ago

scantable plugin to flag negative values in the bandpass (MOPS)

File size: 1.0 KB
Line 
1def _vanvleck_flag(scan, level=5.0):
2 """
3 Find level crossings and flag all data before (and after it).
4 This rectifies the problem caused by the MOPS bandpass being negative.
5 Use levels greater > 0.0 to get a better result.
6 """
7 sel = selector()
8 saved = scan.get_selection()
9 for i in xrange(scan.nrow()):
10 sel.reset()
11 sel.set_scans(scan.getscan(i))
12 sel.set_beams(scan.getbeam(i))
13 sel.set_ifs(scan.getif(i))
14 sel.set_polarisations(scan.getpol(i))
15 sel.set_cycles(scan.getcycle(i))
16 scan.set_selection(sel)
17 spec = scan._getspectrum()
18 msk = mask_not(scan._getmask())
19 n = len(msk)
20 modified = False
21 for j in xrange(n):
22 if j < n and spec[j] <= level:
23 msk[j] = 1
24 modified = True
25 if modified:
26 scan.flag(msk)
27 scan.set_selection(saved)
28
29if not "vanvleck_flag" in dir(scantable):
30 scantable.vanvleck_flag = _vanvleck_flag
31 del _vanvleck_flag
Note: See TracBrowser for help on using the repository browser.