### Makefile for Pyramid domain Data interpolation
###  2008
###
###
# include SEP rules
#include /net/server/xukai/huo/xukai1/C/mpi.Makefile

#include ${SEPINC}/SEP.top
include /net/server/wrk/sep143/xukai2/mpi.Makefile
datadir= ~xukai/huo/Vmodel
BINDIR1= ./bin
GR = Grey  > /dev/null
GRAPH = Graph  > /dev/null



nx=300
ny=300
nz=300
dd=20
hnx=250
hny=150
hnz=60
nxtape=50
nytape=50
nztape=50
nt=1000
dt=0.004
#exec: ${SRCDIR}/XC3dWaveprop.c ${SRCDIR}/XC3dAcouMod.c
#	make ${BINDIR}/XC3dAcouModMain.x
######################################################################################################################
exec: ${OBJDIR}/XC3dWaveprop.o ${OBJDIR}/XC3dAcouMod.o ${OBJDIR}/4t10spacewithbnd1opt2.o
	make ${OBJDIR}/XC3dAcouModMain.o
	make ${BINDIR}/XC3dAcouModMain.x
	make ${BINDIR}/XC3dRndbndVel.x
	make ${BINDIR}/XC3dRndbndVelGrain.x
	make ${BINDIR}/XC3dRndbndVelGrainNew.x
	make ${BINDIR}/XC3dMakesrc.x
	make ${BINDIR}/XFt.x



######################################################################################################################
####### velocity patameters 

######################################################################################################################
######################################################################################################################
######################################################################################################################
## stability limit for 2nd order in time stencil 0.44

vel1.H:
	Vel n1=${nz} n2=${nx} n3=${ny} d1=${dd} d2=${dd} d3=${dd} o1=0 o2=0 o3=0 vc=3000 >$@

src242.H: exec
	${BINDIR}/XC3dMakesrc.x nt=${nt} nx=${nx} ny=${ny} ot=0 dt=${dt} ox=0 oy=0 dx=${dd} dy=${dd} nsx=1 osx=3000 dsx=60 nsy=1 osy=3000 dsy=60 type=2. tdelay=-0.2 fpeak=25 >$@

src6.H: exec
	${BINDIR}/XC3dMakesrc.x nt=${nt} nx=${nx} ny=${ny} ot=0 dt=${dt} ox=0 oy=0 dx=${dd} dy=${dd} nsx=1 osx=3000 dsx=60 nsy=1 osy=3000 dsy=60 type=1. tdelay=-0.2 fpeak=6 >$@

sdep.H:vel1.H
	<vel1.H Window n3=1 >tmp1.H
	Math file1=tmp1.H exp="file1*0+2000" >$@

record1.src%.H : src%.H sdep.H vel1.H
	date        
	<src$*.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record1.src$*.H wvmv=wvmv1.src$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date
   
######################################################################################################################
wvmv%.1t.H:wvmv%.H
	<$< Window3d n3=1 f3=150 n4=1 f4=22 >$@


wvmv1.sum6.H:wvmv1.src6.1t.H
	<$< Scale rscale=16 >$@

wvmv1.sum242.H:wvmv1.src242.1t.H
	<$< Scale rscale=16 >$@
######################################################################################################################

%.kabs.H:%.H
	<$< Pad n1=512 n2=512 | Rtoc | Ft3d center1=1 center2=1 sign1=1 sign2=1 | Cabs | Window f1=181 n1=150 f2=181 n2=150 >$@
######################################################################################################################
######################################################################################################################
###from Makefile2
vel1.rndbnd2.%.H: exec
	<vel1.H ${BINDIR}/XC3dRndbndVel.x nxtape=${nxtape} nytape=${nytape} nztape=${nztape} sigma=60 alpha_mean=0.01 alpha_sigma=0.1 minV=1500 maxV=4200 vel=$@ >/dev/null alpha_maxV=0.0

vel1.rndbnd51.%.H: exec
	<vel1.H ${BINDIR}/XC3dRndbndVel.x nxtape=${nxtape} nytape=${nytape} nztape=${nztape} sigma=60 alpha_mean=0.01 alpha_sigma=0.1 minV=1500 maxV=4200 vel=$@ >/dev/null alpha_maxV=0.0 j_rndbnd=10

