########################################

9000data_hor.H:
	susynlv ref="1:0,.5.;5.,.5." nt=2001 dt=0.002 dxs=0.02 nxs=401 fpeak=20 fxs=0 dxo=0.02 nxo=401 fxo=-4 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > x1.H; echo n2=401 o2=-4 d2=0.02 n3=401 o3=0 d3=0.02 >> x1.H
	susynlv ref="1:0,1.;5.,1." nt=2001 dt=0.002 dxs=0.02 nxs=401 fpeak=20 fxs=0 dxo=0.02 nxo=401 fxo=-4 v00=2.05 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > x2.H; echo n2=401 o2=-4 d2=0.02 n3=401 o3=0 d3=0.02 >> x2.H
	Math file1=x1.H file2=x2.H exp=file1+file2 >> $@
	rm -f junk.*
#
# Window3d <9200data_hor.H n3=101 f3=100 > 92001data_hor.H; echo o3=0 >> 92001data_hor.H
#
9000data_prep.H:
	< 9000data_hor.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=50 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

9000data_source.H:
	< 9000data_hor.H Window3d n3=1 n2=1 min2=0 max1=.6 | Pad n1out=501 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=50 maxsize=1000 > y.H
	< y.H Window3d | Cabs > spec.H
	Math file1=spec.H exp="file1*0" > phse.H; Cmplx spec.H phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

9000_mig_brad.H: # With 401 shots separated 0.01km correct velocity
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=9000data_prep.H D_d=9000data_source.H V=9000vel2_brad.H R=$@ \
	pverb=2  nx_shots=401 ox_shots=0 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-8 maxx=17 & #restart=y

9000_mig_brad_vcorr.H: # With 401 shots separated 0.01km correct velocity
	@date
	Transp plane=12 reshape=1,3 < 9000vel_correct.H > 9000vel_correct_brad.H
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=9000data_prep.H D_d=9000data_source.H V=9000vel_correct_brad.H R=$@ \
	pverb=2  nx_shots=301 ox_shots=1 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-3 maxx=8 & #restart=y

9000_mig_vlomy:
	< 9000data_prep.H ../../ShtPSPI/bin/OWmig2d_par.x sou=9000data_source.H vel=9000vel2.H image=9000_mig_vlomy.H \
	vref=9000vel2_lloyd.H nrefvel=4 nrefvel=1 xmin=-8 xmax=16. ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	noencode=0  node= 12  > /dev/null

9000_mig_vhimy:
	< 9000data_prep.H ../../ShtPSPI/bin/OWmig2d_par.x sou=9000data_source.H vel=9000vel2.2.H image=9000_mig_vhimy.H \
	vref=9000vel2.2_lloyd.H nrefvel=4 nrefvel=1 xmin=-8 xmax=16. ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	noencode=0  node= 12  > /dev/null

#Math file1=9000vel_correct.H exp="2*file1-2" > 9000vel2.2.H
#Transp plane=12 reshape=1,3 < 9000vel2.2.H > 9000vel2.2_transp.H
#Math file1=deltaS.H exp="file1/(-0.0238)*0.021645" > deltaShi.H


9000_mig_vhi.H: # With 401 shots separated 0.01km correct velocity
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=9000data_prep.H D_d=9000data_source.H V=9000vel2.2_transp.H R=$@ \
	pverb=2  nx_shots=401 ox_shots=0 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-8 maxx=17 & #restart=y

9001encz%: # todas as frequencis com a mesma fase que varia em x,z
	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
	< 9000_mig_brad.H Window3d min1=1 max1=4 min3=-.2 n3=21 max5=1.3| Transp plane=13 | ../bin/OW_hgain.x w=4. > 9000_mig_brad_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_brad_transp.H sou_cig=9000_mig_brad_transp.H sou=1source.H \
	vel=9000vel2.H rec_out=9001areal_rec_encz_$*.H sou_out=9001areal_sou_encz_$*.H ngap=40 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node=16 encode=1 encz=1 random=1 nreal=2 eps=1 injection=1 > /dev/null 

9001mig_encz_%.H: # correct velocity areal shot random encode in w
	< 9001areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9001areal_rec_tempz$*.H 
	< 9001areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9001areal_sou_tempz$*.H
	< 9001areal_rec_tempz$*.H ../bin/OWmig2d_areal.x sou=9001areal_sou_tempz$*.H \
	vel=9000vel2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=31 oh=-.3 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 deltaV= ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 9001*tempz$*.H

9002encz%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_brad.H Window3d min1=1 max1=4 min3= -.1 n3= 11 | Transp plane=13 > 9000_mig_brad_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_brad_transp.H sou_cig=9000_mig_brad_transp.H sou=9000data_source.H \
	vel=9000vel2.H rec_out=9002areal_rec_encz_$*.H sou_out=9002areal_sou_encz_$*.H ngap=10 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node=10 encode=1 encz=1 random=1 eps=2 0.75 injection=1 nreal=4 > /dev/null 

9002mig_encz_%.H: # correct velocity areal shot random encode in w
	< 9002areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9002areal_rec_tempz$*.H 
	< 9002areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9002areal_sou_tempz$*.H
	< 9002areal_rec_tempz$*.H ./bin/OWmig2d_areal.x sou=9002areal_sou_tempz$*.H vref=9000vel2_lloyd.H \
	vel=9000vel2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 9002*tempz$*.H

9003encz%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_vhi.H Window3d min1=1 max1=4 min3= -.1 n3= 11 | Transp plane=13 > 9000_mig_vhi_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_vhi_transp.H sou_cig=9000_mig_vhi_transp.H sou=9000data_source.H \
	vel=9000vel2.2.H rec_out=9003areal_rec_encz_$*.H sou_out=9003areal_sou_encz_$*.H ngap=10 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node=10 encode=1 encz=1 random=1 eps=2 0.75 injection=1 nreal=4 > /dev/null 

9003mig_encz_%.H: # correct velocity areal shot random encode in w
	< 9003areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9003areal_rec_tempz$*.H 
	< 9003areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9003areal_sou_tempz$*.H
	< 9003areal_rec_tempz$*.H ./bin/OWmig2d_areal.x sou=9003areal_sou_tempz$*.H vref=9000vel2.2_lloyd.H \
	vel=9000vel2.2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 


#
#
9101encz%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vlo_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_vlo_transp.H sou_cig=9000_mig_vlo_transp.H sou=9000data_source.H \
	vel=9000vel2.H rec_out=9101areal_rec_encz_$*.H sou_out=9101areal_sou_encz_$*.H ngap=10 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node= 12 encode=1 encz=1 random=1 eps=2 0.75 injection=0 nreal=4 > /dev/null 

9101mig_encz_%.H: # correct velocity areal shot random encode in w
	< 9101areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9101areal_rec_tempz$*.H 
	< 9101areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9101areal_sou_tempz$*.H
	< 9101areal_rec_tempz$*.H ./bin/OWmig2d_areal.x sou=9101areal_sou_tempz$*.H vref=9000vel2_lloyd.H \
	vel=9000vel2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 

