##################################################
#
# BP - benchmark
#
#	< shots.small.H Window3d min3=6000 n3=2 400 min2=-8000 max1=9 | Mute vmute=1500 tmute=.5 | Bandpass flo=4 fhi=40 | Interp d1out=0.008 > j1.H; <j1.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \

bpait_recip_left.H:
	< shots.small.H Window3d max3=30000 max1=7.74 min2=-8000 j2=2 | Mute vmute=1500 tmute=.1 > j1.H 
	< j1.H Transp plane=23 | Interp d2out=25 | Transp plane=23 > x.H
	< x.H Create3d keyname1=offset keyname2=sx verb=1 > junk.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=3000 og2=-4000 dg2=12.5 key1=offset ng1=161 og1=-8000 dg1=50 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H
	Interp d2out=25 maxsize=1000 <junk3.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d min2=0.002 > junk4.H
	Merge axis=2 space=n junk31.H junk4.H > junk5.H
	Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H
	Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H
	Sort3d < junk7.H nkeys=2 key2=sx ng2=600 og2=50 dg2=50 key1=offset ng1=641 og1=-8000 dg1=25 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@

bpait.prep.left.H:
	< bpait_recip_left.H Window3d n3=400 | Transf is2d=y wei=y shotmig=y f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Pad beg1=670 end1=670 > $@

bpait_recip.H:
	< shots.small.H Window3d min3=45000 max3=59000 max1=7.74 min2=-8000 | Mute vmute=1500 tmute=.07 | ../Tese/bin/OW_mutmud.x  > j1.H 
	< j1.H Window3d j2=2 | Transp plane=23 | Interp d2out=25 | Transp plane=23 > x.H
	< x.H Create3d keyname1=offset keyname2=sx verb=1 > junk.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=1760 og2=41000 dg2=12.5 key1=offset ng1=161 og1=-8000 dg1=50 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H
	Interp d2out=25 maxsize=1000 <junk3.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d min2=0.002 > junk4.H
	Merge axis=2 space=n junk31.H junk4.H > junk5.H
	Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H
	Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H
	Sort3d < junk7.H nkeys=2 key2=sx ng2=280 og2=45000 dg2=50 key1=offset ng1=641 og1=-8000 dg1=25 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@

bpait.prep.H:
	< shots.small.H Window3d min3=6000 n3=450 min2=-8000 max1=7.74 | Mute vmute=1500 tmute=.07 > j1.H
	<j1.H /net/koko/bob/SEP_BOB/bin/Transf is2d=y wei=y shotmig=y f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=670 end1=670 > $@

bpait.prep.all.H:
	< shots.small.H Window3d min2=-8000 max1=7.74 | Mute vmute=1500 tmute=.07 > j1.H
	<j1.H /net/koko/bob/SEP_BOB/bin/Transf is2d=y wei=y shotmig=y f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=670 end1=670 > $@

bpait.prep.mudmut.H:
	< bpait_recip.H Transf is2d=y wei=y shotmig=y f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Pad beg1=320 end1=320 > $@

bpait.source.H:
	< shots.small.H Window3d n3=1 n2=1 min2=-2000 max1=1.8 | Pad n1out=1290 > x.H
	echo o1=-1.368 >> x.H
	< x.H Transf shotmig=y is2d=y wei=y \
	f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Window3d | Cabs > spec.H; Math file1=spec.H exp="file1*0." > phase.H
	Cmplx spec.H phase.H > j3.H; < j3.H Transp plane=13 > $@

bpait.source1.H:
	Spike n1=1290 k1=646 d1=0.006 o1=-3.87 | Bandpass flo=5 fhi=28 > j1.H
	< j1.H Transf is2d=y wei=y shotmig=y \
	f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > j2.H
	< j2.H Window3d | Cabs > spec.H; Math file1=spec.H exp="file1*0." > phase.H
	Cmplx spec.H phase.H > j3.H; < j3.H Transp plane=13 > $@
!	< j2.H Transp plane=53 > $@

bpait.vel.orig.H:
	Cp /data/bpait/vel.H bpait.vel.orig.H

bpait.vel.edit.H:
	< bpait.vel.orig.H Window3d max1=11000 max2=16000 30000 | ../Tese/bin/OW_editvel.x gradz=0 2.15 v0=2000 vsalt=4510 vh2o=1495 > $@
	< bpait.vel.edit.H Window3d j2=2 | Pad beg2=2000 end2=2000 extend=1  > bpait.vel.edit.pad.H

bpait.vel.edit.transp.H:
	< bpait.vel.edit.pad.H Transp plane=12 reshape=1,3 > $@

bpait.vel.edit.transp.lloyd.H:
	< bpait.vel.edit.transp.H Lloyd_vel nref=4 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

bpait.vel.pad.H:
	< bpait.vel.orig.H Window3d max1=9000 | Pad beg2=2000 end2=2000 extend=1  > $@

bpait.vel.transp.H:
	< bpait.vel.orig.H Window3d max1=9000 | Transp plane=12 reshape=1,3 | Pad beg1=2000 end1=2000 extend=1  > $@

bpait.vel.transp.lloyd.H:
	Lloyd_vel < bpait.vel.transp.H nref=4 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

bpait.vel.edit.lloyd.H:
	< bpait.vel.edit.H Window3d max1=9000 | Transp plane=12 reshape=1,3 | Lloyd_vel nref=4 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

