Changeset 937 for trunk/test/tid.py


Ignore:
Timestamp:
03/27/06 19:36:17 (18 years ago)
Author:
mar637
Message:

asap2 version of test script

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/test/tid.py

    r794 r937  
    44rcParams['verbose'] = 0
    55
     6rcParams['plotter.ganged'] = 0
    67# Don't plot to the screen...
    78del plotter
     
    1213# Create the quotient spectra
    1314data = scantable('data/tid-t002.rpf')
    14 s = data.get_scan('*[^we]')
    15 r = data.get_scan('*[we]')
    16 q = quotient(s,r)
    17 
     15q = data.auto_quotient()
    1816# Set the restfreq for each IF
    1917q.set_restfreqs(freqs= [23694.4700e6,23722.6336e6])
    2018q.set_unit('km/s')
    2119q.set_freqframe('LSRK')
    22 
    2320# Align frequencies
    24 
    25 q.freq_align(perif=True)
    26 
     21qa = q.freq_align(insitu=False)
    2722# Recalculate the az/el
    28 
    29 q.recalc_azel()
    30 
     23qa.recalc_azel()
    3124# Correct for gain curve and opacity
    32 q.gain_el()
    33 q.opacity(0.075)
    34 
     25qa.gain_el()
     26qa.opacity(0.075)
    3527# Average in time
    36 av = average_time(q)
    37 
     28av = qa.average_time()
    3829# Baseline
    3930msk=av.create_mask([-70,-50],[40,60])
    4031av.poly_baseline(msk,1)
    41 
    42 plotter.set_mode('i')
     32plotter.set_mode('i','s')
    4333plotter.plot(av)
    44 plotter.save('output/tid.png')
     34plotter.save('output/tid.png', dpi=80)
    4535
    4636# Do some random processing, just to test these functions
Note: See TracChangeset for help on using the changeset viewer.