__group__ ticket summary component version milestone type severity owner status created _changetime _description _reporter Active Tickets 258 scantable masking bug General 2.0 ASAP 4.0 defect critical Malte Marquarding new 2011-11-11T08:27:17+11:00 2012-02-17T14:05:50+11:00 "The code to handle mask (scantable.create_mask() etc.) has a bug. It causes failure for the data that has more than one IFs whose number of channels are not equal. Problem is that the code always refers 'nChan' field in the header for number of channels. In create_mask(), for example, it will fail (by the exception like 'list index out of range') when IFNO for specified row has different number of channels from 'nChan' field value. " Takeshi Nakazato Active Tickets 194 Create developer doc Documentation 2.0 Unified development task major Malte Marquarding new 2010-08-03T21:37:35+10:00 2010-08-03T21:37:35+10:00 "At a minimum we need: * Scantable schema and description * how to to c++ to python bindings * API docs * overview * note on large data sets * Filler docs (Takeshi?) probably use sphinx for this" Malte Marquarding Active Tickets 279 export to sdfits of data with xpol no longer working c++ 2.0 Unified development defect major Takeshi Nakazato new 2012-09-27T17:24:21+10:00 2012-10-17T11:38:18+11:00 "During the introduction of npol to PKSrec the behaviour seem to have changed. I can no longer export to sdfits. The data is in trunk/test/data/parkes-pol.rpf Could you look into this please. When I tried to change the STwriter to only pass through npol=2 this breaks the MS2writer. " Malte Marquarding Active Tickets 196 Create new Scantable definition wiki 2.0 Unified development enhancement normal Malte Marquarding assigned 2010-08-13T10:20:29+10:00 2011-04-13T18:04:50+10:00 "Make a formal wiki page for the Scantable version 4 design. " Malte Marquarding Active Tickets 206 Tidy up poly_baseline batch mode General 2.0 Unified development defect normal WataruKawasaki new 2010-09-10T10:58:47+10:00 2010-09-14T12:22:27+10:00 "The current interface returns pointers. Return {{{asapfit}}} objects instead. " Malte Marquarding Active Tickets 213 Speed up baseline methods. c++ 2.0 Unified development defect normal WataruKawasaki new 2010-09-28T08:37:27+10:00 2010-09-28T08:43:55+10:00 This is to track the correspondence and development. PLEASE add comments here. Malte Marquarding Active Tickets 224 cmake build for CASA-ASAP General 2.0 Unified development enhancement normal Takeshi Nakazato new 2010-12-02T17:21:37+11:00 2011-05-10T16:12:52+10:00 "Currently, ASAP for CASA is compiled by make system with fixed Makefile. Cross platform build system using cmake must be needed for unified build system on casacore, casa, and asap. Corresponding to CAS-2668." Takeshi Nakazato Active Tickets 235 "manipulation of the on/off header using the ""quotient"" tool in ASAP" General 2.0 Unified development defect normal Takeshi Nakazato new 2011-03-04T13:29:48+11:00 2011-03-04T13:58:23+11:00 "This ticket is raised during single-dish developers meeting held in Tokyo at Jul./Aug. 2010. " Takeshi Nakazato Active Tickets 111 writing to SDFITS fails if data are pol averaged General 2.0 enhancement normal Malte Marquarding assigned 2007-06-27T11:03:31+10:00 2010-08-25T11:30:32+10:00 "using data /DATA/KAPUTAR_3/mopra/MOPSarchive/2007-06-25_2212-M232.rpf with these steps: s=scantable(fileinname) q=s.auto_quotient() tav=q.average_time(weight='tsys') polav=tav.average_pol() polav.set_unit('MHz') plotter.plot(tav) polav.save('NMLTAU2007-06-25_2212','SDFITS') generates: Traceback (most recent call last): File """", line 1, in ? File ""reduction.py"", line 42, in ? polav.save('NMLTAU2007-06-25_2212','SDFITS') File ""/usr/local/lib/python2.3/site-packages/asap/scantable.py"", line 120, in save writer.write(self, name) RuntimeError: poltype = stokes not yet supported in output. The error isn't produced if the average_pol step is omitted." herpin@… Active Tickets 134 scantable constructor doesn't work for multiple asap files python 2.0 enhancement normal Malte Marquarding assigned 2008-05-22T12:12:18+10:00 2010-02-11T12:00:34+11:00 "{{{ scantable(['a.asap', 'b.asap']) }}} doesn't work. Need to handle multiple files via {{{ merge }}} " Malte Marquarding Active Tickets 144 Ability to create empty scantable? General 2.0 enhancement normal Malte Marquarding assigned 2008-09-08T22:26:58+10:00 2010-08-25T11:29:35+10:00 "I've recently written a script to load data stored as a two column text file (velocity and flux density) and use the ASAP Gaussian fitting and it works very nicely. I do this using code like : def gauss_fit_data(x,y,ngauss,outname='none') : # Fit Gaussian profiles to the data f = fitter() f.set_data(x,y) f.set_function(gauss=ngauss) f.fit() .... where x and y are lists containing the data read from the file. I would like to be able to create a scantable with the same data in it using scantable._setspectrum (and scantable._setabcissa, although that doesn't seem to exist?). In particular it would be useful to be able to create/use masks, use the in built statistic etc, but there doesn't seem to be any facility to create an empty scantable? I know that I can fake some of this by loading in an unrelated dataset, but its not very elegant (nor portable)." Simon.Ellingsen@… Active Tickets 186 Move this trac to science redmine General 2.0 task normal Malte Marquarding new 2010-05-24T14:36:30+10:00 2010-05-24T14:36:30+10:00 Malte Marquarding Active Tickets 168 scantable.save doesn't honour frequency conversions for non asap output General 2.0 ASAP 2.4 defect minor Malte Marquarding assigned 2009-06-25T14:41:11+10:00 2010-02-11T11:47:26+11:00 "Frequency conversion is only done on-the-fly for {{{toWorld}}} calls. The underlying coordinate isn't changed from the original though. That means on data export the initial state is saved not taking into account e.g. frame changes. Th current workaround is to use {{{scantable.freq_align}}} before saving. This will regrid the coordinate. We should probably make this the default. However, this won't work for {{{PKSwriter}}} output as this takes information from the {{{STHeader}}}. This needs to be changed and then we can get rid off this class in favour of table keywords. This should become a new ticket and breaks the table format -> major version change." Malte Marquarding Active Tickets 162 Polarisation data extraction General 2.0 ASAP 2.4 defect normal Malte Marquarding reopened 2009-04-30T16:29:45+10:00 2009-05-08T11:23:16+10:00 "I'm trying to extract ascii files with just the polarisation data in them, but the files produced look identical. I was wondering if you might be able to spot anything obvious that I'm missing. Here are the commands I'm using: {{{ plotscans.set_selection() sel = plotscans.get_selection() sel.set_polarisations(['I','Plinear']) plotscans.set_selection(sel) plotscans.save(source+'IPl','ASCII') sel.reset() plotscans.set_selection() sel = plotscans.get_selection() sel.set_polarisations(['RR', 'LL']) plotscans.set_selection(sel) plotscans.save(source+'RRLL','ASCII') }}}" james.green@… Active Tickets 225 Uncalibrated units option General 2.0 ASAP 2.4 enhancement normal Malte Marquarding new 2010-12-07T05:11:07+11:00 2010-12-07T05:11:07+11:00 "It would be nice to have a third option for scantable.set_unit(), which would be something like ""counts"" or None, with a Y-axis label of ""Counts""." kuiper@… Active Tickets 212 On OSX 10.6: IERSpredict not found despite successfully running asap_update_data General 2.0 ASAP 3.0 defect blocker Malte Marquarding assigned 2010-09-23T16:03:42+10:00 2010-09-24T19:36:12+10:00 "Then, when running my reduction scripts, these are the messages, the script fails finally on average_pols (which may or may not be related). The same script runs fine on kaputar. {{{ 010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309) Requested data table IERSpredict cannot be found in the searched directories: 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ ./ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ ./data/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Users/balt/aips++/data//ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Library/Python/2.6/site-packages/asap/data//ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Users/balt/aips++/data//geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Library/Python/2.6/site-packages/asap/data//geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/geodetic/ 2010-09-23 05:50:35 INFO MeasIERS::initMeas(MeasIERS::Files) Cannot read IERS table IERSpredict 2010-09-23 05:50:35 INFO MeasIERS::initMeas(MeasIERS::Files) + Calculations will proceed with lower precision 2010-09-23 05:50:35 INFO MeasTable::dUT1(Double) No requested dUT1 data available from IERS tables. 2010-09-23 05:50:35 INFO MeasTable::dUT1(Double) + Proceeding with probably less precision. 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309) Requested data table TAI_UTC cannot be found in the searched directories: 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ ./ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ ./data/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Users/balt/aips++/data//ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Library/Python/2.6/site-packages/asap/data//ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/ephemerides/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Users/balt/aips++/data//geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /Library/Python/2.6/site-packages/asap/data//geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/geodetic/ 2010-09-23 05:50:35 WARN MeasIERS::fillMeas(MeasIERS::Files, Double) (file measures/Measures/MeasIERS.cc, line 309)+ /local/mar637/pyenv/share/casacore/data/geodetic/ 2010-09-23 05:50:35 SEVERE MeasTable::dUTC(Double) (file measures/Measures/MeasTable.cc, line 6295) Cannot read leap second table TAI_UTC 2010-09-23 05:50:35 WARN CoordinateUtil::makeFrequencyMachine Unable to convert between the input and output SpectralCoordinates 2010-09-23 05:50:35 WARN CoordinateUtil::makeFrequencyMachine+ this probably means one is in the REST frame which requires 2010-09-23 05:50:35 WARN CoordinateUtil::makeFrequencyMachine+ the radial velocity - this is not implemented yet 2010-09-23 05:50:35 SEVERE CoordinateUtil::makeFrequencyMachine A trial conversion failed - something wrong with frequency conversion machine 2010-09-23 05:50:35 SEVERE CoordinateUtil::makeFrequencyMachine A trial conversion failed - something wrong with frequency conversion machine Aligned at reference Epoch 2010/09/21/13:04:07 (UTC) in frame BARY Traceback (most recent call last): File ""./process.py"", line 41, in iav = av.average_pol(weight='tsys') AttributeError: 'NoneType' object has no attribute 'average_pol' }}}" balt.indermuehle@… Active Tickets 181 plotting a saved file with one IF doesn't work General 2.0 ASAP 3.0 defect major Malte Marquarding assigned 2010-03-28T04:12:40+11:00 2011-03-04T14:01:30+11:00 "Code to demo bug: {{{ d = scantable('2010-03-24_1445_316p81-0p06.rpf') m = d.auto_quotient() m.set_freqframe('LSRK') m.set_unit('GHz') m.freq_align() iffreqs = [23.6945e9,24.13939e9] m.set_restfreqs(iffreqs) # Set the (1,1) and (4,4) freFs pdq = m.average_pol(weight='tsys') av = average_time(pdq,weight='tsys',align=False) bsp = av.copy() bsp.save('blspec',overwrite=True) resp=scantable('blspec') plotter.plot(resp) }}} works but {{{ ifno=0 sel=selector() sel.set_ifs(ifno) bsp.set_selection(sel) bsp.save('only1IF',overwrite=True) resp=scantable('only1IF') plotter.plot(resp) }}} dies with the message {{{ /usr/lib64/python2.6/site-packages/asap/asapplotter.py in _plot(self, scan) 765 y = scan._getspectrum(r, polmodes[scan.getpol(r)]) 766 else: --> 767 y = scan._getspectrum(r) 768 m = scan._getmask(r) 769 from matplotlib.numerix import logical_not, logical_and RuntimeError: Table DataManager error: Invalid operation: MSM: no array in row 0 of /local/home/gaf/DATA/Parkes/MMB_NH3/temp14634_856/table.f1 }}}" G.Fuller@… Active Tickets 261 Problems running ASAP General 2.0 ASAP 4.0 defect major Malte Marquarding new 2012-02-06T02:25:31+11:00 2012-02-12T09:48:05+11:00 "Hi, I'm running ubuntu 10.04 Lucid and I'm having a few issues after attempting a manual installation. The installation completed without errors but when I try to run the commands ASAP>s = scantable (""2008-03-12_0932-M999.rpf"") Found ATMOPRA data. Auto averaging integrations Importing 2008-03-12_0932-M999.rpf... ASAP>s.summary() I get the following error: WARN: Requested data table IERSeop97 cannot be found in the searched directories: ./ ./data/ /home/jack/aips++/data//ephemerides/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jack/aips++/data//geodetic/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ Cannot read IERS table IERSeop97 Calculations will proceed with lower precision WARN: Requested data table IERSpredict cannot be found in the searched directories: ./ ./data/ /home/jack/aips++/data//ephemerides/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jack/aips++/data//geodetic/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ Cannot read IERS table IERSpredict Calculations will proceed with lower precision No requested dUT1 data available from IERS tables. Proceeding with probably less precision. WARN: Requested data table TAI_UTC cannot be found in the searched directories: ./ ./data/ /home/jack/aips++/data//ephemerides/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jack/aips++/data//geodetic/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ SEVERE: Cannot read leap second table TAI_UTC WARN: Unable to convert between the input and output SpectralCoordinates this probably means one is in the REST frame which requires the radial velocity - this is not implemented yet WARN: Requested data table TAI_UTC cannot be found in the searched directories: ./ ./data/ /home/jack/aips++/data//ephemerides/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jack/aips++/data//geodetic/ /usr/local/lib/python2.6/dist-packages/asap-4.0.0-py2.6.egg/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ SEVERE: Cannot read leap second table TAI_UTC WARN: Unable to convert between the input and output SpectralCoordinates this probably means one is in the REST frame which requires the radial velocity - this is not implemented yet SEVERE: Couldn't convert frequency frame. Also when I first load ASAP up I get this warning: Loading ASAP... /usr/lib/python2.6/dist-packages/pytz/__init__.py:32: UserWarning: Module matplotlib was already imported from /usr/lib/python2.6/dist-packages/matplotlib/__init__.pyc, but /usr/lib/python2.6/dist-packages is being added to sys.path from pkg_resources import resource_stream Welcome to ASAP v4.0.0 (2011-10-05) - the ATNF Spectral Analysis Package Any help would be greatly appriciated. Kind regards - Jack " py08ja@… Active Tickets 173 NH3 line fitting General 2.0 ASAP 4.0 enhancement normal Malte Marquarding assigned 2009-09-11T08:38:50+10:00 2012-10-15T15:53:03+11:00 "Hi, Following on from ticket:166 it would be nice to have ammonia line fitting in ASAP. I know it is possible to export data to CLASS format, but if I can keep everything within ASAP, so much the better. THe old SPC program did have a NH3fit routine which could be ported to ASAP? Cheers, Stacy." Stacy Mader Active Tickets 233 Test Interactive flagging module, flagplotter python 2.0 ASAP 4.0 enhancement normal Malte Marquarding new 2011-02-25T12:37:06+11:00 2011-05-10T16:23:11+10:00 "Please test the new module, flagplotter, on stand alone ASAP. You can load flagplotter and play with it by 1. building trunk codes, 2. running the attached Python script, 'guiflag.py' to define a function guiflag, and 3. invoking the guiflag function, e.g., guiflag(sdfile='INPUT_FILE_NAME',outfile='OUTPUT_NAME.asap') You should hit return on the terminal to finish flagging and save file. There are a few known issues: a. Buttons in the new toolbar are not nice. It's possible to change the faces to pictures like plotms in CASA. b. Page change and flagging operations are somehow slow. c. Buttons should be disabled during slow operation. Description of the new toolbar buttons: * 'region': add a region to the selections by pressing left mouse button twice which specify start and end channel * 'panel': add panels to the selections by pressing left mouse on the panels. * 'clear': clear all the selections * 'flag': flag the selections * 'unflag': unflag the selections * 'statistics': calculate statistics of the selections (results are output to logger) * 'notation': add notes on the plot * Page '+' and '-': go back and forth pages * 'Quit': finish flagplotter (without saving flags)" Kana Sugimoto Active Tickets 281 Several rms values for one spectrum General 1.0 ASAP 4.0 task normal Malte Marquarding assigned 2012-10-16T02:03:52+11:00 2012-10-17T11:37:26+11:00 "Hi, I'm reducing Mopra data taken in position-switching mode. After removing the bad scans and performing the quotient calculation, I merged the relevant scantables, time-averaged the whole bunch and finally averaged both polarisations. When I plot the resulting 'scantable', I get one single spectrum (as expected). However, when I use the task 'stats' to get the rms depth, it gives me several different values, corresponding to different times. Is this normal? And if so, how should I interpret these values? Thanks, Flor" flor.allaert@… Active Tickets 284 exception raising for bad scans python 2.0 ASAP 4.0 enhancement normal Malte Marquarding new 2012-10-23T03:19:04+11:00 2012-10-23T03:19:04+11:00 "I got this error which I think comes from the reader code SEVERE: ArrayBase::validateConformance shape [1] differs from [2] At least, I can't find it in the asap package and I can't trap it with try: ... except: I suspect there's a bad scan than I need to recognize and skip. I know which scan it is from where the program stopped. However, it's one scan among dozens of scans that day, and a hundred or more days that I'm processing to make an index. I can't put in ad hoc code to skip bad scans. " kuiper@… Active Tickets 286 use logging module instead of print General 2.0 ASAP 4.0 enhancement normal Malte Marquarding new 2012-10-30T05:25:36+11:00 2012-10-30T05:37:25+11:00 Proposal to use 'logging' so unwanted info messages can be turned off. kuiper@… Active Tickets 273 Problem with ASAP on Fedora 17 or iPython > 0.10 python 2.0 Unified development defect blocker Malte Marquarding new 2012-07-11T16:34:09+10:00 2012-07-11T16:34:09+10:00 " I've been trying to get ASAP to work on Fedora 17. I have built it without errors but cannot run it because of changes to ipython (from version 0.11 onwards). I have asap runing ok with Fedora 16 with ipython version 0.10. However when I update ipython in Fedora 16 asap stops working. The output from the Fedora 17 is as follows: running asap results in the following output: [admin@maths14 ~]$ asap Loading ASAP... RUNNING->/usr/bin/ipython-ipythondir /home/admin/.asap usage: ipython [-h] [--profile TERMINALIPYTHONAPP.PROFILE] [-c TERMINALIPYTHONAPP.CODE_TO_RUN] [--pylab [TERMINALIPYTHONAPP.PYLAB]] [--autocall TERMINALINTERACTIVESHELL.AUTOCALL] [--ipython-dir TERMINALIPYTHONAPP.IPYTHON_DIR] [--gui TERMINALIPYTHONAPP.GUI] [--logappend TERMINALINTERACTIVESHELL.LOGAPPEND] [--colors TERMINALINTERACTIVESHELL.COLORS] [--log-level TERMINALIPYTHONAPP.LOG_LEVEL] [--ext TERMINALIPYTHONAPP.EXTRA_EXTENSION] [--cache-size TERMINALINTERACTIVESHELL.CACHE_SIZE] [--logfile TERMINALINTERACTIVESHELL.LOGFILE] [--no-autoindent] [--pprint] [--classic] [--term-title] [--no-confirm-exit] [--autoindent] [--no-term-title] [--deep-reload] [--color-info] [--init] [--pydb] [--no-color-info] [--autoedit-syntax] [--confirm-exit] [--no-autoedit-syntax] [--quick] [--banner] [--automagic] [--no-automagic] [--nosep] [-i] [--quiet] [--no-deep-reload] [--no-pdb] [--debug] [--pdb] [--no-pprint] [--no-banner] ipython: error: argument -i/--i: ignored explicit argument u'pythondir' [admin@maths14 ~]$ Changing the call to ipython and the asap script as follows: #$ip -ipythondir ""${HOME}/.asap"" $* $ip --ipython-dir=""${HOME}/.asap"" $* results in the following message: [admin@maths14 ~]$ asap Loading ASAP... RUNNING->/usr/bin/ipython-ipythondir /home/admin/.asap /usr/lib/python2.7/site-packages/IPython/utils/path.py:417: UserWarning: Found old IPython config file u'/home/admin/.asap/ipy_user_conf.py' (modified by user) warnings.warn(""Found old IPython config file %r (modified by user)""%f) /usr/lib/python2.7/site-packages/IPython/utils/path.py:417: UserWarning: Found old IPython config file u'/home/admin/.asap/ipythonrc' (modified by user) warnings.warn(""Found old IPython config file %r (modified by user)""%f) /usr/lib/python2.7/site-packages/IPython/utils/path.py:429: UserWarning: The IPython configuration system has changed as of 0.11, and these files will be ignored. See http://ipython.github.com/ipython-doc/dev/config for details of the new config system. To start configuring IPython, do `ipython profile create`, and edit `ipython_config.py` in /profile_default. If you need to leave the old config files in place for an older version of IPython and want to suppress this warning message, set `c.InteractiveShellApp.ignore_old_config=True` in the new config. `c.InteractiveShellApp.ignore_old_config=True` in the new config."""""") Python 2.7.3 (default, Apr 30 2012, 21:18:11) Type ""copyright"", ""credits"" or ""license"" for more information. IPython 0.12 -- An enhanced Interactive Python. ? -> Introduction and overview of IPython's features. %quickref -> Quick reference. help -> Python's own help system. object? -> Details about 'object', use 'object??' for extra details. In [1]: It looks like the ipy_user_conf.py file needs to be changed to the new format for ipython versions >=0.11 " Simon Ellingsen Active Tickets 272 Proble building asap-4.0.0 of fedora 17 (x64) General 2.0 Unified development defect major Malte Marquarding new 2012-07-05T13:10:05+10:00 2012-08-24T22:47:02+10:00 "I'm trying to build ASAP-4.0.0 on fedora 17 (x64). I built it on fedora 16 (32bit) and that was ok. See below for the error. Any clues would be appreciated. Thanks Kym [ {{{ root@maths14 asap]# cd asap-4.0.0/ [root@maths14 asap-4.0.0]# scons cfitsioincdir=/usr/include/cfitsio scons: Reading SConscript files ... Restoring previous command-line options from 'options.cache' Checking for C++ library casa_casa... (cached) yes Checking for C header file Python.h... (cached) yes Checking for C library python2.7... (cached) yes Checking for C++ library boost_python... (cached) yes Checking for C++ library pyrap... (cached) no Checking for C header file numpy/config.h... (cached) no Checking for C header file numpy/numpyconfig.h... (cached) yes Checking for C library m... (cached) yes Checking for C library cfitsio... (cached) yes Checking for C library wcs... (cached) yes Checking for C library rpfits... (cached) yes Checking for C library lapack... (cached) yes Checking for C library blas... (cached) yes Checking for C library gfortran... (cached) yes Checking for C++ library stdc++... (cached) yes scons: done reading SConscript files. scons: Building targets ... g++ -o build/STFitter.o -c -O2 -DWCSLIB_GETWCSTAB -fPIC -DHAVE_LIBPYRAP -DAIPS_USENUMPY -I/usr/include/cfitsio -I/usr/local/include/casacore -I/usr/include/python2.7 -I/usr/lib64/python2.7/site-packages/numpy/core/include -Iexternal-alma src/STFitter.cpp src/STFitter.cpp: In member function 'bool asap::Fitter::computeEstimate()': src/STFitter.cpp:121:30: error: request for member 'getAmpl' in 'listGauss.casa::SpectralList::operator[](i)', which is of pointer type 'casa::SpectralElement*' (maybe you meant to use '->' ?) src/STFitter.cpp:122:30: error: request for member 'getCenter' in 'listGauss.casa::SpectralList::operator[](i)', which is of pointer type 'casa::SpectralElement*' (maybe you meant to use '->' ?) src/STFitter.cpp:123:30: error: request for member 'getFWHM' in 'listGauss.casa::SpectralList::operator[](i)', which is of pointer type 'casa::SpectralElement*' (maybe you meant to use '->' ?) scons: *** [build/STFitter.o] Error 1 scons: building terminated because of errors. }}} " anonymous Active Tickets 288 SDFITS export import doesn't keep flags. General 2.0 Unified development defect major Malte Marquarding assigned 2013-03-08T10:38:44+11:00 2013-03-15T13:20:53+11:00 "SDFITS does not preserve flag data Hi there, I am writing out ASAP data in SDFITS format and have noticed that flagged scans or cycles become unflagged when read back into ASAP. ie. the process of writing to SDFITS file appears to delete flag information. If I save in ASAP format, the flags are preserved. But I need to write into SDFITS format because the data then go through livedata/gridzilla. Later. Andrew xxx " Malte Marquarding Active Tickets 294 Show IF and scan number of spectrum on flagplotter (from CAS-5140) python 2.0 Unified development enhancement minor Kana Sugimoto new 2013-07-24T20:00:40+10:00 2013-07-24T20:00:40+10:00 "Reporter: ALMA-CSV team (Shinya Komugi) the task sdflag will benefit greatly from some improvements, mostly owing to the fact that single dish data has many rows compared to the interferometric data. Below are few that are important for efficient data reduction. 1) 2) when using interactive=T, the plots should show what IF number and Scan is shown, instead of just the rownumber in the data table. " Kana Sugimoto Active Tickets 195 Move Lorentzian fitting codes to casacore c++ 2.0 Unified development enhancement normal Malte Marquarding new 2010-08-06T17:07:54+10:00 2010-08-06T17:07:54+10:00 "The following five codes for Lorentzian fitting, namely, Lorentzian1D.h, Lorentzian1D.tcc, Lorentzian1D2.tcc, Lorentzian1DParam.h and Lorentzian1DParam.tcc, are currently placed in the src/ directory of ASAP. Should be moved to casacore (under scimath/Functionals/). " WataruKawasaki Active Tickets 200 Routine for processing cross-scan data General 2.0 Unified development task normal Malte Marquarding new 2010-08-19T16:28:31+10:00 2012-09-04T08:41:11+10:00 "As the DFBs at Parkes do not have a total-power mode, one cannot use 'SPOT' scans (as done with the Multibeam correlator) to obtain flux calibration factors. However we can do 'spider-scans' (see below) and it would be nice to have a function in ASAP which could process a series of 'spider-scans' to extract flux calibration factors and other information. For the Multibeam correlator, total-power SPOT observations produce the following for each polarization: Beam FWHM ampR (amplitude of RA scans) ampD (amplitude of DEC scans) ON-OFF (baseline) The calibration factors for each pol are then obtained using the following formula: N x ((ON-OFF x 2) / (ampR + ampD)) N is the flux of the source in Jy (i.e., obtained from MIRIAD's calplot utility). AZ-EL ""spider"" scans (..., -90, -45, 0, 45, 90, ...) are referenced to an RA/DEC position. The scans are conducted along the Great Circle passing through the source." stacy.mader@… Active Tickets 215 how to get columns not returned by get_column_names() General 2.0 Unified development question normal Malte Marquarding assigned 2010-10-11T13:36:13+11:00 2010-10-12T10:28:17+11:00 When asap reads in one of my FITS files, what does it do with columns that are not returned by scantable.get_column_names()? More specifically, does it keep the data and can I get at them somehow? I'm attaching an example SDFITS file. kuiper@… Active Tickets 216 If your platform isn't listed - RedHat General 2.0 Unified development enhancement normal Malte Marquarding new 2010-10-11T13:53:40+11:00 2010-10-11T13:53:40+11:00 "This is an IBM Blade Server we are using. [sweston@ska-blade1 ~]$ uname -a Linux ska-blade1.aut.ac.nz 2.6.18-164.11.1.el5 #1 SMP Wed Jan 6 13:26:04 EST 2010 x86_64 x86_64 x86_64 GNU/Linux In the list of required packages, yum isn't finding some of them: libstdc++ libgfortran lapack blas wget wcslib Some of the others fail with dependency issues: python-matplotlib-tk Trying to casacore: [sweston@ska-blade1 casacore-1.0.1]$ scons --cfitsio-incdir=/usr/include/cfitsio scons: Reading SConscript files ... Restoring previous command-line options from 'options.cache' Checking for C library dl... yes Checking for C library gfortran... yes Checking for C library blas... no " stuart.weston@… Active Tickets 217 How does auto_quotient(mode='time') recognize scans? General 2.0 Unified development question normal Malte Marquarding assigned 2010-10-19T05:48:52+11:00 2010-10-19T09:36:03+11:00 "What information in an SDFITS record tells ASAP that a scan is an OFF or an ON? The only thing I see in the column names that might do it is SRCTYPE. If that is what it is, what are acceptable values, and what SDFITS column name would map into that in the scantable?" kuiper@… Active Tickets 220 Leap second table TAI_UTC seems out-of-date. General 2.0 Unified development question normal Malte Marquarding new 2010-11-10T10:07:05+11:00 2010-11-10T10:07:05+11:00 I'm trying to read in SDFITS files created from 1989-1995 Tid data and I get this message. Is there a work-around or fix? kuiper@… Active Tickets 221 Retaining scan numbers General 2.0 Unified development question normal Malte Marquarding new 2010-11-10T10:55:47+11:00 2011-01-11T16:13:35+11:00 Can I get ASAP to use the scan numbers that appear in the SCAN column? If not, is the information available somewhere in the scantable? anonymous Active Tickets 222 ERROR: SDFITS file contains invalid beam number. General 2.0 Unified development question normal Malte Marquarding new 2010-11-12T12:39:52+11:00 2012-09-19T09:55:05+10:00 "What does this error mean? From my data: In [6]: for i in range(32): print tbdata[i]['BEAM'], ...: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 From the header: TTYPE16 = 'BEAM ' TFORM16 = 'I ' " kuiper@… Active Tickets 223 polarization numbers General 2.0 Unified development question normal Malte Marquarding new 2010-11-23T09:58:39+11:00 2010-11-23T09:58:39+11:00 "This is a question but it might turn into a defect. How does asap get its polnos from an SDFITS file? This allows me to get at the table header and columns: ASAP>import Data_Reduction.FITS.superasap as S The first is a scantable with the table HDU atached. The procedure groups scans. ASAP>tab_1994_338, scan_groups, aves = \ S.report_Tid_SDFITS(""1994/1994_338.fits"",plot=True) Whar ASAP thinks the polarization codes are: ASAP>tab_1994_338.getpolnos() asap>:(0,) What the polarization codes are in the SDFITS table: ASAP>tab_1994_338.data.field('CRVAL4') asap>: array([-2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2], dtype=int16) ASAP>tab_1994_338.header['CTYPE4'] asap>:'STOKES' " kuiper@… Active Tickets 226 independently scaled subplots General 2.0 Unified development enhancement normal Malte Marquarding new 2010-12-07T05:16:39+11:00 2010-12-07T05:16:39+11:00 For situations in which one merely wants to have an overview of what is in a scantable, and not compare spectra particularly, it would be good to have an option in which the subplots are independently scaled. I believe that is the matplotlib default so it might just involve bypassing some asap code. kuiper@… Active Tickets 227 ignoring channels when scaling a plot General 2.0 Unified development question normal Malte Marquarding new 2010-12-07T08:25:09+11:00 2010-12-07T08:25:09+11:00 How can I tell the plotter to ignore the masked channels in a scan table when scaling the plot? kuiper@… Active Tickets 228 sub-scantables in scantable General 2.0 Unified development defect normal Malte Marquarding new 2010-12-11T08:06:15+11:00 2010-12-15T09:21:07+11:00 "Reading a 50 row SDFITS extension into asap resulted in 45 'scans' but three of them are mini-scantables. However, all the get_[something]() methods return 50 values consistent with the SDFITS original. Only getscannos() gives a list of scans that can be fetched with get_scan() and that summary() will show. Also something funny happens with the scan numbers, which go from 0 to 51 (instead of 49) with gaps not related to the scan grouping into miniscans. I will attach a more detailed discussion and the SDFITS file." kuiper@… Active Tickets 229 Installing from source General 2.0 Unified development question normal Malte Marquarding new 2011-01-11T09:09:47+11:00 2011-01-11T09:09:47+11:00 "We have a IBM Blade Server (RH): Linux ska-blade1.aut.ac.nz 2.6.18-164.11.1.el5 #1 SMP Wed Jan 6 13:26:04 EST 2010 x86_64 x86_64 x86_64 GNU/Linux Trying to follow the instructions to install from source. Big issues on dependencies. At the moment got to casacore, ok this needs blas which in turn has other dependencies (rpmlib libgfortran ....). " stuart.weston@… Active Tickets 234 improve asapplotter by using Matplotlib 1.0.0 functions General 2.0 Unified development defect normal Kana Sugimoto assigned 2011-02-25T12:43:23+11:00 2011-02-25T12:43:44+11:00 "Matplotlib 1.0.0 has improvements in plotting subplots etc. Look into the new APIs and update asapplotter." Kana Sugimoto Active Tickets 255 RuntimeError: Couldn't convert frequency frame General 2.0 Unified development defect normal Malte Marquarding assigned 2011-10-13T15:49:27+11:00 2011-10-17T10:40:52+11:00 "Hello, I am still going not far with ASAP, after some update this is my return: jose@ubuntu:~$ sudo asap_update_data [sudo] password for jose: Checking if an update is required. Data already at latest available version. If you still get errors running asap, please report this. jose@ubuntu:~$ asap Loading ASAP... Welcome to ASAP v4.0.0 (2011-10-05) - the ATNF Spectral Analysis Package Please report any bugs via: http://svn.atnf.csiro.au/trac/asap/simpleticket [IMPORTANT: ASAP is 0-based] Type commands() to get a list of all available ASAP commands. ASAP>s = scantable('2008-03-12_0932-M999.rpf') Found ATMOPRA data. Auto averaging integrations Importing 2008-03-12_0932-M999.rpf... ASAP>print s --------------------------------------------------------------------------- RuntimeError Traceback (most recent call last) /home/jose/ in () /usr/lib/pymodules/python2.7/asap/scantable.pyc in __str__(self) 438 def __str__(self): 439 tempFile = tempfile.NamedTemporaryFile() --> 440 Scantable._summary(self, tempFile.name) 441 tempFile.seek(0) 442 asaplog.clear() RuntimeError: Couldn't convert frequency frame. ASAP>s = scantable('2008-03-12_0932-M999.rpf') WARN: Requested data table IERSeop97 cannot be found in the searched directories: ./ ./data/ /home/jose/aips++/data//ephemerides/ /usr/lib/pymodules/python2.7/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jose/aips++/data//geodetic/ /usr/lib/pymodules/python2.7/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ Cannot read IERS table IERSeop97 Calculations will proceed with lower precision WARN: Requested data table IERSpredict cannot be found in the searched directories: ./ ./data/ /home/jose/aips++/data//ephemerides/ /usr/lib/pymodules/python2.7/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jose/aips++/data//geodetic/ /usr/lib/pymodules/python2.7/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ Cannot read IERS table IERSpredict Calculations will proceed with lower precision No requested dUT1 data available from IERS tables. Proceeding with probably less precision. WARN: Requested data table TAI_UTC cannot be found in the searched directories: ./ ./data/ /home/jose/aips++/data//ephemerides/ /usr/lib/pymodules/python2.7/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jose/aips++/data//geodetic/ /usr/lib/pymodules/python2.7/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ SEVERE: Cannot read leap second table TAI_UTC WARN: Unable to convert between the input and output SpectralCoordinates this probably means one is in the REST frame which requires the radial velocity - this is not implemented yet WARN: Requested data table TAI_UTC cannot be found in the searched directories: ./ ./data/ /home/jose/aips++/data//ephemerides/ /usr/lib/pymodules/python2.7/asap/data//ephemerides/ /usr/local/share/casacore/data/ephemerides/ /usr/local/share/casacore/data/ephemerides/ /home/jose/aips++/data//geodetic/ /usr/lib/pymodules/python2.7/asap/data//geodetic/ /usr/local/share/casacore/data/geodetic/ /usr/local/share/casacore/data/geodetic/ SEVERE: Cannot read leap second table TAI_UTC WARN: Unable to convert between the input and output SpectralCoordinates this probably means one is in the REST frame which requires the radial velocity - this is not implemented yet Found ATMOPRA data. Auto averaging integrations Importing 2008-03-12_0932-M999.rpf... Thanks for any light... ASAP> " jvazquez@… Active Tickets 280 possible incorrect sdfits headers General 2.0 Unified development defect normal Malte Marquarding assigned 2012-10-03T14:24:33+10:00 2012-10-03T14:35:27+10:00 "I am trying to process Parkes mosaicing data partially through ASAP. The standard route is to use livedata, then gridzilla. But ASAP is preferred partially due to the superior baseline fitter (there are other reasons). I have tried two methods to reduce the data, whose steps are outlined below: Method1: 1. Livedata - read in rpfits, flag off-source data, quotient, baseline, write sdfits 2. Gridzilla - read in sdfits, map, write fits data cube Method2: 1. Livedata - read in rpfits, flagg off-source data, write sdfits 2. ASAP - remove RFI, quotient, baseline, write sdfits 3. Gridzilla - read in sdfits, map, write fits data cube The problem is that the resultant cubes for the two methods are offset in velocity by about 1.18km/s. I suspect that the ASAP stage is either incorrectly writing header information or writing it in a format that subsequent gridzilla processing does not recognise. I come to this conclusion partly because in the output sdfits header from ASAP (Method2), one line says: SPECSYS = 'TOPOCENT' / Doppler reference frame (transformed) whilst the sdfits output from livedata (Method1) says: SPECSYS = 'LSRK ' / Doppler reference frame (transformed) Note that manually editing the sdfits file from ASAP to read ""LSRK"" does NOT change the final data cube, and it is still offset in velocity. I have attached the following files: raw rpfits file (2012-05-13_1506-P817_G335.0-1.0_lon_dfb4.rpf) sdfits file made using livedata in step 1 of Method1 (method1.sdfits) sdfits file made using livedata in step 1 of Method2 (2012-05-13_1506-P817_G335.0-1.0_lon_dfb4.sdfits) sdfits file made using ASAP in step 2 of Method2 (1506.sdfits) ASAP script used to produce 1506.sdfits above (red.py) I have included the files in a tar file. Later. Andrew xxx " Andrew Walsh Active Tickets 282 test_nrofiller fails General 2.0 Unified development defect normal Takeshi Nakazato new 2012-10-17T13:04:34+11:00 2012-10-22T12:57:05+11:00 "{{{ mar637@loki:test> nosetests test_nrofiller.py E... ====================================================================== ERROR: test_nrofiller.TestNRO.test_frequency ---------------------------------------------------------------------- Traceback (most recent call last): File ""/usr/lib/pymodules/python2.6/nose/case.py"", line 183, in runTest self.test(*self.arg) File ""/home/mar637/repos/asap/trunk/test/test_nrofiller.py"", line 59, in test_frequency abc=self.st._getabcissa(0) RuntimeError: Couldn't convert frequency frame. ---------------------------------------------------------------------- Ran 4 tests in 0.453s FAILED (errors=1) }}}" Malte Marquarding Active Tickets 285 error in frequency frame conversion, old data, old script, new install. General 2.0 Unified development defect normal Malte Marquarding assigned 2012-10-23T17:27:51+11:00 2012-10-23T17:56:55+11:00 "Welcome to ASAP v4.0.0 (2011-10-05) - the ATNF Spectral Analysis Package ASAP>from asap import * ASAP>file='2008-03-12_0932-M999.rpf' ASAP>s = scantable(file, average=True) Found ATMOPRA data. Auto averaging integrations Importing 2008-03-12_0932-M999.rpf... ASAP>print s ---->print(s) --------------------------------------------------------------------------- RuntimeError Traceback (most recent call last) /Users/erik/ in () /Library/Python/2.6/site-packages/asap-4.0.0-py2.6.egg/asap/scantable.pyc in __str__(self) 438 def __str__(self): 439 tempFile = tempfile.NamedTemporaryFile() --> 440 Scantable._summary(self, tempFile.name) 441 tempFile.seek(0) 442 asaplog.clear() RuntimeError: Couldn't convert frequency frame. " emuller Active Tickets 287 Improve pointing plots General 2.0 Unified development enhancement normal Kana Sugimoto assigned 2012-12-12T13:58:33+11:00 2013-07-24T20:16:06+10:00 "Inprove plots by asapplotter.plotpointings: - enable plotting by src type - better axis labels" Kana Sugimoto Active Tickets 290 Implement sideband separation algorithm in ASAP (from CAS-4141) General 2.0 Unified development defect normal Kana Sugimoto assigned 2013-03-15T11:29:39+11:00 2013-07-24T20:22:48+10:00 "ALMA Band-9 and 10 is observed with double side band receiver. It is necessary to develop a way to supress and/or separate signals from the sidebands. The algorithm to separate sidebands recommended by Commisioning Science Verification team of ALMA is based on Fourier Transform technique (Emerson, Klein, and Haslam 1979, A&A 76, 92)." Kana Sugimoto Active Tickets 291 alignment error for Tid data General 2.0 Unified development defect normal Malte Marquarding assigned 2013-03-22T13:48:37+11:00 2013-07-02T09:15:45+10:00 "The NH3 features of the averaged spectrum after processing don't line up. There is a shift of approx 50km/s. It looks like an issue with frequency alignment of data from several epochs. {{{ The follow is the process before merging a1 = sd.scantable('2011_260_t199/2011-09-17_074510_T199.rpf') # lupus tid 2011 260 some bad scans, but 1,1 and 2,2 present LupIR a2 = sd.scantable('2011_260_t199/2011-09-17_083236_T199.rpf') # lupus tid 2011 260 nothing LupIR a3 = sd.scantable('2011_260_t199/2011-09-17_092613_T199.rpf') # lupus tid 2011 260 mostly nothing with 1 scan with 2,2 LupIR (0-25) and LupINW (26-42) a4 = sd.scantable('2012_087_t199/2012-03-27_194022_T199.rpf') # lupus tid 2012 087 some 1,1 (LupIR, only one scan) a5Scan = sd.scantable('2012_091_t199/2012-03-31_182325_T199.rpf') # lupus tid 2012 091 noisy but some 1,1 transition LupINW (0-8) LupISE (9-34) LupISW (35-84) LupIIIMMS (85-109) LupIIIS (110-134) LupIIIN(135-159) LupIV(160-182) a5 = a5Scan.get_scan(range(0,85)) a6 = sd.scantable('2012_092_t199/2012-04-01_185240_T199.rpf') # lupus tid 2012 092 bad start, nothing det LISW a7 = sd.scantable('2012_092_t199/2012-04-01_203555_T199.rpf') # lupus tid 2012 092 1,1 overall noise is bad LupIR (0,1) LupISW (2 scan) a8Scan = sd.scantable('2012_099_t199/2012-04-08_182004_T199.rpf') # lupus tid 2012 099 good noise lvl, few 1,1 and one 2,2 last scan bit dodgy LupIR (0-1) LupISW (2-7) LupIIIMMS(8-21) a8 = a8Scan.get_scan(range(0,7)) a9 = sd.scantable('2012_172_t199/2012-06-20_095048_T199.rpf') # lupus tid 2012 172 Lupus I NW (0-7) central hyperfine detected 2 scans (LupINW01-02) a10 = sd.scantable('2012_172_t199/2012-06-20_101136_T199.rpf') # lupus tid 2012 172 Lupus I NW (0-17) and SE (18-35), possible central hyperfine detected. a11 = sd.scantable('2012_182_t199/2012-06-30_165737_T199.rpf') # lupus tid 2012 182 no detection covering LupISE (0-6) and SW (7-36) a12 = sd.scantable('2012_209_t199/2012-07-27_150854_T199.rpf') # lupus tid 2012 209 LupI SW only one scan src a13 = sd.scantable('2012_209_t199/2012-07-27_151516_T199.rpf') # lupus tid 2012 209 LupI SW 0-14 nothing a14Scan = sd.scantable('2012_209_t199/2012-07-27_154021_T199.rpf') # lupus tid 2012 209 LupI SW 0-19 Lup3MMS 20-40 (scan 39-33 1,1 dectection) a14 = a14Scan.get_scan(range(0,19)) a15 = sd.scantable('2012_251_t199/2012-09-07_112832_T199.rpf') # lupus tid 2012 251 Lupus IR reasonable 1,1 detection a16 = sd.scantable('2012_251_t199/2012-09-07_114049_T199.rpf') # lupus tid 2012 251 Lupus IR (0,1) detection 1,1 and 2,2 and sw a17 = sd.scantable('2012_251_t199/2012-09-07_115617_T199.rpf') # lupus tid 2012 251 Lupus IR #the arrays you want to process separately and merge arrayScan = [a1,a2,a3,a4,a6,a7,a8,a9,a10,a11,a13,a14,a15,a16,a17] qArray = []#leave this blank for merging the scan tables counter = 0#counter for naming convention for i in arrayScan: qString = 'q'+str(counter)#naming convention for the qArray and merging scan tabled qString = i.auto_quotient()# Make the quotient spectra qString.set_freqframe('LSRK')# Plot/select in velocity qString.set_unit('km/s') # Correct for gain/el effects qString.recalc_azel()# Tid does not write the elevation qString.gain_el() qString.opacity(0.05) qString.freq_align() qArray.append(qString)#need this for merging the scanTables counter = counter + 1 totalScan = sd.merge(qArray) #merging the scan tables print(totalScan) totalScan.set_restfreqs([23694.4700e6]) # 1,1 dectection cropScan = totalScan.average_pol() msk = cropScan.create_mask([-100,-70],[70,100]) cropScan.poly_baseline(mask=msk, order=1) sd.plotter.plot(cropScan) sd.plotter.set_layout(4,4) #change the format of plotter.plot sd.plotter.set_histogram() sd.plotter.plot(cropScan) cropScan.auto_poly_baseline() sd.plotter.plot() }}}" Malte Marquarding Active Tickets 295 asapplotter.plotpointing2 needs improvement (from CAS-5317) python 2.0 Unified development enhancement normal Kana Sugimoto new 2013-07-24T20:05:12+10:00 2013-07-24T20:05:12+10:00 "asapplotter.plotpointing2 plots the pointings in absolute degrees. The precision is not enough to figure out the pointing spacing from the plot. Please add an option to make the plot in *relative arcseconds* " Kana Sugimoto Active Tickets 296 Enable polarization average of spectra without time averating. (from CAS-5194) c++ 2.0 Unified development enhancement normal Kana Sugimoto assigned 2013-07-24T20:08:05+10:00 2013-07-24T20:10:56+10:00 "Reporter: Erik Muller (EA-ARC) polarizations are not averaged correctly together -- currently, they include temporal average by default. " Kana Sugimoto Active Tickets 209 Change FIT PARAMETER Column from Double to Float c++ 2.0 Unified development enhancement normal Malte Marquarding assigned 2010-09-14T12:14:24+10:00 2011-10-12T10:37:04+11:00 "The fitting is done as Float, so no need for conversions. This breaks the Scantable version so should be done when we change the schema." Malte Marquarding