#! /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 - 2}'`
RANGE2=`echo $NZ | awk '{print 1 "-" $1 - 1}'`

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

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

rm -f CONSEC_DIFFS.fits
bin_op JUNK2.fits minus JUNK1.fits -o CONSEC_DIFFS.fits

#Cleanup

rm -f JUNK[12].fits