3d_mkcig.H:
	Wavelet n1=21 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=.04 domain=time |Scale dclip=-1 > j1.H
	Math file1=j1.H exp="file1*-1" > j2.H;echo d1=10 o1=-100 >> j2.H
	./bin/OW_makecig.x > 1.H n2=11 o2=-100 d2=20 ;echo d3=0.001 >> 1.H
	Transp plane=13 < 1.H | Filter filter=j2.H pad=zero | Nmo3d adj=1 vel=5000 maxsize=1000 | Transp plane=45 > $@
	echo d1=10 >> $@

3d_mkcig_flat.H:
	Wavelet n1=21 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=.04 domain=time |Scale dclip=-1 > j1.H
	Math file1=j1.H exp="file1*-1" > j2.H;echo d1=10 o1=-100 >> j2.H
	./bin/OW_makecig.x > 1.H n2=11 o2=-100 d2=20 px=0 py=0 ;echo d3=0.001 >> 1.H
	Transp plane=13 < 1.H | Filter filter=j2.H pad=zero | Nmo3d adj=1 vel=5000 maxsize=1000 | Transp plane=45 > $@
	echo d1=10 >> $@

3d_rot2d:
	< 3d_mkcig.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=0 > temp.H 
	< temp.H Transp plane=23 reshape=1,3,5 maxsize=1000 | ./bin/OW_3dslant_2d.x nodes=14 nax=7 \
	nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 nxw=3 rot=1 rec=0 > 3d_rotsou_2d.H
	< 3d_mkcig.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=1 > temp.H
	< temp.H Transp plane=23 reshape=1,3,5 maxsize=1000 | ./bin/OW_3dslant_2d.x nodes=14 nax=7 \
	nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 nxw=3 rot=1 rec=1 > 3d_rotrec_2d.H

3d_rot2.5d:
	< 3d_mkcig5d.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=0 > temp.H 
	< temp.H Transp plane=23 reshape=1,3,5 maxsize=1000 | ./bin/OW_3dslant_2d.x nodes=14 nax=7 \
	nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 nxw=3 rot=1 rec=0 > 3d_rotsou_2.5d.H
	< 3d_mkcig5d.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=1 > temp.H
	< temp.H Transp plane=23 reshape=1,3,5 maxsize=1000 | ./bin/OW_3dslant_2d.x nodes=14 nax=7 \
	nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 nxw=3 rot=1 rec=1 > 3d_rotrec_2.5d.H

3d_rot2dx:
	< 3d_mkcig.H Window3d n4=1 f4=25 | ./bin/OW_3dslant_2d.x nodes=14 nax=16 nahx=41 oax=-30 oahx=-40 dax=2 dahx=2 \
	nxw=15 rot=1 rec=0 > 3d_rotsou_2dx.H

3d_mkcig5d.H:
	Wavelet n1=21 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=.04 domain=time |Scale dclip=-1 > j1.H
	Math file1=j1.H exp="file1*-1" > j2.H;echo d1=10 o1=-100 >> j2.H
	./bin/OW_makecig.x n5=11 o5=-100 d5=20 n2=11 o2=-100 d2=20 > 1.H ;echo d3=0.001 >> 1.H
	Transp plane=13 < 1.H | Filter filter=j2.H pad=zero | Nmo3d adj=1 vel=5000 maxsize=1000 | \
	Transp plane=24 reshape=1,3,4,5 maxsize=1000 | Nmo3d adj=1 vel=5000 maxsize=1000 | Transp plane=23 reshape=1,3,5 > $@
	echo d1=10 >> $@

3d_rot3d:
	< 3d_mkcig5d.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 oax=-18 dax=1 nahx=7 oahx=-15 dahx=5 nxw=3 \
	nay=7 oay=-18 day=1 nahy=7 oahy=-15 dahy=5 nyw=3 rot=1 rec=0 > 3d_rotsou_3d.H
#	< 3d_mkcig5d.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 oax=-18 dax=1 nahx=31 oahx=-20 dahx=2 nxw=7 \
#	nay=7 oay=-18 day=1 nahy=31 oahy=-20 dahy=2 nyw=7 rot=1 rec=1 > 3d_rotsou_3d.H

