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

9200data_hor.H:
	susynlv ref="1:0,.5.;5.,.5." nt=2001 dt=0.002 dxs=0.01 nxs=301 fpeak=20 fxs=0 dxo=0.01 nxo=401 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > x1.H; echo n2=401 o2=-2 d2=0.01 n3=301 o3=0 d3=0.01 >> x1.H
	susynlv ref="1:0,1.;5.,1." nt=2001 dt=0.002 dxs=0.01 nxs=301 fpeak=20 fxs=0 dxo=0.01 nxo=401 fxo=-2 v00=2.5 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > x2.H; echo n2=401 o2=-2 d2=0.01 n3=301 o3=0 d3=0.01 >> x2.H
	Math file1=x1.H file2=x2.H exp=file1+file2 >> $@

#
# Window3d <9200data_hor.H n3=101 f3=100 > 92001data_hor.H; echo o3=0 >> 92001data_hor.H
#
9200data_prep.H:
	< 9200data_hor.H  /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=32 f_max=40 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

9200_mig_brad.H: # With 401 shots separated 0.01km correct velocity
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=9200data_prep.H D_d=9100data_source.H V=4000velow_brad_pad.H R=$@ \
	pverb=2  nx_shots=301 ox_shots=0 \
	par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-4 maxx=7.2 & #restart=y

9200areal: # correct velocity areal shots
	< 9200_mig_brad.H Window3d min1=1.75 n1=1 | Transp > 9201_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=9201_mig_brad_transp.H sou_cig=9201_mig_brad_transp.H sou=4004data_recip_source_pad.H  vel=4000vel10m.H \
	rec_out=9201areal_rec_comb1_enc0_vel.H sou_out=9201areal_sou_comb1_enc0_vel.H ngap=1 nrefvel=1 vsal=6 \
	padtraces=121 encode=0 combine=1 wwr=9201wwr.H wws=9201wws.H > /dev/null 
#	./bin/OW_arealsht.x rec_cig=9201_mig_brad_transp.H sou_cig=9201_mig_brad_transp.H sou=9100data_source.H vel=4000vel10m.H \

9200mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 
	< 9201areal_rec_comb1_enc0_vel.H Transp plane=13 > 9200areal_rec_temp.H
	< 9201areal_sou_comb1_enc0_vel.H Transp plane=13 > 9200areal_sou_temp.H
	< 9200areal_rec_temp.H ./bin/OWmig2d_areal.x sou=9200areal_sou_temp.H \
	vel=4000vel10m.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=161 oh=-.8 dh=0.01 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws=9200migws.H wr=9200migwr.H > $@ 
	Rm 9200areal*temp.H