Opened 17 years ago
Closed 16 years ago
#125 closed defect (fixed)
Further to last message
Reported by: | Owned by: | Malte Marquarding | |
---|---|---|---|
Priority: | normal | Milestone: | ASAP 2.2 |
Component: | General | Version: | 2.0 |
Severity: | normal | Keywords: | |
Cc: | m.burton@… |
Description (last modified by )
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')
Change History (3)
comment:1 by , 17 years ago
Description: | modified (diff) |
---|---|
Owner: | changed from | to
Reporter: | changed from | to
Status: | new → assigned |
comment:2 by , 17 years ago
Cc: | added |
---|---|
Reporter: | changed from | to
comment:3 by , 16 years ago
Resolution: | → fixed |
---|---|
Status: | assigned → closed |
Note:
See TracTickets
for help on using tickets.
This is broken in the release, but will be fixed in the next release.
Please use this instead
I think this (easier) version might also work.