

#define DATA ../Data

DATA/marm%.slow :  DATA/smoothpar_marm%  DATA/velpar_marm%
	pushd DATA; cake marm%.slow 

DATA/marm%.vel :  DATA/marm%.slow
	pushd DATA; cake marm%.vel 

DATA/marm%.ETT: DATA/ettpar_marm% DATA/marm%.slow
	pushd DATA; cake marm%.ETT

DATA/marm%.TT: DATA/ttpar_marm% DATA/marm%.slow
	pushd DATA; cake marm%.TT

/* can only make these on a CM5 */

DATA/marm%.GR2: DATA/greenpar2_marm% DATA/marm%.vel
	pushd DATA; cake marm%.GR2

DATA/marm%.TT1: DATA/marm%.TT DATA/ttpar_marm%
	pushd DATA; cake marm%.TT1



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

marmslow.v: DATA/truemarm.slow
	Window <DATA/truemarm.slow j1=4 j2=4 >tmp.H
	Byte <tmp.H allpos=y | Ta2vplot label1=Depth label2=Distance title="Marmousi Slowness Model"  >/dev/null out=tmp.v
	vppen tmp.v scale=.5 vpstyle=n rotate=-90 | vppen vpstyle=n align=lb > marmslow.v
	

marm%-vel.v : DATA/marm%.vel
	Window <DATA/marm%.vel j1=3 j2=3 >tmp.H
	Byte <tmp.H allpos=y | Ta2vplot fastplot=15 >/dev/null out=marm%-vel.v

$=================================================
$ Contour plots of traveltimes over velocity field
$=================================================

#define CONTARGS nc=20 c0=.1 dc=.1 plotfat=4 label1=Depth label2=Distance

marm%-ett.v:  DATA/marm%.ETT marm%-vel.v
	Contour <DATA/marm%.ETT CONTARGS title="Eikonal FD Traveltimes" >/dev/null out=tmp.v
	vppen erase=once vpstyle=n marm%-vel.v tmp.v >marm%-ett.v

marm%-mtt.v:  DATA/marm%.TT marm%-vel.v
	Window < DATA/marm%.TT n3=1 | Contour CONTARGS title="Band-Limited Traveltimes" >/dev/null out=tmp.v
	vppen erase=once vpstyle=n marm%-vel.v tmp.v >marm%-mtt.v

marm%-ptt.v:  ../Parax/parax-map3.H marm%-vel.v
	Window < ../Parax/parax-map3.H n3=1 f3=0 | Contour CONTARGS title="Paraxial Ray-Tracing Traveltimes" >/dev/null out=tmp.v
	vppen erase=once vpstyle=n marm%-vel.v tmp.v >marm%-ptt.v

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

marm%-tt.v: marm%-ett.v marm%-mtt.v marm%-ptt.v
	vppen vpstyle=n gridnum=1,3  marm%-ett.v marm%-ptt.v marm%-mtt.v | vppen vpstyle=n yscale=.6 >marm%-tt.v


marm%-tt.v3: marm%-ett.v marm%-mtt.v marm%-ptt.v
	vppen vpstyle=n marm%-ett.v marm%-ptt.v marm%-mtt.v  >marm%-tt.v3

$================================================
$ Time slices of greens func. 
$================================================

two-way-%.v1 : ../FDModel/two-way-%.HH
	Window <../FDModel/two-way-%.HH j1=4 j2=4 >tmp.H
	Byte <tmp.H | Ta2vplot fastplot=15 >/dev/null out=two-way-%.v1
	Rm tmp.H

two-way-%.v : ../FDModel/two-way-%.HH
	Window <../FDModel/two-way-%.HH  j1=2 j2=2 >tmp.H
	Byte <tmp.H pclip=99 | Ta2vplot label1=Depth label2=Distance title="Two-way wave equation"  >/dev/null out=two-way-%.v
	Rm tmp.H

one-way-%.v : DATA/marm2.GR2
	Window <DATA/marm2.GR2 n3=375 n2=1 min2=% >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot transp=n label1=Distance label2=Depth title="Outgoing wave-equation" >/dev/null out=one-way-%.v
	Rm tmp.H

param-rec-%.v : DATA/marm2.TT1
	Window <DATA/marm2.TT1 n2=1 min2=% >tmp.H
	Byte <tmp.H pclip=100 | Ta2vplot label2=Distance label1=Depth title="Parametric reconstruction" >/dev/null out=param-rec-%.v
	Rm tmp.H

green-comp-%.v: two-way-%.v one-way-%.v param-rec-%.v 
	vppen vpstyle=n gridnum=1,3  two-way-%.v one-way-%.v param-rec-%.v | vppen vpstyle=n yscale=.4 >green-comp-%.v

green-comp-%.v3: two-way-%.v one-way-%.v param-rec-%.v 
	vppen vpstyle=n  two-way-%.v one-way-%.v param-rec-%.v >green-comp-%.v3

$================================================
$ Time slices of greens func. with traveltime overlays at .9 sec
$================================================

#define CONTARG8a plotfat=12 dash=2 nc=1 c0=.91 plotcol=7 fastplot=15
#define CONTARG8 label1=Depth label2=Distance  plotcol=0 plotfat=6 dash=2 nc=1 c0=.91

eikonal-tov-.9.v : two-way-.9.v1 DATA/marm2.ETT
	Contour <DATA/marm2.ETT  CONTARG8a >/dev/null out=tmp2.v
	Contour <DATA/marm2.ETT  title="Time slice + Eikonal Traveltime" CONTARG8 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-.9.v1 tmp2.v tmp3.v >eikonal-tov-.9.v
 