9501combx%: # todas as frequencis com a mesma fase que varia em x,z
#	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 |Mute tmute=0.6 vmute=0.0 > 9000_mig_vlo_transp.H
	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vlo_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_vlo_transp.H sou_cig=9000_mig_vlo_transp.H sou=../4004data_recip_source_pad.H \
	vel=9000vel2.H rec_out=9101areal_rec_encz_$*.H sou_out=9101areal_sou_encz_$*.H ngap=21 nrefvel=1 vsal=6 \
	padtraces=41 node=16 encode=1 gold=1 perfect=1 jshift=1 fac=$* encz=1 eps=1 injection=1 maxang=80 > /dev/null 

9501mig_combx_%.H: # correct velocity areal shot random encode in w
	< 9101areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9101areal_rec_tempz$*.H 
	< 9101areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9101areal_sou_tempz$*.H
	Math file1=9000vel2.H exp="file1*0+1.8" > 9000vel1.8.H
	< 9101areal_rec_tempz$*.H ./../bin/OWmig2d_areal.x sou=9101areal_sou_tempz$*.H vref=9000vel2_lloyd.H 9000vel1.8_lloyd.H \
	vel=9000vel2.H 9000vel1.8.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 maxang=80 > $@ 

9601combx%: # todas as frequencis com a mesma fase que varia em x,z
#	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 |Mute tmute=0.6 vmute=0.0 > 9000_mig_vlo_transp.H
	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vlo_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_vlo_transp.H sou_cig=9000_mig_vlo_transp.H sou=../4004data_recip_source_pad.H \
	vel=9000vel2.H rec_out=9601areal_rec_encz_$*.H sou_out=9601areal_sou_encz_$*.H ngap=21 nrefvel=1 vsal=6 \
	padtraces=41 node=16 encode=1 sgnpi=1 encz=1 eps=1 maxang=80 > /dev/null 

9601mig_combx_%.H: # correct velocity areal shot random encode in w
	< 9601areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9601areal_rec_tempz$*.H 
	< 9601areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9601areal_sou_tempz$*.H
	Math file1=9000vel2.H exp="file1*0+1.8" > 9000vel1.8.H
	< 9601areal_rec_tempz$*.H ./../bin/OWmig2d_areal.x sou=9601areal_sou_tempz$*.H vref=9000vel2_lloyd.H 9000vel1.8_lloyd.H \
	vel=9000vel2.H 9000vel1.8.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 node=16 maxang=80 > $@ 

9102encz%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_vhimy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vhi_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_vhi_transp.H sou_cig=9000_mig_vhi_transp.H sou=../4004data_recip_source_pad.H \
	vel=9000vel2.2.H rec_out=9102areal_rec_encz_$*.H sou_out=9102areal_sou_encz_$*.H ngap=20 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node= 12 encode=1 encz=1 gold=1 perfect=1 eps=1 injection=1 nreal=4 > /dev/null 

9102mig_encz_%.H: # correct velocity areal shot random encode in w
	< 9102areal_rec_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9102areal_rec_tempz$*.H 
	< 9102areal_sou_encz_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9102areal_sou_tempz$*.H
	< 9102areal_rec_tempz$*.H ./bin/OWmig2d_areal.x sou=9102areal_sou_tempz$*.H vref=9000vel2.2_lloyd.H \
	vel=9000vel2.2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 

9201enc0%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_vlomy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vlo_transp.H
	../bin/OW_arealsht.x rec_cig=9000_mig_vlo_transp.H sou_cig=9000_mig_vlo_transp.H sou=9000data_source.H \
	vel=9000vel2.H rec_out=9201areal_rec_enc0_$*.H sou_out=9201areal_sou_enc0_$*.H ngap=25 nrefvel=1 vsal=6 \
	padtraces=41 node= 12 encode=0 deltat=0.02 dt=0.008 > /dev/null 

9201mig_enc0_%.H: # correct velocity areal shot random encode in w
	< 9201areal_rec_enc0_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9201areal_rec_temp0$*.H 
	< 9201areal_sou_enc0_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9201areal_sou_temp0$*.H
	< 9201areal_rec_temp0$*.H ../bin/OWmig2d_areal.x sou=9201areal_sou_temp0$*.H vref=9000vel2_lloyd.H \
	vel=9000vel2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=15 oh=-.14 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=0 deltat=0.048 dt=0.012 migration=1 > $@ 

9202enc0%: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_vhimy.H Window3d min1=1 max1=4 | Transp plane=13 > 9000_mig_vhi_transp.H
	../bin/OW_arealsht.x rec_cig=9000_mig_vhi_transp.H sou_cig=9000_mig_vhi_transp.H sou=9000data_source.H \
	vel=9000vel2.2.H rec_out=9202areal_rec_enc0_$*.H sou_out=9202areal_sou_enc0_$*.H ngap=25 nrefvel=1 vsal=6 \
	padtraces=41 node= 12 encode=0 deltat=0.02 dt=0.008 > /dev/null 

9202mig_enc0_%.H: # correct velocity areal shot random encode in w
	< 9202areal_rec_enc0_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9202areal_rec_temp0$*.H 
	< 9202areal_sou_enc0_$*.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 9202areal_sou_temp0$*.H
	< 9202areal_rec_temp0$*.H ../bin/OWmig2d_areal.x sou=9202areal_sou_temp0$*.H vref=9000vel2.2_lloyd.H \
	vel=9000vel2.2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=15 oh=-.14 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=0 deltat=0.048 dt=0.012 migration=1 > $@ 

44413mig_encw.H: # migration
	Math file1=../44413areal_rec_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 | Window3d squeeze=n n4=1 > 44413areal_rec_temp.H 
	Math file1=../44413areal_sou_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 | Window3d squeeze=n n4=1 > 44413areal_sou_temp.H 
	< 44413areal_rec_temp.H ./bin/OWmig2d_areal.x sou=44413areal_sou_temp.H \
	vel=../../Tese/4000vel.H nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 \
	wbott=0 migration=144413Inv vsal=3 node= 8 traces_in_taper=25 > $@ 


xxx.H: # correct velocity areal shot random encode in w
	< ../44413areal_rec_encz_velow1.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 2rec.H 
	< ../44413areal_sou_encz_velow1.H Transp plane=13 > junk.H;Math file1=junk.H exp="file1/100000000" > 1sou.H
	< 2rec.H ./bin/OWmig2d_areal.x sou=1sou.H vref=9000vel2.2_lloyd.H \
	vel=9000vel2.2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=21 oh=-.2 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 

##########################################
#
# Testing inversion
#
9101Inv_encz_2000.H: # correct velocity areal shot random encode in w
#	Math file1=9101areal_rec_encz_2000.H junk.H exp="file1/100000000" | Window3d squeeze=n j4=4 | Transp plane=13 > 9101areal_rec_tempz.H 
#	Math file1=9101areal_sou_encz_2000.H junk.H exp="file1/100000000" | Window3d squeeze=n j4=4 | Transp plane=13 > 9101areal_sou_tempz.H 
	Math file1=9101areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_rec_tempz.H 
	Math file1=9101areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_sou_tempz.H 
	< 9101areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9101areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H grad=9101Grad.H dso=9101DSO.H bkgimg=9101Bkgimg.H\
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node= 10 > $@

