#! /bin/csh -f
#
#jhz 2015-07-01  updated for Miriad SWARM3.0.2
#jhz 2014-04-30 released for Miriad SWARM2.0.0
#SWARM test real data
#corrected data from /sma/SMAusers/taco/SWARMTest24apr14/
#for target: sgrb2n
#
set dt    = 140424 
set hour =  12:18:22 
set sb    = usb
set nn    = 1
set rxid  = 0
set rx    = rx${rxid}
set edge  = 750 
set fname = "$dt""_""$rx".$sb
echo $fname
set fcal  = neptune 
set gcal  = 1924-292 
set bcal  = 3c454.3 
set font  = sgrb2n
set refant  = 2
set avetime = 5
set source  = $font
#set DATPATH = /data/engineering/mir_data/
set DATPATH = /sma/reduction/smamiriad/data/
goto SMALOD         # load the data
#goto INSPECT        # inspect the data
#goto SPLIT          # process standard
#goto SMAFIX         # apply Tsys
#goto QVACK          # elminate first interval after switch source
#goto UVFLAG         # flag select visibilities
#goto BPASS          # do bandpass
#goto BPPLT          # plot bandpass solutions
#goto SMAFLUX        # flux scale
#goto UVLIN          # separate line and contimuum
#goto SMABLFLAG      # flag bad data points
#goto SMAMFCAL       # calculate complex gain from a gain calibrator
#goto SMAGPPLT       # plot gain solutions
#goto SELFCAL        # gain cali
#goto SMASPC         # plot spec
#goto CMAP           # do images
#goto CDISP          # display
#goto LMAP           # do images
#goto LDISP          # display  
#goto ZAP            # clean 
SMALOD:
echo $fname
\rm -fr $fname
smalod in=${DATPATH}/$dt"_"$hour out=$dt rxif=$rxid \
   options=circular sideband=$nn
exit
if($nn == "1" || $nn == "0") then
   echo $nn
   echo "click enter for continue"
   $<
   goto SMALOD
   else
   echo "down!"
   exit
endif
$<
INSPECT:
uvindex  vis=$fname
$<
SPLIT:
#
## if1 = s1,s2,s3 ...,s24
#
\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.s50
uvaver  vis=$fname  out=$fname.s50 select='window(50)'
$<
\rm -rf $fname.s49
uvaver  vis=$fname  out=$fname.s49 select='window(51)'
$<
#
## if2 = s25,s26,s27, ..., s48
#
\:200
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.s49.tsys
smafix vis=$fname.s49 out=$fname.s49.tsys \
        device=/xs xaxis=time \
        yaxis=systemp nxy=2,4   \
        options=tsyscorr 
$<

QVACK:
qvack vis=$fname.s49.tsys mode=source interval=1
$<
UVFLAG:
uvflag vis=$fname.s49.tsys select='ant(7),time(16:15,18:00)' \
    flagval=flag
$<
BPASS:
smamfcal vis=$fname.s49.tsys select='source('$bcal')' edge=$edge,$edge \
     weight=2 refant=5 interval=1000000   
$<
BPPLT:
smagpplt vis=$fname.s49.tsys device=/xs \
        yaxis=amp \
        options=bandpass  \
        nxy=2,4
$<
smagpplt vis=$fname.s49.tsys device=/xs \
        yaxis=phase \
        options=bandpass  \
        nxy=2,4
$<
SMAFLUX:
\rm -r $fname.s49.tsys.bp
uvaver vis=$fname.s49.tsys out=$fname.s49.tsys.bp
$<
smaflux vis=$fname.s49.tsys.bp  \
	select='source('$fcal')' mirhome=$MIR
$<
\rm -r $fname.s49.tsys.bp.a
uvaver vis=$fname.s49.tsys.bp out=$fname.s49.tsys.bp.a
$<
UVLIN:
\rm -r $font.ch0 $font.line
uvlin vis=$fname.s49.tsys.bp.a chans=1000,8000,9500,15000 order=1 \
       mode=chan0 select='source('$font')' \
       out=$font.ch0
uvlin vis=$fname.s49.tsys.bp.a chans=1000,8000,9500,15000 order=1 \
       mode=line select='source('$font')' \
       out=$font.line
