#! /bin/sh

#
suwind itmin=0 itmax=1500 < farfield_gun.su > big_gun.su


# fake data no multiples no noise
sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=.5  dt=.004 vel=1500  |
suaddevent type=nmo nt=1500 t0=.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=.9  dt=.004 vel=1600 |
suaddevent type=nmo nt=1500 t0=1.0  dt=.004 vel=1620 |
suaddevent type=nmo nt=1500 t0=1.1  dt=.004 vel=1650 |
suaddevent type=nmo nt=1500 t0=1.3  dt=.004 vel=1700 |
suaddevent type=nmo nt=1500 t0=1.6  dt=.004 vel=1850 |
suaddevent type=nmo nt=1500 t0=1.9  dt=.004 vel=2000 |
suaddevent type=nmo nt=1500 t0=2.0  dt=.004 vel=2050 amp=2 |
suaddevent type=nmo nt=1500 t0=2.1  dt=.004 vel=2100 |
suaddevent type=nmo nt=1500 t0=2.4  dt=.004 vel=2150 |
suaddevent type=nmo nt=1500 t0=2.5  dt=.004 vel=2200 |
suaddevent type=nmo nt=1500 t0=2.7  dt=.004 vel=2250 |
suaddevent type=nmo nt=1500 t0=3.0  dt=.004 vel=2350 |
suaddevent type=nmo nt=1500 t0=3.1  dt=.004 vel=2400 |
suaddevent type=nmo nt=1500 t0=3.3  dt=.004 vel=2450 |
suaddevent type=nmo nt=1500 t0=3.5  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=3.7  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=4.0  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=4.1  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.3  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.5  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.6  dt=.004 vel=2600 |
suaddevent type=nmo nt=1500 t0=5.0  dt=.004 vel=2600 |
suaddevent type=nmo nt=1500 t0=5.1  dt=.004 vel=2650 |
suaddevent type=nmo nt=1500 t0=5.5  dt=.004 vel=2650 |
sufilter f=0,5,40,50  > fake_no_noise.su

# fake data no multiples
sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=.5  dt=.004 vel=1500  |
suaddevent type=nmo nt=1500 t0=.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=.9  dt=.004 vel=1600 |
suaddevent type=nmo nt=1500 t0=1.0  dt=.004 vel=1620 |
suaddevent type=nmo nt=1500 t0=1.1  dt=.004 vel=1650 |
suaddevent type=nmo nt=1500 t0=1.3  dt=.004 vel=1700 |
suaddevent type=nmo nt=1500 t0=1.6  dt=.004 vel=1850 |
suaddevent type=nmo nt=1500 t0=1.9  dt=.004 vel=2000 |
suaddevent type=nmo nt=1500 t0=2.0  dt=.004 vel=2050 amp=2 |
suaddevent type=nmo nt=1500 t0=2.1  dt=.004 vel=2100 |
suaddevent type=nmo nt=1500 t0=2.4  dt=.004 vel=2150 |
suaddevent type=nmo nt=1500 t0=2.5  dt=.004 vel=2200 |
suaddevent type=nmo nt=1500 t0=2.7  dt=.004 vel=2250 |
suaddevent type=nmo nt=1500 t0=3.0  dt=.004 vel=2350 |
suaddevent type=nmo nt=1500 t0=3.1  dt=.004 vel=2400 |
suaddevent type=nmo nt=1500 t0=3.3  dt=.004 vel=2450 |
suaddevent type=nmo nt=1500 t0=3.5  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=3.7  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=4.0  dt=.004 vel=2500 |
suaddevent type=nmo nt=1500 t0=4.1  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.3  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.5  dt=.004 vel=2550 |
suaddevent type=nmo nt=1500 t0=4.6  dt=.004 vel=2600 |
suaddevent type=nmo nt=1500 t0=5.0  dt=.004 vel=2600 |
suaddevent type=nmo nt=1500 t0=5.1  dt=.004 vel=2650 |
suaddevent type=nmo nt=1500 t0=5.5  dt=.004 vel=2650 |
sufilter f=0,5,40,50 | suaddnoise sn=50   > fake.su