9102Inv_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9102areal_rec_encz_2000.H junk.H exp="file1/100000000" | Window3d squeeze=n n4=10 | Transp plane=13 > 9102areal_rec_tempz.H 
	Math file1=9102areal_sou_encz_2000.H junk.H exp="file1/100000000" | Window3d squeeze=n n4=10 | Transp plane=13 > 9102areal_sou_tempz.H 
	< 9102areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9102areal_sou_tempz.H \
	vel=9000vel2.2.H vref= 9000vel2.2_lloyd.H grad=9102Grad.H dso=9102DSO.H bkgimg=9102Bkgimg.H\
	nrefvel=1 vsal=3 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node= 12 > $@

#9101Inv_encz_2000.H: # correct velocity areal shot random encode in w
#	Math file1=9101areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_rec_tempz.H 
#	Math file1=9101areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_sou_tempz.H 
#	< 9101areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9101areal_sou_tempz.H \
#	vel=9000vel2.H vref=9000vel2_lloyd.H grad=9101Grad.H dso=9101DSO.H bkgimg=9101Bkgimg.H\
#	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
#	traces_in_taper=25 node= 10 > $@
#
#9102Inv_encz_2000.H: # correct velocity areal shot random encode in w
#	Math file1=9102areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_rec_tempz.H 
#	Math file1=9102areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_sou_tempz.H 
#	< 9102areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9102areal_sou_tempz.H \
#	vel=9000vel2.2.H vref= 9000vel2.2_lloyd.H grad=9102Grad.H dso=9102DSO.H bkgimg=9102Bkgimg.H\
#	nrefvel=1 vsal=3 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
#	traces_in_taper=25 node= 12 > $@

9103Inv_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9101areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_rec_tempz.H 
	Math file1=9101areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_sou_tempz.H 
	< 9101areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9101areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H grad=9103Grad.H dso=9103DSO.H bkgimg=9103Bkgimg.H\
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node= 10 impulse=1 > $@

9104Inv_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9102areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_rec_tempz.H 
	Math file1=9102areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_sou_tempz.H 
	< 9102areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9102areal_sou_tempz.H \
	vel=9000vel2.2.H vref= 9000vel2.2_lloyd.H grad=9104Grad.H dso=9104DSO.H bkgimg=9104Bkgimg.H\
	nrefvel=1 vsal=3 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node= 12  impulse=1 > $@

9201Inv_enc0_1000.H: # correct velocity areal shot random encode in w
	Math file1=9201areal_rec_enc0_1000.H junk.H exp="file1/100000000" | Transp plane=13 > 9201areal_rec_temp0.H 
	Math file1=9201areal_sou_enc0_1000.H junk.H exp="file1/100000000" | Transp plane=13 > 9201areal_sou_temp0.H 
	< 9201areal_rec_temp0.H ./bin/OW_lbfgsc_areal.x sou=9201areal_sou_temp0.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H grad=9201Grad.H dso=9201DSO.H bkgimg=9201Bkgimg.H\
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=15 oh=-0.14 -.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 deltat=0.048 dt=0.012 icwind=1 node= 12 > $@

9202Inv_enc0_1000.H: # correct velocity areal shot random encode in w
	Math file1=9202areal_rec_enc0_1000.H junk.H exp="file1/100000000" | Transp plane=13 > 9202areal_rec_temp0.H 
	Math file1=9202areal_sou_enc0_1000.H junk.H exp="file1/100000000" | Transp plane=13 > 9202areal_sou_temp0.H 
	< 9202areal_rec_temp0.H ./bin/OW_lbfgsc_areal.x sou=9202areal_sou_temp0.H \
	vel=9000vel2.2.H vref= 9000vel2.2_lloyd.H grad=9202Grad.H dso=9202DSO.H bkgimg=9202Bkgimg.H\
	nrefvel=1 vsal=3 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=15 oh=-0.14 -.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 deltat=0.048 dt=0.012 icwind=1 node= 12 > $@
##########################################
#
# Testing forward and adjoint
#
9111for_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9101areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_rec_tempz.H 
	Math file1=9101areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9101areal_sou_tempz.H 
	< 9101areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9101areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaS= \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 adj=0  \
	traces_in_taper=25 node= 10 > $@

#9102Inv_encz_2000.H: # correct velocity areal shot random encode in w
#	Math file1=9102areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_rec_tempz.H 
#	Math file1=9102areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9102areal_sou_tempz.H 
#	< 9102areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9102areal_sou_tempz.H \
#	vel=9000vel2.2.H vref= 9000vel2.2_lloyd.H grad=9102Grad.H dso=9102DSO.H bkgimg=9102Bkgimg.H\
#	nrefvel=1 vsal=3 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
#	traces_in_taper=25 node= 12 > $@


9000vel2.2_lloyd.H: # it is slowness !!!!!!
	Lloyd_vel < 9000vel2.2_transp.H nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

9000vel2_lloyd.H: # it is slowness !!!!!!
	Lloyd_vel < 9000vel2_brad.H nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

9000vel1.8_lloyd.H: # it is slowness !!!!!!
	Math file1=9000vel2_brad.H exp="file1*0.+1.8" > j.H
	Lloyd_vel < j.H nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

9005Dimg_encz.5r.H: # correct velocity areal shot random encode in w
	Cat 9001areal_sou_encz_0001.H 9001areal_sou_encz_0011.H 9001areal_sou_encz_0021.H 9001areal_sou_encz_0031.H 9001areal_sou_encz_0041.H axis=4 > junk.H;Math file1=junk.H exp="file1/100000000" | \
	Transp plane=13 > 9001areal_sou_tempz.H #;echo sets next: in="/net/mad/scr1/claudio/9001areal_sou_tempz.H@" >> 9001areal_sou_tempz.H
	Cat 9001areal_rec_encz_0001.H 9001areal_rec_encz_0011.H 9001areal_rec_encz_0021.H 9001areal_rec_encz_0031.H 9001areal_rec_encz_0041.H axis=4 > junk.H;Math file1=junk.H exp="file1/100000000" | \
	Transp plane=13 > 9001areal_rec_tempz.H #;echo sets next: in="/net/mad/scr1/claudio/9001areal_rec_tempz.H@" >> 9001areal_rec_tempz.H
	< 9001areal_rec_tempz.H ./bin/OWmig2d_areal.x sou=9001areal_sou_tempz.H \
	vel=9000vel2.H 9000vel_correct.H vref=9000vel2_lloyd.H deltaS=deltaS.H deltaS1.H ws=wvfld_S.5r.H wr=wvfld_R.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 
	Rm 9001*tempz$*.H

9005Dslow_encz.5r.H: # correct velocity areal shot random encode in w
	< wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=wvfld_S.5r.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9005Dimg_encz.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 
	Rm 9001*tempz$*.H

9006Dimg_encz.4r.H: # correct velocity areal shot random encode in w
	Math file1=9002areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9002areal_rec_tempz.H 
	Math file1=9002areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9002areal_sou_tempz.H 
	< 9002areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9002areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaS=deltaS.H ws= wvfld_S.5r.H wr= wvfld_R.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 
	Rm 9001*tempz$*.H

9006Dslow_encz.4r.H: # correct velocity areal shot random encode in w
	< 9002areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9002areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9006Dimg_encz.4r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node=8 16 > $@ 

