source: trunk/test/tid.py@ 951

Last change on this file since 951 was 950, checked in by mar637, 19 years ago

re-enebled insitu freq_align

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 902 bytes
RevLine 
[746]1#!/usr/bin/env python
2from asap import *
3
4rcParams['verbose'] = 0
5
[937]6rcParams['plotter.ganged'] = 0
[746]7# Don't plot to the screen...
8del plotter
9plotter = asapplotter(False)
10
11print "Test of Tidbinbilla"
12
13# Create the quotient spectra
14data = scantable('data/tid-t002.rpf')
[937]15q = data.auto_quotient()
[752]16# Set the restfreq for each IF
17q.set_restfreqs(freqs= [23694.4700e6,23722.6336e6])
18q.set_unit('km/s')
19q.set_freqframe('LSRK')
[794]20# Align frequencies
[950]21q.freq_align()
[794]22# Recalculate the az/el
[950]23q.recalc_azel()
[794]24# Correct for gain curve and opacity
[950]25q.gain_el()
26q.opacity(0.075)
[752]27# Average in time
[950]28av = q.average_time()
[746]29# Baseline
30msk=av.create_mask([-70,-50],[40,60])
31av.poly_baseline(msk,1)
[937]32plotter.set_mode('i','s')
[746]33plotter.plot(av)
[937]34plotter.save('output/tid.png', dpi=80)
[746]35
[794]36# Do some random processing, just to test these functions
37av.smooth('gauss',5)
38av.scale(1.05)
39av.add(0.05)
[746]40
41print "Tid test finished successfully"
Note: See TracBrowser for help on using the repository browser.