#! /bin/sh

vnmo=1500

interoff=-262
offref=-3237
depthref=1000
pmin=-2000
pmax=2000
dp=8
igopt=2        

pmula=-800
pmulb=47

data=gain.jon=1.cdp.su

# view gather in the radon transform domain
choose=0             # =1 use suradon as a radon domain filter

#sunmo vnmo=$vnmo  <  $data |
#suradon offref=$offref depthref=$depthref  \
#interoff=$interoff pmin=$pmin pmax=$pmax dp=$dp choose=$choose igopt=$igopt  |
#suximage perc=95 title="original" &


# filter data in radon domain
choose=1  # data - multiples
sunmo vnmo=$vnmo  <  $data |
suradon offref=$offref depthref=$depthref  pmula=$pmula pmulb=$pmulb \
interoff=$interoff pmin=$pmin pmax=$pmax dp=$dp choose=$choose igopt=$igopt |
sunmo vnmo=$vnmo invert=1 > radon.$data


choose=2          # capture multiples
sunmo vnmo=$vnmo  <  $data |
suradon offref=$offref depthref=$depthref  pmula=$pmula pmulb=$pmulb \
interoff=$interoff pmin=$pmin pmax=$pmax dp=$dp choose=$choose igopt=$igopt |
sunmo vnmo=$vnmo invert=1 > multiples.$data

exit 0

# view data minus multiples in radon domain
choose=0
sunmo vnmo=$vnmo  <  junk1.su |
suradon offref=$offref depthref=$depthref  \
interoff=$interoff pmin=$pmin pmax=$pmax dp=$dp choose=$choose igopt=$igopt  |
suximage perc=95 title="data minus multiples " &

# view multiples in radon domain
choose=0
sunmo vnmo=$vnmo  <  junk2.su |
suradon offref=$offref depthref=$depthref  \
interoff=$interoff pmin=$pmin pmax=$pmax dp=$dp choose=$choose igopt=$igopt  |
suximage perc=99 title="multiples" &

# view semblance panels
suvelan nv=150 dv=15 fv=1450 < junk1.su |
suximage perc=99 cmap=hsv2 f2=1450 d2=15  title="data minus multiples" legend=1 &
suvelan nv=150 dv=15 fv=1450 < junk2.su |
suximage perc=99 cmap=hsv2 f2=1450 d2=15  title="multiples" legend=1 &

exit 0


exit 0