9006Dslow_encz.DSO.1r.H: # correct velocity areal shot random encode in w
	< 9002areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9002areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9006Bkgimg.H 9006DSO.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 

9007Dimg_encz.4r.H: # correct velocity areal shot random encode in w
	Math file1=9003areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9003areal_rec_tempz.H 
	Math file1=9003areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9003areal_sou_tempz.H 
	< 9003areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9003areal_sou_tempz.H \
	vel=9000vel2.2.H vref=9000vel2.2_lloyd.H deltaS=deltaShi.H  \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 

9007Dslow_encz.4r.H: # correct velocity areal shot random encode in w
	< 9003areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9003areal_sou_tempz.H \
	vel=9000vel2.2.H vref=9000vel2.2_lloyd.H deltaI=9007Dimg_encz.4r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 

9007Dslow_encz.DSO.1r.H: # correct velocity areal shot random encode in w
	< 9003areal_rec_tempz.H ./bin/OWwemva_opers.x sou=9003areal_sou_tempz.H \
	vel=9000vel2.2.H vref=9000vel2.2_lloyd.H deltaI=9006Bkgimg.H 9007DSO.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 

9006Dimg_encz.5r.H: # correct velocity areal shot random encode in w
	Cat 9001areal_sou_encz_0001.H 9001areal_sou_encz_0011.H 9001areal_sou_encz_0021.H 9001areal_sou_encz_0031.H 9001areal_sou_encz_0041.H axis=4 > junk.H;Math file1=junk.H exp="file1/100000000" | \
	Transp plane=13 > 9001areal_sou_tempz.H #;echo sets next: in="/net/mad/scr1/claudio/9001areal_sou_tempz.H@" >> 9001areal_sou_tempz.H
	Cat 9001areal_rec_encz_0001.H 9001areal_rec_encz_0011.H 9001areal_rec_encz_0021.H 9001areal_rec_encz_0031.H 9001areal_rec_encz_0041.H axis=4 > junk.H;Math file1=junk.H exp="file1/100000000" | \
	Transp plane=13 > 9001areal_rec_tempz.H #;echo sets next: in="/net/mad/scr1/claudio/9001areal_rec_tempz.H@" >> 9001areal_rec_tempz.H
	< 9001areal_rec_tempz.H ./bin/OWmig2d_areal.x sou=9001areal_sou_tempz.H \
	vel=9000vel2.H 9000vel_correct.H vref=9000vel2_lloyd.H deltaS=deltaS.H deltaS1.H ws=wvfld_S.5r.H wr=wvfld_R.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=31 oh=-.3 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 
	Rm 9001*tempz$*.H

9006Dslow_encz.5r.H: # correct velocity areal shot random encode in w
	< wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=wvfld_S.5r.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9006Dimg_encz.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 
	Rm 9001*tempz$*.H

9006Dslow_encz.5r.ltd.H: # correct velocity areal shot random encode in w
	< 9006Dimg_encz.5r.H Window3d f2=12 n2=7 > 9006Dimg_encz.5r.ltd.H
	< wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=wvfld_S.5r.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9006Dimg_encz.5r.ltd.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 
	Rm 9001*tempz$*.H

9006Dslow_encz.5r.ltd1.H: # correct velocity areal shot random encode in w
	< 9006Dimg_encz.5r.H Window3d f2=12 n2=7 > 9006Dimg_encz.5r.ltd.H
	< wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=wvfld_S.5r.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9006Dimg_encz.5r.ltd.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.05 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 
	Rm 9001*tempz$*.H

##############################################
#
# Using Gold codes
#
9002agold%: # correct velocity areal shot random encode in w
	< 9000_mig_brad.H Window3d min1=1 max1=4 min3= -.1 n3= 11 | Transp plane=13 |Mute tmute=0.6 vmute=0 > 9000_mig_brad_transp.H
	/net/koko/wrk1/sep136/claudio1/Bin/OW_arealsht.x rec_cig=9000_mig_brad_transp.H sou_cig=9000_mig_brad_transp.H \
	sou=../4004data_recip_source_pad.H 9000data_source.H vel=9000vel2.H rec_out=9002areal_rec_goldt.H sou_out=9002areal_sou_goldt.H ngap=17 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node= 10 jshift=2 fac=$* \
	> /dev/null

9002bgold%: # correct velocity areal shot random encode in w
	< 9000_mig_brad.H Window3d min1=1 max1=4 min3= -.1 n3= 11 | Transp plane=13 > 9000_mig_brad_transpb.H
	Math file1=9000_mig_brad_transp.H file2=9000_mig_brad_transpb.H exp="file2-file1" > 9000_mig_brad_transpm.H
	/net/koko/wrk1/sep136/claudio1/Bin/OW_arealsht.x rec_cig=9000_mig_brad_transpm.H sou_cig=9000_mig_brad_transpm.H \
	sou=../4004data_recip_source_pad.H 9000data_source.H vel=9000vel2.H rec_out=9002areal_rec_goldb.H sou_out=9002areal_sou_goldb.H ngap=17 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node= 10 jshift=2 fac=$* \
	> /dev/null
	Cat 9002areal_sou_goldt.H 9002areal_sou_goldb.H axis=4> 9002areal_sou_gold.H
	Cat 9002areal_rec_goldt.H 9002areal_rec_goldb.H axis=4 > 9002areal_rec_gold.H

9002mig_gold.H:
	< 9002areal_rec_gold.H Transp plane=13 > 9001areal_rec_tempg.H
	< 9002areal_sou_gold.H Transp plane=13 > 9001areal_sou_tempg.H
	< 9001areal_rec_tempg.H ./bin/OWmig2d_areal.x sou=9001areal_sou_tempg.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 migration=1 > $@ 
	Rm 9001*tempg$*.H

9002Dimg_gold.H: # computing the perturbed image (forward modeling)
	< 9002areal_rec_gold.H Transp plane=13 > 9001areal_rec_tempg.H
	< 9002areal_sou_gold.H Transp plane=13 > 9001areal_sou_tempg.H
	< 9001areal_rec_tempg.H ./bin/OWwemva_opers.x sou=9001areal_sou_tempg.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaS=deltaS.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 adj=0 vsal=3 node= 8 \
	traces_in_taper=25 > $@ 
	
9002Dslow_gold.H: # computing the slowness perturbation (adjoint)
	< 9001areal_rec_tempg.H ./bin/OWwemva_opers.x sou=9001areal_sou_tempg.H \
	vel=9000vel2.H vref= deltaI=9002Dimg_gold.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 wbott=0 adj=1 vsal=3 node= 8\
	traces_in_taper=25 > $@ 

##############################################
#
# Shot profile
#
MIG_PARA = offset=y ntpx=25 nws=193 oper="phsft" nxpad_beg=50 nxpad_end=50 image_xmin=1 image_xmax=4 # "foufd" 
OFF_PARA = nhx=11 ohx=-.1
BSC_PARA = typbsc=2 nopbsc=3
BINDIR = /net/koko/tang/bin/LINUX

