| 1 | asap0000
|
|---|
| 2 | ========
|
|---|
| 3 | Hi
|
|---|
| 4 |
|
|---|
| 5 | Jim Lovell has discovered a strange corruption bug in asap.
|
|---|
| 6 |
|
|---|
| 7 | When we run:
|
|---|
| 8 |
|
|---|
| 9 | scans = scantable('2005-03-09_172811_t002.rpf')
|
|---|
| 10 | scans.summary(verbose=True)
|
|---|
| 11 |
|
|---|
| 12 | repeatably the source position become corrupted (generally it seems to be
|
|---|
| 13 | after the third "scantable" construction.
|
|---|
| 14 |
|
|---|
| 15 | The above data is in /DATA/VULCAN_1/cphillip/vulcan/data/asap, but I think
|
|---|
| 16 | any data will do.
|
|---|
| 17 |
|
|---|
| 18 | Cheers
|
|---|
| 19 | Chris
|
|---|
| 20 |
|
|---|
| 21 |
|
|---|
| 22 | ---------------------------------------------------------
|
|---|
| 23 | **OK
|
|---|
| 24 |
|
|---|
| 25 | ASAP> scans.summary(verbose=True)
|
|---|
| 26 |
|
|---|
| 27 | --------------------------------------------------------------------------------
|
|---|
| 28 | Scan Table Summary
|
|---|
| 29 | --------------------------------------------------------------------------------
|
|---|
| 30 | Beams: 1
|
|---|
| 31 | IFs: 2
|
|---|
| 32 | Polarisations: 1
|
|---|
| 33 | Channels: 4096
|
|---|
| 34 |
|
|---|
| 35 | Observer: Auto_Spec
|
|---|
| 36 | Obs Date: 2005/03/09/17:29:22
|
|---|
| 37 | Project: NGC6334Fw
|
|---|
| 38 | Obs. Type:
|
|---|
| 39 | Antenna Name: DSS-43
|
|---|
| 40 | Flux Unit: K
|
|---|
| 41 | Rest Freqs: [2.39639e+10] [Hz]
|
|---|
| 42 | Abcissa: Channel
|
|---|
| 43 | Cursor: Beam[0] IF[0] Pol[0]
|
|---|
| 44 |
|
|---|
| 45 | Scan Source Position (J2000) Time Integration
|
|---|
| 46 | FreqIDs
|
|---|
| 47 | --------------------------------------------------------------------------------
|
|---|
| 48 | 0 NGC6334Fw 17:18:23.4 -35.47.01.3 17:30:37 1 x 02:31.0
|
|---|
| 49 | [0, 1] [0]
|
|---|
| 50 | 1 NGC6334F 17:20:53.4 -35.47.01.3 17:33:43 1 x 02:33.9
|
|---|
| 51 | [0, 1] [0]
|
|---|
| 52 | 2 NGC6334F 17:20:53.4 -35.47.01.3 17:36:34 1 x 02:33.9
|
|---|
| 53 | [2, 3] [0]
|
|---|
| 54 | 3 NGC6334Fe 17:23:23.4 -35.47.01.3 17:39:40 1 x 02:33.9
|
|---|
| 55 | [2, 3] [0]
|
|---|
| 56 | 4 NGC6334Fw 17:18:23.4 -35.47.01.3 17:42:11 1 x 01:22.9
|
|---|
| 57 | [4, 5] [0]
|
|---|
| 58 |
|
|---|
| 59 | Table contains 5 integration(s) in 5 scan(s).
|
|---|
| 60 |
|
|---|
| 61 |
|
|---|
| 62 | FreqID Frame RefFreq(Hz) RefPix Increment(Hz)
|
|---|
| 63 | --------------------------------------------------------------------------------
|
|---|
| 64 | 0 TOPO 2.396822e+10 2048 3906.25
|
|---|
| 65 | 1 TOPO 2.4143763e+10 2048 3906.25
|
|---|
| 66 | 2 TOPO 2.3968219e+10 2048 3906.25
|
|---|
| 67 | 3 TOPO 2.4143763e+10 2048 3906.25
|
|---|
| 68 | 4 TOPO 2.3968219e+10 2048 3906.25
|
|---|
| 69 | 5 TOPO 2.4143762e+10 2048 3906.25
|
|---|
| 70 | --------------------------------------------------------------------------------
|
|---|
| 71 |
|
|---|
| 72 |
|
|---|
| 73 | **BAD
|
|---|
| 74 |
|
|---|
| 75 | ASAP> scans = scantable('2005-03-09_172811_t002.rpf')
|
|---|
| 76 | Importing data...
|
|---|
| 77 | Auto averaging integrations...
|
|---|
| 78 | Weighting type selected : None
|
|---|
| 79 |
|
|---|
| 80 | ASAP> scans.summary(verbose=True)
|
|---|
| 81 |
|
|---|
| 82 | --------------------------------------------------------------------------------
|
|---|
| 83 | Scan Table Summary
|
|---|
| 84 | --------------------------------------------------------------------------------
|
|---|
| 85 | Beams: 1
|
|---|
| 86 | IFs: 2
|
|---|
| 87 | Polarisations: 1
|
|---|
| 88 | Channels: 4096
|
|---|
| 89 |
|
|---|
| 90 | Observer: Auto_Spec
|
|---|
| 91 | Obs Date: 2005/03/09/17:29:22
|
|---|
| 92 | Project: NGC6334Fw
|
|---|
| 93 | Obs. Type:
|
|---|
| 94 | Antenna Name: DSS-43
|
|---|
| 95 | Flux Unit: K
|
|---|
| 96 | Rest Freqs: [2.39639e+10] [Hz]
|
|---|
| 97 | Abcissa: Channel
|
|---|
| 98 | Cursor: Beam[0] IF[0] Pol[0]
|
|---|
| 99 |
|
|---|
| 100 | Scan Source Position (J2000) Time Integration
|
|---|
| 101 | FreqIDs
|
|---|
| 102 | --------------------------------------------------------------------------------
|
|---|
| 103 | 0 NGC6334Fw 17:47:49.9 +16.00.27.6 17:30:37 1 x 02:31.0
|
|---|
| 104 | [0, 1] [0]
|
|---|
| 105 | 1 NGC6334F 17:03:08.2 -44.33.03.4 17:33:43 1 x 02:33.9
|
|---|
| 106 | [0, 1] [0]
|
|---|
| 107 | 2 NGC6334F 03:46:52.4 +56.43.05.6 17:36:34 1 x 02:33.9
|
|---|
| 108 | [2, 3] [0]
|
|---|
| 109 | 3 NGC6334Fe 17:47:49.9 +16.00.27.6 17:39:40 1 x 02:33.9
|
|---|
| 110 | [2, 3] [0]
|
|---|
| 111 | 4 NGC6334Fw 17:03:08.2 -44.33.03.4 17:42:11 1 x 01:22.9
|
|---|
| 112 | [4, 5] [0]
|
|---|
| 113 |
|
|---|
| 114 | Table contains 5 integration(s) in 5 scan(s).
|
|---|
| 115 |
|
|---|
| 116 |
|
|---|
| 117 | FreqID Frame RefFreq(Hz) RefPix Increment(Hz)
|
|---|
| 118 | --------------------------------------------------------------------------------
|
|---|
| 119 | 0 TOPO 2.396822e+10 2048 3906.25
|
|---|
| 120 | 1 TOPO 2.4143763e+10 2048 3906.25
|
|---|
| 121 | 2 TOPO 2.3968219e+10 2048 3906.25
|
|---|
| 122 | 3 TOPO 2.4143763e+10 2048 3906.25
|
|---|
| 123 | 4 TOPO 2.3968219e+10 2048 3906.25
|
|---|
| 124 | 5 TOPO 2.4143762e+10 2048 3906.25
|
|---|
| 125 | --------------------------------------------------------------------------------
|
|---|
| 126 |
|
|---|
| 127 |
|
|---|
| 128 | Seems to be a good old rpfits (fortran) memory problem.
|
|---|
| 129 |
|
|---|
| 130 | repeatedly opening an rpfits file with the read buggers up the refbeam field.
|
|---|
| 131 | as the summary function looks at the refbeam to determine the positions it
|
|---|
| 132 | fails as the reader returns a "random" value.
|
|---|
| 133 |
|
|---|
| 134 | I have a chat to Mark.
|
|---|
| 135 |
|
|---|
| 136 | Malte.
|
|---|
| 137 |
|
|---|
| 138 | Mark,
|
|---|
| 139 |
|
|---|
| 140 | The attached data has only one beam.
|
|---|
| 141 | This only occurs on this data set.
|
|---|
| 142 | If I convert it to sdfits the problem is solved as you don't write this column
|
|---|
| 143 | out unless it is MX data.
|
|---|
| 144 | For now I have put a kludge in, which checks nBeam >1 and only then gets the
|
|---|
| 145 | refbeam from the data.
|
|---|
