1001data_hor.H:
	susynlv ref="1:0,.5.;10.,.5" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=20 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > $@; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

1002data_hor.H:
	susynlv ref="1:0,0.;10.,2.5" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=20 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk2.su
	< junk2.su Su2sep > junk2.H
	<junk2.H Infill3d |Window3d j1=4 > $@; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

1001data_recip_prep.H:
	< 1001data_hor.H Window3d |  /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=48 maxsize=1000 > j1.H
	< j1.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@

1002data_recip_prep.H:
	< 1002data_hor.H Window3d |  /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=48 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@

1000data_recip_source.H:
	Window3d < 1001data_hor.H n3=1 f3=100 n2=1 min2=1 > x.H
	< x.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=48 maxsize=1000 > y.H
	< y.H Transp plane=35 | Cabs | Stack |Window3d > 1000_spec.H;
	Math file1=1000_spec.H exp="file1*0" > 1000_phse.H; Cmplx 1000_spec.H 1000_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

1001_mig_brad.H: # With 201 shots separated 0.02km correct velocity
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=1001data_recip_prep.H D_d=1000data_recip_source.H V=4000vel_brad_pad.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=0 \
	par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y

1002_mig_brad.H: # With 201 shots separated 0.02km correct velocity
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=1002data_recip_prep.H D_d=1000data_recip_source.H V=4000vel_brad_pad.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=0 \
	par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y

1001areal: # correct velocity areal shots
	< 1001_mig_brad.H Window3d min1=0 max1=4 max5=1. | Transp plane=13 > 1001_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=1001_mig_brad_transp.H sou_cig=1001_mig_brad_transp.H sou=1000data_recip_source.H vel=4000vel.H \
	rec_out=1001areal_rec_comb1_enc0_vel.H sou_out=1001areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 wwr=1.H wws=2.H> /dev/null 

1001mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15
	< 1001areal_rec_comb1_enc0_vel.H Transp plane=13 > 1001areal_rec_temp.H
	< 1001areal_sou_comb1_enc0_vel.H Transp plane=13 > 1001areal_sou_temp.H
	< 1001areal_rec_temp.H ./bin/OWmig2d_areal.x sou=1001areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws=1001_GS.H wr=1001_GR.H > $@ 
	Rm 1001areal*temp.H

1002areal: # correct velocity areal shots
	< 1002_mig_brad.H Window3d min1=0 max1=4 max5=1. | Transp plane=13 > 1002_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=1002_mig_brad_transp.H sou_cig=1002_mig_brad_transp.H sou=1000data_recip_source.H vel=4000vel.H \
	rec_out=1002areal_rec_comb1_enc0_vel.H sou_out=1002areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 wwr=1.H wws=2.H> /dev/null 

1002mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15
	< 1002areal_rec_comb1_enc0_vel.H Transp plane=13 > 1002areal_rec_temp.H
	< 1002areal_sou_comb1_enc0_vel.H Transp plane=13 > 1002areal_sou_temp.H
	< 1002areal_rec_temp.H ./bin/OWmig2d_areal.x sou=1002areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws=1002_GS.H wr=1002_GR.H > $@ 
	Rm 1002areal*temp.H

1003mig_combine.H:
	Math file1=1001mig_comb1_enc0_vel.H file2=1002mig_comb1_enc0_vel.H exp=file1+file2 > $@

#
#
#
1011data_hor.H:
	susynlv ref="1:0,.5.;10.,.5" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=32 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > $@; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> $@
	susynlv ref="1:0,.75;10.,.75" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=32 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d | Window3d j1=4 > $@; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

1011data_recip_prep.H:
	< 1011data_hor.H Window3d |  /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=48 f_max=60 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@

1011data_recip_source.H:
	< junk2.H Infill3d  | Window3d j1=4 > junk3.H; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> junk3.H
	Window3d < junk3.H n3=1 f3=100 n2=1 min2=1 > x.H
	< x.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=2 f_min1=6 f_max1=48 f_max=60 maxsize=1000 > y.H
	< y.H Transp plane=35 | Cabs | Stack |Window3d > 1000_spec.H;
	Math file1=1000_spec.H exp="file1*0" > 1000_phse.H; Cmplx 1000_spec.H 1000_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

1011_mig_brad.H: # With 201 shots separated 0.02km correct velocity
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=1011data_recip_prep.H D_d=1011data_recip_source.H V=4000vtemp_0.02.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=0 \
	par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y

1011areal: # correct velocity areal shots
	< 1011_mig_brad.H Window3d min1=0 max1=4 max5=1. | Transp plane=13 > 1011_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=1011_mig_brad_transp.H sou_cig=1011_mig_brad_transp.H sou=1011data_recip_source.H vel=4000velow.H \
	rec_out=1011areal_rec_comb1_enc0_vel.H sou_out=1011areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 wwr=1.H wws=2.H> /dev/null 

1011mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15
	< 1011areal_rec_comb1_enc0_vel.H Transp plane=13 > 1011areal_rec_temp.H
	< 1011areal_sou_comb1_enc0_vel.H Transp plane=13 > 1011areal_sou_temp.H
	< 1011areal_rec_temp.H ./bin/OWmig2d_areal.x sou=1011areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws=1011_GS.H wr=1011_GR.H > $@ 
	Rm 1011areal*temp.H

1012areal: # correct velocity areal shots
	< 1011_mig_brad.H Window3d min1=2 n1=1 max5=1. | Transp > 1011_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=1011_mig_brad_transp.H sou_cig=1011_mig_brad_transp.H sou=1011data_recip_source.H vel=4000velow.H \
	rec_out=1012areal_rec_comb1_enc0_vel.H sou_out=1012areal_sou_comb1_enc0_vel.H ngap=1 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 wwr=1.H wws=2.H> /dev/null 

1012mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15
	< 1012areal_rec_comb1_enc0_vel.H Transp plane=13 > 1012areal_rec_temp.H
	< 1012areal_sou_comb1_enc0_vel.H Transp plane=13 > 1012areal_sou_temp.H
	< 1012areal_rec_temp.H ./bin/OWmig2d_areal.x sou=1012areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws=1012_GS.H wr=1012_GR.H > $@ 
	Rm 1012areal*temp.H

#
#
#