#! /bin/sh

data="case_1_1.su"

# general parameters for processing data
perc=97			#clip value for display of traces
wagc=0.01		#AGC window length
f1=200,300,1400,2000 	#bandpass filter frequencies (Hz)
f2=200,300,1000,1400
num_cdp=1		#number of cdps (seconds in horizontal direction) to stack
#tmin=0.00		#minimum time to display
tmax=0.211              #maximum time to display
time_interval=0.0166666667
time_interval=0.00833333333
dx=0.8333


ve=8.0
title="Case Inlet, Day 1, Line 1: sparker"


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 |

#suwind j=3 |

sumute key=tracl ntaper=16 xmute=1 tmute=0.008 |
sufilter f=$f1 |
supef minlag=0.0024 maxlag=0.020 pnoise=0.001 mincorr=0.010 maxcor = 0.06 |
sugain tpow=1.0 pbal=1 |
#suacor ntout=500 sym=0 |
#supef minlag=0.0084 maxlag=0.0254 pnoise=0.001 mincorr=0.01 maxcor = 0.2 |
#sufilter f=$f1 |
#sugain tpow=-0.9 |

#migrate the data
suchw key1=gx key2=cdp |
sushw key=cdp a=1 b=1 |
sushw key=offset a=0 |
sustolt cdpmin=1 cdpmax=6224 tmig=0.0 vmig=1500 smig=0.6 fmax=1500 \
        vscale=0.9 lstaper=40 lbtaper=200 dxcdp=$dx |
#sushw key=cdp a=093252 b=1.6666666666 |
suchw key1=cdp key2=gx |

suchw key1=dt key2=dt b=0.75 >case_1_1_filt.su
#suchw key1=dt key2=dt b=1 >case_1_1_filt.su
#suximage perc=98 <case_1_1_filt.su &




# first panel of data  *************************************************
starttime=201100
endtime=202100
lefttic=2011
righttic=2021

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image1.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

# second panel of data  ******************
starttime=202100
endtime=203100
lefttic=2021
righttic=2031

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image2.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=203100
endtime=204100
lefttic=2031
righttic=2041

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image3.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=204100
endtime=205100
lefttic=2041
righttic=2051

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image4.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=205100
endtime=205959
lefttic=2051
righttic=2100

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image5.ps\
	perc=$perc\
	n2tic=2\
	d2=$time_interval\
	f2=$lefttic\
	x2beg=$lefttic\
	x2end=$righttic\
	nbpi=600\
	wbox=6.00\
	hbox=9.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=205100
endtime=210044
lefttic=2051
righttic=2101

suwind <case_1_1_filt.su key=cdp min=$starttime max=$endtime |
suchw key1=cdp key2=cdp d=100 |
##suflip flip=2 |
supsimage >c_1_1_image6.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