#! /bin/sh

echo "Bandpass Filter Test" 

indata=ozdata.25.sw
rm -f tmp*
tpow=2.0

#------------------------------------------------
# Show the original Shotplot and Spectrum first..
#------------------------------------------------

sugain <$indata tpow=$tpow >tmp0 
   
suxwigb <tmp0 perc=90 xbox=10 ybox=10 wbox=400 hbox=600 \
        label1="Traveltime [s]" label2="Offset [m]"\
        title="Original Shot Gather" verbose=0 key=offset &

suwind <tmp0 key=offset min=-50 max=50 |
sustack key=dt | 
sugabor fmax=125 band=6 |
suximage xbox=420 ybox=10 wbox=400 hbox=600 bclip=2e05 wclip=0 \
         label1="Traveltime [s]" label2="Frequency [Hz]" \
         title="Frequency Spectrum (Near Offset)" \
         grid1=dot grid2=dot cmap=hsv2 legend=1 \
         verbose=0 &

#------------------------------------------------
# Bandpass Filter Test...
#------------------------------------------------

ok=false
while [ $ok = false ]
do
  rm -f tmp1
  if [ $selection ]
  then
    echo "Bandpass Filter Test"
    echo "Press A to add a filter"
    echo "Press S to start over"
    >/dev/tty
    read choice1
      case $choice1 in
      [sS]) cp tmp0 tmp1
            echo "-> Using original data";;
      [aA]) cp tmp2 tmp1
            echo "-> Using filtered data";;
      esac
  else
    cp tmp0 tmp1
  fi
  echo "Select a filterband [Hz]:"
  echo "Input: a,b,c,d  -  a:=lowcut    b:=lowpass"
  echo "                   c:=highpass  d:=highcut"
  >/dev/tty 
  read band
  sufilter <tmp1 f=$band amps=0,1,1,0 >tmp2
 
#------------------------------------------------
# Plot the filtered data...
#------------------------------------------------
   
  suxwigb <tmp2 xbox=10 ybox=10 wbox=400 hbox=600 \
          label1="Traveltime [s]" label2="Offset [m]"\
          title="Filtered Data" verbose=0 perc=90 \
          key=offset &

  suwind <tmp2 key=offset min=-50 max=50 |
  sustack key=dt | 
  sugabor fmax=125 band=6 |
  suximage xbox=420 ybox=10 wbox=400 hbox=600 \
           labeil1="Traveltime [s]" label2="Frequency [Hz]" \
           title="Frequency Spectrum (Near Offset)" \
           grid1=dot grid2=dot legend=1 bclip=2e05 wclip=0 \
           cmap=hsv2 verbose=0 &

  sufilter <tmp2 f=$band amps=1,0,0,1 |
  suxwigb  xbox=830 ybox=10 wbox=400 hbox=600 perc=90\
           label1="Traveltime [s]" label2="Offset [m]"\
           title="Rejected Data" verbose=0 perc=90 \
           key=offset &

#------------------------------------------------
# Go back or exit...
#------------------------------------------------

  echo "Press 1 for more filter testing"
  echo "Press 2 for PS output of current file and EXIT"
  >/dev/tty
  read selection
    case $selection in
	1) clear
           ok=false;;
	2) pause exit
	   ok=true;;
     esac
  done

rm -f tmp*

exit