
all&: polar-greens.v rect-greens.v tsl-greens.v zsl-greens.v const-mod.v adj-mig.v ata-mig.v const-mod.v3 adj-mig.v3 ata-mig.v3 tadj-mig.v tadj-mig.v3 const-modt.v const-modt.v3

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

../Data/%.R41: ../Data/par_% ../Data/%.vel
        pushd ../Data; cake  %.R41

../Data/%.PR41: ../Data/par_% ../Data/%.vel
        pushd ../Data; cake  %.PR41

../Data/%.%1.mod: ../Data/%.GR ../Data/%.scat ../Data/modpar_%1_%
        pushd ../Data; cake  %.%1.mod

../Data/%.%1.modt: ../Data/%.GR ../Data/%.scat ../Data/modpar_%1_%
        pushd ../Data; cake  %.%1.modt

../Data/%.%1.modmig: ../Data/%.GR ../Data/%.%1.mod ../Data/migpar_%1_%
        pushd ../Data; cake  %.%1.modmig

../Data/%.%1.modadj: ../Data/%.GR ../Data/%.%1.mod ../Data/migpar_%1_%
        pushd ../Data; cake  %.%1.modadj

../Data/%.%1.modtadj: ../Data/%.GR ../Data/%.%1.modt ../Data/migpar_%1_%
        pushd ../Data; cake  %.%1.modtadj

../Data/%.scat: ../Data/scatpar_% if exist ../Data/scatpar_%
        pushd ../Data; cake  %.scat

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

../Data/%.GR: ../Data/greenpar_%
	 pushd ../Data; cake ../Data/%.GR

/*../Data/% :  if not {{ pushd ../Data ; cake -q % }}*/

/*../Data/% :  if not ok ../Data/%*/
/*        pushd ../Data; cake  %*/


#define PLOTARGS labelsz=12 titlesz=15
#define TWOPLOTS xscale=.52 txsquare=y 
#define ONEPLOTS xscale=.52 txsquare=y 


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

rect-greens.v: ../Data/data4a.R41 
	Window < ../Data/data4a.R41 n3=1 f3=2 >junk.H
	Byte pclip=100 <junk.H | Ta2vplot title=Amplitude head=/dev/null out=junk1.v PLOTARGS
	Window < ../Data/data4a.R41 n3=1 f3=3 >junk.H
	Byte <junk.H | Ta2vplot title=Phase head=/dev/null out=junk2a.v PLOTARGS
	Contour <junk.H nc=10 fastplot=20 head=/dev/null out=junk2c.v transp=y
	vp_Overlay junk2a.v junk2c.v >junk2.v
	vppen gridnum=2,1 junk1.v junk2.v vpstyle=n | vppen yscale=2 vpstyle=n >rect-greens.v
	/bin/rm junk1.v junk2.v junk2a.v junk2c.v
	Rm junk.H

polar-greens.v: ../Data/data4a.PR41
	Window < ../Data/data4a.PR41 n3=1 f3=2 >junk.H
	Byte pclip=100  <junk.H | Ta2vplot title=Amplitude head=/dev/null out=junk1.v PLOTARGS
	Window < ../Data/data4a.PR41 n3=1 f3=3 >junk.H
	Byte  <junk.H | Ta2vplot title=Phase head=/dev/null out=junk2a.v PLOTARGS
	Contour <junk.H nc=10 fastplot=20 head=/dev/null out=junk2c.v transp=y
	vp_Overlay junk2a.v junk2c.v >junk2.v
	vppen gridnum=2,1 junk1.v junk2.v vpstyle=n | vppen yscale=2 vpstyle=n >polar-greens.v
	/bin/rm junk1.v junk2.v junk2a.v junk2c.v
	Rm junk.H

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

data4.scat.v: ../Data/data4.scat scatter_label
	Byte <../Data/data4.scat pclip=100 | Ta2vplot title="Scatterer field" head=/dev/null out=stdout PLOTARGS transp=0  | vp_annotate  text=scatter_label vpstyle=n batch=y >data4.scat.v

data4.1.mod.v: ../Data/data4.1.mod
	Byte  <../Data/data4.1.mod pclip=100 | Ta2vplot title="Offset=0 modeling" head=/dev/null out=data4.1.mod.v PLOTARGS transp=n 

data4.2.mod.v: ../Data/data4.2.mod
	Byte  <../Data/data4.2.mod pclip=100 | Ta2vplot title="Offset=375m modeling" head=/dev/null out=data4.2.mod.v PLOTARGS transp=n

data4.3.mod.v: ../Data/data4.3.mod
	Byte  <../Data/data4.3.mod pclip=100 | Ta2vplot title="Offset=750m modeling" head=/dev/null out=data4.3.mod.v PLOTARGS transp=n

const-mod.v3: data4.scat.v data4.1.mod.v data4.2.mod.v data4.3.mod.v
	vppen vpstyle=n data4.scat.v data4.1.mod.v data4.2.mod.v data4.3.mod.v >const-mod.v3