3d_mkcig5d_2dips.H:
	Wavelet n1=21 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=.04 domain=time |Scale dclip=-1 > j1.H
	Math file1=j1.H exp="file1*-1" > j2.H;echo d1=10 o1=-100 >> j2.H
	./bin/OW_makecig.x n5=21 o5=-200 d5=20 n2=21 o2=-200 d2=20 px=15 py=7 phx=7 phy=-7 > 1.H ;echo d3=0.001 >> 1.H
	Transp plane=13 < 1.H | Filter filter=j2.H pad=zero | Nmo3d adj=1 vel=5000 maxsize=1000 | \
	Transp plane=24 reshape=1,3,4,5 maxsize=1000 | Nmo3d adj=1 vel=6000 maxsize=1000 | Transp plane=23 reshape=1,3,5 > $@
	echo d1=10 >> $@

3d_rot2.5dips:
	< 3d_mkcig5d_2dips.H ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=0 > temp.H 
	< temp.H Transp plane=23 reshape=1,3,5 maxsize=1000 | ./bin/OW_3dslant_2d.x nodes=14 nax=7 nahx=7 oax=-18 oahx=-15 dax=1 dahx=5 \
	nxw=3 rot=1 rec=0 | Transp plane=23 reshape=1,3,5 maxsize=1000 > 3d_rotsou_2.5dips.H

3dperm10:
	Transp plane=15 maxsize=1000 < 3d_rotrec_2d.H | Transp plane=24 maxsize=1000 > 3dr.H ; echo n6=1 o4=0 d4=20 >> 3dr.H
	Transp plane=15 maxsize=1000 < 3d_rotsou_2d.H | Transp plane=24 maxsize=1000 > 3ds.H ; echo n6=1 o4=0 d4=20 >> 3ds.H
	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Window3d f5=1 | Cabs | Smooth rect1=11 > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+3000" > 3d.vel3000.H
	./bin/OW3d_arealsht.x rec_cig=3dr.H sou_cig=3ds.H sou=11111.H vel=3d.vel3000.H \
	rec_out=3d.crec10.perm.H sou_out=3d.csou10.perm.H extrapolator=pspi lloyd=0 nrefvel=1 vsal=7000 \
	padtraces_x=25 padtraces_y=0 25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=1 \
	random=1 nreal=1 ngapx=25 ngapy=25 zmin=0 zmax=1000 node=14 zzz=zzz.H \
	> /dev/null 

3dperm10ffd:
#	Transp plane=15 maxsize=1000 < 3d_rotrec_2d.H | Transp plane=24 maxsize=1000 > 3dr.H ; echo n6=1 o4=0 d4=20 >> 3dr.H
#	Transp plane=15 maxsize=1000 < 3d_rotsou_2d.H | Transp plane=24 maxsize=1000 > 3ds.H ; echo n6=1 o4=0 d4=20 >> 3ds.H
#	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
#	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Window3d f5=1 | Cabs | Smooth rect1=11 > j1.H
#	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+3000" > 3d.vel3000.H
	./bin/OW3d_arealsht.x rec_cig=3dr.H sou_cig=3ds.H sou=11111.H vel=3d.vel3000.H \
	rec_out=3d.crec10.perm.ffd.H sou_out=3d.csou10.perm.ffd.H extrapolator=ffd \
	padtraces_x=25 padtraces_y=0 25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=1 \
	random=1 nreal=1 ngapx=25 ngapy=25 zmin=0 zmax=1000 node=14 \
	> /dev/null 

