#! /bin/bash

#Compute the egression power of DOPP_X_NU_5-7mHz.fits for a large pupil,
#and make a temporal movie of it.

#FREQ_MIN='0.0025'
#FREQ_MAX='0.0045'
#FREQ_MIN='0.0050'
#FREQ_MAX='0.0070'
#FREQ_MIN='0.0070'
#FREQ_MAX='0.0090'
FREQ_MIN='0.0090'
FREQ_MAX='0.0110'

#...

rm -f BANDPASS.fits
get_zpass DOPP_X_NU.fits $FREQ_MIN $FREQ_MAX BANDPASS.fits

NX=`print_head BANDPASS.fits | grep NAXIS1 | awk '{print $3;}'`
NY=`print_head BANDPASS.fits | grep NAXIS2 | awk '{print $3;}'`
CX=`echo $NX | awk '{print 0.5*$1;}'`
CY=`echo $NY | awk '{print 0.5*$1;}'`

#rm -f PUPIL.fits
#ann_mask BANDPASS.fits $CX $CY 0.010 0.12 PUPIL.fits

SUFFIX=`echo $FREQ_MIN $FREQ_MAX | awk '{
    printf("%.1f-%.1fmHz.fits", 1000*$1, 1000*$2);
}'`

#rm -f EGRESS_X_NU_$SUFFIX
#regress_plane_eikonal BANDPASS.fits 0.0 1 egr EGRESS_X_NU_$SUFFIX \
#-p PUPIL.fits -dispers 2.0

#rm -f EGRESS_PWR_X_T_$SUFFIX
#nu2t_pwr EGRESS_X_NU_$SUFFIX EGRESS_PWR_X_T_$SUFFIX -samp 45.0

rm -f DOPP_PWR_X_T_9.0-11.0mHz.fits
nu2t_pwr BANDPASS.fits DOPP_PWR_X_T_9.0-11.0mHz.fits -samp 45.0

#Cleanup

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