const-mod.v: data4.scat.v data4.1.mod.v data4.2.mod.v data4.3.mod.v
	vppen gridnum=2,2 vpstyle=n data4.scat.v data4.1.mod.v data4.2.mod.v data4.3.mod.v | vppen vpstyle=n xscale=.5 >const-mod.v

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

data4.1.modt.v: ../Data/data4.1.modt
	Byte  <../Data/data4.1.modt pclip=100 | Ta2vplot title="Offset=0 modeling + Tpow" head=/dev/null out=data4.1.modt.v PLOTARGS transp=n 

data4.2.modt.v: ../Data/data4.2.modt
	Byte  <../Data/data4.2.modt pclip=100 | Ta2vplot title="Offset=375m modeling + Tpow" head=/dev/null out=data4.2.modt.v PLOTARGS transp=n

data4.3.modt.v: ../Data/data4.3.modt
	Byte  <../Data/data4.3.modt pclip=100 | Ta2vplot title="Offset=750m modeling + Tpow" head=/dev/null out=data4.3.modt.v PLOTARGS transp=n

const-modt.v3: data4.scat.v data4.1.modt.v data4.2.modt.v data4.3.modt.v
	vppen vpstyle=n data4.scat.v data4.1.modt.v data4.2.modt.v data4.3.modt.v >const-modt.v3

const-modt.v: data4.scat.v data4.1.modt.v data4.2.modt.v data4.3.modt.v
	vppen gridnum=2,2 vpstyle=n data4.scat.v data4.1.modt.v data4.2.modt.v data4.3.modt.v | vppen vpstyle=n xscale=.5 >const-modt.v

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

data4.1.modmig.v: ../Data/data4.1.modmig scatter_label
	Byte  <../Data/data4.1.modmig pclip=100 gpow=.5 | Ta2vplot title="Offset=0 migration" head=/dev/null out=stdout PLOTARGS transp=n   | vp_annotate text=scatter_label vpstyle=n batch=y >data4.1.modmig.v

data4.2.modmig.v: ../Data/data4.2.modmig scatter_label
	Byte  <../Data/data4.2.modmig pclip=100 gpow=.5 | Ta2vplot title="Offset=375m migration" head=/dev/null out=stdout PLOTARGS transp=n  | vp_annotate text=scatter_label vpstyle=n batch=y >data4.2.modmig.v

data4.3.modmig.v: ../Data/data4.3.modmig scatter_label
	Byte  <../Data/data4.3.modmig pclip=100 gpow=.5 | Ta2vplot title="Offset=750m migration" head=/dev/null out=stdout PLOTARGS transp=n  | vp_annotate text=scatter_label vpstyle=n batch=y >data4.3.modmig.v

ata-mig.v3: data4.scat.v data4.1.modmig.v data4.2.modmig.v data4.3.modmig.v
	vppen vpstyle=n data4.scat.v data4.1.modmig.v data4.2.modmig.v data4.3.modmig.v >ata-mig.v3

ata-mig.v: data4.scat.v data4.1.modmig.v data4.2.modmig.v data4.3.modmig.v
	vppen gridnum=2,2 vpstyle=n data4.scat.v data4.1.modmig.v data4.2.modmig.v data4.3.modmig.v | vppen vpstyle=n xscale=.5 >ata-mig.v

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

data4.1.modtadj.v: ../Data/data4.1.modtadj scatter_label
	Byte  <../Data/data4.1.modtadj pclip=100 gpow=.5 | Ta2vplot title="Offset=0, Tpow + adjoint" head=/dev/null out=stdout PLOTARGS transp=n  | vp_annotate text=scatter_label vpstyle=n batch=y >data4.1.modtadj.v

data4.2.modtadj.v: ../Data/data4.2.modtadj scatter_label
	Byte  <../Data/data4.2.modtadj pclip=100 gpow=.5 | Ta2vplot title="Offset=375m, Tpow + adjoint" head=/dev/null out=stdout PLOTARGS transp=n | vp_annotate text=scatter_label vpstyle=n batch=y >data4.2.modtadj.v

data4.3.modtadj.v: ../Data/data4.3.modtadj scatter_label
	Byte  <../Data/data4.3.modtadj pclip=100 gpow=.5 | Ta2vplot title="Offset=750m, Tpow + adjoint" head=/dev/null out=stdout PLOTARGS transp=n | vp_annotate text=scatter_label vpstyle=n batch=y >data4.3.modtadj.v

tadj-mig.v3: data4.scat.v data4.1.modtadj.v data4.2.modtadj.v data4.3.modtadj.v
	vppen vpstyle=n data4.scat.v data4.1.modtadj.v data4.2.modtadj.v data4.3.modtadj.v >tadj-mig.v3

tadj-mig.v: data4.scat.v data4.1.modtadj.v data4.2.modtadj.v data4.3.modtadj.v
	vppen gridnum=2,2 vpstyle=n data4.scat.v data4.1.modtadj.v data4.2.modtadj.v data4.3.modtadj.v | vppen vpstyle=n xscale=.5 >tadj-mig.v

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