9002Dimg_shotprof.H:
	Window3d max4=5.4 squeeze=n < 9000data_prep.H > 1.H
	Math file1=9000vel2_brad.H exp="1/file1" | Pad beg1=100 extend=1 > 9000slow2_brad.H
	Transp plane=12 reshape=1,3 < deltaS.H | Pad beg1=100 extend=1 > deltaS_transp.H
	${BINDIR}/wemva3d_baspix.x csou=9000data_source.H crec=1.H bslw=9000slow2_brad.H \
                                   dimg=$@ dslw=deltaS_transp.H adj=n add=n image_xmin=1 image_xmax=4 \
                                   comsou=n ${MIG_PARA} ${BSC_PARA} ${OFF_PARA} radian=n

#	/net/koko/claudio/Research/Others/yaxun/wemva3d_baspix.x csou=9000data_source.H crec=1.H \
#	bslw=9000slo2_brad.H nhx=11 ohx=-.1 \
#	dslw=deltaS_brad.H dimg=$@ oper=phsft nws=193 adj=n ntpx=15 radian=n offset=y comsou=n \
#	image_xmin=1 image_xmax=4

9002Dslow_shotprof.H:
	${BINDIR}/wemva3d_baspix.x csou=9000data_source.H crec=1.H bslw=9000slow2_brad.H \
                                   dslw=$@ dimg=9002Dimg_shotprof.H adj=y add=n\
                                   comsou=n ${MIG_PARA} ${BSC_PARA} ${OFF_PARA} radian=n

#	< 9002Dimg_shotprof1.H Window3d min1=1 max1=4 squeeze=n | Pad beg1=61 n1out=273 > 2.H
#	/net/koko/claudio/Research/Others/yaxun/wemva3d_baspix.x csou=9000data_source.H crec=1.H \
#	bslw=9000slo2_brad.H dimg=9002Dimg_shotprof.H dslw=$@ oper=phsft nws=193 adj=y ntpx=15 radian=n offset=y comsou=n

##############################################
#
#44413areal_sou_encz_velow1000.H
#
4000velow_lloyd.H: # it is slowness !!!!!!
	< ../../Tese/4000velow.H Transp plane=12 reshape=1,3 > 4000velow_transp.H
	Lloyd_vel < 4000velow_transp.H nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

4000deltaS.H:
	Math file1=../../Tese/4000velow.H exp="file1*0-0.2" > $@

444131velow_lloyd.H: # it is slowness !!!!!!
	< ../../Tese/44413vlow1.9.H Transp plane=12 reshape=1,3 > 44413vlow1.9_transp.H
	Lloyd_vel < 44413vlow1.9_transp.H nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

444131deltaS.H:
	Math file1=../../Tese/4000velow.H exp="file1*0-0.1" > $@

44413Dimg_encz.5r.H: # correct velocity areal shot random encode in w
#	Cat ../../Tese/44413areal_sou_encz_velow1000.H ../../Tese/44413areal_sou_encz_velow2000.H \
#	../../Tese/44413areal_sou_encz_velow3000.H ../../Tese/44413areal_sou_encz_velow4000.H \
#	../../Tese/44413areal_sou_encz_velow5000.H axis=4 > junk.H 
#	Math file1=junk.H exp="file1/100000000/100000000" | Transp plane=13 > 44413areal_sou_5r.H 
#	Cat ../../Tese/44413areal_rec_encz_velow1000.H ../../Tese/44413areal_rec_encz_velow2000.H \
#	../../Tese/44413areal_rec_encz_velow3000.H ../../Tese/44413areal_rec_encz_velow4000.H \
#	../../Tese/44413areal_rec_encz_velow5000.H axis=4 > junk.H 
#	Math file1=junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_5r.H 
	Math file1=../44413areal_rec_encz_velow105.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_5r.H 
	Math file1=../44413areal_sou_encz_velow105.H exp="file1/100000000" | Transp plane=13 > 44413areal_sou_5r.H 
	< 44413areal_rec_5r.H ./bin/OWmig2d_areal.x sou=44413areal_sou_5r.H \
	vel=../../Tese/4000velow.H vref=4000velow_lloyd.H deltaS=4000deltaS.H ws=4000wvfld_S.5r.H wr=4000wvfld_R.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 

44413Dslow_encz.5r.H: # correct velocity areal shot random encode in w
	< 4000wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=4000wvfld_S.5r.H \
	vel=../../Tese/4000velow.H vref=4000velow_lloyd.H deltaI=44413Dimg_encz.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 

444131Dimg_encz.5r.H: # correct velocity areal shot random encode in w
#	Cat ../../Tese/44413areal_sou_encz_velow1000.H ../../Tese/44413areal_sou_encz_velow2000.H \
#	../../Tese/44413areal_sou_encz_velow3000.H ../../Tese/44413areal_sou_encz_velow4000.H \
#	../../Tese/44413areal_sou_encz_velow5000.H axis=4 > junk.H 
#	Math file1=junk.H exp="file1/100000000/100000000" | Transp plane=13 > 44413areal_sou_5r.H 
#	Cat ../../Tese/44413areal_rec_encz_velow1000.H ../../Tese/44413areal_rec_encz_velow2000.H \
#	../../Tese/44413areal_rec_encz_velow3000.H ../../Tese/44413areal_rec_encz_velow4000.H \
#	../../Tese/44413areal_rec_encz_velow5000.H axis=4 > junk.H 
#	Math file1=junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_5r.H 
	Math file1=../444131areal_rec_encz_velow105.H exp="file1/100000000" | Transp plane=13 > 444131areal_rec_5r.H 
	Math file1=../444131areal_sou_encz_velow105.H exp="file1/100000000" | Transp plane=13 > 444131areal_sou_5r.H 
	< 444131areal_rec_5r.H ./bin/OWmig2d_areal.x sou=444131areal_sou_5r.H \
	vel=../../Tese/44413vlow1.9.H vref=444131velow_lloyd.H deltaS=444131deltaS.H ws=444131wvfld_S.5r.H wr=444131wvfld_R.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 every=0 encode=1 node= 8 > $@ 

444131Dslow_encz.5r.H: # correct velocity areal shot random encode in w
	< 444131wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=444131wvfld_S.5r.H \
	vel=../../Tese/44413vlow1.9.H vref=444131velow_lloyd.H deltaI=444131Dimg_encz.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 every=0 encode=1 node= 16 > $@ 

#
#
##############################################
#
dottst%:
	< 9001areal_rec_encz_$*.H Window3d n4=1 squeeze=n | Transp plane=13 | Scale scale_to=1 > 9001areal_rec_tempz$*.H
	< 9001areal_sou_encz_$*.H Window3d n4=1 squeeze=n | Transp plane=13 | Scale scale_to=1 > 9001areal_sou_tempz$*.H
	< 9001areal_rec_tempz$*.H ./bin/OWwemva_opers.x sou=9001areal_sou_tempz$*.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 dot=1 \
	traces_in_taper=25 node= 1 > /dev/null
	Rm 9001*tempz$*.H

dotimg:
	< 9001areal_rec_encz_1000.H Window3d n4=1 squeeze=n | Transp plane=13 | Scale scale_to=1 > 9001areal_rec_tempz1000.H
	< 9001areal_sou_encz_1000.H Window3d n4=1 squeeze=n | Transp plane=13 | Scale scale_to=1 > 9001areal_sou_tempz1000.H
	< 9001areal_rec_tempz1000.H ./bin/OWwemva_opers.x sou=9001areal_sou_tempz1000.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 dot=1 dotimg=1 \
	traces_in_taper=25 node= 1 > /dev/null
	Rm 9001*tempz*.H

