#define DEBUG
#define FIGDIR ../Fig
#define CMLIB -v -lcmssl -lcmtiming 
#define XCMLIB -lcmfortgraph_slice -lcmsr_slice -lX11
#define XCMLIBOLD -lcmfortgraph -lcmsr -lX11
#define BINDIR /r3/martin/bin/MTYPE
#define SEQ -S0-1
#include <SEP.defs>

all&: vert2c.sin vert2c.imp

clean&: 
	-Rm *seis* 
	-Rm *snap*
	-Rm wavelet2
	-Rm wavelet
	-Rm scatter.h
	-Rm j*

test&:  wavelet2 scatter.h  
	time cmattach SEQ BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=scatter.h par=parfile  <wavelet2 >junk

/* 2 component modeling of drill bit source in an elastic medium */
/* vert  =  vertical force                   */

vert2c.imp&:  wavelet2 scatter.h  
	time cmattach SEQ BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=scatter.h par=parfile  <wavelet2 >junk xseis=xseis.imp zseis=zseis.imp xsnap=xsnap.imp zsnap=zsnap.imp


vert2c.sin&:   wavelet scatter.h 
	time cmattach SEQ BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=scatter.h par=parfile  <wavelet >junk xseis=xseis.sin zseis=zseis.sin xsnap=xsnap.sin zsnap=zsnap.sin


wavelet:
	Sinusoid nt=2000 dt=0.001 omega=10 >junk
	Spike n1=2000 n2=1 n3=1 d1=0.001 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk junk2 space=n axis=2 |Scale dscale=1000000. >wavelet

wavelet2:
	Wavelet n1=128 d1=0.001 domain=time wavelet=ricker1 tdelay=0.050 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Spike n1=128 n2=1 n3=1 d1=0.001 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk junk2 space=n axis=2 >wavelet2

homoiso.h:
	Geninput nx=208 nz=208 dx=10. dz=10.  symmetry=2 vp=3.0 vs=1.8  rho=2.3   >homoiso.h

borehole.h: homoiso.h
	Geninput nx=1 nz=60 dx=10. dz=10.  symmetry=2 vp=3.0 vs=0  rho=2.3 |Transp plane=12  >junk1
	echo "n3=5 n2=1" >>junk1
	Window n1=148 n2=1 <homoiso.h >junk2
	echo "n3=5 n2=1" >>junk2
	Merge junk1 junk2 axis=1 space=n >borehole.h

background.h: homoiso.h borehole.h
	Window f2=0 n2=59 <homoiso.h >junk1
	Window f2=60      <homoiso.h >junk2
	Merge junk1 borehole.h junk2 axis=2 space=n > background.h

scatter.h: background.h
	Spike n1=208 n2=208 k1=120,120,60  k2=60,120,160 nsp=3 d1=0.01 d2=0.01 >junk
	Triconv rect1=4 rect2=4 <junk |Scale dscale=64 >junk2
	Cp <junk2 >junk3
	Merge junk2 junk3 axis=3 space=n >junk4
	Merge junk2 junk3 junk4 axis=3 space=n >junk5
	Merge junk2 junk5 axis=3 space=n >junk6
	Add junk6 background.h add=1.,0. mode=product >scatter.h


Sinusoid: Sinusoid.o
	FC FCOPTS Sinusoid.o SEPLIB SYSLIB -o BINDIR/Sinusoid

clean&:
	-RM_CMD *.x *.o  a.out core junk*


xsnap.imp.T:
	Taplot gainpanel=every <xsnap.imp >xsnap.imp.T

zsnap.imp.T:
	Taplot gainpanel=every <zsnap.imp >zsnap.imp.T

xsnap.sin.T:
	Taplot gainpanel=every <xsnap.sin >xsnap.sin.T

zsnap.sin.T:
	Taplot gainpanel=every <zsnap.sin >zsnap.sin.T

%.T:	
	Transp plane=12 <% |Reverse which=2 >junk
	Taplot <junk gainpanel=every  >%.T



figures&: xseis.imp.T zseis.imp.T xseis.sin.T zseis.sin.T
	Ta2vplot < xseis.imp.T title="Imp Z-Src in Borehole   x-rec" head=/dev/null out=FIGDIR/xseis.imp.v 
	Ta2vplot < zseis.imp.T title="Imp Z-Src in Borehole   z-rec" head=/dev/null out=FIGDIR/zseis.imp.v
	Ta2vplot < xseis.sin.T title="Sin Z-Src in Borehole   x-rec" head=/dev/null out=FIGDIR/xseis.sin.v 
	Ta2vplot < zseis.sin.T title="Sin Z-Src in Borehole   z-rec" head=/dev/null out=FIGDIR/zseis.sin.v
	Window n3=1 f3=3 <scatter.h >junk
	Taplot <junk pclip=100. |Ta2vplot title="Drill bit source modeling" label1="depth [km]" label2="distance [km]" head=/dev/null out=FIGDIR/scatter.v
	Window n2=1 <wavelet2 >junk
	Graph <junk title="Imp Wavelet" |Vppen vpstyle=n scale=0.5 head=/dev/null out=FIGDIR/wav.imp.v
	Window n2=1 <wavelet >junk
	Graph <junk title="The Pump: Sin Wavelet"   |Vppen vpstyle=n scale=0.5 head=/dev/null out=FIGDIR/wav.sin.v



movies&: xsnap.imp.T zsnap.imp.T xsnap.sin.T zsnap.sin.T
	Ta2vplot < xsnap.imp.T title="Imp Z-Src in Borehole   x-rec" head=/dev/null out=FIGDIR/xseis.imp.v3
	Ta2vplot < zsnap.imp.T title="Imp Z-Src in Borehole   z-rec" head=/dev/null out=FIGDIR/zseis.imp.v3
	Ta2vplot < xsnap.sin.T title="Sin Z-Src in Borehole   x-rec" head=/dev/null out=FIGDIR/xseis.sin.v3 
	Ta2vplot < zsnap.sin.T title="Sin Z-Src in Borehole   z-rec" head=/dev/null out=FIGDIR/zseis.sin.v3
	

#include <SEP.obj.rules>
#include <SEP.prog.rules>
