#############################################################################
#              Script for the imaging of Orion 9 point mosaic
##############################################################################
## Initialize
## Include relevant tools if not already done:
##
include 'vlafiller.g';
include 'ms.g';
include 'autoflag.g';
include 'autoflag.g';
include 'flagger.g';
include 'msplot.g';
include 'calibrater.g';
include 'viewer.g';
include 'imager.g';
include 'image.g';
include 'imagepol.g';
include 'interactivemask.g'

# get a summary:
dc.summary('dvd-orion.ms');

#
########################################################################
# Here to the MEM Clean
#######################################################################
#
imgr:=imager('dvd-orion.ms');
imgr.setdata(mode='none',fieldid=[3:11],spwid=[1,2]);
imgr.setimage(nx=300,ny=300,cellx='2.0arcsec',
     celly='2.0arcsec',stokes='I',fieldid=7,spwid=[1,2]);

# imgr.weight(type='uniform',mosaic='T');
# uniform weighting also was not optimum, I got a better image using: 
imgr.weight(type="briggs" , rmode="norm" , robust=-1, mosaic=T);
imgr.setvp(dovp=T,usedefaultvp=T,dosquint=F);
imgr.setmfcontrol(cyclefactor=3.0);

# let's give it a first try without mask
# 1.
imgr.mem(algorithm='mfentropy',niter=10,sigma='10mJy',
	targetflux='100Jy', displayprogress=T,fixed=F,
	complist='',prior='',
        model='warmels-orion-vlamem1.mod',
	image='warmels-orion-vlamem1.img',
        residual='warmels-orion-vlamem1.res);
dv.gui();
# 
# Hmm, not realy what we wanted but now we can create a first mask
#
# imask:=interactivemask(refimage='warmels-orion-vlamem1.img',
#        mask='warmels-orion-vlamem1.msk');
# imask.start(); 
#
# Let's have a second MEM round and see what comes out
# 2.
imgr.mem(algorithm='mfentropy',niter=10,sigma='10mJy',
	targetflux='100Jy', displayprogress=T,fixed=F,
	complist='',prior='',
        model='warmels-orion-vlamem2.mod',
	image='warmels-orion-vlamem2.img',
        residual='warmels-orion-vlamem2.res',
	mask='warmels-orion-vlamem1.msk');
#
# Much better already. Now we can create a second mask
#
# imask:=interactivemask(refimage='orion-vlamem-rein2.img',
#        mask='orion-vlamem-rein2.msk');
# imask.start(); 
# 3.
imgr.mem(algorithm='mfentropy',niter=10,sigma='10mJy',
	targetflux='100Jy', displayprogress=T,fixed=F,
	complist='',prior='',
        model='warmels-orion-vlamem3.mod',
	image='warmels-orion-vlamem3.img',
        residual='warmels-orion-vlamem3.res',
	mask='warmels-orion-vlamem2.msk');
#
# imask:=interactivemask(refimage='orion-vlamem-rein3.img',
#        mask='orion-vlamem-rein3.msk');
# imask.start(); 
#
# 4.
imgr.mem(algorithm='mfentropy',niter=40,sigma='5mJy',
        targetflux='100Jy', displayprogress=T,fixed=F,
        complist='',prior='',
        model='warmels-orion-vlamem4.mod',
	image='warmels-orion-vlamem4.img',
        residual='warmels-orion-vlamem4.res',
	mask='warmels-orion-vlamem3.msk');
#
# Now do the image feathering
imgr.setvp(dovp=T,usedefaultvp=T);
imgr.feather(image='warmels-orion.feathered4.img',
     highres='warmels-orion-vlamem4.img', lowres='dvd-orion.gbt.img');
#
# That's it.
#