Changeset 937 for trunk/test/parkes-pol.py
- Timestamp:
- 03/27/06 19:36:17 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/test/parkes-pol.py
r744 r937 3 3 4 4 rcParams['verbose'] = 0 5 rcParams['plotter.ganged'] = 0 5 6 6 7 # Don't plot to the screen... … … 16 17 data_1665.set_freqframe('LSRK') 17 18 18 # Make a copy for the 1667 transition19 data_1667 = data_1665.copy() 20 data_1667.set_restfreqs(lines = ['OH1667'])21 data_1667.set_unit('km/s')22 data_166 7.set_freqframe('LSRK')19 # Look at the first scan 20 from asap._asap import selector 21 selection = selector() 22 selection._setscans([0]) 23 data_1665._setselection(selection) 23 24 25 d1_5 = data_1665.copy() 26 d1_7 = data_1665.copy() 24 27 25 # Look at the first scan 26 d1_5 = data_1665.get_scan(0) 27 d1_7 = data_1667.get_scan(0) 28 d1_7.set_restfreqs([1667.0],'MHz') 29 # merge the two scans back together into a new scantable 30 plotscans = merge(d1_5,d1_7) 31 print plotscans.summary() 32 del d1_5,d1_7 28 33 29 34 # Baseline 30 msk = d1_5.create_mask([-30,-25],[-5,0]) 31 d1_5.poly_baseline(msk,1) 32 msk = d1_7.create_mask([-30,-25],[-5,0]) 33 d1_7.poly_baseline(msk,1) 34 35 msk = plotscans.create_mask([-30,-25],[-5,0]) 36 plotscans.poly_baseline(msk,1) 35 37 # Plot the results 36 plotter.plot(d1_5,d1_7)37 38 plotter.set_mode('p','s') 38 39 plotter.set_layout(2,1) 39 40 plotter.set_range(-30,0) 40 plotter.set_cursor(pol=['I','Q','U','V']) 41 plotter.set_cursor(pol=['I','Plinear']) 42 plotter.set_cursor(pol=['RR','LL']) 43 plotter.save('output/parkes_rrll.png') 41 selection._reset() 42 selection._setpolstrings(['I','Q', 'U', 'V']) 43 #plotter.set_selection(selection) 44 #selection._setpolstrings(['I','Plinear']) 45 #plotter.set_selection(selection) 46 #selection._setpolstrings(['RR','LL']) 47 plotter.plot(plotscans) 48 plotter.set_selection(selection) 49 plotter.save('/nfs/wwwpeople/mmarquar/parkes.png', dpi=80) 50 #plotter.save('output/parkes_rrll.png') 44 51 45 52 print "Parkes-Pol Test successful"
Note: See TracChangeset
for help on using the changeset viewer.