9005Dslow_encz.5ra.H: # correct velocity areal shot random encode in w
	< wvfld_R.5r.H ./bin/OWmig2d_areal.x sou=wvfld_S.5r.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9005Dimg_encz.5r.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=0 \
	traces_in_taper=0 every=0 encode=1 node= 8 > $@ 
	Rm 9001*tempz$*.H

9005Dimg_encz_%.H: # correct velocity areal shot random encode in w
	< 9001areal_rec_encz_$*.H Transp plane=13 | Scale scale_to=1 > 9001areal_rec_tempz$*.H
	< 9001areal_sou_encz_$*.H Transp plane=13 | Scale scale_to=1 > 9001areal_sou_tempz$*.H
	< 9001areal_rec_tempz$*.H ./bin/OWmig2d_areal.x sou=9001areal_sou_tempz$*.H \
	vel=9000vel2.H 9000vel_correct.H vref=9000vel2_lloyd.H deltaS=deltaS.H deltaS1.H ws= wvfld_S.H wr= wvfld_R.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=1 31 oh=0 -.3 dh=0.02 wbott=0 adj=0 \
	traces_in_taper=25 > $@ 
	Rm 9001*tempz$*.H

9005Dslow_encz_%.H: # correct velocity areal shot random encode in w
	< wvfld_R.H ./bin/OWmig2d_areal.x sou=wvfld_S.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=9005Dimg_encz_$*.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=25 node= 8 ws= wvfld_S1.H wr= wvfld_R1.H Dimg= DI.H> $@ 
	Rm 9001*tempz$*.H

9005Dslow_encz_tang_mydimg.H: # correct velocity areal shot random encode in w
	< 9005Dimg_encz_tang.H Window3d n3=375 | Transp plane=12 reshape=1,3 > xxx.H
	< wvfld_R.H ./bin/OWmig2d_areal.x sou=wvfld_S.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H deltaI=xxx.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 wbott=0 adj=1 \
	traces_in_taper=0 25 node= 8 ws=wvfld_S1.H wr=wvfld_R1.H Dimg= DI.H> $@ 
	Rm 9001*tempz$*.H
##########################################
#
# Testing inversion
#
9006Inv_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9002areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9002areal_rec_tempz.H 
	Math file1=9002areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9002areal_sou_tempz.H 
	< 9002areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9002areal_sou_tempz.H \
	vel=9000vel2.H vref=9000vel2_lloyd.H grad=9006Grad.H dso=9006DSO.H bkgimg=9006Bkgimg.H\
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node=12 > $@

9007Inv_encz_2000.H: # correct velocity areal shot random encode in w
	Math file1=9003areal_rec_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9003areal_rec_tempz.H 
	Math file1=9003areal_sou_encz_2000.H junk.H exp="file1/100000000" | Transp plane=13 > 9003areal_sou_tempz.H 
	< 9003areal_rec_tempz.H ./bin/OW_lbfgsc_areal.x sou=9003areal_sou_tempz.H \
	vel=9000vel2.2.H vref=9000vel2.2_lloyd.H grad=9007Grad.H dso=9007DSO.H bkgimg=9007Bkgimg.H\
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node= 12 > $@

44413Dimg_encz.H: # computing the perturbed image (forward modeling)
	Math file1=../44413areal_rec_encw_velow1.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_temp.H 
	Math file1=../44413areal_sou_encw_velow1.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_sou_temp.H 
	Math file1=../../Tese/44413vlow1.9.H exp="1/2-1/file1" > 44413deltaS.H
	< 44413areal_rec_temp.H ./bin/OWwemva_opers.x sou=44413areal_sou_temp.H \
	vel=../../Tese/44413vlow1.9.H vref= deltaS=44413deltaS.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 adj=0 vsal=3 node=8 \
	traces_in_taper=25 > $@ 
	
44413Dslow_encz.H: # computing the slowness perturbation (adjoint)
	< 44413areal_rec_temp.H ./bin/OWwemva_opers.x sou=44413areal_sou_temp.H \
	vel=../../Tese/44413vlow1.9.H vref= deltaI=44413Dimg_encz.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 wbott=0 adj=1 vsal=3 node=8\
	traces_in_taper=25 > $@ 

44413Dslow_encz_a.H: # computing the slowness perturbation (adjoint) with DSO output
	Window3d n4=1 < 44413DSO_a.H > 44413DSO_aa.H
	< 44413areal_rec_temp.H ./bin/OWwemva_opers.x sou=44413areal_sou_temp.H \
	vel=../../Tese/44413vlow1.9.H vref= deltaI=44413DSO_aa.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 wbott=0 adj=1 vsal=3 node=8\
	traces_in_taper=25 > $@ 

44414Dimg_encz.H: # computing the perturbed image (forward modeling)
	Math file1=../44413areal_rec_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_rec_temp.H 
	Math file1=../44413areal_sou_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_sou_temp.H 
	Math file1=../../Tese/junkvelhi.H exp="1/2-1/file1" > 44414deltaS.H
	< 44414areal_rec_temp.H ./bin/OWwemva_opers.x sou=44414areal_sou_temp.H \
	vel=../../Tese/junkvelhi.H vref= deltaS=44414deltaS.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 adj=0 vsal=3 node=8 \
	traces_in_taper=25 > $@ 
	
44414Dslow_encz.H: # computing the slowness perturbation (adjoint)
	< 44414areal_rec_temp.H ./bin/OWwemva_opers.x sou=44414areal_sou_temp.H \
	vel=../../Tese/junkvelhi.H vref= deltaI=44414Dimg_encz.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 wbott=0 adj=1 vsal=3 node=8\
	traces_in_taper=25 > $@ 

44414Dslow_encz_a.H: # computing the slowness perturbation (adjoint) with DSO output
	Window3d n4=1 < 44414DSO_a.H > 44414DSO_aa.H
	< 44414areal_rec_temp.H ./bin/OWwemva_opers.x sou=44414areal_sou_temp.H \
	vel=../../Tese/junkvelhi.H vref= deltaI=44414DSO_aa.H \
	nrefvel=1 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 wbott=0 adj=1 vsal=3 node=8\
	traces_in_taper=25 > $@ 

44413Inv_encw.H: # correct velocity areal shot random encode in w
	Math file1=../44413areal_rec_encw_velow1.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_temp.H 
	Math file1=../44413areal_sou_encw_velow1.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_sou_temp.H 
	< 44413areal_rec_temp.H ./bin/OW_lbfgsc_areal.x sou=44413areal_sou_temp.H \
	vel=../../Tese/44413vlow1.9.H vref= 444131velow_lloyd.H grad=44413Grad.H Dimage=44413DSO.H bkgimg=44413Bkgimg.H\
	nrefvel=2 vsal=4 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1 resmig=1 \
	traces_in_taper=25 node= 12 > $@