bpait.vel.mud.H:
	Window3d min2=45000 max2=55000 < bpait.vel.orig.H max1=2200 min1=1500 | Clip clip=1745 | Smooth rect1=5 rect2=7 > mask1.H
	< bpait.vel.orig.H Window3d min2=37000 | ../Tese/OW_editvel_mud1.x mask1=mask1.H maskout=maskout.H> xxx.H 
	Window3d min2=50000 max2=52000 < xxx.H max1=2000 min1=1700 | Clip clip=2060 | Smooth rect1=7 rect2=7 > mask1.H
	< xxx.H ../Tese/OW_editvel_mud1.x mask1=mask1.H maskout=maskout.H | Window3d j2=2 | Pad beg2=1000 end2=1000 extend=1 > $@
	< $@ Transp plane=12 reshape=1,3 > bpait.vel.mud.shtp.H

bpait.vrms.orig.H:
	< bpait.vel.orig.H Window3d min2=40000 | ../Tese/OW_vint2rms.x | Pad extend=1 beg2=3000 end2=3000 > $@

bpait.vrms.mud.H:
	< bpait.vel.mud.H Window3d min1=40000 max1=65000 | Transp | ../Tese/OW_vint2rms.x | Pad extend=1 beg2=3000 end2=3000 > $@

bpait.vel.mud.lloyd.H:
	Transp plane=12 reshape=1,3 < bpait.vel.mud.H | \
	Lloyd_vel nref=4 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

bpait.vel.mud.lloyd1v.H:
	< bpait.vel.orig.H Window3d min2=40000 | ../Tese/OW_editvel_mud.x vh2o=1486 | Transp plane=12 reshape=1,3 | \
	Lloyd_vel nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

bpait.mig.H: # With 401 shots separated 0.01km correct velocity
	Window3d < bpait.prep.left.H squeeze=n j4=2 > x.H
	nohup /net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=x.H D_d=bpait.source.H V=bpait.vel.edit.transp.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=50 \
	par=bpait.gosis.P taper=y taper_pcent=2 marine=y minx=-24750. maxx=46000. & #restart=y

bpait.mig1.H: # With 401 shots separated 0.01km correct velocity
	Window3d < bpait.prep.left.H squeeze=n f4=1 j4=2 > x1.H
	/net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=x1.H D_d=bpait.source.H V=bpait.vel.edit.transp.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=100 \
	par=bpait.gosis.P taper=y taper_pcent=2 marine=y minx=-24750. maxx=46000. & #restart=y

bpait.ang.left.H:
#	Math file1=bpait.mig.H file2=bpait.mig1.H exp="file1+file2"|Window3d min1=0 max1=16000 |Transp plane=13|Scale scale_to=1 \
#	|Interp d2out=25 > bpait.mig.left.H;
	< bpait.mig.left.H Window3d min1=4000 | ../Tese/bin/OW_3dslant.x oax=-20 oah=-40 dax=2 dah=2 nx=7 nodes=10 teste=1 > $@

bpait.mig.left.rot:
	< bpait.mig.left.H Window3d min1=0 4000 max1=9000 min3=6000 max3=9000 | ../Tese/bin/OW_3dslant.x oax=-30 oah=-30 dax=1.5 dah=1.5 nx=15 nodes=14 rot= 1 > bpait.mig.left.rot.sd.H
#	< bpait.mig.left.H Window3d min1=4000 | ../Tese/bin/OW_3dslant.x oax=-30 oah=-30 dax=1.5 dah=1.5 nx=15 nodes=14 rot=1 rec=1 > bpait.mig.left.rot.rd.H

bpait.mig.all.H: # With 401 shots separated 0.01km correct velocity
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=bpait.prep.all.H D_d=bpait.source.H V=bpait.vel.transp.H R=$@ \
	pverb=2  nx_shots=1348 ox_shots=50 \
	par=bpait.gosis.P taper=y taper_pcent=2 marine=y minx=-17000. maxx=78000. & #restart=y

bpait.mig.mud.H: # With 401 shots separated 0.01km correct velocity
#	Window3d < bpait.prep.all.H squeeze=n min4=45000 max1=3000 > xxx.H
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=bpait.prep.mudmut.H D_d=bpait.source.H V=bpait.vel.mud.shtp.H R=$@ \
	pverb=2  nx_shots=280 ox_shots=45000 \
	par=bpait.gomud.P taper=y taper_pcent=2 marine=y minx=28000. maxx=71000. & #restart=y

bpait.mig.mud.vorig.H:
#	Window3d < bpait.prep.all.H squeeze=n min4=45000 max1=3000 > xxx.H
	Window3d < bpait.vel.transp.H squeeze=n max3=5000 j1=2 > xxx.H
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=bpait.prep.mudmut.H D_d=bpait.source.H V=xxx.H bpait.vel.transp.H R=$@ \
	pverb=2  nx_shots=280 ox_shots=45000 \
	par=bpait.gomud.P taper=y taper_pcent=2 marine=y minx=28000. maxx=71000. & #restart=y

bpait.mig.mud.short.H: # With 401 shots separated 0.01km correct velocity
	Window3d < bpait.prep.all.H squeeze=n min1=-1000 max1=1000 min4=54000 | Pad beg1=200 > yyy.H
	@date
	/net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=bpait.prep.mudmut.H D_d=bpait.source.H V=bpait.vel.mud.H R=$@ \
	pverb=2  nx_shots=269 ox_shots=54000 \
	par=bpait.gomud1.P taper=y taper_pcent=2 marine=y minx=40000. maxx=70000. & #restart=y

