#! /bin/sh
## next steps generates synthetic vibe sweep
#suvibro dt=0.002 sweep=1 f1=15 f2=120 tv=12 phz=90 t1=0.1 t2=0.1 taper=1 >synsweep.su

## next 2 steps convert data from seg2 into su format
#seg2segy 2001.dat 510
#segyread tape=2001.sgy verbose=1 endian=0 conv=0 | segyclean > 2001.su
## next 2 steps do same as above for day2 shots
#seg2segy 2511.dat 496
#segyread tape=2511.sgy verbose=1 endian=0 conv=0 | segyclean > 2511.su
## next 2 steps do same as above for day3 shots # below should be 455
#seg2segy 3007.dat 461
#segyread tape=3007.sgy verbose=1 endian=0 conv=0 | segyclean > 3007.su
## next 2 steps do same as above for day4 shots
#seg2segy 3462.dat 501
#segyread tape=3462.sgy verbose=1 endian=0 conv=0 | segyclean > 3462.su
## next 2 steps do same as above for day4 last few shots
#seg2segy 3963.dat 26
#segyread tape=3963.sgy verbose=1 endian=0 conv=0 | segyclean > 3963.su
## next steps do vibroseis correlation
#sugain agc=1 wagc=0.4 < 3963.su | suxcor sufile=synsweep.su vibroseis=1000 > 3963.corr

## next steps generates geometry
### another example for brute geometry is in file stack_tlp or geom_tlp (real geom)
#sugethw < 1112.su output=geom key=fldr,tracf,offset,cdp,sx,sy,gx,gy > hdrfile
##edit hdrfile
#a2b--this did not work with excel output because of dos vs unit characters(?)
### OR
#suchw < 3963.corr key1=cdp,offset key2=gx,gx key3=sx,sx b=1,5 c=1,-5 >3963.geom

## next steps do sort, agc, bp, nmo, stack
vnmo=1500,1800,2400
tnmo=0,.4,1.0
#suwind < data.geom key=tracf min=1 max=144 tmax=2.0 dt=.002 > datatmp.geom
#susort < datatmp.geom cdp offset > data.cdp
sushw <data.cdp key=sstat a=-180 | sustatic sign=-1 hdrs=1 | 
sugain agc=1 wagc=0.4 | 
#sugain < data.cdp agc=1 wagc=0.4 | 
sufilter f=15,25,100,120 | sufilter f=55,57,63,65 amps=1,0,0,1 |
sunmo vnmo=$vnmo tnmo=$tnmo | sustack key=cdp verbose=1 > brtstk4
#suximage