Further to last message

[m.burton]

Malte

Further to my last message, I noticed that there was a directory problem (asap didn't seem to be able to cope with a space in a directory name), so I removed the space (note this wasn't a problem in the old version of asap on my laptop). When I fix this I then get the following error message. Note that reading in the files is taking an absolute age. Again much longer then before??

Michael

ASAP>execfile('071001_asap.txt') Auto averaging integrations Importing 2007-10-01_082248_T040.rpf... Importing 2007-10-01_091752_T040.rpf... Importing 2007-10-01_094543_T040.rpf... Importing 2007-10-01_114242_T040.rpf... Importing 2007-10-01_114813_T040.rpf... Importing 2007-10-01_123418_T040.rpf... Importing 2007-10-01_125836_T040.rpf... Couldn't find any match. Couldn't find any match.

---

Traceback (most recent call last):

File "<ipython console>", line 1, in ? File "071001_asap.txt", line 33, in ?

q = quotient(src,ref,preserve=False)

File "/Library/Frameworks/Python.framework/Versions/2.4/lib/python2.4/site-packages/asap/asapmath.py", line 110, in quotient

s = scantable(stm._quotient(source, reference, preserve))

ArgumentError: Python argument types in

stmath._quotient(stmath, NoneType, NoneType, bool)

did not match C++ signature:

_quotient(asap::STMathWrapper {lvalue}, asap::ScantableWrapper, asap::ScantableWrapper, bool)

*

This is my script:

# Reduction of Tid 34m data for 2007-10-01, Day of Year 274 # Manually adjust relevant source to be reduced (via s=sXY) # And also change source name (3 times: once in plot header, twice in output files) # files = ['2007-10-01_082248_T040.rpf', '2007-10-01_091752_T040.rpf', '2007-10-01_094543_T040.rpf', '2007-10-01_114242_T040.rpf', ' 2007-10-01_114813_T040.rpf', '2007-10-01_123418_T040.rpf', '2007-10-01_125836_T040.rpf']

d = scantable(files)

s1=d.get_scan(range(0,16)) s2=d.get_scan(range(16,32)) s3=d.get_scan(range(32,48)) s4=d.get_scan(range(48,64)) s5=d.get_scan(range(66,84)) s6=d.get_scan(range(84,100)) s7=d.get_scan(range(100,116)) s8=d.get_scan(range(116,132)) s9=d.get_scan(range(132,141))

# change this line manually s=s9

# sel = selector() # select the second IF for this day ONLY sel.set_ifs([1]) s.set_selection(sel) # Scaling factor to get Tsys approx in K scaled = s.scale(3.33,tsys=True,insitu=False) # Separate the on and off-source scans ref = scaled.get_scan('*[ew]') src = scaled.get_scan('*[^{ew]') # take the quotient q = quotient(src,ref,preserve=False) q.set_freqframe('LSRK') q.set_unit('km/s') # calculate antenna Az, El q.recalc_azel() # gain-elevation correction for 34m at Ka-band q.gain_el([0.1757578928574,0.000981543992646,-1.040739732816e-05]) # approx opacity on day 191 #q.opacity(0.05) q.freq_align() # HC5N rest freq q.set_restfreqs(freqs=31951.777e6) q.set_fluxunit('Jy') a = q.average_time() plotter.set_mode('i','t') # change plot header manually plus two output file names #plotter.set_title('HC5N G0.84+0.1 DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G0.50+0.1 DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G5.89-0.3 DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G14.99-0. DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G16.59-0. DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G22.36+0. DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G2.54+0.2 DSS-34 2007-10-01',fontsize=18) #plotter.set_title('HC5N G28.31-0. DSS-34 2007-10-01',fontsize=18) plotter.set_title('HC5N G28.15-0. DSS-34 2007-10-01',fontsize=18)}

plotter.set_range(xstart=-100,xend=150,ystart=-2,yend=4)

# remove residual baseline #msk = a.create_mask([-100,-10],[30,230],[310,420]) #a.poly_baseline(msk,order=2) a.auto_poly_baseline(edge=[20],order=5, threshold=8, insitu=True, plot=False) # smooth smoo=a.smooth(kernel='hanning',insitu=False) plotter.plot(smoo)

a.save('9','ASCII','overwrite') plotter.save('9.png')