#bpait.mig.mud.H: # With 401 shots separated 0.01km correct velocity
#	Window3d < bpait.prep.all.H squeeze=n min4=45000 max1=3000 > xxx.H
#	@date
#	nohup Parallel command=/net/koko/claudio/Research/Others/ShotMig/ShotMig \
#	files="xxx.H,bpait.source.H,bpait.vel.transp.H,$@" \
#	tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" nblock=449 axis=4,4,4,6 \
#	file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 device=eth0 \
#	pverb=2  nx_shots=449 ox_shots=45000 \
#	par=bpait.gomud.P taper=y taper_pcent=2 marine=y minx=28000. maxx=71000. & #restart=y

bpait.mig.cluster.H: # With 401 shots separated 0.01km correct velocity
	@date
	Parallel command=/net/koko/claudio/Research/Others/ShotMig/bin/LINUX32/ShotMigNew \
	files="bpait.prep.all.H,bpait.source.H,bpait.vel.transp.H,$@" \
	tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" nblock=1348 axis=4,4,4,6 \
	file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 device=eth0 \
	pverb=2  nx_shots=1348 ox_shots=50 \
	par=bpait.gosis.P taper=y taper_pcent=2 marine=y minx=-17000. maxx=78000. & #restart=y

bpait.mig.mud.mute.window.H:
#	< bpait.mig.mud.mute.H Window3d min1=45000 n1=881 | Transp plane=13 | Interp d2out=12.5 maxsize=2000 > 1.H
	< 1.H ../Tese/OW_windowing.x hor1=yel3dmud hor2=red3dmud wavz=9 pickr=pickmud.H > $@

bpait_rot_mud:
	< bpait.mig.mud.H Window3d max5=4000 min1=45000 max1=60000 | Transp plane=12 reshape=2,3,6 > bpait.mig.temp.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.sou.H par=bpait_rot.mud.P positive_rot=0
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.rec.H par=bpait_rot.mud.P positive_rot=1
	Rm bpait.mig.temp.H
sd

bpait_rot:
	< bpait.mig.H Window3d min1=5000 max1=15000 | Transp plane=12 reshape=2,3 | Interp d3out=12.5 maxsize=2000 > bpait.mig.temp.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.sou.H par=bpait_rot.P positive_rot=0
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.rec.H par=bpait_rot.P positive_rot=1
	Rm bpait.mig.temp.H

#########################################
#
# encoding all cig's at once - conv.rand.
# when using horizons min1=0 max1=21000
# for initial migration
#
bpait00: # correct velocity areal shot random encode in w
	< bpait.rot.rec.H Window3d min2=5000 max2=15000 | Transp plane=23 > junk1.H
	< bpait.rot.sou.H Window3d min2=5000 max2=15000 | Transp plane=23 > junk2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec.H sou_out=bpait.sou.H \
	padtraces=101 encode=1 random=1 encw=1 eps=1 injection=1 nreal=10 ngap=1 zmin=5000 zmax=8990 node=16 \
	> /dev/null 

bpait.mig.00.H: # correct velocity areal shot random encode in w
	< bpait.rec.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=1 61 oh=0 -.2 60 dh=12.5 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01: # correct velocity areal shot random encode in w
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec1.H sou_out=bpait.sou1.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal.P > /dev/null 
#oh= -750. dh=12.5
bpait.mig.01.H: # no reflection on h
	< bpait.rec1.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait.mig.02.H: # reflection on h source
	< bpait.rec1.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01a: # reflection in h receiver
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out= bpait.rec1a.H sou_out= bpait.sou1a.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal.P pickr=bpait.pickr.H picks=bpait.picks.H > /dev/null 

bpait.mig.03.H: # correct velocity areal shot random encode in w
	< bpait.rec1a.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1a.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01b: # reflection in h both
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec1b.H sou_out=bpait.sou1b.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal.P pickr=bpait.pickr.H picks=bpait.picks.H > /dev/null 

bpait.mig.03b.H: # correct velocity areal shot random encode in w
	< bpait.rec1b.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1b.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01c: # no rotation and reflection in h source
	< bpait.mig.H Window3d min1=5000 max1=15000 | Transp plane=13 | Interp d2out=12.5 maxsize=2000 > junk0.H
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk0.H sou_cig=junk0.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec1c.H sou_out=bpait.sou1c.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal1c.P pickr=bpait.pickr1c.H picks=bpait.picks1c.H > /dev/null 

bpait.mig.03c.H: # correct velocity areal shot random encode in w
	< bpait.rec1c.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1c.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01d: # rotation, w/no reflection in h source window 13
	< bpait.rot.rec.H Transp plane=23 > junk1.H
	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec1d.H sou_out=bpait.sou1d.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal1d.P pickr=bpait.pickr1d.H picks=bpait.picks1d.H > /dev/null 

bpait.mig.03d.H: # correct velocity areal shot random encode in w
	< bpait.rec1d.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1d.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait01e: # rotation, w/no reflection in h source window 7
	< bpait.rot.rec.H Transp plane=23 > junk1.H
	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec1e.H sou_out=bpait.sou1e.H soucig=bpait.souw.H \
	reccig=bpait.recw.H par=bpait_areal1d.P pickr=bpait.pickr1d.H picks=bpait.picks1d.H > /dev/null 