data4.1.modadj.v: ../Data/data4.1.modadj scatter_label
	Byte  <../Data/data4.1.modadj pclip=100 gpow=.5 | Ta2vplot title="Offset=0 adjoint modeling" head=/dev/null out=stdout PLOTARGS transp=n  | vp_annotate text=scatter_label vpstyle=n batch=y >data4.1.modadj.v

data4.2.modadj.v: ../Data/data4.2.modadj scatter_label
	Byte  <../Data/data4.2.modadj pclip=100 gpow=.5 | Ta2vplot title="Offset=375m adjoint modeling" head=/dev/null out=stdout PLOTARGS transp=n | vp_annotate text=scatter_label vpstyle=n batch=y >data4.2.modadj.v

data4.3.modadj.v: ../Data/data4.3.modadj scatter_label
	Byte  <../Data/data4.3.modadj pclip=100 gpow=.5 | Ta2vplot title="Offset=750m adjoint modeling" head=/dev/null out=stdout PLOTARGS transp=n | vp_annotate text=scatter_label vpstyle=n batch=y >data4.3.modadj.v

adj-mig.v3: data4.scat.v data4.1.modadj.v data4.2.modadj.v data4.3.modadj.v
	vppen vpstyle=n data4.scat.v data4.1.modadj.v data4.2.modadj.v data4.3.modadj.v >adj-mig.v3

adj-mig.v: data4.scat.v data4.1.modadj.v data4.2.modadj.v data4.3.modadj.v
	vppen gridnum=2,2 vpstyle=n data4.scat.v data4.1.modadj.v data4.2.modadj.v data4.3.modadj.v | vppen vpstyle=n xscale=.5 >adj-mig.v

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

data4.scat.1: ../Data/data4.scat
	Window <../Data/data4.scat n1=1 f1=127 >data4.scat.1

data4.1.modadj.1: ../Data/data4.1.modadj
	Window  <../Data/data4.1.modadj n1=1 f1=126 >data4.1.modadj.1

data4.2.modadj.1: ../Data/data4.2.modadj
	Window  <../Data/data4.2.modadj n1=1 f1=126 >data4.2.modadj.1

data4.3.modadj.1: ../Data/data4.3.modadj
	Window  <../Data/data4.3.modadj n1=1 f1=126 >data4.3.modadj.1

adj-mig2.v: data4.scat.1 data4.1.modadj.1 data4.2.modadj.1 data4.3.modadj.1
	Merge data4.scat.1 data4.1.modadj.1 data4.2.modadj.1 data4.3.modadj.1 axis=2 space=no | Dots >/dev/null out=adj-mig2.v labels=Input:near:mid:far title="Adjoint modeling"

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

data4.1.modtadj.1: ../Data/data4.1.modtadj
	Window  <../Data/data4.1.modtadj n1=1 f1=126 >data4.1.modtadj.1

data4.2.modtadj.1: ../Data/data4.2.modtadj
	Window  <../Data/data4.2.modtadj n1=1 f1=126 >data4.2.modtadj.1

data4.3.modtadj.1: ../Data/data4.3.modtadj
	Window  <../Data/data4.3.modtadj n1=1 f1=126 >data4.3.modtadj.1

tadj-mig2.v: data4.scat.1 data4.1.modtadj.1 data4.2.modtadj.1 data4.3.modtadj.1
	Merge data4.scat.1 data4.1.modtadj.1 data4.2.modtadj.1 data4.3.modtadj.1 axis=2 space=no | Dots >/dev/null out=tadj-mig2.v labels=Input:near:mid:far title="Tpow + Adjoint modeling"

$=============================
data4.1.modmig.1: ../Data/data4.1.modmig
	Window  <../Data/data4.1.modmig n1=1 f1=126 >data4.1.modmig.1

data4.2.modmig.1: ../Data/data4.2.modmig
	Window  <../Data/data4.2.modmig n1=1 f1=126 >data4.2.modmig.1

data4.3.modmig.1: ../Data/data4.3.modmig
	Window  <../Data/data4.3.modmig n1=1 f1=126 >data4.3.modmig.1

ata-mig2.v: data4.scat.1 data4.1.modmig.1 data4.2.modmig.1 data4.3.modmig.1
	Merge data4.scat.1 data4.1.modmig.1 data4.2.modmig.1 data4.3.modmig.1 axis=2 space=no | Dots >/dev/null out=ata-mig2.v labels=Input:near:mid:far title="Migration/Inversion"

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



clean&:
	-/bin/rm junk*v
	-Rm junk.H
	-Rm data4*.1
	-rm data4.*.mod.v
	-rm data4.*.modt.v
	-rm data4.*.modadj.v
	-rm data4.*.modtadj.v
	-rm data4.*.modmig.v
	-rm data4.scat.v

