#define DEBUG
#define CMLIB -v -lcmssl -lcmtiming 
#define XCMLIB -lcmfortgraph_slice -lcmsr_slice -lX11
#define XCMLIBOLD -lcmfortgraph -lcmsr -lX11
#define BINDIR /r3/martin/bin/MTYPE
#define CURDIR [[pwd]]
#define SEQ -S1

#include <SEP.defs>


preoas&:  wavelet.oas
	Geninput symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 nx=208 nz=208 dx=40. dz=40. | Transp plane=12 >modelinput.oas

/* make a two component source */
wavelet.oas:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Spike n1=128 n2=1 n3=1 d1=0.004 >junk2
	Merge junk junk2 space=n axis=2 >wavelet.oas


prehan&:  wavelet.han
	Geninput symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 nx=208 nz=208 dx=40. dz=40. |Transp plane=12>modelinput.han

wavelet.han:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=25. | Scale dscale=1. >wavelet.han



boundarytest&: prehan
	time BINDIR/Ti2d  par=parfile.han uzsnap=reczsnap.han uxsnap=recxsnap.han <wavelet.han >jrec.han
	time BINDIR/Ti2dnonrec  par=parfile.han uzsnap=nonreczsnap.han uxsnap=nonrecxsnap.h <wavelet.han >jnonrec.han
	Add  recxsnap.han nonrecxsnap.han mode=subtract >jaddx
	Add  reczsnap.han nonreczsnap.han mode=subtract >jaddz
	


runhan&: 
	time BINDIR/Ti2d  par=parfile.han <wavelet.han >junk.han

run&: 
	time cmattach   BINDIR/Ultimod par=parfile.oas <wavelet.oas >junk

run2&: 
	time cmattach    BINDIR/Ultimod par=parfile.oas <wavelet.oas >junk

runTi2d&:
	time Ti2d par=parfile.oas >junk <wavelet.oas

spike.h:
	Spike n1=128 d1=0.001 nsp=1 k1=1 >spike.h

#define SIZE  nx=208 nz=50 dx=20. dz=20.

wavelet.layer.oas:
	Wavelet n1=128 d1=0.0015 domain=time wavelet=ricker2 tdelay=0.065 fund=15 fhigh=30. | Scale dscale=-1000000. >wavelet.layer.oas

wavelet.isolayer.oas:
	Wavelet n1=128 d1=0.002 domain=time wavelet=ricker1 tdelay=0.055 fund=15 fhigh=30. | Scale dscale=1. >wavelet.isolayer.oas



homoiso.h:
	Geninput nx=208 nz=208 dx=40. dz=40.  symmetry=2 vp=4.0 vs=2.3  rho=2.5   >homoiso.h


isolayer.h:
	Geninput nx=208 nz=70 dx=20. dz=20.  symmetry=2 vp=4.0 vs=2.3  rho=2.5   >junk0
	Geninput nx=208 nz=138 dx=20. dz=20.  symmetry=2 vp=4.5 vs=3.2  rho=2.0   >junk1
	Merge junk0 junk1  space=n axis=2 |Transp plane=12 >isolayer.h



layer.h: 
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 vp=4.0 vs=3.3  rho=2.5   >junk1
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 >junk0
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 vp=3.0 vs=2.0  rho=2.5     >junk3
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 >junk2
	Geninput nx=208 nz=8 dx=40. dz=15.  symmetry=2 vp=4.0 vs=2.3  rho=2.5     >junk4
	Merge junk0 junk1 junk2 junk3 junk4 space=n axis=2 |Transp plane=12 >layer.h

step.h: layer.h 
	Window <layer.h n2=130 >junkleft
	Window <layer.h n2=78 f1=20 >junkb
	Window <layer.h n2=78 n1=20 >junkc
	Merge junkb junkc space=n axis=1 >junkright
	Merge junkleft junkright space=n axis=2 >step.h
	-RM_CMD junk*  

small.h: step.h
	Window n1=80 n2=80 <step.h >small.h

/* smallmode&: small.h wavelet.layer.oas */

smallmod&:
        echo "d1=0.004" >>wavelet.layer.oas
        echo "d1=0.04 d2=0.04" >>small.h
        time cmattach  BINDIR/Ultimod modeling=1 moduli=small.h par=parfile.small.oas <wavelet.layer.oas >junk

smallmig&:
	time cmattach   BINDIR/Ultimod modeling=0 moduli=small.h par=parfile.small.oas <wavelet.layer.oas >junk




semmod&: step.h wavelet.layer.oas
	echo "d1=0.004" >>wavelet.layer.oas
	echo "d1=0.04 d2=0.04" >>step.h
	time cmattach -p 8k   BINDIR/Ultimod modeling=1 moduli=step.h par=parfile.test.oas <wavelet.layer.oas >junk

semmig&:
	time cmattach   BINDIR/Ultimod modeling=0 moduli=step.h par=parfile.test.oas <wavelet.layer.oas >junk


modroll&: step.h wavelet.layer.oas
	time cmattach -p 8k BINDIR/Ultimod modeling=1 moduli=step.h par=parfile.roll.oas <wavelet.layer.oas >junk

