0014data:
	Wavelet n1=501 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1. domain=time > j2.H;echo o1=-1. >> j2.H
	<j2.H Transf f_min=0 f_min1=4 f_max1=36 f_max=42 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13|Window3d squeeze=n > 0015data_source.H
	Vel xa1=-10.2 xb1=9.8 z1=0 alfa1=10 vr1=2.8 vc=3 const1=1 n1=401 o1=0 d1=0.005 n2=1001 \
	o2=-10.2 d2=0.02 n3=1 o3=0 d3=1 | Pad extend=1 beg2=1000 end2=1000 > 0015refl.H
	Spike n2=1001 o2=-10.2 d2=0.02 n3=1 o3=0 d3=1 n1=401 o1=0 d1=0.005 nsp=1 k1=1 k2=0 k3=1 > 1.H
	Math file1=1.H exp="file1*0+3" | Pad extend=1 beg2=1000 end2=1000 > 2.H
	../ShtPSPI/bin/OWmod2d_par.x sou=0015data_source.H vel=2.H \
	vreflect=0015refl.H shots=0015data_born.H \
	vsal=4 nrefvel=1 traces_in_taper=140 padtraces=60 xmin=-14 xmax=14 \
	nrec=201 orec=-2 drec=.02 nsht=75 osht=-6 dsht=0.08 zmin=0 zmax=2 node=16 > /dev/null
#	Spike n1=501 d1=0.004 o1=-1. n2=1 d2=1 o2=0 n3=1 o3=0 d3=1 k1=251 nsp=1 | Bandpass flo=0 fhi=30 > j2.H
#	<j2.H Transf f_min=0 f_min1=4 f_max1=30 f_max=38 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs|Smooth rect1=11 > j2.H
#	Math file1=j2.H exp="file1*0">j3.H;Cmplx j2.H j3.H|Transp plane=13 > 0015data_source.H

0015data_hor.H:
	susynlv ref="1:-8.,0.;8.,2.82" nt=500 dt=0.004 dxs=0.08 nxs=101 fpeak=15 fxs=-4. dxo=0.02 nxo=201 fxo=-2. v00=3 > junk1.su
	< junk1.su Su2sep > junk1.H
	< junk1.H Infill3d > $@; echo n2=201 o2=-2. d2=0.02 n3=101 o3=-4 d3=0.08 >> $@

0015data_prep.H:
	< 0015data_hor.H Transf shotmig=y is2d=y wei=y f_min=0 f_min1=4 f_max1=24 f_max=30 maxsize=1000 > j2.H
	< j2.H Pad beg1=100 end1=100 > $@

0015data_source.H:
	Spike n1=501 d1=0.004 o1=-1. n2=1 d2=1 o2=0 n3=1 o3=0 d3=1 k1=251 nsp=1 | Bandpass flo=0 fhi=30 > j2.H
	<j2.H Transf f_min=0 f_min1=4 f_max1=30 f_max=38 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs|Smooth rect1=11 > j2.H
	Math file1=j2.H exp="file1*0">j3.H;Cmplx j2.H j3.H|Transp plane=13 > 0015data_source.H

0016mig:
	Spike n2=1501 o2=-15. d2=0.02 n3=1 o3=0 d3=1 n1=401 o1=0 d1=0.005 nsp=1 k1=1 k2=0 k3=1 > 1.H
	Math file1=1.H exp="file1*0+1/2.7"|Transp plane=12 reshape=1,3 | Pad beg1=1000 end1=1000 extend=1 > 0015slow.H
	../Tese/bin/OWmig2d_areal.x sou=0015data_source.H rec=0015data_born.H  \
	vel=0015slow.H bkgimg=0015dana.bkim.H \
	zmin=0 zmax=2 nh=51 oh=-0.5 migration=1 areal=0 \
	padtraces=0 traces_in_taper=25 node=14 vsal=4 nrefvel=1 > /dev/null

0017rotsd:
	Window3d <0015dana.bkim.H min1=-6 n1=301 |Transp plane=13| ./bin/OW_3dslant.x nax=5 21 nah=51 oax=8 -20 oah=-50 dax=1 2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 > 0016.rot.sou.H
	Window3d <0015dana.bkim.H min1=-6 n1=301 |Transp plane=13| ./bin/OW_3dslant.x nax=5 21 nah=51 oax=8 -20 oah=-50 dax=1 2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 rec=1 > 0016.rot.rec.H
#
0018perm:
#	Window3d < 0016.rot.sou.H n3=1 min3=-3 squeeze=n | Pad beg3=200 end3=200|Transp plane=13 > s1.H
#	Window3d < 0016.rot.rec.H n3=1 min3=-3 squeeze=n |Pad beg3=200 end3=200|Transp plane=13 > r1.H
	< 0016.rot.sou.H Window3d min3=-5 n3=201 squeeze=n | Transp plane=13 > s1.H
	< 0016.rot.rec.H Window3d min3=-5 n3=201 squeeze=n | Transp plane=13 > r1.H
	Gpow gpow=-1 < 0015slow.H | Window3d | Transp > 1.H
	../Tese/bin/OW_arealsht.x rec_cig=r1.H sou_cig=s1.H sou=0015data_source.H vel=1.H \
	rec_out=0018.crec.perm.H sou_out=0018.csou.perm.H nrefvel=1 vsal=4 \
	padtraces=100 traces_in_taper=60 encode=1 encz=1 eps=0 random=1 injection=1 nreal=1 \
	ngap=51 20 zmin=0 zmax=2 node=14 > /dev/null 
#	ws=0018ws.H wr=0018wr.H \

0019mig:
	Spike n2=1501 o2=-15. d2=0.02 n3=1 o3=0 d3=1 n1=401 o1=0 d1=0.005 nsp=1 k1=1 k2=0 k3=1 > 1.H
	Math file1=1.H exp="file1*0+1/2.7"|Transp plane=12 reshape=1,3 | Pad beg1=1000 end1=1000 extend=1 > 0015slow.H
	Transp plane=13 < 0018.csou.perm.H > s1.H
	Transp plane=13 < 0018.crec.perm.H > r1.H
	../Tese/bin/OWmig2d_areal.x sou=s1.H rec=r1.H  \
	vel=0015slow.H bkgimg=0019.bkim.encrsxz.H \
	zmin=0 zmax=2 nh=21 oh=-0.2 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=4 nrefvel=1 > /dev/null