#!/bin/bash -l
#SBATCH --ntasks=1
#SBATCH --ntasks-per-node=1
#SBATCH --time=12:00:00
#SBATCH --job-name=mssplit
#SBATCH --no-requeue
#SBATCH --export=NONE

# Runs mssplit for 30 beams and 54 channel averaging as an 
# array job on the first argument. 
# One output measurement set is produced per beam
# Selects only the scan where the selected beam is centered on 1934.
# Parsets for the mssplit command are generated in the current directory
# Run with the following command to do beams 0 to 29
# 
# sbatch --array 0-29 mssplit.job msname

beam=`echo ${SLURM_ARRAY_TASK_ID} | awk '{printf "%02d",$1}'`
msin=mssplit_${SLURM_JOB_ID}_b${beam}.in
ms=$1

cat > $msin <<EOF

# Input measurement set
vis         = $ms

# Output measurement set
outputvis   = beam${beam}.ms

# Beam and Scan selection
beams       = [${SLURM_ARRAY_TASK_ID}]
scans       = [${SLURM_ARRAY_TASK_ID}]

# The input channel range, indexing is one-based.
channel     = 1-2592

# Defines the number of channel to average to form the one output channel
width       = 54


EOF

aprun -n 1 -N 1 mssplit -c $msin