44413Inv_gencw.H: # correct velocity areal shot random encode in w
#	Math file1=../44413areal_rec_gencw_velow1a.H junk.H exp="file1/100000000" | Transp plane=13 | Window3d j4=5 squeeze=n > 44413areal_rec_temp.H 
#	Math file1=../44413areal_sou_gencw_velow1a.H junk.H exp="file1/100000000" | Transp plane=13 | Window3d j4=5 squeeze=n > 44413areal_sou_temp.H 
	Math file1=../44413areal_rec_gencw_velow1a.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_rec_temp.H 
	Math file1=../44413areal_sou_gencw_velow1a.H junk.H exp="file1/100000000" | Transp plane=13 > 44413areal_sou_temp.H 
	< 44413areal_rec_temp.H ./bin/OW_lbfgsc_areal.x sou=44413areal_sou_temp.H \
	vel=../../Tese/44413vlow1.9.H vref= grad=44413Grad_a.H Dimage=44413DSO_a.H bkgimg=44413Bkgimg_a.H Velout=44413Velout_a.H \
	nrefvel=2 vsal=4 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1 dso=1 resmig_lin= 1 rho=0.95 \
	traces_in_taper=25 node= 8 > $@

44414Inv_gencw.H: # correct velocity areal shot random encode in w
	Math file1=../44414areal_rec_gencw_vehi1a.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_rec_temp.H 
	Math file1=../44414areal_sou_gencw_vehi1a.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_sou_temp.H 
	< 44414areal_rec_temp.H ./bin/OW_lbfgsc_areal.x sou=44414areal_sou_temp.H \
	vel=../../Tese/junkvelhi.H vref= grad=44414Grad_a.H dso=44414DSO_a.H bkgimg=44414Bkgimg_a.H Velout=44414Velout_a.H \
	nrefvel=2 vsal=4 xmin=-1.22 xmax=6.2 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1 dso=1 rho=1.1 \
	traces_in_taper=25 node= 8 > $@

44414Inv_encw.H: # correct velocity areal shot random encode in w
	Math file1=../44413areal_rec_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_rec_temp.H 
	Math file1=../44413areal_sou_encw_vehi1.H junk.H exp="file1/100000000" | Transp plane=13 > 44414areal_sou_temp.H 
	< 44414areal_rec_temp.H ./bin/OW_lbfgsc_areal.x sou=44414areal_sou_temp.H \
	vel=../../Tese/junkvelhi.H vref= 444131velow_lloyd.H grad=44414Grad.H dso=44414DSO.H bkgimg=44414Bkgimg.H\
	nrefvel=2 vsal=4 xmin=-4000 xmax=3000 zmin=0 zmax=1 nh=11 oh=-.1 dh=0.02 wbott=0 niter=2 inversion=1  \
	traces_in_taper=25 node=8 > $@

concept.inv.slow.H: # correct velocity areal shot random encode in w
#	Window3d < concept.crec.H n4=1 squeeze=n > xxx.H
	< concept.crec.H ./bin/OW_lbfgsc_areal.x sou=concept.csou.H ws= concept.S.H wr= concept.R.H \
	vel=concept.bvel.slow.H vref= 444131velow_lloyd.H grad=concept.Grad.slow.H dso=concept.DSO.slow.H bkgimg=concept.Bkgimg.slow.H\
	nrefvel=1 vsal=2000 zmin=0 zmax=2000 nh=33 oh=-160 dh=10 wbott=0 niter=1 inversion=1  \
	traces_in_taper=25 node=16 > $@

concept.inv.fast.H: # correct velocity areal shot random encode in w
#	Window3d < concept.crec.H n4=1 squeeze=n > xxx.H
	< concept.crec.H ./bin/OW_lbfgsc_areal.x sou=concept.csou.H ws= concept.S.H wr= concept.R.H \
	vel=concept.bvel.fast.H vref= 444131velow_lloyd.H grad=concept.Grad.fast.H dso=concept.DSO.fast.H bkgimg=concept.Bkgimg.fast.H\
	nrefvel=1 vsal=3000 zmin=0 zmax=2000 nh=33 oh=-160 dh=10 wbott=0 niter=1 inversion=1  \
	traces_in_taper=25 node=12 > $@

wem.inv.slow.H: # correct velocity areal shot random encode in w
	< wem.crec.H ./bin/OW_lbfgsc_areal.x sou=wem.csou.H ws= \
	vel=wem.vmod.slow.H grad=wem.Grad.slow.H dso=wem.DSO.slow.H bkgimg=wem.Bkgimg.slow.H\
	nrefvel=1 vsal=2000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 niter=1 inversion=1  \
	traces_in_taper=25 node=8 16 > $@

wem.inv.fast.H: # correct velocity areal shot random encode in w
	< wem.crec.H ./bin/OW_lbfgsc_areal.x sou=wem.csou.H \
	vel=wem.vmod.fast.H grad=wem.Grad.fast.H dso=wem.DSO.fast.H bkgimg=wem.Bkgimg.fast.H\
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 niter=1 inversion=1  \
	traces_in_taper=25 node=8 16 > $@

wem.vel.slow.lloyd.H:
	Transp plane=12 reshape=1,3 < wem.vmod.slow.H > junk.H
	< junk.H Lloyd_vel nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

wem.vel.fast.lloyd.H:
	Transp plane=12 reshape=1,3 < wem.vmod.fast.H > junk.H
	< junk.H Lloyd_vel nref=1 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

wem.dimg.slow.oper.H:
	Pad < wem.csou.H beg1=200 end1=200 > wem.csou.pad.H; echo o1=-2000 d1=10 n4=1 o4=-600 d4=200 >> wem.csou.pad.H
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.slow.H vref=wem.vel.slow.lloyd.H deltaS=wem.dslow.slow.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=0  \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

wem.dimg.fast.oper.H:
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.fast.H vref=wem.vel.fast.lloyd.H deltaS=wem.dslow.fast.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=0  \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

wem.dslo.slow.oper.H:
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.slow.H vref=wem.vel.slow.lloyd.H deltaI=wem.dimg.slow.oper.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=1 \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

wem.dslo.fast.oper.H:
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.fast.H vref=wem.vel.fast.lloyd.H deltaI=wem.dimg.fast.oper.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=1 \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

wem.dimg.slow1.oper.H:
	Pad < wem.csou.H beg1=200 end1=200 > wem.csou.pad.H; echo o1=-2000 d1=10 n4=1 o4=-600 d4=200 >> wem.csou.pad.H
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.slow.H vref=wem.vel.slow.lloyd.H deltaS=wem.dslow.slow.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=0  \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

wem.dslo.slow1.oper.H:
	< wem.crec.H ./bin/OWwemva_opers.x sou=wem.csou.pad.H \
	vel=wem.vmod.slow.H vref=wem.vel.slow.lloyd.H deltaI=wem.dimg.slow1.oper.H \
	nrefvel=1 vsal=3000 zmin=0 zmax=1500 nh=33 oh=-160 dh=10 wbott=0 adj=1 \
	xmin=-3000 xmax=3000 zmin=0 traces_in_taper=25 node=12 > $@

##########################################
#
# Bill Symes' model
#
7000data.H:
	< /net/koko/tang/Research/wemva/tripl.crec.coarse.H Pad beg1=200 end1=200 > $@ #; echo o3=4.83164 d3=0.1666 >> $@

7000velow_pady.H:
	</net/koko/tang/Research/wemva/tripl.vwrg.H Pad extend=1 beg1=200 end1=200 > $@

