#!/bin/csh

set inp = rejected.su

# f-k filtering  
suwind < $inp tmax=3.0 key=gelev | \
#sufilter f=10,20,100,120 amps=0,1,1,0|\
sudipfilt dx=20 dt=0.008 slopes=0.0,0.0002,0.0005 amps=1,0,0 > $inp.fkfilt.su

#Adding first arrival times
sushw key=d1 a=1000 < $inp | suxwigb xcur=.5 mpicks=pickfile_20_sources_B1_milisec.txt
a2b < pickfile_20_sources_B1_milisec.txt n1=2 > rejected_fa.bin 
b2a < rejected_fa.bin > tmp.ascii 
sushw < $inp.fkfilt.su infile=rejected_fa.bin key=tstat,tracf |\
sustatic hdrs=1 sign=-1 > $inp.fkfilt.flat.su 

#Generating corridor stack
sustack < $inp.fkfilt.flat.su key=cdp repeat=1 nrepeat=10 > $inp.stacked.su