\rm -r $gcal.ch0
uvlin vis=$fname.s49.tsys.bp.a chans=1000,9499,9500,15000 order=1 \
       mode=chan0 select='source('$gcal')' \
       out=$gcal.ch0
$<
SMABLFLAG:
#
## flag bad data points
#
smablflag vis=$gcal.ch0 axis=time,ampl device=/xs
smablflag vis=$font.ch0 axis=time,ampl device=/xs
$<
SELFCAL:
selfcal vis=$gcal.ch0 interval=0.5 options=phase  
$<
gpcopy vis=$gcal.ch0 out=$font.ch0 mode=copy
gpcopy vis=$gcal.ch0 out=$font.line mode=copy
puthd in=$font.ch0/interval value=100
puthd in=$font.line/interval value=100
\rm -r $font.ch0.p $font.line.p $gcal.ch0.p
uvaver vis=$font.ch0 out=$font.ch0.p
uvaver vis=$font.line out=$font.line.p
uvaver vis=$gcal.ch0 out=$gcal.ch0.p
$<
SMAMFCAL:
smamfcal vis=$gcal.ch0.p interval=3 refant=1 options=nopassol \
 flux=3.55 select='source('$gcal')'
$<
SMAGPPLT:
smagpplt vis=$gcal.ch0.p  device=/xs \
        yaxis=amp \
        options=gain  \
        nxy=2,4
$<
$<
smagpplt vis=$gcal.ch0.p  device=/xs \
        yaxis=phase \
        options=gain \
        nxy=2,4
$<
puthd in=$font.ch0.p/interval value=1000 type=double
puthd in=$font.line.p/interval value=1000 type=double
gpcopy vis=$gcal.ch0.p out=$font.ch0.p mode=copy
gpcopy vis=$gcal.ch0.p out=$font.line.p mode=copy
$<
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  \
        axis=freq,both device=/xs nxy=1,1 hann=15 
$<
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
$<
CMAP:
\rm -r $font.map $font.beam
invert vis=$font.ch0.p map=$font.map \
     beam=$font.beam imsize=512,512 cell=0.3 sup=0 
$<
\rm -r $font.icmp
clean map=$font.map beam=$font.beam out=$font.icmp \
        gain=0.08 cutoff=0.01 niters=25000 \
        region=arcsec,'boxes(-5,-5,5,5)'
$<
\rm -r $font.icln
restor model=$font.icmp beam=$font.beam map=$font.map out=$font.icln
$<
CDISP:
cgdisp in=$font.icln type=c region=arcsec,'boxes(-25,-25,25,25)' \
    xybin=1,1 device=/xs nxy=1,1 options=beambr,trlab,3val \
        labtyp=hms,dms cols1=7 csize=0.5,0.5,0.5 \
        slev=a,0.2 levs1=-8,8,16,32,64,128,256
$<
cgdisp in=$font.icln type=c region=arcsec,'boxes(-15,-15,15,15)' \
    xybin=1,1 device=/xs nxy=1,1 options=beambr,trlab,3val \
        labtyp=arcsec,arcsec cols1=7 csize=0.5,0.5,0.5 \
        slev=a,0.2 levs1=-8,8,16,32,64,128,256

$<
LMAP:
\rm -r $font.map $font.beam 
invert vis=$font.line.p map=$font.map \
     beam=$font.beam imsize=512,512 cell=0.3 sup=0 \
     line=vel,100,-100,5
$<
\rm -r $font.icmp
clean map=$font.map beam=$font.beam out=$font.icmp \
	gain=0.08 cutoff=0.02 niters=25000 \
	region=arcsec,'boxes(-5,-5,5,5)(1,100)'
$<
\rm -r $font.icln
restor model=$font.icmp beam=$font.beam map=$font.map out=$font.icln
$<
LDISP:
cgdisp in=$font.icln type=c region=arcsec,'boxes(-25,-25,25,25)(1,100)' \
    xybin=1,1 device=/xs nxy=10,10 options=beambr,trlab,3val \
        labtyp=arcsec,arcsec cols1=7 csize=0.5,0.5,0.5 \
        slev=a,0.3 levs1=-4,4,8,16,32,64,128,192,256
exit
ZAP:
\rm -r $fname.s49.tsys.bp $fname.s49.tsys.bp.a $fname.s49.tsys $fname.s49 $gcal.ch0