Custom Query (241 matches)

ASAP closes and reports a segmentation fault after I issue a time average command. I've listed here the process I'm using to get there, and the result from ASAP:

a=scantable('2006-10-04_1924-M165.rpf')
sel=selector()

execfile("vanvleckflag.py")
vanvleck_flag(a,20)

r = a.get_scan('*_R')
s = a.get_scan('*[^_R]')
q = quotient(s,r)
q.set_unit()

sel.reset()
mask=q.create_mask([1500,6000])

q.set_unit('GHz')

q.poly_baseline(order=0,mask=mask)
q.flag(mask)
q.flag(mask_not(mask))

q.set_selection()
sel = selector()
sel.set_ifs(3)
q.set_selection(sel)
tav=q.average_time()

/usr/local/bin/asap: line 54: 24596 Segmentation fault $ip -ipythondir "${HOME}/.asap" -p 'asap' $*

Can the following gain curve be added to ASAP for Parkes (K-band). Questions regarding the polynomial form should be directed to John Reynolds

The current curve for K-band (22.235GHz) is:

     Gain = 0.210528
               +0.296759E-01 * El
               -0.293124E-03 * El^2
               +0.264295E-06 * El^3

where El is elevation in degrees. The curve has a peak normalised to 1.00 (at an El of close to 55 degrees), and has approximate values of 0.845 at 30 degrees El, and 0.70 at 90 degrees El (useful as a sanity check).

Need a function to retrieve integration times, similar to scantable.get_time().