#ALMA moment mapping script created by George J. Bendo for the ALMA Data
#Processing Workshop at UCLan (23-24 January 2019).
#
#This is a custom-written script created for use with the NGC 3256 data imaged
#during the workshop.  The script cound be executed using the command
#execfile('momentmap.py') at the command line, or the individual lines could
#be copied and pasted into a CASA session.



#Create the moment 0 map, which is based on the integrated line intensity.
os.system('rm -rf ngc3256_sci.spw1.cube.I.manual.image.mom0')
immoments(imagename='ngc3256_sci.spw1.cube.I.manual.image',
  moments=[0],
  axis='spectral',
  region='box[[70pix,70pix],[170pix,170pix]]',
  chans='85~110',
  outfile='ngc3256_sci.spw1.cube.I.manual.image.mom0')

#Create the moment 1 map, which is equivalent to the mean velocity.
os.system('rm -rf ngc3256_sci.spw1.cube.I.manual.image.mom1')
immoments(imagename='ngc3256_sci.spw1.cube.I.manual.image',
  moments=[1],
  axis='spectral',
  region='box[[70pix,70pix],[170pix,170pix]]',
  chans='85~110',
  excludepix=[1.5e-3],
  outfile='ngc3256_sci.spw1.cube.I.manual.image.mom1')

#Create the moment 2 map, which is equivalent to the velocity dispersion.
os.system('rm -rf ngc3256_sci.spw1.cube.I.manual.image.mom2')
immoments(imagename='ngc3256_sci.spw1.cube.I.manual.image',
  moments=[2],
  axis='spectral',
  region='box[[70pix,70pix],[170pix,170pix]]',
  chans='85~110',
  excludepix=[1.5e-3],
  outfile='ngc3256_sci.spw1.cube.I.manual.image.mom2')