#
    # Set the various paraments of imaging
    #
    FITSFILE := '/home/rohir2/jmcmulli/EVLATST1/Data/Shk116/1046P59.ABCD34.FITS';
    MS       := '1034.34.ms';          # Name of the MS
    ALGO     := 'wfclark';             # The algorithm to use for
    IMSIZE   := [7200,7200];           # The image size
    CELLSIZE := ['0.5arcsec','0.5arcsec']; # The cellsize
    SPWID    := [1,2];                 # Spectral windows to image
    FIELDID  := 1;                     # Field IDs to image
    NFACETS  := 2;                     # No. of facets
    WPROJPLANES:=128;                  # No. of gridding convolution functions used in 
                                       #     w-project fourier transform
    STOKES   := 'I';                   # Stokes value
    MODEL_IMAGE_NAME:='1046.im';       # Name of the image files.
                                       # The restored Clean-ed image
                                       # and the residual imgaes are
                                       # stored in images with
                                       # ``.clean'' and ``.res''
                                       # appended to this name.
    NITER    := 10000;                 # No. of Clean components
    THRESHOLD:= '100e-3mJy'            # The Cleaning threshold.  The
                                       # iterations are stopped when
                                       # the magnitude of the
                                       # strongest components is lower
                                       # than this value.
    DOINTERACTIVE:= F;                 # Set it to T to do interactive
                                       # box setting in-between iterations.

# ----

# Import
    include 'ms.g';
    myms:=fitstoms(fitsfile=FITSFILE,msfile=MS);
    myms.done();

# ----

# Image

    include 'imager.g';
    im:=imager(MS);
    im.setdata(spwid=SPWID,fieldid=FIELDID,mode='none');
    im.setimage(nx=IMSIZE[1],ny=IMSIZE[2],
                cellx=CELLSIZE[1],celly=CELLSIZE[2],
                facets=NFACETS,stokes=stokes,
                spwid=SPWID,fieldid=FIELDID);
     im.weight(type="briggs");

#    Options
#    im.weight(type="briggs" , rmode="norm", robust=0.5,
                      fieldofview='3arcmin')
#    im.filter(bmaj='1.8arcsec', bmin='1.8arcsec');

    im.setoptions(ftmachine='wproject',wprojplanes=WPROJPLANES);
    im.clean(algorithm=ALGO,
             niter=NITER,threshold=THRESHOLD,
             interactive=DOINTERACTIVE,
             model=MODEL_IMAGE_NAME,
             image=spaste(MODEL_IMAGE_NAME,'.clean'),
             residual=spaste(MODEL_IMAGE_NAME,'.res'));