Changeset 937 for trunk/test/tid.py
- Timestamp:
- 03/27/06 19:36:17 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/test/tid.py
r794 r937 4 4 rcParams['verbose'] = 0 5 5 6 rcParams['plotter.ganged'] = 0 6 7 # Don't plot to the screen... 7 8 del plotter … … 12 13 # Create the quotient spectra 13 14 data = scantable('data/tid-t002.rpf') 14 s = data.get_scan('*[^we]') 15 r = data.get_scan('*[we]') 16 q = quotient(s,r) 17 15 q = data.auto_quotient() 18 16 # Set the restfreq for each IF 19 17 q.set_restfreqs(freqs= [23694.4700e6,23722.6336e6]) 20 18 q.set_unit('km/s') 21 19 q.set_freqframe('LSRK') 22 23 20 # Align frequencies 24 25 q.freq_align(perif=True) 26 21 qa = q.freq_align(insitu=False) 27 22 # Recalculate the az/el 28 29 q.recalc_azel() 30 23 qa.recalc_azel() 31 24 # Correct for gain curve and opacity 32 q.gain_el() 33 q.opacity(0.075) 34 25 qa.gain_el() 26 qa.opacity(0.075) 35 27 # Average in time 36 av = average_time(q) 37 28 av = qa.average_time() 38 29 # Baseline 39 30 msk=av.create_mask([-70,-50],[40,60]) 40 31 av.poly_baseline(msk,1) 41 42 plotter.set_mode('i') 32 plotter.set_mode('i','s') 43 33 plotter.plot(av) 44 plotter.save('output/tid.png' )34 plotter.save('output/tid.png', dpi=80) 45 35 46 36 # Do some random processing, just to test these functions
Note: See TracChangeset
for help on using the changeset viewer.