bpait.mig.03e.H: # correct velocity areal shot random encode in w
	< bpait.rec1e.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.sou1e.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=61 jh=2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait02: # correct velocity areal shot random encode in w
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec2.H sou_out=bpait.sou2.H \
	par=bpait_areal2.P > /dev/null 

bpait.mig.04.H: # correct velocity areal shot random encode in w
	< bpait.rec2.H Transp plane=13 > bpait.rec.temp2.H
	< bpait.sou2.H Transp plane=13 > bpait.sou.temp2.H
	< bpait.rec.temp2.H ~claudio/research.ln/Prst_ExpRefl/Tese/InvVel/bin/OWmig2d_areal.x sou=bpait.sou.temp2.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.transp.lloyd.H zmin=5000 zmax=8990. nh=1 21 jh=1 2 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait03: # 
#	< bpait.rot.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.sou.H Transp plane=23 > junk2.H
	/net/koko/claudio/Research/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source.H \
	vel=bpait.vel.edit.pad.H vref=bpait.vel.edit.lloyd.H rec_out=bpait.rec3.H sou_out=bpait.sou3.H \
	par=bpait_areal3.P > /dev/null 

###############################################
#
# Mud area (right part of the model)
#
bpait.ang.mud.H:
	Window3d < bpait.mig.mud.H min1=45000 max1=57000 max5=4000 | Transp plane=13 | Transp plane=34 \
	| Scale scale_to=1 | Pad beg2=120 end2=120 end1=100 > xx.H
	< xx.H OFF2ANG amax=45 eps=2 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=4000 > $@

bpait.ang.pick.H:
	Window3d n3=474 min2= -35 n2 =219 < bpait.ang.mud.H > xxx.H; ../Tese/bin/OW_autopick.x \
	cig=xxx.H hor1=red1d.txt hor2= blu1d.txt mask=maskpick.H nwind=25 > $@
	Math file1=xxx.H file2=maskpick.H exp="file1*file2" > temp2.H
	< temp2.H ../Tese/bin/off2ang_sd.x hmin=-500 hmax=500 dhx=25 adj=0 style=1 > bpait.mig.mud.mask.H

#bpait.mig.mud.window.H:
#	Window3d < bpait.mig.mud.H min1=45000 max1=58000 max5=5000 f3=14 n3=12 | Transp plane=13 > 1.H
#	< 1.H ../Tese/bin/OW_windowing.x hor1=HH1 hor2=HH2 wavz=15 pickr= pickmud.H mask=maskwind.H > 2.H
#	Smooth rect1=31 < maskwind.H > 3.H; Math file1=2.H file2=3.H exp="file1*file2" > $@

bpait.mig.mud.window%.H:
	Window3d < bpait.mig.mud.H min1=48000 max1=55000 max5=5000 f3=14 n3=12 | Transp plane=13 > 1.H
	< 1.H ../Tese/bin/OW_windowing.x hor1=HH1 hor2=HH2 wavz=$* pickr=pickmud.H mask=maskwind.H > $@ #2.H
#	Smooth rect1=31 < maskwind.H > 3.H; Math file1=2.H file2=3.H exp="file1*file2" > $@

bpait.mig.mud.gaus11.H:
#	Window3d < bpait.mig.mud.H min1=45000 max1=58000 max5=5000 f3=14 n3=12 | Transp plane=13 > 1.H
	Window3d < bpait.mig.mud.pp11.H > 1.H
	< 1.H ~/research.ln/lin1/lin1/Prog/Gausfil.x \
	sigma10=5. sigma20=5. > $@

bpait.mig.mud.picks11.H:
#	Window3d < bpait.mig.mud.H min1=45000 max1=58000 max5=5000 f3=14 n3=12 | Transp plane=13 > 1.H
	Window3d < bpait.mig.mud.pp11.H > 1.H
	< 1.H ~/research.ln/lin1/lin1/Prog/Nonlpick.x \
	qpick=bpait.mig.mud.qpick11.H pick=bpait.mig.mud.pick11.H > $@

bpait.mig.mud.pck1.window%.H:
	Window3d < bpait.mig.mud.H min1=45000 j1=2 n1=295 max5=4000 min3=-200 n3=17 | Transp plane=13 > 1.H
	< 1.H ../Tese/bin/OW_windowing.x hor1=HH2 red1.out yel1.out hor2= red1.out wavz=$* pickr= pickmud.H mask= maskwind.H > $@ #2.H
#	Smooth rect1=31 < maskwind.H > 3.H; Math file1=2.H file2=3.H exp="file1*file2" > $@

#####
#
# Seminar 090130
bpait.mig.mud.pck2.window%.H:
	Window3d < bpait.mig.mud.H min1=48000 n1=281 max5=4000 min3=-200 n3=17 | Transp plane=13 > 1.H
	< 1.H ../Tese/bin/OW_windowing.x hor1=yel1.out hor2=red1.out wavz=$* pickr= pickmud.H mask=maskwind.H > 2.H
	Smooth rect1=31 < maskwind.H > 3.H; Math file1=2.H file2=3.H exp="file1*file2" > $@

bpait.mig.dips.H:
	< bpait.mig.mud.pck2.window21.H Window3d min2=0 n2=1 | Dip method=sergey_s nl_iter=30 niter=30 eps=0.1 > j1.H
	Smooth rect1=5 < j1.H > j2.H ; Math file1=j2.H exp="file1*20/0.8" | Clip chop=greater clip=20 > $@

