# Finite-difference

moduli.H: vel.H
	Add mode=product vel.H vel.H      > c11.H
	< vel.H Add scale=0. | Add add=1. > rho.H
	Merge axis=3 space=n c11.H rho.H  >$@
	echo 'd1=0.01 d2=0.0125' >>$@
	Rm c11.H rho.H
wave.H:
	Wavelet wavelet=ricker1 fund=55 n1=64 d1=0.001 tdelay=0.03 >wave.H


#Sadly, Ultimod is no more available in seplib!!!!!!
fdseis%: ./par/ulti%.P wave.H moduli.H overmovie.x layers.v
	< wave.H Ultimod par=./par/ulti.P par=./par/ulti$*.P moduli=moduli.H >junk
	< zseis Transp plane=12 | Window j2=2  | Window min1=1. | Halfint inv=1 |\
	$(grey) title="Shot $*" label1='time (s)' label2='x (km)' out=${RESDIR}/fdseis$*.v pclip=100
	< zsnap Window j1=3 j2=3 | Taplot gainpanel=all pclip=99 >zsnap.T
	overmovie.x zsnap.T ${RESDIR}/fdsnap$*.v3 35
	< zsnap Window n3=1 f3=20 | $(grey) $(ratio) out=junk1.v pclip=100
	vp_Overlay junk1.v layers.v > ${RESDIR}/fdsnap$*.v

fdmseis%:
	Fdmod intag=vel1.H no_stdout=0 par=./par/fdmsnap$*.P hsfile=fdmsism$*.H outtag=fdmsnap$*.H > /dev/null 
# tmax=6 nxs=1 oxs=6000 dxs=25 nzs= 1 ozs= 0 2750 dzs= 10 jt=10 bx=100 ex=100 verb=0
#	< wave.H Ultimod par=./par/ulti.P par=./par/ulti$*.P moduli=moduli.H >junk
#	< zseis Transp plane=12 | Window j2=2  | Window min1=1. | Halfint inv=1 |\
#	$(grey) title="Shot $*" label1='time (s)' label2='x (km)' out=${RESDIR}/fdseis$*.v pclip=100
#	< zsnap Window j1=3 j2=3 | Taplot gainpanel=all pclip=99 >zsnap.T
#	overmovie.x zsnap.T ${RESDIR}/fdsnap$*.v3 35
#	< zsnap Window n3=1 f3=20 | $(grey) $(ratio) out=junk1.v pclip=100
#	vp_Overlay junk1.v layers.v > ${RESDIR}/fdsnap$*.v
#