Let’s assume you just got your data from the telescope; now you need to look at the data, but don’t
want to read the
SPECX manual or this one. Then follow these steps, depending on what you
want to do. You can learn about the intricacies later. Occasional margin notes, like the one
here(3), indicate sections where additional information can be found in the main part of this
SPECX is simple. All the necessary environment variables should have been setup correctly as
part of your Starlink login. You can start
SPECX from any directory in which you have write access but
for now it’s easier if you start it up in the directory that contains your data. In order to start
% sign represents the shell prompt for input; its actual form may be something Unix-y like
3:45|iiwi~|; in this case giving the time and a reminder of the name of the workstation you are
SPECX will next dump out part of a long introductory text. You might want to read this; if so use the
space bar to jump a page at a time. If not, use
Control-c to jump to the end. At the end, the system
responds with the
SPECXcomplains that it cannot find your data, this is likely if you started
SPECXin a directory that does not contain your data, you should exit
SPECX(using the command exit), and before restarting up
or wherever you happen to have stored your data.
SPECX will then look for data first in
your current directory but then in your data directory (as defined by the DATADIR
s-t-d (system default is xw for a Sun workstation such as IIWI)
s-p-d t (set to plot on terminal screen)
s-h-d (default is Postscript-landscape – shown as ps_l; but note alternatives)
The default state of
SPECX is to look for spectra taken with the
DAS; if that’s what you want, there is no
need to do anything to get it. However, if you have (old) data taken with the
AOSC, say, a number of
parameters have to be changed. See Section 3.3.1 for this and other additional information on setting
To read in a spectrum (GSD format)(3.4) type (e.g. spectrum 131):
See Section 3.4 for more information, especially if you want to look at map data.
Plot on screen:
DAS spectra do not look very pretty (cf. Figure 8) in most cases at this point until you
have done the step described in section 2.5.
SPECX under Unix also has a bad habit of leaving the
cursor on the plot window.
Change plot scales(3.4.1) with:
and answer the questions.
The x-axis default is velocity; if you want to change it use(3.16.1):
and complete the dialog. The basic options are channels, velocity and frequency scales.
To average(3.11.1) three spectra (131, 132 and 133) type:
r-g-d 132; ave
r-g-d 133; ave
Plot result, with
Incidentally, note that two or more commands may be placed on one line, separated by semicolons, as shown in this example.
As described elsewhere,
DAS spectra usually (3.10) consist of several overlapping (in frequency/velocity
space) sections (or ‘sub-bands’). Except for single sub-band (125-MHz) spectra it is necessary to
combine these sub-bands to produce a nice spectrum. This is done with the command das-merge;
This command removes the large data values at the ends of the subbands, and combines the subbands into one continuous spectum. It is sufficient to take the defaults for most purposes; however, see Section 3.10 for caveats etc. The final spectrum should be free of evidence of sub-band edges, except in pathological cases.
das-mergeoperation on individual spectra before averaging; the average can be formed first and das-merge applied to the result.
before attempting the das-merge command.
The (3.18) simplest way is to type
to send the last plot you made on the screen. The Unix version of
SPECX creates a plot file called
specx_pgplot.ps if you have asked for Postscript files to be made (these are the required type for
our normal printer). This plot file will be overwritten every time you issue a new plot
To send any other plot directly to the (laser) printer, first change plot device to hardcopy:
then open a new plotfile:
and either close the plotfile;
or open a new plotfile;
or reset the plot device:
all of which create the plotfile specx_pgplot.ps. To plot it on the default (local) printer use
$ lp specx_pgplot.ps
and if you like this plot enough you can save the Postscript file with
$ mv specx_pgplot.ps yourplotfilename.ps
Data may be written to files, which can contain many spectra, or saved temporarily in storage registers, each of which contains only one spectrum.
To write data to a file, first open one(3.9.1):
SPECX the file name etc4.
Each file can be addressed by its name, or a number. A file name has by default the filetype .sdf,
which will be added to whatever name you give. Several files may be open at once. To discover a file’s
number (you will need this) type:
Then, to write the current data (that seen with the n command) to (say) file 2 type:
The data are stored in the next available space.
To list the contents of a file type:
and answer the questions.
To save the current data in a temporary storage register(3.9.2), type, e.g.:
To read data from file you must have read (R), or read-write (RW), access. The default state of newly-created files is write-only. l-o-f will tell you the current access. To change the latter enter(3.9.1):
and answer the questions.
To read a spectrum (say, number 16, from file 2):
rea-s 2 16
This places the result in the x-register, from where it is immediately available for plotting or whatever you might have in mind.
To get data from a temporary storage area(3.9.2) (say, 3):
One usually finds that one’s spectral line lies atop some kind of baseline curvature; to display the line to better advantage one fits a linear, polynomial or other function to the line-free regions of the spectrum. It is possible to specify the “fit regions” either by typing in the numbers, or by using the cursor with the plot. The latter is where the interactive mode comes in handy. First type:
to remove a linear baseline from the current spectrum. Define two fitting regions by placing the cursor in the respective positions, and typing l, r, and a (left, right, accept, region) for each. The x-register will contain the result.
To remove a polynomial baseline, type(3.7.2)
and define a number of regions to which the baseline is to be fitted. If you use the cursor as before, exit the
plot screen with an ‘e’ (for ‘exit’). If you supply the fit ranges in non-interactive mode terminate the entries with
you will be prompted for the order of polynomial to be fitted. The model fit curve is in the x-register and the original
data in the y-register.6
Then to display the fit superposed on the original before subtraction, type, say
over 1 5
The overlay function expects two numbers – the line width and colour.
Generally you will be confident of the results of your fit, so you can skip the above step. To see the final baseline-subtracted result type
Map files(3.19) (‘datacubes’ to the initiated) consist of a potentially large number of related spectra in
one place under a common header. The spectra are taken on a regular spatial grid of points. This,
and a few other things, you must tell
SPECX before anything else. First create a map file
Give it a name (whatever name you give is automatically appended with the string _map.sdf), tell the system the grid spacing of the observations, and a few other things (answer the questions). Be sure to set the map size to be adequate for your needs, or else you will have to recreate it later. Only one map file may be open at one time.
To place the current spectrum in the map file, use the add-to-map command(3.19.2):
To overwrite an existing spectrum in the map, first use
to set map access (or not, as the case may be).
Several techniques are useful. First, typing(3.19.5)
will display a set of postage-stamp spectra on the page. You provide a number of parameters defining the plot.
to obtain a contour plot.
Or, better yet (provided you are working on a Sun workstation or something similar), type(3.19.7)
to obtain a gray-scale plot. This can be turned into various colour schemes, with a number of options. I suggest you experiment; using the interactive mode is a good way to do this – type ‘h’ within the plot to see the options available – or in non-interactive mode use the command
Finally, either using
s-p-d h and sending one of the above commands, or typing
see-map, will save
the image to a file (e.g. specx_pgplot.ps) suitable for output to the printer.
Any command may be aborted while it awaits input by typing
Data may be ‘binned’ (averaged) over n channels by typing
and smoothed using
Several spectra may be plotted on the same frame to good effect by offsetting vertically with respect to one another using
and then using the
One spectrum may be subtracted from another by using
or divided by another (channel-by-channel) by typing:
This just gets you started. By now you should have learned enough for simple operations. However,
like as not you will want more, such as how to build procedures for repetitive operations. So, you will
have to read the full introduction to
SPECX in this cookbook, or Rachael Padman’s full
SPECX manual. Also, experimentation is a good thing; the worst that can happen is a core dump (a
large file called
core will be created in your directory). You should delete this file if this
4The VMS version of
SPECX requires the number of spectra expected in addition; Unix allows the files to be expanded
SPECX prompts you for an ‘EOF’ (end-of-file); what the Unix-types don’t tell you is that this is a
Control-d; it is
Control-z under VMS.
SPECX uses an HP calculator-like stack; the x-register is on the top, and is the one which is always plotted.