bpaitrot_mud%:
	< bpait.mig.mud.pck2.window$*.H Transp plane=23 > bpait.mig.temp.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.pck2.sou$*.H \
	par=bpait_rot.mud.P positive_rot=0
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.pck2.rec$*.H \
	par=bpait_rot.mud.P positive_rot=1
	Rm bpait.mig.temp.H

bpaitmud%:
	Window3d < bpaitmig_rot_rd$*.H max1=3000 > j1.H
	Window3d < bpaitmig_rot_sd$*.H max1=3000 > j2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x \
	rec_cig=j1.H sou_cig=j2.H sou=bpait.source1.H bpait.source.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H weight= weight1.H \
	rec_out=bpait.areal.rec.mud$*.H sou_out=bpait.areal.sou.mud$*.H \
	padtraces=101 encode=0 1 random=1 encz=1 eps=1 injection=1 nreal=1 0 ngap=40 50 1 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig.%.mud.H:
	< bpait.areal.rec.mud$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mud$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=13 oh=-150 1 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpaitrotspd_%:
	<bpait.mig.mud.pck2.window$*.H Window3d max3=53000 max1=3000 | \
	../Tese/bin/OW_3dslant.x nax=31 nah=51 oax=-30 oah=-50 dax=2 dah=2 \
	nxw=21 rot=1 upshft=0 nodes=16 zmin=0 > j1sd$*.H
	Bandpass flo=0.01 fhi=0.2 < j1sd$*.H > bpaitmig_rot_sd$*.H
	<bpait.mig.mud.pck2.window$*.H Window3d max3=53000 max1=3000 | \
	../Tese/bin/OW_3dslant.x nax=31 nah=51 oax=-30 oah=-50 dax=2 dah=2 \
	nxw=21 rot=1 upshft=0 rec=1 nodes=16 zmin=0 > j1rd$*.H
	Bandpass flo=0.01 fhi=0.2 < j1rd$*.H > bpaitmig_rot_rd$*.H

1bpaitmudspd_%:
#	Wavelet n1=750 dt=0.006 wavelet=ricker2 fund=10 | Transf f_min=1 f_min1=4 f_max1=32 f_max=40 \
#	is2d=y wei=y shotmig=y maxsize=1000 > 1source.H
	Spike n1=751 n2=1 n3=1 d1=0.004 d2=1 d3=1 o1=-1.5 o2=0 o3=0 nsp=1 k1=376 | Bandpass flo=4 fhi=48 \
	| Transf f_min=1 f_min1=4 f_max1=40 f_max=48 is2d=y wei=y shotmig=y maxsize=1000 > 1source.H
	Window3d < bpaitmig_rot_rd$*.H max1= 3000 | ../Tese/bin/OW_hgain.x w=1. > j1.H
	Window3d < bpaitmig_rot_sd$*.H max1= 3000 | ../Tese/bin/OW_hgain.x w=1. > j2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x \
	rec_cig=j1.H sou_cig=j2.H sou=1source.H ../Tese/40011data_source.H bpait.source1.H bpait.source.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H weight= weight1.H \
	rec_out=bpait.areal.rec.mudspd$*.H sou_out=bpait.areal.sou.mudspd$*.H \
	padtraces=101 encode=1 encz=1 eps=1 injection=1 nreal=1 ngap=40 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig.%.mudspd.H:
	< bpait.areal.rec.mudspd$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mudspd$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=17 oh=-200 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait.mig1.%.mudspd.H:
	< bpait.areal.rec.mudspd$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mudspd$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=17 oh=-200 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpaitrotspdsh_%:
	Window3d max2=53000 max1=3000 < bpait.mig.dips.H > 1dips.H
	<bpait.mig.mud.pck2.window$*.H Window3d max3=53000 max1=3000 | \
	../Tese/bin/OW_3dslant.x nax=31 nah=51 oax=-30 oah=-50 dax=2 dah=2 \
	nxw=21 rot=1 upshft=1 nodes=16 zmin=0 dips=1dips.H nxchunk=2 > j1sdsh$*.H
	Bandpass flo=0.01 fhi=0.2 < j1sdsh$*.H > bpaitmig_rot_sdsh$*.H
	<bpait.mig.mud.pck2.window$*.H Window3d max3=53000 max1=3000 | \
	../Tese/bin/OW_3dslant.x nax=31 nah=51 oax=-30 oah=-50 dax=2 dah=2 \
	nxw=21 rot=1 upshft=1 rec=1 nodes=16 zmin=0 dips=1dips.H nxchunk=2 > j1rdsh$*.H
	Bandpass flo=0.01 fhi=0.2 < j1rdsh$*.H > bpaitmig_rot_rdsh$*.H

2bpaitmudspdsh_%:
#	Wavelet n1=750 dt=0.006 wavelet=ricker2 fund=10 | Transf f_min=1 f_min1=4 f_max1=32 f_max=40 \
#	is2d=y wei=y shotmig=y maxsize=1000 > 1source.H
	Spike n1=751 n2=1 n3=1 d1=0.004 d2=1 d3=1 o1=-1.5 o2=0 o3=0 nsp=1 k1=376 | Bandpass flo=4 fhi=48 \
	| Transf f_min=1 f_min1=4 f_max1=40 f_max=48 is2d=y wei=y shotmig=y maxsize=1000 > 1source.H
	Window3d < bpaitmig_rot_rdsh$*.H max1= 3000 > j1.H
	Window3d < bpaitmig_rot_sdsh$*.H max1= 3000 > j2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x \
	rec_cig=j1.H sou_cig=j2.H sou=1source.H bpait.source1.H bpait.source.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H weight= weight1.H \
	rec_out=bpait.areal.rec.mudspdsh$*.H sou_out=bpait.areal.sou.mudspdsh$*.H \
	padtraces=101 encode=0 1 random=1 encz=1 eps=1 injection=1 nreal=1 0 ngap=40 50 1 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig.%.mudspdsh.H:
	< bpait.areal.rec.mudspdsh$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mudspdsh$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=17 oh=-200 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait.mig1.%.mudspdsh.H:
	< bpait.areal.rec.mudspdsh$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mudspdsh$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=17 oh=-200 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

