summaryrefslogtreecommitdiffstats
path: root/ast/ast_tester/a20070718_00010_02_cube.ast
blob: f0059388a1321d3ef12477c45f6e4299f0e67ae0 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
 Begin FrameSet 	# Set of inter-related coordinate systems
#   Title = "3-d compound coordinate system" 	# Title of coordinate system
#   Naxes = 3 	# Number of coordinate axes
#   Domain = "SKY-DSBSPECTRUM" 	# Coordinate system domain
#   Epoch = 2007.54208230851 	# Julian epoch of observation
#   Lbl1 = "Right ascension" 	# Label for axis 1
#   Lbl2 = "Declination" 	# Label for axis 2
#   Lbl3 = "Radio velocity (USB)" 	# Label for axis 3
#   System = "Compound" 	# Coordinate system type
#   Uni1 = "hh:mm:ss.s" 	# Units for axis 1
#   Uni2 = "ddd:mm:ss" 	# Units for axis 2
#   Uni3 = "km/s" 	# Units for axis 3
#   Dir1 = 0 	# Plot axis 1 in reverse direction
 IsA Frame 	# Coordinate system description
    Nframe = 4 	# Number of Frames in FrameSet
#   Base = 1 	# Index of base Frame
    Currnt = 4 	# Index of current Frame
    Nnode = 5 	# Number of nodes in FrameSet
    Nod1 = 3 	# Frame 1 is associated with node 3
    Nod2 = 4 	# Frame 2 is associated with node 4
    Nod3 = 5 	# Frame 3 is associated with node 5
    Nod4 = 2 	# Frame 4 is associated with node 2
    Lnk2 = 1 	# Node 2 is derived from node 1
    Lnk3 = 1 	# Node 3 is derived from node 1
    Lnk4 = 1 	# Node 4 is derived from node 1
    Lnk5 = 1 	# Node 5 is derived from node 1
    Frm1 = 	# Frame number 1
       Begin Frame 	# Coordinate system description
          Title = "Data grid indices; first pixel at (1,1,1)" 	# Title of coordinate system
          Naxes = 3 	# Number of coordinate axes
          Domain = "GRID" 	# Coordinate system domain
#         Lbl1 = "Data grid index 1" 	# Label for axis 1
#         Lbl2 = "Data grid index 2" 	# Label for axis 2
#         Lbl3 = "Data grid index 3" 	# Label for axis 3
#         Uni1 = "pixel" 	# Units for axis 1
#         Uni2 = "pixel" 	# Units for axis 2
#         Uni3 = "pixel" 	# Units for axis 3
          Ax1 = 	# Axis number 1
             Begin Axis 	# Coordinate axis
                Label = "Data grid index 1" 	# Axis Label
                Symbol = "g1" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax2 = 	# Axis number 2
             Begin Axis 	# Coordinate axis
                Label = "Data grid index 2" 	# Axis Label
                Symbol = "g2" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax3 = 	# Axis number 3
             Begin Axis 	# Coordinate axis
                Label = "Data grid index 3" 	# Axis Label
                Symbol = "g3" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
       End Frame
    Frm2 = 	# Frame number 2
       Begin Frame 	# Coordinate system description
          Title = "Pixel coordinates; first pixel at (-8.5,-8.5,-2048.5)" 	# Title of coordinate system
          Naxes = 3 	# Number of coordinate axes
          Domain = "PIXEL" 	# Coordinate system domain
#         Lbl1 = "Pixel coordinate 1" 	# Label for axis 1
#         Lbl2 = "Pixel coordinate 2" 	# Label for axis 2
#         Lbl3 = "Pixel coordinate 3" 	# Label for axis 3
#         Uni1 = "pixel" 	# Units for axis 1
#         Uni2 = "pixel" 	# Units for axis 2
#         Uni3 = "pixel" 	# Units for axis 3
          Ax1 = 	# Axis number 1
             Begin Axis 	# Coordinate axis
                Label = "Pixel coordinate 1" 	# Axis Label
                Symbol = "p1" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax2 = 	# Axis number 2
             Begin Axis 	# Coordinate axis
                Label = "Pixel coordinate 2" 	# Axis Label
                Symbol = "p2" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax3 = 	# Axis number 3
             Begin Axis 	# Coordinate axis
                Label = "Pixel coordinate 3" 	# Axis Label
                Symbol = "p3" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
       End Frame
    Frm3 = 	# Frame number 3
       Begin Frame 	# Coordinate system description
          Title = "Axis coordinates; first pixel at (-8.5,-8.5,-2048.5)" 	# Title of coordinate system
          Naxes = 3 	# Number of coordinate axes
          Domain = "AXIS" 	# Coordinate system domain