vel1.rndbnd52.%.H: exec
	<vel1.H ${BINDIR}/XC3dRndbndVelGrainNew.x nxtape=${nxtape} nytape=${nytape} nztape=${nztape} sigma=60 alpha_mean=0.01 alpha_sigma=0.1 minV=1500 maxV=4200 vel=$@ >/dev/null alpha_maxV=0.0 j_rndbnd=10

vel1.rndbnd53.%.H: exec
	<vel1.H ${BINDIR}/XC3dRndbndVelGrainNew1.x nxtape=${nxtape} nytape=${nytape} nztape=${nztape} sigma=60 alpha_mean=0.01 alpha_sigma=0.1 minV=800 maxV=3300 vel=$@ >/dev/null alpha_maxV=0.0 j_rndbnd=10

record2.src6.%.H : src6.H sdep.H vel1.rndbnd2.%.H
	date
	<src6.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd2.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record2.src6.$*.H wvmv=wvmv2.src6.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

record51.src6.%.H : src6.H sdep.H vel1.rndbnd51.%.H
	date
	<src6.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd51.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record51.src6.$*.H wvmv=wvmv51.src6.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

record52.src6.%.H : src6.H sdep.H vel1.rndbnd52.%.H
	date
	<src6.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd52.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record52.src6.$*.H wvmv=wvmv52.src6.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

record53.src6.%.H : src6.H sdep.H vel1.rndbnd53.%.H
	date
	<src6.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd53.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record53.src6.$*.H wvmv=wvmv53.src6.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date
######################################################################################################################
######################################################################################################################
record2.src242.%.H : src242.H sdep.H vel1.rndbnd2.%.H
	date
	<src242.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd2.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record2.src242.$*.H wvmv=wvmv2.src242.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

record51.src242.%.H : src242.H sdep.H vel1.rndbnd51.%.H
	date
	<src242.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd51.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record51.src242.$*.H wvmv=wvmv51.src242.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

record52.src242.%.H : src242.H sdep.H vel1.rndbnd52.%.H
	date
	<src242.H ${BINDIR}/XC3dAcouModMain.x vel=vel1.rndbnd52.$*.H nxtape=${nxtape} nytape=${nytape} nztape=${nztape} data=record52.src242.$*.H wvmv=wvmv52.src242.$*.H needwvmv=1 jtmv=40 rndbnd=2 sdep=sdep.H
	date

######################################################################################################################a
rndbnd%:
	make vel1.rndbnd$*.1.H
	make vel1.rndbnd$*.2.H
	make vel1.rndbnd$*.3.H
	make vel1.rndbnd$*.4.H
	make vel1.rndbnd$*.5.H
	make vel1.rndbnd$*.6.H
	make vel1.rndbnd$*.7.H
	make vel1.rndbnd$*.8.H
	make vel1.rndbnd$*.9.H
	make vel1.rndbnd$*.10.H
	make vel1.rndbnd$*.11.H
	make vel1.rndbnd$*.12.H
	make vel1.rndbnd$*.13.H
	make vel1.rndbnd$*.14.H
	make vel1.rndbnd$*.15.H
	make vel1.rndbnd$*.16.H

6rcd%:
	make record$*.src6.1.H
	make record$*.src6.2.H
	make record$*.src6.3.H
	make record$*.src6.4.H
	make record$*.src6.5.H
	make record$*.src6.6.H
	make record$*.src6.7.H
	make record$*.src6.8.H
	make record$*.src6.9.H
	make record$*.src6.10.H
	make record$*.src6.11.H
	make record$*.src6.12.H
	make record$*.src6.13.H
	make record$*.src6.14.H
	make record$*.src6.15.H
	make record$*.src6.16.H

242rcd%:
	make record$*.src242.1.H
	make record$*.src242.2.H
	make record$*.src242.3.H
	make record$*.src242.4.H
	make record$*.src242.5.H
	make record$*.src242.6.H
	make record$*.src242.7.H
	make record$*.src242.8.H
	make record$*.src242.9.H
	make record$*.src242.10.H
	make record$*.src242.11.H
	make record$*.src242.12.H
	make record$*.src242.13.H
	make record$*.src242.14.H
	make record$*.src242.15.H
	make record$*.src242.16.H