###############################
#
# 2 horizons
#
# Refer to 
# #####
# #
# # Seminar 090130
# bpait.mig.mud.pck2.window%.H:
# to check the windowing flow
#
bpait2mud%:
	Spike n1=751 n2=1 n3=1 d1=0.004 d2=1 d3=1 o1=-1.5 o2=0 o3=0 nsp=1 k1=376 | Bandpass flo=4 fhi=52 \
	| Transf f_min=1 f_min1=4 f_max1=52 f_max=60 is2d=y wei=y shotmig=y maxsize=1000 > 1source.H
	Window3d < bpait.rot.mud.pck2.rec$*.H max1=3000  | Transp plane=23 | ../Tese/bin/OW_hgain.x w=1. > j1.H
	Window3d < bpait.rot.mud.pck2.sou$*.H max1=3000  | Transp plane=23 | ../Tese/bin/OW_hgain.x w=1. > j2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x \
	rec_cig=j1.H sou_cig=j2.H sou=1source.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H \
	rec_out=bpait.areal.rec.rot2h.mud$*.H sou_out=bpait.areal.sou.rot2h.mud$*.H \
	padtraces=101 encode=1 random=1 encz=1 eps=1 injection=1 nreal=2 ngap=40 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig2.%.mud.H:
	< bpait.areal.rec.rot2h.mud$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.rot2h.mud$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=17 oh=-200 1 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait.Dimag2h.mud.%.H:
	< bpait.areal.rec.rot2h.mud$*.H Transp plane=13 > r1.H
	< bpait.areal.sou.rot2h.mud$*.H Transp plane=13 > s1.H
	< r1.H ../Tese/bin/OWwemva_opers.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H deltaS=bpait.Dslow.mud.H \
	xmin=zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 adj=0 \
	traces_in_taper=101 every=0 encode=1 node=8 > $@ 

bpait.Dslow2h.mud.%.H:
	< bpait.areal.rec.rot2h.mud$*.H Transp plane=13 > r1.H
	< bpait.areal.sou.rot2h.mud$*.H Transp plane=13 > s1.H
	< bpait.Dimag2h.mud.$*.H Scale scale_to=1 > j1.H
	< r1.H ../Tese/bin/OWwemva_opers.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H deltaI=j1.H \
	zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 adj=1 \
	traces_in_taper=101 every=0 encode=1 node=8 > $@ 

slown%:
	< bpait.areal.rec.rot2h.mud$*.H Window3d squeeze=n f1= 25 n1= 16 96 |Transp plane=13 > r1.H
	< bpait.areal.sou.rot2h.mud$*.H Window3d squeeze=n f1= 25 n1= 16 96 |Transp plane=13 > s1.H
	< r1.H ../Tese/bin/OW_lbfgsc_areal.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H grad=bpait.Grad.mud.rot2h.slow$*.H \
	Dimage=bpait.DSO.mud.rot2h.slow$*.H bkgimg=bpait.Bkgimg.mud.rot2h.slow$*.H Velout=bpait.Velout.mud.rot2h.slow$*.H \
	zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 niter=5 dso=1 inversion=1 damp=1000 horv0= dv0.out\
	traces_in_taper=51 padtraces=51 node=16 bounds=1 vmin=1400 vmax=5000 vsal=5000 nrefvel=3 xmin=48000 xmax=55000 \
	> bpait.inv.mud.rot2h.slow%.H
# 	dxk= .2 nxk= 150 dzk= .012 nzk= 200 >
#
###############################


###############################################
#
# Inversion
#
#bpait.Drho.mud.H:
#	Window3d < bpait.vel.orig.H min1=48000 max1=8990 min2=2987.5 n2=1123 > j1.H
#	Window3d < bpait.vel.edit.H min1=48000 max1=8990 min2=2987.5 n2=1123 > j2.H
#	Math file1=j1.H file2=j2.H exp="file2/file1-1" > $@

bpait.Dslow.mud.H:
	< bpait.vel.orig.H Pad extend=1 end2=3000 | Window3d j2=2 min2=12000 max2=92450 > j1.H
	Window3d < bpait.vel.mud.H > j2.H
	Math file1=j1.H file2=j2.H exp="file2-file1" > $@

bpait.Dimag.mud.%.H:
	Math file1=bpait.areal.rec.mudspd$*.H exp="file1/190.595" > j1.H
	< j1.H Window3d squeeze=n f1= 50 n1= 5 |Transp plane=13 > r1.H
	Math file1=bpait.areal.sou.mudspd$*.H exp="file1/189.365" > j1.H
	< j1.H Window3d squeeze=n f1= 50 n1= 5 |Transp plane=13 > s1.H
	< r1.H ../Tese/bin/OWwemva_opers.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H deltaS=bpait.Dslow.mud.H \
	xmin=zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 adj=0 \
	traces_in_taper=101 every=0 encode=1 node=16 > $@ 