3dperm01ffd:
#	Transp plane=15 maxsize=1000 < 3d_rotrec_2d.H | Transp plane=24 maxsize=1000 | Pad beg1=150 end1=150 beg3=75 end3=75 > 3dr.H ; echo n6=1 o4=0 d4=20 >> 3dr.H
#	Transp plane=15 maxsize=1000 < 3d_rotsou_2d.H | Transp plane=24 maxsize=1000 | Pad beg1=150 end1=150 beg3=75 end3=75 > 3ds.H ; echo n6=1 o4=0 d4=20 >> 3ds.H
#	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
#	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Window3d f5=1 | Cabs | Smooth rect1=11 > j1.H
#	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="1/(file1*0+3000)" > 3d.vel3000.H
	./bin/OW3d_arealsht.x rec_cig=3dr.H sou_cig=3ds.H sou=11111.H vel=3d.vel3000.H \
	rec_out=3d.crec01.perm.ffd.H sou_out=3d.crec01.perm.ffd.H extrapolator=pspi nrefvel=1 vsal=7000 ffd \
	padtraces_x=25 padtraces_y=25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=0 \
	random=1 nreal=1 ngapx=25 ngapy=25 zmin=0 zmax=1000 node=14 nzblk=101 \
	> /dev/null 

3dperm50:
#	Transp plane=15 maxsize=1000 < 3d_rotrec_2.5d.H | Transp plane=24 maxsize=1000 | Window3d squeeze=n min5=500 > 3dr.H ; echo n6=1 >> 3dr.H
#	Transp plane=15 maxsize=1000 < 3d_rotsou_2.5d.H | Transp plane=24 maxsize=1000 | Window3d squeeze=n min5=500 > 3ds.H ; echo n6=1 >> 3ds.H
#	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
#	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Window3d f5=1 | Cabs | Smooth rect1=11 > j1.H
#	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+3000" > 3d.vel3000.H
#	Transp plane=15 maxsize=1000 < 3d_rotrec_2.5d.H | Transp plane=24 maxsize=1000|\
#	Pad beg1=75 end1=75 beg3=75 end3=75 > 3dr.H ; echo n6=1 >> 3dr.H
#	Transp plane=15 maxsize=1000 < 3d_rotsou_2.5d.H | Transp plane=24 maxsize=1000|\
#	Pad beg1=75 end1=75 beg3=75 end3=75 > 3ds.H ; echo n6=1 >> 3ds.H
	Rm 3d.c*50a.perm.H
	./bin/OW3d_arealsht.x rec_cig=3dr.H sou_cig=3ds.H sou=11111.H vel=3d.vel3000.H \
	rec_out=3d.crec50a.perm.H sou_out=3d.csou50a.perm.H extrapolator=ffd \
	padtraces_x=25 padtraces_y=25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=1 \
	random=1 injection=1 nreal=2 ngapx=25 ngapy=25 zmin=0 500 zmax=1000 node=14 nzblk=30 nwblk=30 \
	> /dev/null 

3dperm50flat:
#	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=18 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
#	Transp plane=21 reshape=1,5 maxsize=1000 < 3d_mkcig_flat.H | Transp plane=12 maxsize=1000 | Transp plane=34 maxsize=1000 | \
#	Window3d min1=0 min3=0 n1=1 n3=1 squeeze=n | Pad beg1=180 end1=180 beg3=175 end3=175 > 3d_flat.H ; echo n6=1 >> 3d_flat.H
	Transp plane=21 reshape=1,5 maxsize=1000 < 3d_mkcig_flat.H | Transp plane=12 maxsize=1000 | Transp plane=34 maxsize=1000 | \
	Pad beg1=25 end1=25 beg3=25 end3=25 | Window3d squeeze=n min5=300 max5=660 > 3d_flat.H ; echo n6=1 >> 3d_flat.H
	Spike n1=351 d1=0.004 o1=-0.7 n2=1 d2=1 o2=0 nsp=1 k1=176 | Bandpass flo=2 |\
	Transf f_min=0 f_min1=4 f_max1=52 f_max=60 wei=y is2d=y | Window3d | Cabs | Smooth rect1=11 > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=501 d2=20 o2=0 n1=501 d1=20 o1=0 > j1.H