migroll&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod modeling=0 moduli=step.h par=parfile.roll.oas <wavelet.layer.oas >junk

migrolliso&: homoiso.h wavelet.layer.oas
	time cmattach   BINDIR/Ultimod  datapath=/scr2/martin/ modeling=0 moduli=homoiso.h par=parfile.rolliso.oas  <wavelet.layer.oas >junk

modrollX&: step.h wavelet.layer.oas
	time cmattach -p 8k BINDIR/Ultimod modeling=1 moduli=step.h par=parfile.rollX.oas <wavelet.layer.oas >junk

migrollX&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod modeling=0 moduli=step.h par=parfile.rollX.oas <wavelet.layer.oas >junk



modrollfree&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod datapath=/scr2/martin/ modeling=1 moduli=step.h par=parfile.rollfree.oas <wavelet.layer.oas >junk

migrollfree&: step.h wavelet.layer.oas
	time cmattach   BINDIR/Ultimod datapath=/scr2/martin/ modeling=0 moduli=step.h par=parfile.rollfree.oas <wavelet.layer.oas >junk

modrollXfree&: step.h wavelet.layer.oas
	time cmattach   BINDIR/Ultimod datapath=/scr2/martin/ modeling=1 moduli=step.h par=parfile.rollXfree.oas <wavelet.layer.oas >junk

migrollXfree&: step.h wavelet.layer.oas
	time cmattach   BINDIR/Ultimod datapath=/scr2/martin/ modeling=0 moduli=step.h par=parfile.rollXfree.oas <wavelet.layer.oas >junk


testmig&: step.h wavelet.layer.oas
	Window f3=8 n3=1 <Zxseis.roll.oas >Zxseis.testmig.oas
	Window f3=8 n3=1 <Zzseis.roll.oas >Zzseis.testmig.oas
	time cmattach -p 8k  BINDIR/Ultimod modeling=0 moduli=step.h par=parfile.testmig.oas <wavelet.layer.oas >junk


modstep&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfile.step.oas <wavelet.layer.oas >junk

migstep&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfile.step.oas <wavelet.layer.oas >junk

modstep1&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfile1.step.oas <wavelet.layer.oas >junk

migstep1&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfile1.step.oas <wavelet.layer.oas >junk

modstep2&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfile2.step.oas <wavelet.layer.oas >junk

migstep2&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfile2.step.oas <wavelet.layer.oas >junk


modstepX&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfileX.step.oas <wavelet.layer.oas >junk

migstepX&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfileX.step.oas <wavelet.layer.oas >junk

modstepX1&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfileX1.step.oas <wavelet.layer.oas >junk

migstepX1&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfileX1.step.oas <wavelet.layer.oas >junk

modstepX2&: step.h wavelet.layer.oas
	time cmattach  BINDIR/Ultimod-0.7 modeling=1 moduli=step.h par=parfileX2.step.oas <wavelet.layer.oas >junk

migstepX2&: step.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=step.h par=parfileX2.step.oas <wavelet.layer.oas >junk



modlayer&: layer.h wavelet.layer.oas
	time cmattach   BINDIR/Ultimod-0.7 modeling=1 moduli=layer.h par=parfile.layer.oas <wavelet.layer.oas >junk

miglayer&: layer.h wavelet.layer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0 moduli=layer.h par=parfile.layer.oas <wavelet.layer.oas >junk


modisosplit&: homoiso.h wavelet.oas
	time cmattach BINDIR/Ultisplit  modeling=1 moduli=homoiso.h par=parfile.homoiso.oas <wavelet.oas >junk

modiso&: homoiso.h wavelet.oas
	time cmattach BINDIR/Ultimod  type=3 src_type=0 ncomp=3 modeling=1 moduli=homoiso.h par=parfile.homoiso.oas <wavelet.oas >junk

migiso&: homoiso.h wavelet.oas
	time cmattach  BINDIR/Ultimod modeling=0 moduli=homoiso.h par=parfile.homoiso.oas <wavelet.oas >junk

verysmall&: homoiso.h
	Window n1=80 n2=80 <homoiso.h >vsmall.h
	time cmattach  BINDIR/testjunk modeling=0 moduli=vsmall.h par=parfile.homoiso.oas <wavelet.oas >junk
	

testhan&: homoiso.h wavelet.layer.oas
	Lin2d  centered=1 moduli=homoiso.h par=parfile.homoiso.han <wavelet.layer.oas >junk



modisolayer&:  isolayer.h wavelet.isolayer.oas
	time cmattach -p 8k   BINDIR/Ultimod-0.7 modeling=1 moduli=isolayer.h par=parfile.isolayer.oas <wavelet.isolayer.oas >junk

migisolayer&:  isolayer.h wavelet.isolayer.oas
	time cmattach -p 8k  BINDIR/Ultimod-0.7 modeling=0  moduli=isolayer.h par=parfile.isolayer.oas <wavelet.isolayer.oas >junk




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

%.T:	
	Taplot <% gainpanel=every  >%.T

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