bpait.Dslow1.mud.%.H:
	Math file1=bpait.areal.rec.mudspd$*.H exp="file1/190.595" > j1.H
	< j1.H Window3d squeeze=n f1= 50 n1= 5 |Transp plane=13 > r1.H
	Math file1=bpait.areal.sou.mudspd$*.H exp="file1/189.365" > j1.H
	< j1.H Window3d squeeze=n f1= 50 n1= 5 |Transp plane=13 > s1.H
	< bpait.Dimag.mud.$*.H Scale scale_to=1 > j1.H
	< r1.H ../Tese/bin/OWwemva_opers.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H deltaI=j1.H \
	zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 adj=1 \
	traces_in_taper=101 every=0 encode=1 node=16 > $@ 

bpait.inv.mud.spdsh.slow%.H: #
	Transp < bpait.areal.rec.mudspdsh$*.H plane=13 > r1.H
	Transp < bpait.areal.sou.mudspdsh$*.H plane=13 > s1.H
	< r1.H ../Tese/bin/OW_lbfgsc_areal.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H grad=bpait.Grad.mud.spdsh.slow$*.H \
	Dimage=bpait.DSO.mud.spdsh.slow$*.H bkgimg=bpait.Bkgimg.mud.spdsh.slow$*.H Velout=bpait.Velout.mud.spdsh.slow$*.H\
	zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 niter=1 dso=1 inversion=1 vmin=1400 vmax=5500 \
	traces_in_taper=0 node=16 > $@

bpait.inv.mud.spd.slow%.H: #
	< bpait.areal.rec.mudspd$*.H Window3d squeeze=n f1=50 n1=5 |Transp plane=13 > r1.H
	< bpait.areal.sou.mudspd$*.H Window3d squeeze=n f1=50 n1=5 |Transp plane=13 > s1.H
	< r1.H ../Tese/bin/OW_lbfgsc_areal.x sou=s1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H grad=bpait.Grad.mud.spd.slow$*.H \
	Dimage=bpait.DSO.mud.spd.slow$*.H bkgimg=bpait.Bkgimg.mud.spd.slow$*.H Velout=bpait.Velout.mud.spd.slow$*.H\
	zmin=0 zmax=3000 nh=17 oh=-200 jh=1 dh=25 wbott=0 niter=1 dso=1 inversion=1 vmin=1400 vmax=5500 \
	traces_in_taper=101 node=16 > $@


bpait00mud%:
#	< bpait.rot.mud.rec.H Transp plane=23 > junk1.H
#	< bpait.rot.mud.sou.H Transp plane=23 > junk2.H
	Window3d < bpait.mig.mud.window$*.H max1=4000 squeeze=n > 1.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x \
	rec_cig=bpait.mig.mud.window$*.H sou_cig=bpait.mig.mud.window$*.H sou=bpait.source1.H bpait.source.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H rec_out=bpait.areal.rec.mud$*.H sou_out=bpait.areal.sou.mud$*.H \
	padtraces=101 encode=0 random=1 encw=1 eps=1 injection=1 nreal=1 0 ngap=30 50 1 zmin=0 zmax= 3000 node=16 \
	> /dev/null 

bpait.mig.00.mud%.H:
#	< bpait.areal.rec.mud.H Window3d n4=1 squeeze=n | Transp plane=13 > bpait.rec.temp1.H
#	< bpait.areal.sou.mud.H Window3d n4=1 squeeze=n | Transp plane=13 > bpait.sou.temp1.H
	< bpait.areal.rec.mud$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.sou.mud$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax= 3000. nh=12 oh=-150 1 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H

bpait_rot_mud%:
	< bpait.mig.mud.window$*.H Transp plane=23 > bpait.mig.temp.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.sou$*.H \
	par=bpait_rot.mud.P positive_rot=0
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.rec$*.H \
	par=bpait_rot.mud.P positive_rot=1
	Rm bpait.mig.temp.H

bpait01mudrot%:
	< bpait.rot.mud.rec$*.H Transp plane=23 > junk1.H
	< bpait.rot.mud.sou$*.H Transp plane=23 > junk2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H rec_out=bpait.areal.rot.mud.rec$*.H sou_out=bpait.areal.rot.mud.sou$*.H \
	padtraces=101 encode=0 1 random=1 encw=1 eps=1 injection=1 nreal=1 0 ngap=30 41 50 1 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig.01.mud.rot%.H:
	< bpait.areal.rot.mud.rec$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.rot.mud.sou$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=12 oh=-150 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=14 > $@ 
	Rm bpait.*.temp1.H


bpait_spd_mud%:
	< bpait.mig.mud.window$*.H Transp plane=23 > bpait.mig.temp.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.sou$*.H \
	par=bpait_rot.mud.P positive_rot=0
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=bpait.mig.temp.H Cig_tilt=bpait.rot.mud.rec$*.H \
	par=bpait_rot.mud.P positive_rot=1
	Rm bpait.mig.temp.H

bpait01mudspd%:
	< bpait.rot.mud.rec$*.H Transp plane=23 > junk1.H
	< bpait.rot.mud.sou$*.H Transp plane=23 > junk2.H
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk2.H sou=bpait.source1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H rec_out=bpait.areal.rot.mud.rec$*.H sou_out=bpait.areal.rot.mud.sou$*.H \
	padtraces=101 encode=0 1 random=1 encw=1 eps=1 injection=1 nreal=1 0 ngap=30 41 50 1 zmin=0 zmax=3000 node=16 \
	> /dev/null 

