| 1 | from asap import *
|
|---|
| 2 |
|
|---|
| 3 | # Load in the data
|
|---|
| 4 |
|
|---|
| 5 | #scans=scantable("/nfs/PKCCC3_4/2009-06-20_0527_ON2.rpf")
|
|---|
| 6 | scans=scantable("/nfs/PKCCC3_4/2009-06-21_0322_ON2.rpf")
|
|---|
| 7 |
|
|---|
| 8 | # Set some parameters
|
|---|
| 9 |
|
|---|
| 10 | scans.set_freqframe('LSRK')
|
|---|
| 11 | scans.set_doppler('RADIO')
|
|---|
| 12 | scans.set_unit('km/s')
|
|---|
| 13 | #scans.set_unit('GHz')
|
|---|
| 14 |
|
|---|
| 15 | #H20
|
|---|
| 16 | #scans.set_restfreqs(freqs=22235.08e6)
|
|---|
| 17 |
|
|---|
| 18 | #NH3
|
|---|
| 19 | scans.set_restfreqs(freqs=23694.496e6)
|
|---|
| 20 |
|
|---|
| 21 | # Do the freq switching
|
|---|
| 22 |
|
|---|
| 23 | scans.freq_switch()
|
|---|
| 24 |
|
|---|
| 25 | # Average the polarisations together
|
|---|
| 26 | iav=scans.average_pol()
|
|---|
| 27 |
|
|---|
| 28 | #plotter.set_range(-60,60)
|
|---|
| 29 | plotter.set_range(-100,300)
|
|---|
| 30 | #iav.auto_poly_baseline()
|
|---|
| 31 |
|
|---|
| 32 | plotter.plot(iav)
|
|---|
| 33 |
|
|---|
| 34 | #outname=str(fnames)
|
|---|
| 35 | #outname1=outname[32:44]
|
|---|
| 36 |
|
|---|
| 37 | #plotter.save('save_images/'+outname1+'.ps')
|
|---|
| 38 |
|
|---|
| 39 | plotter.save('save_images/2009-06-21_0322_ON2_fqcsw.ps')
|
|---|
| 40 |
|
|---|
| 41 | #iav = iav.save('IRDC316.76-1','CLASS','overwrite')
|
|---|
| 42 |
|
|---|
| 43 | #iav.save('raw_ascii/2009-06-21_0322_ON2_fqcsw_','ascii')
|
|---|
| 44 |
|
|---|
| 45 | del iav
|
|---|
| 46 | del scans |
|---|