#	Math file1=j1.H exp="file1*0+1/3000" | Pad beg1=500 beg2=500 > j11.H
#	Spike n3=201 d3=10 o3=0 n2=500 d2=20 o2=-10000 n1=501 d1=20 o1=0 > j2.H
#	Math file1=j2.H exp="file1*0+1/3100" | Pad beg1=500 end2=501 > j22.H
#	Spike n3=201 d3=10 o3=0 n2=500 d2=20 o2=-10000 n1=500 d1=20 o1=-10000 > j3.H
#	Math file1=j3.H exp="file1*0+1/3200" | Pad end1=501 end2=501 > j33.H
#	Spike n3=201 d3=10 o3=0 n2=501 d2=20 o2=0 n1=500 d1=20 o1=-10000 > j4.H
#	Math file1=j4.H exp="file1*0+1/3300" | Pad end1=501 beg2=500 > j44.H
#	Math file1=j11.H file2=j22.H file3=j33.H file4=j44.H exp="file1+file2+file3+file4" > j1234.H 
	Rm 3d.c*50flat.perm.H
	./bin/OW3d_arealsht.x rec_cig=3d_flat.H sou_cig=3d_flat.H sou=11111.H vel=j1234.H \
	rec_out=3d.crec50flat.perm.H sou_out=3d.csou50flat.perm.H extrapolator=ffd \
	padtraces_x=25 padtraces_y=25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=1 shape= 0 \
	random=1 injection=1 nreal=1 ngapx=10 ngapy=10 zmin=300 zmax=650 node=14 nzblk=30 nwblk=30 > /dev/null 

3dpermig50flat:
	Rm 3d.bkim.perm.50flat.H
	Window3d < 3d.csou50flat.perm.H squeeze=n n4= 1 n3= 14 f3= 40 > s1.H
	Window3d < 3d.crec50flat.perm.H squeeze=n n4= 1 n3= 14 f3= 40 > r1.H
	./bin/OWmig3d_areal.x sou=s1.H rec=r1.H \
	vel=j1234.H bkgimg=3d.bkim.perm.50flat.H \
	zmin=300 zmax=650 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=ffd nwblk= 1 > /dev/null

3dpermig51flat:
	Rm 3d.bkim.perm.51flat.H
#	Window3d < 3d.csou50flat.perm.H squeeze=n n4= 1 n3= 14 f3= 40 > s1.H
#	Window3d < 3d.crec50flat.perm.H squeeze=n n4= 1 n3= 14 f3= 40 > r1.H
	./bin/OWmig3d_areal.x sou=s1.H rec=r1.H \
	vel=j1234.H bkgimg=3d.bkim.perm.51flat.H \
	zmin=300 zmax=650 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=ffd nwblk= 1 > /dev/null

3dpermh0:
	Transp plane=15 maxsize=1000 < 3d_rotrec_2d.H | Transp plane=24 maxsize=1000 | Window3d min2=0 n2=1 squeeze=n | Pad beg2=5 end2=5 > 3dr.H ; echo n6=1 o4=0 d4=20 >> 3dr.H
	Transp plane=15 maxsize=1000 < 3d_rotsou_2d.H | Transp plane=24 maxsize=1000 | Window3d min2=0 n2=1 squeeze=n | Pad beg2=5 end2=5 > 3ds.H ; echo n6=1 o4=0 d4=20 >> 3ds.H
	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Cabs | Smooth rect1=11 > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+3000" > 3d.vel3000.H
	./bin/OW3d_arealsht.x rec_cig=3dr.H sou_cig=3ds.H sou=11111.H vel=3d.vel3000.H \
	rec_out=3d.crec25_h0.perm.H sou_out=3d.csou25_h0.perm.H extrapolator=pspi lloyd=0 nrefvel=1 vsal=7000 \
	padtraces_x=25 padtraces_y=25 traces_in_taper_x=25 traces_in_taper_y=25 encode=1 encz=1 eps=1 \
	random=1 nreal=1 ngapx=25 ngapy=25 zmin=0 zmax=1000 node=14 > /dev/null 

3dpermigh0:
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+1/3000" > 3d.slo3000.H
	./bin/OWmig3d_areal.x sou=3d.csou25_h0.perm.H rec=3d.crec25_h0.perm.H \
	vel=3d.slo3000.H bkgimg=3d.bkim.perm.25h0.H \
	lloyd=0 vsal=7000 nrefvel=1 zmin=0 zmax=1000 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 extrapolator=pspi > /dev/null

