include 'ms.g'

myms1:=ms('phIII_64ant.ms, readonly=F);

#Output average 60 channel from spw 3 (which has 64 channels) 
# to 1 output channel 
myms1.split(outputms='Cont3mm.ms', fieldids=[3], spwids=[3], 
		  nchan=1, start=3, step=60, whichcol='CORRECTED_DATA');

#select 48 channel no averaging from spw 7 (which has 256 channels)

myms1.split(outputms='Line3mm.ms', fieldids=[3], spwids=[7], nchan=48, 
 start=105, step=1, whichcol='CORRECTED_DATA');

# select multiple spw and average

myms1.split(outputms='Cont1mm.ms', fieldids=[3], spwids=[11,12,15,16,19,2], 
            nchan=1, start=4, step=58, whichcol='CORRECTED_DATA');


# Similar thing could be achieved this way with different selections
# for each spw if needed....


myms1.split(outputms='Cont1mmb.ms', fieldids=[3], spwids=[11,12,15,16,19,2], 
            nchan=[1,1,1,1,1,1], start=[4,4,4,4,4,4], 
            step=[58,58,58,58,58,58], whichcol='CORRECTED_DATA');




myms1.done()