#! /bin/bash

#Extract the 5--7-mHz interval of the Doppler spectrum 
#(DOPP_X_NU_5.0-7.0mHz.fits), and make a temporal movie of it.

INPUT='DOPP.fits'

rm -f J01_TMP01.fits
switch_y_z $INPUT J01_TMP01.fits
rm -f J01_TMP02.fits
switch_x_y J01_TMP01.fits J01_TMP02.fits

rm -f BAD_PIXELS.txt
bc_filter J01_TMP02.fits 100 -hp -b BAD_PIXELS.txt \
    -t 10000.0 #-e MISSING_IMAGES_DOPP.txt 

rm -f J01_TMP03.fits
rfourier_x J01_TMP02.fits J01_TMP03.fits -p 20

rm -f J01_TMP04.fits
switch_x_y J01_TMP03.fits J01_TMP04.fits
rm -f DOPP_X_NU.fits
switch_y_z J01_TMP04.fits DOPP_X_NU.fits

#Cleanup

rm -f J01_TMP0[1-4].fits
#rm -f BAD_PIXELS.txt