3dpermig10:
	./bin/OWmig3d_areal.x sou=3d.csou10.perm.H rec=3d.crec10.perm.H \
	vel=3d.slo3000.H bkgimg=3d.bkim.perm.10.H \
	lloyd=0 vsal=7000 nrefvel=1 zmin=0 zmax=1000 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 extrapolator=pspi > /dev/null

3dpermig10hy:
	./bin/OWmig3d_areal.x sou=3d.csou10.perm.H rec=3d.crec10.perm.H \
	vel=3d.slo3000.H bkgimg=3d.bkim.perm_hy.10.H \
	lloyd=0 vsal=7000 nrefvel=1 zmin=0 zmax=1000 nhx=11 ohx=-100 jhx=1 nhy=11 ohy=-100 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 extrapolator=pspi > /dev/null

3dpermig50:
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+1/3000" > 3d.slo3000.H
	Window3d < 3d.csou50.perm.H squeeze=n n4=1 50 f4=625 > s1.H
	Window3d < 3d.crec50.perm.H squeeze=n n4=1 50 f4=625 > r1.H
	./bin/OWmig3d_areal.x sou=s1.H 3d.csou50.perm.H rec=r1.H 3d.crec50.perm.H \
	vel=3d.slo3000.H bkgimg=3d.bkim.perm.50.H \
	lloyd=1 nrefvel=1 zmin=0 zmax=1000 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 extrapolator=pspi > /dev/null

3ddottst:
#	Window3d < 3d.csou10.perm.H squeeze=n n4=1 f4=312 > s1.H
#	Window3d < 3d.crec10.perm.H squeeze=n n4=1 f4=312 > r1.H
#	Math file1=3d.slo3000.H exp="-file1*0.1" > 3d.dslo300.H
	../Tese/bin/OW3dwemva_opers.x sou=s1.H rec=r1.H \
	vel=3d.slo3000.H deltaI=square5.Dimag.renorm.H dot=1 \
	lloyd=1 nrefvel=1 zmin=0 zmax=1000 nhx=11 ohx=-100 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 extrapolator=ffd born2=0 renorm=0 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 &

# Test Parallel
#	nohup Parallel command=../scr1/S_G_cam.x \
#	               files="data.$*.H,../elf.L7f.sinit.grid20.cam.H,../elf.L7f.ref_vel.grid20.H,../$@" \
#	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
#		     axis="6,4,4,6" file_type="COPY,COPY,COPY,COPY" nblock=4 device=eth0 \
#		     pverb=2 mach_mfile=mfile par=../sisfull.par

3dperm_par:
#	Transp plane=15 maxsize=1000 < 3d_rotrec_2d.H | Transp plane=24 maxsize=1000 > 3dr.H ; echo n6=1 o4=0 d4=20 >> 3dr.H
#	Transp plane=15 maxsize=1000 < 3d_rotsou_2d.H | Transp plane=24 maxsize=1000 > 3ds.H ; echo n6=1 o4=0 d4=20 >> 3ds.H
	Wavelet n1=351 d1=0.004 wavelet=ricker2 phase=0 fund=24 order=8 tdelay=.7 domain=time |Scale dclip=-1 | \
	Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y | Window3d f5=1 | Cabs | Smooth rect1=11 > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > 11111.H
#	Spike n3=201 d3=10 o3=0 n2=1001 d2=20 o2=-10000 n1=1001 d1=20 o1=-10000 > j1.H
#	Math file1=j1.H exp="file1*0+3000" > 3d.vel3000.H
	nohup Parallel command=./bin/OW3d_arealsht.x \
	files="3drr.H,3dss.H,11111.H,3d.vel3000.H,3d.crec10par.perm.H,3d.csou10par.perm.H" \
	tags="rec_cig,sou_cig,sou,vel,rec_out,sou_out" \
	usage="INPUT,INPUT,INPUT,INPUT,OUTPUT,OUTPUT" \
	axis="7,7,4,4,5,5" file_type="DISTRIBUTE,DISTRIBUTE,COPY,COPY,DISTRIBUTE,DISTRIBUTE" nblock=625 device=eth0 par=3dperm10_par.P \
	pverb=2 mach_mfile=mfile