# AIPS++ path used to combine BIMA data of NGC 4826

# 1. concatenate data in new directory

shell ('rm -r ngc4826.combo.ms');
shell ('cp -r ../98apr16/ngc4826.split.ms ngc4826.combo.ms');
shell ('cp -r ../98apr16/ngc4826.split.ms 98apr16.ngc4826.ms');
shell ('cp -r ../98apr22/ngc4826.split.ms 98apr22.ngc4826.ms');

include 'ms.g';
mset:=ms(filename='ngc4826.combo.ms',readonly=F);
mset.concatenate(msfile='98apr22.ngc4826.ms',freqtol='0.1MHz',dirtol='0.1arcsec');
mset.summary();
mset.done();

# 2. Make an image of the target fields
include 'imager.g';
imgr:=imager(filename='ngc4826.combo.ms');
imgr.setdata(fieldid=[1:7],spwid=[1:3],nchan=120,start=47,step=1);
imgr.setimage(nx=256,ny=256,cellx='1arcsec',celly='1arcsec',stokes='I',mode='channel',nchan=30,start=47,step=4,fieldid=1,spwid=[1:3]);

imgr.weight(mode='natural',mosaic=T); # this step did not work
imgr.setvp(dovp=T);
imgr.setoptions(padding=1.0,ftmachine='mosaic');
imgr.clean(algorithm='mfclark',niter=300,gain=0.1,model='ngc4826.cln',residual='ngc4826.resid',image='ngc4826.cm');

# 3. Get the zeroth and first moments
include 'image.g';
im:=image(infile='ngc4826.cm');
im_mom0:=im.moments(outfile='ngc4826.mom0.cm',moments=0,axis=4,mask='indexin(4,[5:26])',includepix=[0.070,1000.0]);
im_mom1:=im.moments(outfile='ngc4826.mom1.cm',moments=1,axis=4,mask='indexin(4,[5:26])',includepix=[0.28,1000.0]);
im_mom0.done();
im_mom1.done();
im.done();

# 4. Check out the moment maps!
include 'viewer.g';
dv.gui();