6wvmv%:
	make wvmv$*.src6.1.1t.H
	make wvmv$*.src6.2.1t.H
	make wvmv$*.src6.3.1t.H
	make wvmv$*.src6.4.1t.H
	make wvmv$*.src6.5.1t.H
	make wvmv$*.src6.6.1t.H
	make wvmv$*.src6.7.1t.H
	make wvmv$*.src6.8.1t.H
	make wvmv$*.src6.9.1t.H
	make wvmv$*.src6.10.1t.H
	make wvmv$*.src6.11.1t.H
	make wvmv$*.src6.12.1t.H
	make wvmv$*.src6.13.1t.H
	make wvmv$*.src6.14.1t.H
	make wvmv$*.src6.15.1t.H
	make wvmv$*.src6.16.1t.H

242wvmv%:
	make wvmv$*.src242.1.1t.H
	make wvmv$*.src242.2.1t.H
	make wvmv$*.src242.3.1t.H
	make wvmv$*.src242.4.1t.H
	make wvmv$*.src242.5.1t.H
	make wvmv$*.src242.6.1t.H
	make wvmv$*.src242.7.1t.H
	make wvmv$*.src242.8.1t.H
	make wvmv$*.src242.9.1t.H
	make wvmv$*.src242.10.1t.H
	make wvmv$*.src242.11.1t.H
	make wvmv$*.src242.12.1t.H
	make wvmv$*.src242.13.1t.H
	make wvmv$*.src242.14.1t.H
	make wvmv$*.src242.15.1t.H
	make wvmv$*.src242.16.1t.H


runmultiple1:
	make record1.src6.H
	make wvmv1.src6.1t.H
	make wvmv1.src6.1t.kabs.H
	make wvmv1.sum6.H
	make wvmv1.sum6.kabs.H
	make record1.src242.H
	make wvmv1.src242.1t.H
	make wvmv1.sum242.H
	make wvmv1.sum242.kabs.H

runmultiple2:
	make rndbnd2
	make 6rcd2 >& log62
	make 242rcd2 >& log2422
	make 6wvmv2 >& logmv62
	make 242wvmv2 >& logmv2422
	make wvmv2.sum6.H
	make wvmv2.sum6.kabs.H
	make wvmv2.src6.1.1t.kabs.H
	make wvmv2.sum242.H

runmultiple51:
	make rndbnd51
	make 6rcd51 >& log651
	make 242rcd51 >& log24251
	make 6wvmv51 >& logmv651
	make 242wvmv51 >& logmv24251
	make wvmv51.sum6.H
	make wvmv51.src6.1.1t.kabs.H
	make wvmv51.sum6.kabs.H
	make wvmv51.sum242.H

runmultiple52:
	make rndbnd52
	make 6rcd52 >& log652
	make 242rcd52 >& log24252
	make 6wvmv52 >& logmv652
	make 242wvmv52 >& logmv24252
	make wvmv52.sum6.H
	make wvmv52.src6.1.1t.kabs.H
	make wvmv52.sum6.kabs.H
	make wvmv52.sum242.H

runmultiple53:
	make rndbnd53
	make 6rcd53 >& log653
	make 6wvmv53 >& logmv653
	make wvmv53.sum6.H
	make wvmv53.sum6.kabs.H


wvmv51.sum%.H:
	Add wvmv51.src$*.1.1t.H wvmv51.src$*.2.1t.H wvmv51.src$*.3.1t.H wvmv51.src$*.4.1t.H >tmp1.H
	Add wvmv51.src$*.5.1t.H wvmv51.src$*.6.1t.H wvmv51.src$*.7.1t.H wvmv51.src$*.8.1t.H >tmp2.H
	Add wvmv51.src$*.9.1t.H wvmv51.src$*.10.1t.H wvmv51.src$*.11.1t.H wvmv51.src$*.12.1t.H >tmp3.H
	Add wvmv51.src$*.13.1t.H wvmv51.src$*.14.1t.H wvmv51.src$*.15.1t.H wvmv51.src$*.16.1t.H >tmp4.H
	Add tmp1.H tmp2.H tmp3.H tmp4.H >$@

wvmv52.sum%.H:
	Add wvmv52.src$*.1.1t.H wvmv52.src$*.2.1t.H wvmv52.src$*.3.1t.H wvmv52.src$*.4.1t.H >tmp1.H
	Add wvmv52.src$*.5.1t.H wvmv52.src$*.6.1t.H wvmv52.src$*.7.1t.H wvmv52.src$*.8.1t.H >tmp2.H
	Add wvmv52.src$*.9.1t.H wvmv52.src$*.10.1t.H wvmv52.src$*.11.1t.H wvmv52.src$*.12.1t.H >tmp3.H
	Add wvmv52.src$*.13.1t.H wvmv52.src$*.14.1t.H wvmv52.src$*.15.1t.H wvmv52.src$*.16.1t.H >tmp4.H
	Add tmp1.H tmp2.H tmp3.H tmp4.H >$@

