#! /bin/sh

data="olympia_line2.su"

# general parameters for processing data
perc=97			#clip value for display of traces
wagc=0.01		#AGC window length
f1=50,100,1100,2000 	#bandpass filter frequencies (Hz)
f2=40,80,800,1400
num_cdp=1		#number of cdps (seconds in horizontal direction) to stack
#tmin=0.00		#minimum time to display
tmax=0.201		#maximum time to display
time_interval=0.0166666667
time_interval=0.00833333333


ve=8.0
title="Totten line 2 sparker, decon"


suwind <$data tmax=$tmax |
sushw key=delrt a=0 |
sushw key=tstat a=8 |
sustatic hdrs=1 |
#sustack key=cdp |
#sushw key=dt a=250 |
#suresamp nt=801 |


#sumute key=tracl ntaper=100 xmute=1 tmute=0.006 |
sugain tpow=0.8 |
#supef minlag=0.010 maxlag=0.0604 pnoise=0.001 mincorr=0.01 maxcor = 0.2 |
sufilter f=$f1 |
sugain pbal=1 |
supef minlag=0.0016 maxlag=0.03 pnoise=0.001 mincorr=0.01 maxcor = 0.2 |
#sufilter f=$f2 |
#sugain tpow=-0.9 |

#migrate the data
#sushw key=cdp a=1 b=1 |
#sustolt cdpmin=1 cdpmax=4448 tmig=0.0 vmig=1500 smig=1.0 fmax=1500 \
#        vscale=1.0 ltaper=40 lbtaper=200 dxcdp=1.67 |
#sushw key=cdp a=093252 b=1.6666666666 |

suchw key1=dt key2=dt b=0.75 >olympia_line2_filt.su






# first panel of data  ******************
starttime=225200
endtime=230000
lefttic=2252
righttic=2300

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_1.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=4.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1

# second panel of data  ******************
starttime=230000
endtime=232000
lefttic=2300
righttic=2320

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_2.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=10.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1
	

# third panel of data  ******************
starttime=232000
endtime=234000
lefttic=2320
righttic=2340

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_3.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=10.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1


# fourth panel of data  ******************
starttime=234000
endtime=240000
lefttic=2340
righttic=2400

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_4.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=10.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1


# fifth panel of data  ******************
starttime=000000
endtime=002000
lefttic=00
righttic=20

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_5.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=10.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1


# sixth panel of data  ******************
starttime=002000
endtime=004000
lefttic=20
righttic=40

suwind <olympia_line2_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
suwind tmax=0.143 |
#suflip flip=2 |
supsimage >totten2_6.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=10.0\
	ybox=0.8\
	xbox=1.5\
	label1="approximate depth (km)"\
	label2="$title \n      v.e. ~ $ve"\
	titlefont=Helvetica-Bold\
	titlesize=16\
	labelsize=12\
	style=normal\
	interp=1
#	xcur=2.0 wt=0 va=1


# seventh panel of data  ******************
#starttime=143200
#endtime=145200
#lefttic=1452
#righttic=1432

#suwind <olympia_line1_filt.su key=cdp min=$starttime max=$endtime |
#suchw key1=cdp key2=cdp d=100 |
#suwind tmax=0.143 |
#suflip flip=2 |
#supsimage >totten1_7.ps\
#	perc=$perc\
#	n2tic=2\
#	d2=$time_interval\
#	f2=$lefttic\
#	x2beg=$lefttic\
#	x2end=$righttic\
#	nbpi=600\
#	wbox=6.00\
#	hbox=10.0\
#	ybox=0.8\
#	xbox=1.5\
#	label1="approximate depth (km)"\
#	label2="$title \n      v.e. ~ $ve"\
#	titlefont=Helvetica-Bold\
#	titlesize=16\
#	labelsize=12\
#	style=normal\
#	interp=1
##	xcur=2.0 wt=0 va=1