# Start with .FITS file corresponding to UV data with flanking fields removed

fitsname := '/home/imager-c/wbrisken/tst2.0/1046P59.ABCD34.FITS'
msname   := 'shk116.ms'

filldata := function()
{
	include 'ms.g'

	ms := fitstoms(msname, fitsname)
	ms.done()
}

smalldirtyimage := function()
{
	include 'imager.g'

	imgr := imager(msname)
	imgr.setdata(fieldid=1, spwid=[1:2])
	imgr.setimage(nx=1000, ny=1000, cellx='10arcsec', celly='10arcsec',
		stokes='I', mode='mfs', spwid=[1,2], fieldid=1)
	imgr.uvrange(0, 10000.0)
	imgr.setoptions(ftmachine=wproject,wprojplanes=256)
	imgr.weight(type='briggs', rmode='norm', robust=-2)
	imgr.makeimage(type='corrected', image='smalldirty.im')
	imgr.done()
}


hugerimage1 := function()
{
	include 'imager.g'

	imgr := imager(msname)
	imgr.setdata(fieldid=1, spwid=[1])
	imgr.setimage(nx=7500, ny=6912, cellx='0.36arcsec', celly='0.36arcsec',
		stokes='I', mode='mfs', spwid=[1], fieldid=1)
	imgr.setoptions(ftmachine='wproject',wprojplanes=256)
	imgr.weight(type='briggs', rmode='norm', robust=0.4)
	imgr.clean(algorithm='wfhogbom', 
	    model='huger.1.model',
	    residual='huger.1.residual',
	    image='huger.1.im',
	    gain=0.2,
	    niter=120000)
	imgr.done()
}

hugerimage12 := function()
{
	include 'imager.g'

	imgr := imager(msname)
	imgr.setdata(fieldid=1, spwid=[1,2])
	imgr.setimage(nx=7500, ny=6912, cellx='0.36arcsec', celly='0.36arcsec',
		stokes='I', mode='mfs', spwid=[1,2], fieldid=1)
	imgr.setoptions(ftmachine='wproject',wprojplanes=256)
	imgr.weight(type='briggs', rmode='norm', robust=0.4)
	imgr.clean(algorithm='wfhogbom', 
	    model='huger.12.model',
	    residual='huger.12.residual',
	    image='huger.12.im',
	    gain=0.2,
	    niter=120000)
	imgr.done()
}

hugerimage2 := function()
{
	include 'imager.g'

	imgr := imager(msname)
	imgr.setdata(fieldid=1, spwid=[2])
	imgr.setimage(nx=7500, ny=6912, cellx='0.36arcsec', celly='0.36arcsec',
		stokes='I', mode='mfs', spwid=[2], fieldid=1)
	imgr.setoptions(ftmachine='wproject',wprojplanes=256)
	imgr.weight(type='briggs', rmode='norm', robust=0.4)
	imgr.clean(algorithm='wfhogbom', 
	    model='huger.2.model',
	    residual='huger.2.residual',
	    image='huger.2.im',
	    gain=0.2,
	    niter=120000)
	imgr.done()
}

hugeraddifs := function()
{
	include 'image.g'

	im := imagecalc(outfile='huger.im',pixels='0.5*(huger.1.im+huger.2.im)')
	im.done()
}

fullredux := function()
{
	filldata()
	hugerimage1()
	hugerimage2()
	hugeraddifs()
}