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

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

9001encz%: # 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_encz_$*.H sou_out=9001areal_sou_encz_$*.H ngap=17 rndcte=$* nrefvel=1 vsal=6 \
	padtraces=41 node= 1 encode=1 encz=1 random=1 eps=0.75 combine=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=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 9001*tempz$*.H

9000vel2_lloyd.H: # it is slowness !!!!!!
	Lloyd_vel < 9000vel2_brad.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.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

##############################################
#
#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 > $@ 

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

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

##########################################
#
# <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