##################################################
## Recipe for aips++ for the calibration and mapping
## of an HI spectral line observation.
##
## This example uses the test dataset AH582.960603.LINE,
## which is a VLA D-array observation of the well known
## merging galaxy pair NGC 4038/9 "The Antennae" (see
## Hibbard et al., 2001, AJ, 122, 2969). 
##
## We assume that you have downloaded the test dataset 
## AH582.960603.LINE and put it in your working directory,
## and started up aips++
## Some details of observations:
## 127 channel datacube. 
## A single IF
## Two calibrators: 1331+305 flux calibrator FIELD_ID=1, 
##                  1130-148 phase calibrator FIELD_ID=2
## Source is in FIELD_ID=3
## The names and FIELD_IDs are printed out to the logger when 
## the ms is created, and can be reviewed using the ms.summary tool.
##
## JHIBBARD 1/18/02
## Last tested on aips++ v1.7 build 355 on 2002MAY1
##################################################

#########################################
## User Specified Parameters
## Note: Here I just put parameters that are used several places
##       in this script. Parameters that are used just for e.g.
##       calibration or mapping are put at the top of those
##       separate scripts. 

FITSFILENAME:='/export/data_1/jhibbard/DATA/n4038/AH582/AH582.960603.LINE';

BASENAME		:= 'n4038test';
PHASECAL		:= '1130-148';
PHASECAL_FIELDID	:= 2;
SOURCE_FIELDID		:= 3;

#########################################
## Data Reduction Control Parameters
## These specify which parts of calibration/mapping will
## be turned on, for times when I restart in the middle.
##

READ_DATA := 1;
DO_FLAG	  := 1;
DO_CAL    := 1;
DO_SELFCAL:= 0;      # this script does not exist yet
MAP_CAL	  := 0;
MAP_CH0	  := 0;
MAP_CUBE  := 0;

#########################################
##
## Get data into ms
##
######################################### 

dowait:=T;

if(READ_DATA){
  include 'ms.g'
  myms:=fitstoms(fitsfile=FITSFILENAME, msfile=spaste(BASENAME,'.ms'));
  myms.done();
}

## flagging data:
if(DO_FLAG){
  include 'specline_flag.g'
}

## Gain and Bandpass Calibration:
if(DO_CAL){
  include 'specline_calib.g'
}

## Map Phase calibrator to make sure everything is ok:
if(MAP_CAL){
  include 'specline_mapcal.g'
}

## Map "channel 0":
if(MAP_CH0){
  include 'specline_mapch0.g'
}

## Use channel 0 to get better phase corrections
## (ie, selfcal) 
if(DO_SELFCAL){
  include 'specline_selfcal.g'
}

## Make a continuum subtracted spectral line datacube
if(MAP_CUBE){
  include 'specline_mapcube.g';
  include 'specline_subcont.g';
}

## THINGS TO DO:
## Clean the continuum subtracted line cube
##   include 'specline_cleancube.g';
## Make better moment maps
##  include 'specline_analysis.g';
##   edit specline_cleancube.g to use mask
## Make some line profile plots and save spectra
## Make KNTR channel maps