Last change
on this file since 748 was
746,
checked in by phi196, 19 years ago
|
Test Tid data
|
-
Property svn:eol-style set to
native
-
Property svn:executable set to
*
-
Property svn:keywords set to
Author Date Id Revision
|
File size:
799 bytes
|
Line | |
---|
1 | #!/usr/bin/env python |
---|
2 | from asap import * |
---|
3 | |
---|
4 | rcParams['verbose'] = 0 |
---|
5 | |
---|
6 | # Don't plot to the screen... |
---|
7 | del plotter |
---|
8 | plotter = asapplotter(False) |
---|
9 | |
---|
10 | print "Test of Tidbinbilla" |
---|
11 | |
---|
12 | # Create the quotient spectra |
---|
13 | data = scantable('data/tid-t002.rpf') |
---|
14 | s = data.get_scan('*[^we]') |
---|
15 | r = data.get_scan('*[we]') |
---|
16 | q = quotient(s,r) |
---|
17 | av = average_time(q) |
---|
18 | |
---|
19 | # Set the restfreq for each IF |
---|
20 | av.set_restfreqs(freqs= [23694.4700e6,23722.6336e6]) |
---|
21 | av.set_unit('km/s') |
---|
22 | av.set_freqframe('LSRK') |
---|
23 | av.smooth('gauss',5) |
---|
24 | av.scale(1.05) |
---|
25 | av.add(0.05) |
---|
26 | |
---|
27 | # Baseline |
---|
28 | msk=av.create_mask([-70,-50],[40,60]) |
---|
29 | av.poly_baseline(msk,1) |
---|
30 | |
---|
31 | plotter.plot(av) |
---|
32 | plotter.set_mode('i') |
---|
33 | plotter.save('output/tid.png') |
---|
34 | |
---|
35 | # These are currently broken as Tid does not set the elevation correctly |
---|
36 | av.gain_el() |
---|
37 | av.opacity(0.075) |
---|
38 | |
---|
39 | print "Tid test finished successfully" |
---|
Note: See
TracBrowser
for help on using the repository browser.