
#define DATA ../Data

DATA/%.PGR2 : DATA/greenpar2_% DATA/%.vel
	pushd DATA; cake %.PGR2

DATA/%.GR2 : DATA/greenpar2_% DATA/%.vel
	pushd DATA; cake %.GR2

DATA/%.PR44 : DATA/par_% DATA/%.vel
	 pushd DATA; cake %.PR44

DATA/%.R43 : DATA/par_% DATA/%.vel
	pushd DATA; cake %.R43

DATA/%.TT : DATA/ttpar_% DATA/%.vel
	 pushd DATA; cake %.TT

DATA/%.TT1 : DATA/%.TT DATA/ttpar_%
	 pushd DATA; cake %.TT1

DATA/%.vel : DATA/velpar_% 
	 pushd DATA; cake %.vel

$================================================

tsl-greens.v : ../Data/data4a.R43
	Window <../Data/data4a.R43 n2=1 f2=64 >junk.H
	Byte <junk.H pclip=100 | Ta2vplot head=/dev/null out=tsl-greens.v title="Time slice" PLOTARGS transp=0 wheretitle=t wherexlabel=b
	Rm junk.H

zsl-greens.v : ../Data/data4a.R43
	Window <../Data/data4a.R43 n3=1 f3=15 >junk.H
	Byte <junk.H pclip=100 | Ta2vplot head=/dev/null out=zsl-greens.v title="Depth slice" PLOTARGS transp=0 wheretitle=t wherexlabel=b
	Rm junk.H

$================================================

twophase.v: if not exist twophase.v
	Spike n1=128 nsp=2 k1=11,21 mag=1,.9 | Ftplot plot=upha >/dev/null title="Two events" out=tmp.v
	vppen big=n tmp.v vpstyle=n align=lb xsize=4 ysize=3 >twophase.v
	rm tmp.v

blob64.v:  DATA/data6.PGR2
	Window <../Data/data6.PGR2  n3=1 f3=91 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot >/dev/null out=blob64.v title="64 frequencies" label1="Angle" label2="Time"
	Rm tmp.H

blob32.v:  DATA/data6.PR44
	Window <../Data/data6.PR44  n3=1 f3=100 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot >/dev/null out=blob32.v title="32 frequencies" label1="Angle" label2="Time"
	Rm tmp.H

blob8.v:  DATA/data6a.PR44
	Window <../Data/data6a.PR44  n3=1 f3=100 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot >/dev/null out=blob8.v title="8 frequencies" label1="Angle" label2="Time"
	Rm tmp.H

blob4.v:  DATA/data6b.PR44
	Window <../Data/data6b.PR44  n3=1 f3=100 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot >/dev/null label1="Angle" label2="Time" out=blob4.v title="4 frequencies"
	Rm tmp.H

4blob.v : blob64.v blob32.v blob8.v blob4.v
	vppen vpstyle=n gridnum=2,2 blob64.v blob32.v blob8.v blob4.v | vppen xscale=.6 txsquare=y vpstyle=n >4blob.v

4blob.v3 : blob64.v blob32.v blob8.v blob4.v
	vppen vpstyle=n blob64.v blob32.v blob8.v blob4.v >4blob.v3

$================================================

blob-tt.v: DATA/data6.TT
	Window <DATA/data6.TT n3=1 f3=1 | Contour nc=12 >/dev/null out=blob-tt.v title="Maximum energy traveltimes" label1=Depth label2=Distance

blob-amp.v: DATA/data6.TT blob-vels.v
	Window <DATA/data6.TT n3=1 f3=0 | Byte clip=1.2 | Ta2vplot allpos=y >/dev/null out=tmp.v title="Maximum energy amplitudes" label1=Depth label2=Distance
	vppen erase=once vpstyle=n tmp.v blob-vels.v >blob-amp.v

blob-phase.v: DATA/data6.TT blob-vels.v
	Window <DATA/data6.TT n3=1 f3=2 | Byte clip=3.14 |  Ta2vplot >/dev/null out=tmp.v title="Maximum energy phase" label1=Depth label2=Distance	
	vppen erase=once vpstyle=n tmp.v blob-vels.v >blob-phase.v

blob-vels.v: DATA/data6.vel
	Contour <DATA/data6.vel nc=3 >/dev/null out=blob-vels.v fastplot=10 title="Blob, velocity field" label1=Depth label2=Distance


blob-ttamp.v: blob-tt.v blob-amp.v blob-phase.v 
	vppen vpstyle=n gridnum=3,1 blob-amp.v blob-tt.v blob-phase.v | vppen align=lb vpstyle=n  xscale=.8 yscale=1.2 txsquare=y  >blob-ttamp.v
	
blob-ttamp.v3:  blob-tt.v blob-amp.v blob-phase.v
	vppen vpstyle=n blob-tt.v blob-amp.v blob-phase.v  >blob-ttamp.v3


$================================================

blob-compt.v: all-t.v max-t.v
	vppen vpstyle=n gridnum=2,1 all-t.v max-t.v |  vppen align=lb vpstyle=n  xscale=.8  txsquare=y  >blob-compt.v

blob-compt.v3: all-t.v max-t.v
	vppen vpstyle=n all-t.v max-t.v  >blob-compt.v3

all-t.v: DATA/data6.GR2
	Window <DATA/data6.GR2  n2=1 f2=90 | Transp >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot label1=Depth label2=Distance title="Full outgoing wavefield" >/dev/null out=all-t.v
	Rm tmp.H

max-t.v: DATA/data6.GR2
	Window <DATA/data6.TT1  n2=1 f2=90 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot label2=Distance label1=Depth title="Maximum energy wavefield" >/dev/null out=max-t.v
	Rm tmp.H

blob-compz.v: all-z.v max-z.v
	vppen vpstyle=n gridnum=2,1 all-z.v max-z.v |  vppen align=lb vpstyle=n  xscale=.8  txsquare=y  >blob-compz.v

blob-compz.v3: all-z.v max-z.v
	vppen vpstyle=n  all-z.v max-z.v >blob-compz.v3

all-z.v: DATA/data6.GR2
	Window <DATA/data6.GR2  n3=1 f3=110  | Transp >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot label1=Time label2=Distance title="Full outgoing wavefield" >/dev/null out=all-z.v
	Rm tmp.H

max-z.v: DATA/data6.TT1
	Window <DATA/data6.TT1  n1=1 f1=110 >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot label2=Distance label1=Time title="Maximum energy wavefield" >/dev/null out=max-z.v
	Rm tmp.H

$==================================================================

blob-over.v: all-t.v DATA/data6.TT
	Window <DATA/data6.TT n3=1 f3=1 | Contour nc=12 c0=.02 dc=.02 >/dev/null out=tmp.v fastplot=10 plotfat=10 
	vppen erase=once vpstyle=n all-t.v tmp.v >blob-over.v
	rm tmp.v

clean&:
	-Rm tmp.H
	-rm tmp*.v
	-rm all-t.v all-z.v  blob-amp.v blob-tt.v 
	-rm blob-phase.v blob-vels.v max-t.v max-z.v 
	-rm blob4.v blob8.v blob32.v blob64.v

