#! /bin/csh -f
#
#SWARM test real data (for Miriad version SWARM1.0.2 and higher)
#JHZ: 2014-04-16
#JHZ: 2014-11-14
#
set dt    = 140411 
set hour =  13:48:03
set sb    = usb
set nn    = 1
set rxid  = 0 
set rx    = rx${rxid}
set edge  = 1010 
set fname = "$dt""_""$rx".$sb
echo $fname
set fcal  = venus 
set gcal  = MWC-Good 
set bcal  = 3c454.3 
set font  = MWC-Good 
set refant  = 2
set avetime = 5
set source  = $font
set DATPATH = /data/engineering/mir_data/
 goto SMALOD         # load the data
#goto INSPECT        # inspect the data
#goto SPLIT          # split the spectral chunks 
#goto SMAFIX         # apply Tsys
#goto BPASS          # do bandpass
#goto SMAFLUX        # flux scale
#goto UVLIN          # separate line and contimuum
#goto SELFCAL        # gain cali
#goto SMASPC         # plot spec
#goto MAP            # do images
#goto DISP           # display  
#goto ZAP            # clean 
SMALOD:
echo $fname
\rm -fr $fname
smalod in=${DATPATH}/$dt"_"$hour out=$dt rxif=$rxid \
   options=circular,cdebug sideband=$nn
$<
INSPECT:
uvindex  vis=$fname
$<
SPLIT:
\rm -rf $fname.if1
uvaver  vis=$fname  out=$fname.if1 select='window(2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25)'
$<
\rm -rf $fname.sw1
uvaver  vis=$fname  out=$fname.sw1 select='window(50)'
$<
\rm -rf $fname.sw2
uvaver  vis=$fname  out=$fname.sw2 select='window(51)'
$<
\rm -rf $fname.if2
uvaver  vis=$fname  out=$fname.if2 select='window(26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49)'
$<
SMAFIX:
\rm -r $fname.sw2.tsys
smafix vis=$fname.sw2 out=$fname.sw2.tsys \
        device=/xs xaxis=time \
        yaxis=systemp nxy=2,4   \
        options=tsyscorr 
$<
BPASS:
smamfcal vis=$fname.sw2.tsys select='source('$bcal')' edge=$edge,$edge \
     weight=2 refant=5 interval=1000000   
$<
SMAFLUX:
\rm -r $fname.sw2.tsys.bp
uvaver vis=$fname.sw2.tsys out=$fname.sw2.tsys.bp
smaflux vis=$fname.sw2.tsys.bp  \
	select='source('$fcal')' mirhome=$MIR
$<
\rm -r $fname.sw2.tsys.bp.a
uvaver vis=$fname.sw2.tsys.bp out=$fname.sw2.tsys.bp.a
$<
UVLIN:
\rm -r $font.ch0 $font.line
uvlin vis=$fname.sw2.tsys.bp.a chans=1000,8000,9500,15000 order=1 \
       mode=chan0 select='source('$font')' \
       out=$font.ch0
uvlin vis=$fname.sw2.tsys.bp.a chans=1000,8000,9500,15000 order=1 \
       mode=line select='source('$font')' \
       out=$font.line
$<
SELFCAL:
selfcal vis=$font.ch0 interval=0.1 options=phase select='source('$font')' 
$<
gpcopy vis=$font.ch0 out=$font.line mode=copy
\rm -r $font.ch0.p $font.line.p
uvaver vis=$font.ch0 out=$font.ch0.p
uvaver vis=$font.line out=$font.line.p
$<
SMASPC:
smauvspec vis=$font.line.p \
        select='source('$font')'  \
        interval=10000  options=nobase,avall \
        axis=chan,both device=/xs nxy=1,1 hann=13
$<
smauvspec vis=$font.line.p \
        select='source('$font')'  \
        interval=10000  options=nobase,avall \
        axis=freq,both device=/xs nxy=1,1 hann=13 
$<
puthd in=$font.line.p/restfreq value=231.9009
smauvspec vis=$font.line.p \
        select='source('$font')'  \
        interval=10000  options=nobase,avall \
        axis=vel,both device=/xs nxy=1,1 hann=13
$<
MAP:
\rm -r mwc349.map mwc349.beam 
invert vis=$font.line.p map=mwc349.map \
     beam=mwc349.beam imsize=512,512 cell=0.2 sup=0 \
     line=vel,100,-100,5
$<
\rm -r mwc349.icmp
clean map=mwc349.map beam=mwc349.beam out=mwc349.icmp \
	gain=0.08 cutoff=0.02 niters=25000 \
	region=arcsec,'boxes(-0.5,-0.5,0.5,0.5)(1,100)'
$<
\rm -r mwc349.icln
restor model=mwc349.icmp beam=mwc349.beam map=mwc349.map out=mwc349.icln
$<
cgdisp in=mwc349.icln type=c region=arcsec,'boxes(-25,-25,25,25)(1,100)' \
    xybin=1,1 device=/xs nxy=7,7 options=beambr,trlab,3val \
        labtyp=arcsec,arcsec cols1=7 csize=0.5,0.5,0.5 \
        slev=a,5 levs1=-35,8,16,32,64,128,256
$<
DISP:
cgdisp in=mwc349.icln type=c region=arcsec,'boxes(-3,-3,3,3)(30,54)' \
    xybin=1,1 device=/xs nxy=5,5 options=beambr,trlab,3val \
        labtyp=arcsec,arcsec cols1=7 csize=0.5,0.5,0.5 \
	slev=a,40 levs1=-35,8,10,13,17,22,28,35
exit
ZAP:
\rm -r $fname.sw2.tsys.bp $fname.sw2.tsys $fname.sw2
\rm -r mwc349.icmp mwc349.map
Update: JHZ 2014-04-21