wvmv53.sum%.H:
	Add wvmv53.src$*.1.1t.H wvmv53.src$*.2.1t.H wvmv53.src$*.3.1t.H wvmv53.src$*.4.1t.H >tmp1.H
	Add wvmv53.src$*.5.1t.H wvmv53.src$*.6.1t.H wvmv53.src$*.7.1t.H wvmv53.src$*.8.1t.H >tmp2.H
	Add wvmv53.src$*.9.1t.H wvmv53.src$*.10.1t.H wvmv53.src$*.11.1t.H wvmv53.src$*.12.1t.H >tmp3.H
	Add wvmv53.src$*.13.1t.H wvmv53.src$*.14.1t.H wvmv53.src$*.15.1t.H wvmv53.src$*.16.1t.H >tmp4.H
	Add tmp1.H tmp2.H tmp3.H tmp4.H >$@

wvmv2.sum%.H:
	Add wvmv2.src$*.1.1t.H  wvmv2.src$*.2.1t.H wvmv2.src$*.3.1t.H wvmv2.src$*.4.1t.H >tmp1.H
	Add wvmv2.src$*.5.1t.H  wvmv2.src$*.6.1t.H wvmv2.src$*.7.1t.H wvmv2.src$*.8.1t.H >tmp2.H
	Add wvmv2.src$*.9.1t.H  wvmv2.src$*.10.1t.H wvmv2.src$*.11.1t.H wvmv2.src$*.12.1t.H >tmp3.H
	Add wvmv2.src$*.13.1t.H wvmv2.src$*.14.1t.H wvmv2.src$*.15.1t.H wvmv2.src$*.16.1t.H >tmp4.H
	Add tmp1.H tmp2.H tmp3.H tmp4.H >$@

run:
	make runmultiple1
	make runmultiple2
	make runmultiple51
	make runmultiple52

######################################################################################################################
######################################################################################################################
velview.v:
	 < vel1.H Window n3=1 f3=150 | Scale rscale=0.001 | Grey wantscalebar=1 label1='z(m)' label2='x(m)' title='a ' barlabel='V(km/s)' out=va.v labelsz=9 titlesz=12 >/dev/null bias=2.85
	< vel1.rndbnd2.1.H Window n3=1 f3=150 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=1 label1='z(m)' label2='x(m)' title='b ' barlabel='V(km/s)' out=vb.v labelsz=9 titlesz=12  >/dev/null
	< vel1.rndbnd51.1.H Window n3=1 f3=150 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=1 label1='z(m)' label2='x(m)' title='c ' barlabel='V(km/s)' out=vc.v labelsz=9 titlesz=12 >/dev/null
	< vel1.rndbnd52.1.H Window n3=1 f3=150 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=1 label1='z(m)' label2='x(m)' title='d ' barlabel='V(km/s)' out=vd.v labelsz=9 titlesz=12 >/dev/null
	vppen < va.v vpstyle=n xshift=0 yshift=10 > 1.v ; 
	vppen < vb.v vpstyle=n xshift=13 yshift=10 > 2.v ; 
	vppen < vc.v vpstyle=n xshift=0 yshift=0 > 3.v ; 
	vppen < vd.v vpstyle=n xshift=13 yshift=0 > 4.v ; 
	vppen 1.v 2.v 3.v 4.v erase=o scale=0.5 > $@ 

wvmvsix.v:
	< wvmv1.src6.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=23.1375 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=10.2675 yshift=23.1375 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=0 yshift=15.425 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=15.425 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=0 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=10.2675 yshift=7.7125 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=0 yshift=0 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=10.2675 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1

wvmvsixhoriz.v:
	< wvmv1.src6.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src6.1.1t.H Scale rscale=16 | Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum6.H Grey clip=0.025 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=7.7125 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=0 yshift=0 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=10.2675 yshift=7.7125 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=0 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=20.535 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=20.535 yshift=0 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=30.8025 yshift=7.7125 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=30.8025 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1.horiz

wvmvtwentyfive.v:
	< wvmv1.src242.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=23.1375 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=10.2675 yshift=23.1375 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=0 yshift=15.425 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=15.425 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=0 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=10.2675 yshift=7.7125 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=0 yshift=0 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=10.2675 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1