#         Lbl1 = "Axis 1" 	# Label for axis 1
#         Lbl2 = "Axis 2" 	# Label for axis 2
#         Lbl3 = "Axis 3" 	# Label for axis 3
#         Uni1 = "pixel" 	# Units for axis 1
#         Uni2 = "pixel" 	# Units for axis 2
#         Uni3 = "pixel" 	# Units for axis 3
          Ax1 = 	# Axis number 1
             Begin Axis 	# Coordinate axis
                Label = "Axis 1" 	# Axis Label
                Symbol = "a1" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax2 = 	# Axis number 2
             Begin Axis 	# Coordinate axis
                Label = "Axis 2" 	# Axis Label
                Symbol = "a2" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
          Ax3 = 	# Axis number 3
             Begin Axis 	# Coordinate axis
                Label = "Axis 3" 	# Axis Label
                Symbol = "a3" 	# Axis symbol
                Unit = "pixel" 	# Axis units
                Format = "%3.1f" 	# Format specifier
             End Axis
       End Frame
    Frm4 = 	# Frame number 4
       Begin CmpFrame 	# Compound coordinate system description
#         Title = "3-d compound coordinate system" 	# Title of coordinate system
#         Naxes = 3 	# Number of coordinate axes
#         Domain = "SKY-DSBSPECTRUM" 	# Coordinate system domain
#         Lbl1 = "Right ascension" 	# Label for axis 1
#         Lbl2 = "Declination" 	# Label for axis 2
#         Lbl3 = "Radio velocity (USB)" 	# Label for axis 3
#         Uni1 = "hh:mm:ss.s" 	# Units for axis 1
#         Uni2 = "ddd:mm:ss" 	# Units for axis 2
#         Uni3 = "km/s" 	# Units for axis 3
#         Dir1 = 0 	# Plot axis 1 in reverse direction
       IsA Frame 	# Coordinate system description
          FrameA = 	# First component Frame
             Begin SkyFrame 	# Description of celestial coordinate system
                Naxes = 2 	# Number of coordinate axes
                Epoch = 2007.54208230851 	# Julian epoch of observation
                System = "FK5" 	# Coordinate system type
                ObsLat = 0.346026050148997 	# Observers geodetic latitude (rads)
                ObsLon = -2.71363306946838 	# Observers geodetic longitude (rads)
                Dut1 = -0.162963400559962 	# UT1-UTC in seconds
                Ax1 = 	# Axis number 1
                   Begin SkyAxis 	# Celestial coordinate axis
                   End SkyAxis
                Ax2 = 	# Axis number 2
                   Begin SkyAxis 	# Celestial coordinate axis
                   End SkyAxis
             IsA Frame 	# Coordinate system description
                SRef1 = 4.90555239050187 	# Ref. pos. RA 18:44:16.2
                SRef2 = -0.0698132003068136 	# Ref. pos. Dec -4:00:00
             End SkyFrame
          FrameB = 	# Second component Frame
             Begin DSBSpecFrame 	# Dual sideband spectral axis
                Naxes = 1 	# Number of coordinate axes
                Epoch = 2007.54208230851 	# Julian epoch of observation
                System = "VRAD" 	# Coordinate system type
                ObsLat = 0.346026069000144 	# Observers geodetic latitude (rads)
                ObsLon = -2.71363307300091 	# Observers geodetic longitude (rads)
                Dut1 = -0.162963400559962 	# UT1-UTC in seconds
                Ax1 = 	# Axis number 1
                   Begin Axis 	# Coordinate axis
                   End Axis
             IsA Frame 	# Coordinate system description
                SoR = "LSRK" 	# Standard of rest
                RefRA = 4.90553899950038 	# Reference RA (rads, FK5 J2000)
                RefDec = -0.0698131700797732 	# Reference Dec (rads, FK5 J2000)
                RstFrq = 372672509000 	# Rest frequency (Hz)
                SrcVel = 87012.623732646 	# Source velocity (m/s)
                SrcVRF = "LSRK" 	# Source velocity rest frame
                UFreq = "GHz" 	# Preferred units for frequency
             IsA SpecFrame 	# Description of spectral coordinate system
                DSBCen = 372576971505.836 	# Central frequency (Hz topo)
                IF = -5251168918.31241 	# Intermediate frequency (Hz)
                SideBn = "USB" 	# Represents upper sideband
             End DSBSpecFrame
       End CmpFrame
    Map2 = 	# Mapping between nodes 1 and 2
       Begin CmpMap 	# Compound Mapping
          Nin = 3 	# Number of input coordinates
          IsSimp = 1 	# Mapping has been simplified
       IsA Mapping 	# Mapping between coordinate systems
          MapA = 	# First component Mapping
             Begin CmpMap 	# Compound Mapping
                Nin = 3 	# Number of input coordinates
                Invert = 0 	# Mapping not inverted
             IsA Mapping 	# Mapping between coordinate systems
                Series = 0 	# Component Mappings applied in parallel
                InvA = 1 	# First Mapping used in inverse direction
                MapA = 	# First component Mapping
                   Begin WinMap 	# Map one window on to another
                      Nin = 2 	# Number of input coordinates
                      Invert = 0 	# Mapping not inverted
                   IsA Mapping 	# Mapping between coordinate systems
                      Sft1 = 10 	# Shift for axis 1
                      Sft2 = 10.4938775599003 	# Shift for axis 2
                   End WinMap
                MapB = 	# Second component Mapping
                   Begin CmpMap 	# Compound Mapping
                      Nin = 1 	# Number of input coordinates
                      IsSimp = 1 	# Mapping has been simplified
                   IsA Mapping 	# Mapping between coordinate systems
                      InvA = 1 	# First Mapping used in inverse direction
                      MapA = 	# First component Mapping
                         Begin WinMap 	# Map one window on to another
                            Nin = 1 	# Number of input coordinates
                            Invert = 0 	# Mapping not inverted
                         IsA Mapping 	# Mapping between coordinate systems
                            Sft1 = 6105939.046167 	# Shift for axis 1
                            Scl1 = -1.63834485003565e-05 	# Scale factor for axis 1
                         End WinMap
                      MapB = 	# Second component Mapping
                         Begin SpecMap 	# Conversion between spectral coordinate systems
                            Nin = 1 	# Number of input coordinates
                            Invert = 0 	# Mapping not inverted
                         IsA Mapping 	# Mapping between coordinate systems
                            Nspec = 2 	# Number of conversion steps
                            Spec1 = "FRTOVL" 	# Convert frequency to rel. velocity
                            Spec1a = 372672509000 	# Rest frequency (Hz)
                            Spec2 = "VLTOVR" 	# Convert relativistic to radio velocity
                         End SpecMap
                   End CmpMap
             End CmpMap
          MapB = 	# Second component Mapping
             Begin CmpMap 	# Compound Mapping
                Nin = 3 	# Number of input coordinates
             IsA Mapping 	# Mapping between coordinate systems
                Series = 0 	# Component Mappings applied in parallel
                MapA = 	# First component Mapping
                   Begin CmpMap 	# Compound Mapping
                      Nin = 2 	# Number of input coordinates
                      Invert = 0 	# Mapping not inverted
                   IsA Mapping 	# Mapping between coordinate systems
                      MapA = 	# First component Mapping
                         Begin MatrixMap 	# Matrix transformation
                            Nin = 2 	# Number of input coordinates
                            Invert = 1 	# Mapping inverted
                         IsA Mapping 	# Mapping between coordinate systems
                            M0 = -2.90888208665722e-05 	# Forward matrix value
                            M1 = 2.90888208665722e-05 	# Forward matrix value
                            Form = "Diagonal" 	# Matrix storage form
                         End MatrixMap
                      MapB = 	# Second component Mapping
                         Begin CmpMap 	# Compound Mapping
                            Nin = 2 	# Number of input coordinates
                         IsA Mapping 	# Mapping between coordinate systems
                            InvA = 1 	# First Mapping used in inverse direction
                            MapA = 	# First component Mapping
                               Begin WcsMap 	# FITS-WCS sky projection
                                  Nin = 2 	# Number of input coordinates
                                  Invert = 0 	# Mapping not inverted
                               IsA Mapping 	# Mapping between coordinate systems
                                  Type = "TAN" 	# Gnomonic projection
                               End WcsMap
                            MapB = 	# Second component Mapping
                               Begin CmpMap 	# Compound Mapping
                                  Nin = 2 	# Number of input coordinates
                               IsA Mapping 	# Mapping between coordinate systems
                                  InvA = 1 	# First Mapping used in inverse direction
                                  MapA = 	# First component Mapping
                                     Begin SphMap 	# Cartesian to Spherical mapping
                                        Nin = 3 	# Number of input coordinates
                                        Nout = 2 	# Number of output coordinates
                                        Invert = 0 	# Mapping not inverted
                                     IsA Mapping 	# Mapping between coordinate systems
                                        UntRd = 1 	# All input vectors have unit length
                                        PlrLg = 0 	# Polar longitude (rad.s)
                                     End SphMap
                                  MapB = 	# Second component Mapping
                                     Begin CmpMap 	# Compound Mapping
                                        Nin = 3 	# Number of input coordinates
                                        Nout = 2 	# Number of output coordinates
                                     IsA Mapping 	# Mapping between coordinate systems
                                        MapA = 	# First component Mapping
                                           Begin MatrixMap 	# Matrix transformation
                                              Nin = 3 	# Number of input coordinates
                                              Invert = 1 	# Mapping inverted
                                           IsA Mapping 	# Mapping between coordinate systems
                                              M0 = -0.0133908234539359 	# Forward matrix value
                                              M1 = 0.981401710741326 	# Forward matrix value
                                              M2 = 0.191497697117288 	# Forward matrix value
                                              M3 = 0.068459122667767 	# Forward matrix value
                                              M4 = 0.191965314976427 	# Forward matrix value
                                              M5 = -0.979011065499037 	# Forward matrix value
                                              M6 = -0.997564050259824 	# Forward matrix value
                                              M7 = 0 	# Forward matrix value
                                              M8 = -0.0697564737441253 	# Forward matrix value
                                              Form = "Full" 	# Matrix storage form
                                           End MatrixMap
                                        MapB = 	# Second component Mapping
                                           Begin SphMap 	# Cartesian to Spherical mapping
                                              Nin = 3 	# Number of input coordinates
                                              Nout = 2 	# Number of output coordinates
                                              Invert = 1 	# Mapping inverted
                                           IsA Mapping 	# Mapping between coordinate systems
                                              UntRd = 1 	# All input vectors have unit length
                                              PlrLg = 4.90555329534916 	# Polar longitude (rad.s)
                                           End SphMap
                                     End CmpMap
                               End CmpMap
                         End CmpMap
                   End CmpMap
                MapB = 	# Second component Mapping
                   Begin ZoomMap 	# Zoom about the origin
                      Nin = 1 	# Number of input coordinates
                   IsA Mapping 	# Mapping between coordinate systems
                      Zoom = 0.001 	# Zoom factor
                   End ZoomMap
             End CmpMap
       End CmpMap
    Map3 = 	# Mapping between nodes 1 and 3
       Begin UnitMap 	# Unit (null) Mapping
          Nin = 3 	# Number of input coordinates
          IsSimp = 1 	# Mapping has been simplified
       IsA Mapping 	# Mapping between coordinate systems
       End UnitMap
    Map4 = 	# Mapping between nodes 1 and 4
       Begin WinMap 	# Map one window on to another
          Nin = 3 	# Number of input coordinates
          IsSimp = 1 	# Mapping has been simplified
       IsA Mapping 	# Mapping between coordinate systems
          Sft1 = -9.5 	# Shift for axis 1
          Sft2 = -9.5 	# Shift for axis 2
          Sft3 = -2049.5 	# Shift for axis 3
       End WinMap
    Map5 = 	# Mapping between nodes 1 and 5
       Begin WinMap 	# Map one window on to another
          Nin = 3 	# Number of input coordinates
          IsSimp = 1 	# Mapping has been simplified
       IsA Mapping 	# Mapping between coordinate systems
          Sft1 = -9.5 	# Shift for axis 1
          Sft2 = -9.5 	# Shift for axis 2
          Sft3 = -2049.5 	# Shift for axis 3
       End WinMap
 End FrameSet