paraxial-tov-.9.v : ../Parax/parax-map3.H two-way-.9.v1 
	Window n3=1 <../Parax/parax-map3.H | Contour CONTARG8a >/dev/null out=tmp2.v
	Window n3=1 <../Parax/parax-map3.H | Contour title="Time slice + Paraxial Traveltime" CONTARG8 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-.9.v1 tmp2.v tmp3.v >paraxial-tov-.9.v
 
bandlim-tov-.9.v : two-way-.9.v1 DATA/marm2.TT
	Window <DATA/marm2.TT n3=1 f3=0| Contour CONTARG8a >/dev/null out=tmp2.v
	Window <DATA/marm2.TT n3=1 f3=0| Contour title="Time slice + Band-limited Traveltime" CONTARG8 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-.9.v1 tmp2.v tmp3.v >bandlim-tov-.9.v

$================================================
$ Time slices of greens func. with traveltime overlays at 1.1 sec
$================================================

#define CONTARG10a plotfat=12 dash=2 nc=1 c0=1.115 plotcol=7 fastplot=15
#define CONTARG10 label1=Depth label2=Distance plotcol=0 plotfat=6 dash=2 nc=1 c0=1.115


eikonal-tov-1.1.v : two-way-1.1.v1 DATA/marm2.ETT
	Contour <DATA/marm2.ETT  CONTARG10a >/dev/null out=tmp2.v
	Contour <DATA/marm2.ETT  title="Time slice + Eikonal Traveltime" CONTARG10 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-1.1.v1 tmp2.v tmp3.v >eikonal-tov-1.1.v
 
paraxial-tov-1.1.v : two-way-1.1.v1 ../Parax/parax-map3.H
	Window n3=1 <../Parax/parax-map3.H | Contour CONTARG10a >/dev/null out=tmp2.v
	Window n3=1 <../Parax/parax-map3.H | Contour title="Time slice + Paraxial Traveltime" CONTARG10 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-1.1.v1 tmp2.v tmp3.v >paraxial-tov-1.1.v
 
bandlim-tov-1.1.v : two-way-1.1.v1 DATA/marm2.TT
	Window <DATA/marm2.TT n3=1 f3=0| Contour CONTARG10a >/dev/null out=tmp2.v
	Window <DATA/marm2.TT n3=1 f3=0| Contour title="Time slice + Band-limited Traveltime" CONTARG10 >/dev/null out=tmp3.v
	vppen erase=once vpstyle=n two-way-1.1.v1 tmp2.v tmp3.v >bandlim-tov-1.1.v

$================================================
$  composites of time slice overlays
$================================================

all-tov%.v : eikonal-tov%.v paraxial-tov%.v bandlim-tov%.v
	vppen vpstyle=n gridnum=1,3   eikonal-tov%.v paraxial-tov%.v bandlim-tov%.v | vppen vpstyle=n yscale=.4 >all-tov%.v

all-tov%.v3 : eikonal-tov%.v paraxial-tov%.v bandlim-tov%.v
	vppen vpstyle=n eikonal-tov%.v paraxial-tov%.v bandlim-tov%.v >all-tov%.v3 


$================================
$ Amplitude and phase plots
$================================

band-amp.v:  DATA/marm2.TT 
	Window < DATA/marm2.TT n3=1 f3=1 >tmp.H
	Byte <tmp.H allpos=y  | Ta2vplot label1=Depth labe2=Distance title="Band-limited, amplitude" >/dev/null out=band-amp.v

band-amp2.v:  DATA/marm2.TT 
	Window < DATA/marm2.TT j1=2 j2=3 n3=1 f3=1 >tmp.H
	Byte <tmp.H allpos=y  | Ta2vplot label1=Depth labe2=Distance title="Band-limited, amplitude" >/dev/null out=band-amp2.v

parax-amp.v:  ../Parax/parax-map3.H
	Window < ../Parax/parax-map3.H n3=1 f3=1 >tmp.H
	Byte <tmp.H allpos=y  | Ta2vplot label1=Depth labe2=Distance title="Paraxial ray-tracing, amplitude" >/dev/null out=parax-amp.v

comp-amp.v: band-amp.v parax-amp.v
	vppen band-amp.v parax-amp.v vpstyle=n gridnum=1,2 | vppen vpstyle=n yscale=.7 >comp-amp.v

comp-amp.v3: band-amp.v parax-amp.v
	vppen band-amp.v parax-amp.v vpstyle=n >comp-amp.v3


band-phase.v:  DATA/marm2.TT
	Window < DATA/marm2.TT n3=1 f3=2 >tmp.H
	Byte <tmp.H  clip=3.2 | Ta2vplot label1=Depth labe2=Distance title="Band-limited, phase" >/dev/null out=band-phase.v

parax-phase.v:   ../Parax/parax-map3.H
	Window < ../Parax/parax-map3.H n3=1 f3=2 >tmp.H
	Byte <tmp.H   | Ta2vplot label1=Depth labe2=Distance title="Paraxial ray-tracing, phase" >/dev/null out=parax-phase.v

comp-phase.v: band-phase.v parax-phase.v
	vppen band-phase.v parax-phase.v vpstyle=n gridnum=1,2 | vppen vpstyle=n yscale=.7 >comp-phase.v

comp-phase.v3: band-phase.v parax-phase.v
	vppen band-phase.v parax-phase.v vpstyle=n >comp-phase.v3



clean&:
	-rm tmp*.v
	-Rm tmp*.H
	-rm Wmod_win_band.v band-amp* band-phase* parax-amp* parax-phase*
	-rm bandlim-tov* eikonal-tov* two-way* one-way* param-rec*
	-rm marm2-vel.v marm2-ett.v marm2-mtt.v marm2-ptt.v
	-rm paraxial-tov*.v
	
