exit
#########################################################################
#User may edit the input variable for OPTION1
#########################################################################
BOTH:
echo "for imaging and plotting with calibrated data of" $InputCaledDATA
echo 'source =' $source
set OK=1
#continuum image
CIMAGE:
set imap = ${source}.imap
set beam = ${source}.beam
set icmp = ${source}.icmp
set icln = ${source}.icln
echo $imap $beam $icmp
#setup image size and cell size
set ims=512
set cel=0.2
set startLineCh = 1
set totalNumbCh = 3072
echo $imap $beam $icmp
\rm -r $imap $beam $icmp $icln
invert vis=$InputCaledDATA map=${source}.imap beam=${source}.beam \
    imsize=$ims cell=$cel sup=0 line=channel,1,$startLineCh,$totalNumbCh \
    select="source("$source")" options=nochk
clean  map=${source}.imap beam=${source}.beam out=${source}.icmp \
       gain=0.1 cutoff=0 niters=2500 region='quarter(1,1)'
restor model=${source}.icmp beam=${source}.beam map=${source}.imap out=${source}.icln
cgdisp in=${source}.icln type=pix region='quarter(1,1)' options=beambr,trlab,3val \
        labtyp=arcsec,arcsec xybin=1,1 device=/xs range=0,0,lin,1
if ($OK == 1 ) then
echo "Do you want to plot the spectrum?  Hit the enter key if yes."
$<
goto PLOTSPC
endif
exit
PLOTSPC:
smauvspec vis=$InputCaledDATA select="source("$source")" options=nobase,avall \
	interval=100000 axis=chan,both device=/xs nxy=1,1
if ($OK == 1 ) then
echo "Do you want to make image cube and momemt0?  Hit the enter key if yes."
$<
goto LCUBE 
endif
exit
LCUBE:
\rm -r $InputCaledDATA.lin
uvlin vis=$InputCaledDATA out=$InputCaledDATA.lin mode=line options=nowin \
	order=1 chans=$lchn1,$lchn2,$rchn1,$rchn2
$<
set imap = ${source}.l.imap
set beam = ${source}.l.beam
set icmp = ${source}.l.icmp
set icln = ${source}.l.icln
#
#setup image size and cell size
set ims=512
set cel=0.2
set startLineCh = $startCH 
set totalNumbCh = $endCH-$startCH 
invert vis=$InputCaledDATA.lin map=$imap beam=$beam out=$icmp \
	imsize=$ims cell=$cel robust=0 options=nochk \
        line=chan,$totalNumbCh,$startCH,1
exit