Custom Query (241 matches)

Running skydip with plot=True produces the following:

taus=skydip("2012-02-11_0500-P000_1568_SKYDIP.rpf",plot=True)

IndexError? Traceback (most recent call last)

/Users/mad134/Home/pyspot/20120211/<ipython console> in <module>()

/Library/Python?/2.6/site-packages/asap-4.0.0-py2.6.egg/asap/opacity.pyc in skydip(data, averagepol, tsky, plot, temperature, pressure, humidity)

197 n = len(airms) 198 for i in range(n):

--> 199 pylab.plot(airms[i], tsyss[i], 'o', color=colors[i])

200 pylab.plot(airms[i], fits[i], '-', color=colors[i]) 201 pylab.figtext(0.7,0.3-(i/30.0),

IndexError?: list index out of range

Using ASAP v4.0.0 (2011-10-05), Mac OSX 10.6.8.

After upgrading to ASAP 4.0, running a script reported the following:

Leopard libedit detected. Found AT???? data. Importing 20111118/2011-11-18_0338-P000_SPOT_3100.rpf...

TypeError? Traceback (most recent call last)

/Users/mad134/Home/pyspot/proc.py in <module>()

89 errmsg('Will only process the first four scans.') 90

---> 91 if(nu/1.0e9 >= 12.0):

92 Diam = 54.5 # Parkes aperture [m] 93 else:

TypeError?: unsupported operand type(s) for /: 'list' and 'float'

The parameter nu is defined as nu = sc.get_restfreqs()[0]. ASAP reports this as:

In [2]: scan.get_restfreqs()[0] Out[2]: [3100000000.0]

So nu is now returned as a list instead of a float as per the previous version of ASAP. Is this a bug or by design in the new version?

Hi Malte,

So I've just tried a few tests in asap and I am fairly confident that when you use get_scan the resultant data file forgets the parallactify option - I tried two simple methods, first setting parallactify on the initial file:

ASAP>data_1665 = scantable('2005-10-27_0154-P484.rpf')
Cross polarization present
Found ATPKSHOH data.
Auto averaging integrations
Importing 2005-10-27_0154-P484.rpf...

ASAP>data_1665.parallactify(True)

ASAP>sourcedata=data_1665.get_scan(55)

ASAP>sourcedata.set_unit('km/s')

ASAP>dmsk2=sourcedata.create_mask([-30,-20],[20,30])
The current mask window unit is km/s

ASAP>sourcedata.poly_baseline(dmsk2,1)

ASAP>plotter.plot(sourcedata)
A new scantable is set to the plotter. The masks and data selections are reset.

ASAP>stokesp = sourcedata.convert_pol('stokes')

ASAP>plotter.plot(stokesp)
A new scantable is set to the plotter. The masks and data selections are reset.

ASAP>sourcedata.parallactify(True)

ASAP>stokesp = sourcedata.convert_pol('stokes')

ASAP>plotter.plot(stokesp)
A new scantable is set to the plotter. The masks and data selections are reset.

*This gives the wrong stokes q and u. The second has parallactify after the get_scan and gives the correct q and u:

ASAP>data_1665 = scantable('2005-10-27_0154-P484.rpf')
Cross polarization present
Found ATPKSHOH data.
Auto averaging integrations
Importing 2005-10-27_0154-P484.rpf...

ASAP>sourcedata=data_1665.get_scan(55)

ASAP>sourcedata.parallactify(True)

ASAP>sourcedata.set_unit('km/s')

ASAP>dmsk2=sourcedata.create_mask([-30,-20],[20,30])
The current mask window unit is km/s

ASAP>sourcedata.poly_baseline(dmsk2,1)

ASAP>plotter.plot(sourcedata)
A new scantable is set to the plotter. The masks and data selections are reset.

ASAP>stokesp = sourcedata.convert_pol('stokes')

ASAP>plotter.plot(stokesp)
A new scantable is set to the plotter. The masks and data selections are reset.




Cheers,

Jimi