----------------------------------------------------------------------------

REDUCING MOPRA OTF MAPS WITH LIVEDATA/GRIDZILLA

1. Make sure you are running the WEEKLY aips++ (check ~/.login.packages).

% which livedata
/nfs/aips++/weekly/sun4sol_gnu/bin/livedata

2. Insert the following lines into your ~/.Xdefaults if necessary:

! Karma resources.
*displayAxisLabelsToggle.state:         True
*OverlayMenu.createAtInit:              True

3. Reload with "xrdb ~/.Xdefaults".

4. Start livedata on a fast machine, like athena.
   (http://www.atnf.csiro.au/internal/computing/hostinfo/epping.html).

% livedata

5. Set "read directory" to where the RPFITS files are, "write directory" to 
where the output SDFITS files should go.

6. Activate the buttons "bandpass calibration", "monitor output", and 
"write data".  Activate "statistics" if desired.

7. Under "bandpass calibration" select the following:
    - spectral smoothing: none
    - method: referenced (*this is the important change*)
    - estimator: mean
    - cross-beam: NOT selected
    - spectral baseline fit: 1 or 2 (use 1 if your line is broad)
    - velocity frame: your preference (LSRK or BARY)
    - preserve continuum: NOT selected
    - rescale frequency axis: selected
The bottom "Validity checking" section is not relevant and ignored for this
method.

8. Under "data monitor" select the following:
    - for dual pol: IF(s) 1,1 or BOTH,BOTH; Pol A,B
    - for dual freq: IF(s) 1,2; Pol A
    - set the channel range for display purposes if desired

9. Select the desired RPFITS files and "Queue selection".

10. Hit the "Start" button at the top of the LiveData window. The baselined
spectra should appear on the Karma windows.  The Tsys should appear on the
multibeam stats plotter if you have activated the "statistics" button (it 
will fluctuate if the 2 channels have very different values).  

11. To view the individual spectra on the Karma windows, select in the
"View" window "Profile Mode: Line, Profile Axis: X".  You can "Pause" the
process in order to view spectra before they scroll off the pane.  By
default 100 spectra are shown.

12. In another window (if available), start gridzilla (in a directory with 
an existing gridzilla.lastexit file if possible).

13. Make the datacube:
    - Frequency or velocity range: give values encompassing emission
    - Rest frequency in MHz, if using velocity range
    - IF, Pol: select IF and pol(s) of interest (e.g. IF 1, Pol A&B)
    - Spectral: selected, continuum and baseline unselected
    - Search path: include directory where SDFITS files reside
    - Image geometry: SIN proj, give ref point if tiling many maps
    - Image centre & extent: autosize OK, pixel scale <= 0.2 arcmin
    - Gridding parameters: MEAN statistic, beam wtg 1, beam FWHM 0.5 arcmin, 
smoothing kernel top-hat FWHM 0.5', cutoff radius 0.25'.
    - Output FITS file basename: your choice
    - IEEE floating format

14. Hit "GO" button on lower right.  Output log appears on glish/tk window.

15. Two images are produced, the actual spectral line data cube and a
"spectracounts" image that indicates how many spectra are used at each grid
point.  The units of the cube are given as Jy/beam but are actually antenna 
temperature in K.