#!/bin/csh -f
#jhz 2016-05-06 continuum imaging              
echo "                                                                      "
echo "######################################################################"
echo "specially for the project 2015B-S044 - a narrow line system           "
echo "######################################################################"
set maxloop = 1
goto CONTMAP 
set uvdat=151226_rx0.lsb.ch0  
blflag vis=$uvdat axis=time,amp device=/xs options=nobase 
$<
set uvdat=151226_rx0.usb.ch0
blflag vis=$uvdat axis=time,amp device=/xs options=nobase
$<
CONTMAP:
set uvdat=151226_rx0.usb.ch0,151226_rx0.lsb.ch0
\rm -r target.map target.beam target.icmp
echo 'uvdat=' $uvdat
set opt = systemp
echo $opt
set i=0
while ($i != $maxloop)
invert vis=$uvdat map=target.map \
        beam=target.beam   \
        imsize=512 cell=.25 sup=0 options=$opt
cgdisp in=target.map type=p region=arcsec,'boxes(-40,-40,40,40)' \
        xybin=1 device=/xs nxy=1 labtyp=arcsec,arcsec \
        range=0,0,lin,2 \
        options=beambr,wedge,trlab csize=1,1,0.8 \
        slev=a,0.0025 
#clean map
clean map=target.map beam=target.beam \
        out=target.icmp gain=0.05 \
        cutoff=0.0015 niters=20000 region=arcsec,'box(-3,-3,3,3)'
\rm -r target.icln
#restore the image with a gaussian beam
restor model=target.icmp beam=target.beam map=target.map \
        out=target.icln
cgdisp in=target.icln type=p region=arcsec,'boxes(-40,-40,40,40)' \
        xybin=1 device=/xs nxy=1 labtyp=arcsec,arcsec \
        range=0,0,lin,2 \
        options=beambr,wedge,trlab csize=1,1,0.8 \
        slev=a,0.0025 
#do statistics in an area off the source
imstat in=target.icln region='boxes(30,30,60,60)' \
        axes=ra,dec
#do statistics in an area on the source
imstat in=target.icln region=quarter \
        axes=ra,dec
#make a contore map
#contours
cgdisp in=target.icln type=c region=arcsec,'boxes(-20,-20,20,20)' \
        xybin=1 device=/xs nxy=1 labtyp=arcsec,arcsec \
        range=0,0,lin,2 \
        options=beambr,wedge,trlab csize=1,1,0.8 \
        slev=a,0.0013 \
        levs1=-3,3,4,6,9,13,18,24,31,40,50 
echo "plot color"
$<
cgdisp in=target.icln type=pixel region=arcsec,'boxes(-20,-20,20,20)' \
        xybin=2 device=/xs \
        nxy=1 labtyp=hms,dms options=beambr,wedge,trlab \
        labtyp=arcsec,arcsec cols1=7 \
        range=0,0,lin,2 csize=1,1,0.8
if($i < $maxloop ) then
selfcal vis=151226_rx0.lsb.ch0 refant=1 options=pha \
interval=1.5 model=target.icmp clip=0.006
gpcopy vis=151226_rx0.lsb.ch0 out=151226.SiO54.uv mode=copy
selfcal vis=151226_rx0.usb.ch0 refant=1 options=pha \
interval=1.5 model=target.icmp clip=0.006
gpcopy vis=151226_rx0.usb.ch0 out=151226.CO21.uv mode=copy
endif
#
DOSTATS:
#do statistics in an area off the source
imstat in=target.icln region='boxes(30,30,60,60)' \
        axes=ra,dec
#do statistics in an area on the source
imstat in=target.icln region=quarter \
        axes=ra,dec
        @ i++
end
echo 'done with continuum mapping!' 
exit