wvmvtwentyfivehoriz.v:
	< wvmv1.src242.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src242.1.1t.H Scale rscale=16 | Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum242.H Grey clip=0.0125 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=7.7125 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=0 yshift=0 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=10.2675 yshift=7.7125 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=0 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=20.535 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=20.535 yshift=0 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=30.8025 yshift=7.7125 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=30.8025 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1.horiz

wvmvsixkabs.v:
	< wvmv1.src6.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey color=j newclip=1 bclip=0 eclip=0.1  label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=23.1375 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=10.2675 yshift=23.1375 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=0 yshift=15.425 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=15.425 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=0 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=10.2675 yshift=7.7125 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=0 yshift=0 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=10.2675 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1

wvmvsixkabshoriz.v:
	 < wvmv1.src6.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1.v >/dev/null
	< wvmv2.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2.v >/dev/null
	< wvmv51.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3.v >/dev/null
	< wvmv52.src6.1.1t.kabs.H Scale rscale=16 | Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4.v >/dev/null
	< wvmv1.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv1s.v >/dev/null
	< wvmv2.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv2s.v >/dev/null
	< wvmv51.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1 label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv3s.v >/dev/null
	< wvmv52.sum6.kabs.H Window f1=25 n1=100 f2=25 n2=100 | Grey clip=0.1  label1=' ' label2=' ' title=' ' wantaxis=no out=wvmv4s.v >/dev/null
	vppen < wvmv1.v vpstyle=n xshift=0 yshift=7.7125 > 1.v ; 
	vppen < wvmv1s.v vpstyle=n xshift=0 yshift=0 > 1s.v ; 
	vppen < wvmv2.v vpstyle=n xshift=10.2675 yshift=7.7125 > 2.v ; 
	vppen < wvmv2s.v vpstyle=n xshift=10.2675 yshift=0 > 2s.v ; 
	vppen < wvmv3.v vpstyle=n xshift=20.535 yshift=7.7125 > 3.v ; 
	vppen < wvmv3s.v vpstyle=n xshift=20.535 yshift=0 > 3s.v ; 
	vppen < wvmv4.v vpstyle=n xshift=30.8025 yshift=7.7125 > 4.v ; 
	vppen < wvmv4s.v vpstyle=n xshift=30.8025 yshift=0 > 4s.v ; 
	vppen 1.v 1s.v 2.v 2s.v 3.v 3s.v 4.v 4s.v erase=o scale=0.3 >tmp.v 
	vp_annotate batch=y <tmp.v >$@ text=my.txt1.horiz

srctwentyfive.v:
	<src242.H Window f2=150 n2=1 f3=150 n3=1 f1=35 n1=30 | Graph plotfat=15 out=1.v label1='t(s) ' label2='amplitude ' title='a' >/dev/null labelsz=12 wheretitle=t titlesz=15 titlefat=4 labelfat=4 
	<src242.H Window f2=150 n2=1 f3=150 n3=1 | ${BINDIR}/XTransf.x | Cabs | Window n1=300 | Graph plotfat=15 out=2.v label1='f(Hz) ' label2='amplitude ' title='b' max2=4.2 >/dev/null labelsz=12 wheretitle=t titlesz=15 titlefat=4 labelfat=4 
	vp_SideBySideAniso 1.v 2.v >$@ 