# fake data waterbottom multiples only
sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=1.0  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=1.1  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=1.2  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=1.3  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=1.4  dt=.004 vel=1600 |
sugain scale=-1 > waterbottom0.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=1.5  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=1.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=1.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=1.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=1.9  dt=.004 vel=1600  > waterbottom1.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=2.0  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=2.1  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=2.2  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=2.3  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=2.4  dt=.004 vel=1600 |
sugain scale=-1 > waterbottom2.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=2.5  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=2.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=2.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=2.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=2.9  dt=.004 vel=1600  > waterbottom3.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=3.0  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=3.1  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=3.2  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=3.3  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=3.4  dt=.004 vel=1600 |
sugain scale=-1 > waterbottom4.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=3.5  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=3.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=3.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=3.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=3.9  dt=.004 vel=1600  > waterbottom5.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=4.0  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=4.1  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=4.2  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=4.3  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=4.4  dt=.004 vel=1600 |
sugain scale=-1 > waterbottom6.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=4.5  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=4.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=4.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=4.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=4.9  dt=.004 vel=1600 > waterbottom7.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=5.0  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=5.1  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=5.2  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=5.3  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=5.4  dt=.004 vel=1600 |
sugain scale=-1 > waterbottom8.su

sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=5.5  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=5.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=5.7  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=5.8  dt=.004 vel=1580 |
suaddevent type=nmo nt=1500 t0=5.9  dt=.004 vel=1600  > waterbottom9.su

susum waterbottom0.su waterbottom1.su > temp1.su
susum waterbottom2.su temp1.su > temp.su
susum waterbottom3.su temp.su  > temp1.su
susum waterbottom4.su temp1.su > temp.su
susum waterbottom5.su temp.su > temp1.su
susum waterbottom6.su temp1.su > temp.su
susum waterbottom7.su temp.su > temp1.su
susum waterbottom8.su temp1.su > temp.su
susum waterbottom9.su temp.su > waterbottom.su

# pegleg multiples
sunull nt=1500 ntr=60 dt=.004 |
sushw key=gx,offset a=3300,-3237 b=25,50 c=1,1 |
suchw key1=sx key2=gx key3=offset a=0 b=1 c=1 |
suaddevent type=nmo nt=1500 t0=1.6  dt=.004 vel=1600 |
suaddevent type=nmo nt=1500 t0=2.1  dt=.004 vel=1550 |
suaddevent type=nmo nt=1500 t0=2.6  dt=.004 vel=1525 |
suaddevent type=nmo nt=1500 t0=3.1  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=3.6  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=4.1  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=4.6  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=5.1  dt=.004 vel=1500 |
suaddevent type=nmo nt=1500 t0=5.6  dt=.004 vel=1500 |
sufilter f=0,5,40,50 | suaddnoise sn=50 > pegleg1.1sec.su


suaddevent < pegleg1.1sec.su type=nmo nt=1500 t0=2.5  dt=.004 vel=2000 |
suaddevent type=nmo nt=1500 t0=3.0  dt=.004 vel=1950 |
suaddevent type=nmo nt=1500 t0=3.5  dt=.004 vel=1900 |
suaddevent type=nmo nt=1500 t0=4.0  dt=.004 vel=1850 |
suaddevent type=nmo nt=1500 t0=4.5  dt=.004 vel=1800 |
suaddevent type=nmo nt=1500 t0=5.0  dt=.004 vel=1750 |
suaddevent type=nmo nt=1500 t0=5.5  dt=.004 vel=1700  \
> pegleg1.1+2.0sec.su


suaddevent <  pegleg1.1+2.0sec.su type=nmo nt=1500 t0=3.8  dt=.004 vel=2400 |
suaddevent type=nmo nt=1500 t0=4.3  dt=.004 vel=2350 |
suaddevent type=nmo nt=1500 t0=4.8  dt=.004 vel=2300 |
suaddevent type=nmo nt=1500 t0=5.1  dt=.004 vel=2250 |
suaddevent type=nmo nt=1500 t0=5.6  dt=.004 vel=2200 \
> pegleg1.1+2.0+3.3sec.su

cp pegleg1.1+2.0+3.3sec.su pegleg.su


susum fake.su waterbottom.su > fake+water.su

susum  fake+water.su pegleg.su > fake+water+pegleg.su

suconv sufile=big_gun.su < fake+water+pegleg.su |
suwind itmin=0 itmax=1499  > fake+water+pegleg_wavelet.su

exit 0
