#############################################################################
#              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');
#
############################################################################
# Directly convert the single-dish GBT image to a model and use this a 
# a starting model for the deconvolution of the high resolution data.
############################################################################
#
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);
imgr.setmfcontrol(stoplargenegatives=-1,cyclefactor=3, 
    cyclespeedup=500);

# setup the scales sizes for multiscale Clean
imgr.setscales(scalemethod='uservector',uservector=[0,3,10,30]);
#
# set the multifield parameters
#
#
# Make a model from the GBT single disk observation
imgr.makemodelfromsd(sdimage='dvd-orion.gbt.im',
    modelimage='warmels-orion-gbt.mod1',maskimage='');
# 
# Now the multiclean clean setup and the model defined. Use the same mask as was
# used for the MEM clean
imgr.clean(algorithm='mfmultiscale',gain=0.2,niter=10000,threshold='0.01Jy',
    model='warmels-orion-gbt.mod1',mask='warmels-orion-vlamem3.msk', 
    image='warmels-orion-combined1.img',residual='warmels-orion-combined1.res');
#
#imgr.done();
#imgr:=image('warmels-orion-combined.img');
#imgr.stats();
#imgr.done()
#
# We are now done with the first round of combined deconvolution
# Let's see if we need a second round with adjusted parammeters
###
imgr.makemodelfromsd(sdimage='dvd-orion.gbt.im',
    modelimage='warmels-orion-gbt.mod2',maskimage='');
#
# let's run multiscale clean again and now more interactively
imgr.clean(algorithm='mfmultiscale',gain=0.2,niter=3000,threshold='0.01Jy',
    model='warmels-orion-gbt.mod2',mask='warmels-orion-vlamem3.msk', 
    image='warmels-orion-combined2.img',residual='warmels-orion-combined2.res',
    interactive=T,npercycle=500);
#
#imgr.done();
#imgr:=image('warmels-orion-combined.img');
#imgr.stats();
#imgr.done()
#
###
imgr.makemodelfromsd(sdimage='dvd-orion.gbt.im',
    modelimage='warmels-orion-gbt.mod3',maskimage='');
#
# let's run clean again and now more interactively
imgr.clean(algorithm='mfmultiscale',gain=0.1,niter=15000,threshold='0.01Jy',
    model='warmels-orion-gbt.mod3',mask='warmels-orion-vlamem3.msk', 
    image='warmels-orion-combined3.img',residual='warmels-orion-combined3.res');
#
#imgr.done();
#imgr:=image('warmels-orion-combined.img');
#imgr.stats();
#imgr.done()
#
###
#imgr.makemodelfromsd(sdimage='dvd-orion.gbt.im',
#    modelimage='warmels-orion-gbt.mod4',maskimage='');
#
#imgr.clean(algorithm='mfmultiscale',gain=0.2,niter=5000,threshold='0.01Jy',
#    model='warmels-orion-gbt.mod4',mask='warmels-orion-vlamem3.msk', 
#    image='warmels-orion-combined4.img',residual='warmels-combined4.res',
#    interactive=T,npercycle=500);