#! /bin/bash

#Create a datacube that is the consecutive differences in the datacube
#named "$INPUT".

INPUT='DOPP.fits'

NZ=`print_head $INPUT | grep NAXIS3 | awk '{print $3;}'`

RANGE1=`echo $NZ | awk '{print 0 "-" $1 - 3}'`
RANGE2=`echo $NZ | awk '{print 1 "-" $1 - 2}'`
RANGE3=`echo $NZ | awk '{print 2 "-" $1 - 1}'`

rm -f JUNK1.fits
restack $INPUT $RANGE1 JUNK1.fits

rm -f JUNK2.fits
restack $INPUT $RANGE2 JUNK2.fits

rm -f JUNK3.fits
restack $INPUT $RANGE3 JUNK3.fits

rm -f SUM13.fits
bin_op JUNK1.fits plus JUNK3.fits -o SUM13.fits

rm -f HALF_SUM13.fits
mon_op SUM13.fits times 0.5 -o HALF_SUM13.fits

rm -f XIENT.fits
bin_op JUNK2.fits minus HALF_SUM13.fits -o XIENT.fits

#Cleanup

rm -f JUNK[123].fits SUM13.fits HALF_SUM13.fits