bpait.mig.01.mud.spd%.H:
	< bpait.areal.rot.mud.rec$*.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.rot.mud.sou$*.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=3000. nh=12 oh=-150 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=14 > $@ 
	Rm bpait.*.temp1.H

#
#
###############################################
###############################################
#
# Testing rotation with slant - stack

bpait.rot_sd.H:
	Window3d < bpait.mig.mud.window35.H n1=650 | ../Tese/bin/OW_3dslant.x oax=-20 oah=-40 dax=1 dah=1 nx=7 rot=1 > $@

bpait.rot_rd.H:
	Window3d < bpait.mig.mud.window35.H n1=650 | ../Tese/bin/OW_3dslant.x oax=-20 oah=-40 dax=1 dah=1 nx=7 rot=1 rec=1 > $@

bpaitspd: # correct velocity areal shot random encode in w
	~claudio/research.ln/Prst_ExpRefl/Tese/bin/OW_arealsht.x rec_cig=bpait.rot_rd.H sou_cig=bpait.rot_sd.H sou=bpait.source1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H rec_out=bpait.areal.spd.mud.rec.H sou_out=bpait.areal.spd.mud.sou.H \
	padtraces=101 encode=0 random=1 encw=1 eps=1 injection=1 nreal=1 ngap=25 zmin=0 zmax=4000 node=16 \
	> /dev/null 

bpait.mig.spd.H: # correct velocity areal shot random encode in w
	< bpait.areal.spd.mud.rec.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.spd.mud.sou.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=4000. nh=12 oh=-150 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=14 > $@ 
	Rm bpait.*.temp1.H

bpait.mig.spd1.H: # correct velocity areal shot random encode in w
	< bpait.areal.spd.mud.rec.H Transp plane=13 > bpait.rec.temp1.H
	< bpait.areal.spd.mud.sou.H Transp plane=13 > bpait.sou.temp1.H
	< bpait.rec.temp1.H ~claudio/research.ln/Prst_ExpRefl/Tese/bin/OWmig2d_areal.x sou=bpait.sou.temp1.H \
	vel=bpait.vel.mud.H vref=bpait.vel.mud.lloyd.H zmin=0 zmax=4000. nh=12 oh=-150 jh=1 wbott=0 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 > $@ 
	Rm bpait.*.temp1.H
###############################################
#
# Inversion
#
bpaitmigvlow:
	Spike n1=1290 d1=0.006 o1=-3.87 nsp=1 k1=646 | Bandpass flo=4 fhi=30 \
	| Transf is2d=y wei=y shotmig=y f_min=1 f_min1=4 f_max1=28 f_max=34 maxsize=1000 > bpait.1source.H
	Window3d < bpait.prep.left.H j4=2 squeeze=n min1=-14000 n1=1121 > bpaitxxx.H
	< bpaitxxx.H ../ShtPSPI/bin/OWmig2d_par.x sou=bpait.1source.H vel=bpait.vel.edit.pad.H \
	vref=bpait.vel.edit.transp.lloyd.H image=bpait.mig.left.vlow.H xmin=-16000 xmax=34100 zmin=5000 zmax=10000 \
	encode=0 node=16 nh=65 oh=-800. dh=25. > /dev/null

#	Window3d < bpait.prep.left.H squeeze=n j4=2 > x.H

bpait.mig.left1.H: # With 401 shots separated 0.01km correct velocity
	nohup /net/koko/claudio/Research/Others/ShotMig/bin/LINUX64/ShotMig \
	D=bpaitxxx.H D_d=bpait.1source.H V=bpait.vel.edit.transp.H R=$@ \
	pverb=2  nx_shots=1 200 ox_shots=5000 \
	par=bpait.gosis.P taper=y taper_pcent=2 marine=y minx=-24750. maxx=46000. & #restart=y

bpait.Drho.H:
	Window3d < bpait.vel.orig.H min1=48000 max1=8990 min2=2987.5 n2=1123 > j1.H
	Window3d < bpait.vel.edit.H min1=48000 max1=8990 min2=2987.5 n2=1123 > j2.H
	Math file1=j1.H file2=j2.H exp="file2/file1-1" > $@

bpait.Dslow.H:
	Window3d < bpait.vel.orig.H min1= 5000 max1=8990 min2=2987.5 n2=1123 > j1.H
	Window3d < bpait.vel.edit.H min1= 5000 max1=8990 min2=2987.5 n2=1123 > j2.H
	Math file1=j1.H file2=j2.H exp="file2-file1" > $@

bpait.inv.slow%.H: #
	Transp < bpait.areal.rec.mud$*.H plane=13 > r1.H
	Transp < bpait.areal.sou.mud$*.H plane=13 > s1.H
	< r1.H ../Tese/bin/OW_lbfgsc_areal.x sou=s1.H \
	vel=bpait.vel.mud.H grad=bpait.Grad.slow$*.H \
	Dimage=bpait.DSO.slow$*.H bkgimg=bpait.Bkgimg.slow$*.H Velout=bpait.Velout.slow$*.H\
	zmin=0 zmax=3000 nh=12 oh=-150 jh=1 dh=25 wbott=0 niter=1 dso=1 inversion=1 \
	traces_in_taper=51 padtraces=51 node=16 bounds=1 vmin=2.5 vmax=3.5 vsal=5 nrefvel=3 > $@