7000source.H:
	< /net/koko/tang/Research/wemva/tripl.csou.H Window3d squeeze=n > $@ #; echo o3=4.83164 d3=0.1666 >> $@

7000_mig_shpt.H: # With 401 shots separated 0.01km wrong velocity
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=7000data.H D_d=7000source.H V=7000velow_pad.H R=$@ \
	pverb=2  nx_shots=401 ox_shots=-2000 \
	par=go_sis_symes.P taper=y taper_pcent=2 offset=n rollalong=y minx=-8000 maxx=8000 & #restart=y

7000_mig_shpty.H: # With 401 shots separated 0.01km wrong velocity
	/net/koko/tang/bin/LINUX/wem3d.x crec=7000data.H csou=7000source.H vmod=7000velow_pady.H imag=$@ comsou=n wnttap=y nvrf=1 pdip= 45 ntpx=20 nws=181 oper="phspi" \
	radian=y wntimg=y wnthes=n nhx=81 ohx=-400


##########################################
#
# <deltaS.H Transp plane=12 reshape=1,3 > deltaS_brad.H
#
wemva.forward.shtp:
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9000data_source.H crec=9000data_prep.H bvel=9000vel2_brad.H \
	dslw=deltaS_brad.H dimg=9005Dimg_shpt_tang.H oper=phsft nws=193 adj=n ntpx=15 radian=n offset=y comsou=n \
	image_xmin=0.08 image_xmax=4.92

wemva.adjoint.shtp:
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9000data_source.H crec=9000data_prep.H bvel=9000vel2_brad.H \
	dslw=deltaS_inverted_brad_shpt.H dimg=9005Dimg_shpt_tang.H oper=phsft nws=193 adj=y ntpx=15 radian=n offset=y comsou=n \
	image_xmin=0.08 image_xmax=4.92

wemva.forward.enc0:
	< 9001areal_rec_enc0.H Transp plane=13 > 9001areal_rec_temp0.H; echo o5=0. >> 9001areal_rec_temp0.H
	< 9001areal_sou_enc0.H Transp plane=13 > 9001areal_sou_temp0.H; echo o5=0. >> 9001areal_sou_temp0.H
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9001areal_sou_temp0.H crec=9001areal_rec_temp0.H bvel=9000vel2_brad.H \
	dslw=deltaS_brad.H dimg=9005Dimg_enc0_tang.H oper=phsft nws=193 adj=n ntpx=15 radian=n offset=n comsou=y image_xmin=0.08 image_xmax=4.92

wemva.adjoint.enc0:
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9001areal_sou_temp0.H crec=9001areal_rec_temp0.H bvel=9000vel2_brad.H \
	dslw=deltaS_inverted_brad0.H dimg=9005Dimg_enc0_tang.H oper=phsft nws=193 adj=y ntpx=15 radian=n offset=n comsou=y image_xmin=0.08 image_xmax=4.92

wemva.forward:
	< 9001areal_rec_encz_0001.H Transp plane=13 > 9001areal_rec_tempz.H; echo o5=0. >> 9001areal_rec_tempz.H
	< 9001areal_sou_encz_0001.H Transp plane=13 > 9001areal_sou_tempz.H; echo o5=0. >> 9001areal_sou_tempz.H
	<deltaS.H Transp plane=12 reshape=1,3 > deltaS_brad.H
	<deltaS1.H Transp plane=12 reshape=1,3 > deltaS1_brad.H
	Math file1=9000vel2_brad.H exp="1/file1" > 9000slo2_brad.H
	/net/koko/claudio/Research/Others/yaxun/wemva3d_baspix.x csou=9001areal_sou_tempz.H crec=9001areal_rec_tempz.H \
	bslw=9000slo2_brad.H \
	dslw=deltaS_brad.H dimg=9005Dimg_encz_tang.H oper=phsft nws=193 adj=n ntpx=15 radian=n offset=n comsou=y \
	image_xmin=0.08 image_xmax=4.92

wemva.adjoint:
	/net/koko/claudio/Research/Others/yaxun/wemva3d_baspix.x csou=9001areal_sou_tempz.H crec=9001areal_rec_tempz.H \
	bslw=9000slo2_brad.H \
	dslw=deltaS_inverted_brad.H dimg=9005Dimg_encz_tang.H oper=phsft nws=193 adj=y ntpx=15 radian=n offset=n comsou=y \
	image_xmin=0.08 image_xmax=4.92

wemva.forward.5r:
	Cat 9001areal_rec_encz_0001.H 9001areal_rec_encz_0011.H 9001areal_rec_encz_0021.H \
	9001areal_rec_encz_0031.H 9001areal_rec_encz_0041.H axis=4 | Transp plane=13 > 9001areal_rec_tempz.5r.H
	echo o5=0. >> 9001areal_rec_tempz.5r.H
	Cat 9001areal_sou_encz_0001.H 9001areal_sou_encz_0011.H 9001areal_sou_encz_0021.H \
	9001areal_sou_encz_0031.H 9001areal_sou_encz_0041.H axis=4 | Transp plane=13 > 9001areal_sou_tempz.5r.H
	echo o5=0. >> 9001areal_sou_tempz.5r.H
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9001areal_sou_tempz.5r.H crec=9001areal_rec_tempz.5r.H bvel=9000vel2_brad.H \
	dslw=deltaS_brad.H dimg=9005Dimg_encz_tang.5r.H oper=phsft nws=193 adj=n ntpx=15 radian=n offset=n comsou=y image_xmin=0.08 image_xmax=4.92

wemva.adjoint.5r:
	/net/koko/tang/bin/LINUX/wemva3d.x csou=9001areal_sou_tempz.5r.H crec=9001areal_rec_tempz.5r.H bvel=9000vel2_brad.H \
	dslw=deltaS_inverted_brad.5r.H dimg=9005Dimg_encz_tang.5r.H oper=phsft nws=193 adj=y ntpx=15 radian=n offset=n comsou=y image_xmin=0.08 image_xmax=4.92

9001enc0: # todas as frequencis com a mesma fase que varia em x,z
	< 9000_mig_brad.H Window3d min1=1 max1=4 min3=-.1 n3=11 | Transp plane=13 > 9000_mig_brad_transp.H
	./../bin/OW_arealsht.x rec_cig=9000_mig_brad_transp.H sou_cig=9000_mig_brad_transp.H sou=9000data_source.H \
	vel=9000vel2.H rec_out=9001areal_rec_enc0.H sou_out=9001areal_sou_enc0.H ngap=17 rndcte= nrefvel=1 vsal=6 \
	padtraces=41 node= 1 encode=0 combine=1> /dev/null 

9001mig_enc0.H: # correct velocity areal shot random encode in w
	< 9001areal_rec_enc0.H Transp plane=13 > 9001areal_rec_tempz0.H
	< 9001areal_sou_enc0.H Transp plane=13 > 9001areal_sou_tempz0.H
	< 9001areal_rec_tempz0.H ./../bin/OWmig2d_areal.x sou=9001areal_sou_tempz0.H \
	vel=9000vel2.H nrefvel=1 xmin=-1.22 xmax=6.2 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=11 oh=-.1 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 deltaV= ws= xxx_GS.H wr= xxx_GR.H > $@ 
	Rm 9000*tempz0.H