#
    # Set the identifier for this run
    #
    PREFIX   := 'test9.';         # This will prepend all output filenames
    #
    # Set the various paraments of imaging
    #

    MS       := 'axaf_avg.ms';         # Name of the MS
    MS       := '/home/imager-b/smyers/2005-06-22/axaf_avg.ms';
    FIELDID  := 1;                     # Field IDs to image
    STOKES   := 'I';                   # Stokes value
    BASE_IMAGE_NAME:='axaf.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.
    DOINTERACTIVE:= F;                 # Set it to T to do interactive
                                       # box setting in-between iterations.
    #
    # Set imaging/deconvolution control parms
    #

    ALGO     := 'wfclark';                 # The algorithm to use for
    IMSIZE   := [3000,3000];               # The image size
    CELLSIZE := ['2.5arcsec','2.5arcsec']; # The cellsize
    SPWID    := [1];                       # Spectral windows to image
    NFACETS  := 1;                         # No. of facets
    WPLANES  := 128;                       # No. of w-projection planes
    NITER    := 50000;                     # No. of Clean components
    THRESHOLD:= '70e-3mJy'                 # The Cleaning threshold.
    ROBUST   := 1;                         # Briggs robustness parameter
    GAIN     := 0.10;                      # Clean gain
    CYCFAC   := 0.10;                      # Major cycle speedup factor

    MODEL_IMAGE_NAME := spaste(PREFIX,BASE_IMAGE_NAME);
    CLEAN_IMAGE_NAME := spaste(MODEL_IMAGE_NAME,'.clean');
    RESID_IMAGE_NAME := spaste(MODEL_IMAGE_NAME,'.res');

# 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",robust=ROBUST);


    im.setoptions(ftmachine='wproject',wprojplanes=WPLANES);

    im.setmfcontrol(cyclefactor=CYCFAC);

    im.clean(algorithm=ALGO, gain=GAIN,
             niter=NITER,threshold=THRESHOLD,
             interactive=DOINTERACTIVE,
             model=MODEL_IMAGE_NAME,
             image=CLEAN_IMAGE_NAME,
             residual=RESID_IMAGE_NAME);

# Statistics (box near center)

    include 'image.g';
    img := image(CLEAN_IMAGE_NAME);
    pfull := drm.box([1,1,1], [IMSIZE,1]);
    img.statistics(statsout=astats, region=pfull);

    print CLEAN_IMAGE_NAME, " : ",
            " Max= ", sprintf("%7.1f", 1.e3*astats.max)," mJy";

    pregion := drm.box([1377,1408,1], [1542,1554,1]); # imsiz 3000
#    pregion := drm.box([2377,2408,1], [2542,2554,1]); # imsiz 5000
    img.statistics(statsout=pstats, region=pregion);

    print CLEAN_IMAGE_NAME, " : ",
            " RMS= ", sprintf("%7.1f", 1.e6*pstats.rms)," uJy";