velfullrndbnd.v:
	<vel1.rndbnd2.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=0 label1='z(m)' label2='x(m)' title='a' barlabel='V(km/s)' out=va.v labelsz=9 titlesz=12 >/dev/null
	<vel1.rndbnd51.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=0 label1='z(m)' label2='x(m)' title='b' barlabel='V(km/s)' out=vb.v labelsz=9 titlesz=12 >/dev/null
	<vel1.rndbnd52.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Grey bias=2.85 wantscalebar=0 label1='z(m)' label2='x(m)' title='c' barlabel='V(km/s)' out=vc.v labelsz=9 titlesz=12 >/dev/null
	 < vel1.rndbnd2.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Rtoc | ${BINDIR}/XFt.x center=1 sign=1 | Transp plane=12 | ${BINDIR}/XFt.x sign=1 center=1 | Cabs | Grey clip=2500 wantscalebar=0 label1='kz(1/m)' label2='kx(1/m)' title='d' out=vakabs.v labelsz=9 titlesz=12 >/dev/null
	< vel1.rndbnd51.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Rtoc | ${BINDIR}/XFt.x center=1 sign=1 | Transp plane=12 | ${BINDIR}/XFt.x sign=1 center=1 | Cabs | Grey clip=2500 wantscalebar=0 label1='kz(1/m)' label2='kx(1/m)' title='e' out=vbkabs.v labelsz=9 titlesz=12 >/dev/null
	< vel1.rndbnd52.1.H Window n3=1 f3=15 | Scale rscale=0.001 | Rtoc | ${BINDIR}/XFt.x center=1 sign=1 | Transp plane=12 | ${BINDIR}/XFt.x sign=1 center=1 | Cabs | Grey clip=2500 wantscalebar=0 label1='kz(1/m)' label2='kx(1/m)' title='f' out=vckabs.v labelsz=9 titlesz=12 >/dev/null
	Cat vel1.rndbnd2.1.H vel1.rndbnd51.1.H vel1.rndbnd52.1.H axis=4 | Window n3=1 f3=15 | Window | Scale rscale=0.001 | Rtoc | ${BINDIR}/XFt.x center=1 sign=1 | Transp plane=12 | ${BINDIR}/XFt.x sign=1 center=1 | Cabs | Stack normalize=y | Graph plotfat=15 label2='amp' label1='kz(1/m)' title=' ' legend=1 legendloc=tr out=vkabskzstk.v titlefat=5 labelfat=5 curvelabel="20m cubic:200m cubic:200m random" >/dev/null dash=1,3,2
	Cat vel1.rndbnd2.1.H vel1.rndbnd51.1.H vel1.rndbnd52.1.H axis=4 | Window n3=1 f3=15 | Window | Scale rscale=0.001 | Rtoc | ${BINDIR}/XFt.x center=1 sign=1 | Transp plane=12 | ${BINDIR}/XFt.x sign=1 center=1 | Cabs | Transp plane=12 | Stack normalize=y | Graph plotfat=15 label2='amp' label1='kx(1/m)' title=' ' legend=1 legendloc=tr out=vkabskxstk.v titlefat=5 labelfat=5 curvelabel="20m cubic:200m cubic:200m random" >/dev/null dash=1,3,2
	vppen < va.v vpstyle=n xshift=0 yshift=10 > 1.v ; 
	vppen < vb.v vpstyle=n xshift=12 yshift=10 > 2.v ; 
	vppen < vc.v vpstyle=n xshift=24 yshift=10 > 3.v ; 
	vppen < vakabs.v vpstyle=n xshift=0 yshift=0 > 4.v ; 
	vppen < vbkabs.v vpstyle=n xshift=12 yshift=0 > 5.v ; 
	vppen < vckabs.v vpstyle=n xshift=24 yshift=0 > 6.v ; 
	vppen < vkabskxstk.v vpstyle=n xshift=0 yshift=0 > 7.v ; 
	vppen < vkabskzstk.v vpstyle=n xshift=12 yshift=0 > 8.v ; 
	vppen 1.v 2.v 3.v 4.v 5.v 6.v erase=o scale=0.3 > $@ 
	vppen 7.v 8.v erase=o scale=0.3 >vkabsstk.v

################################################################################
./Fig/wvmv%.pdf:wvmv%.v
	pstexpen wvmv$*.v wvmv$*.ps 
	epstopdf wvmv$*.ps 
	cp wvmv$*.pdf $@

./Fig/src%.pdf:src%.v
	pstexpen src$*.v src$*.ps fat=1.5
	epstopdf src$*.ps 
	cp src$*.pdf $@

./Fig/v%.pdf:v%.v
	pstexpen v$*.v v$*.ps
	epstopdf v$*.ps 
	cp v$*.pdf $@

######################################################################################################################
ER:
	make exec
	make vel1.H
	make src242.H
	make rndbnd2
	make rndbnd51
	make rndbnd52
	make ./Fig/velview.pdf
	make ./Fig/velfullrndbnd.pdf
	make ./Fig/srctwentyfive.pdf
	make ./Fig/vkabsstk.pdf

CR:
	make exec
	make run	
	make ./Fig/wvmvsixhoriz.pdf
	make ./Fig/wvmvsixkabshoriz.pdf
	make ./Fig/wvmvtwentyfivehoriz.pdf
#####################################3
clean:
	Rm *.H 
	rm Bin/XC* -rf
	rm Bin/XFt.x -rf
	rm Obj/* -rf	 
	rm Fig/* -rf

#include ${SEPINC}/SEP.bottom
