4102areal: # correct velocity areal shot random encode in w
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=4102areal_synth_rec.H sou_out=4102areal_synth_sou.H ngap=41 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 

4130arealmig_vel.H: # correct velocity areal shot random encode in w
	< 4102areal_synth_rec.H Transp plane=13 > 4102areal_synth_rec_pad.H
	< 4102areal_synth_sou.H Transp plane=13 > 4102areal_synth_sou_pad.H
	< 4102areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4102areal_synth_sou_pad.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 std=2 dt=0.008 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ 

4103areal: # correct velocity areal shot random encode in z
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=4103areal_synth_rec.H sou_out=4103areal_synth_sou.H ngap=41 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 combine=1 > /dev/null 

4140arealmig_vel.H: # correct velocity areal shot random encode in w
	< 4103areal_synth_rec.H Transp plane=13 > 4103areal_synth_rec_pad.H
	< 4103areal_synth_sou.H Transp plane=13 > 4103areal_synth_sou_pad.H
	< 4103areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4103areal_synth_sou_pad.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 std=2 dt=0.008 encode=1 encw=1 encz=0 ws= xxx_GS.H wr= xxx_GR.H > $@ 

4104areal: # correct velocity areal shot random encode in z
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=4104areal_synth_rec.H sou_out=4104areal_synth_sou.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 combine=1> /dev/null 

4150arealmig_vel.H: # correct velocity areal shot random encode in w
	< 4104areal_synth_rec.H Transp plane=13 > 4104areal_synth_rec_pad.H
	< 4104areal_synth_sou.H Transp plane=13 > 4104areal_synth_sou_pad.H
	< 4104areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4104areal_synth_sou_pad.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 std=2 dt=0.008 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ 

41041areal: # correct velocity areal shot no encoding ngap=11
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=41041areal_synth_rec.H sou_out=41041areal_synth_sou.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 encw= 0 encz= 1 combine=1> /dev/null 

41501arealmig_vel.H: # correct velocity areal shot no random encode in w
	< 41041areal_synth_rec.H Transp plane=13 > 41041areal_synth_rec_pad.H
	< 41041areal_synth_sou.H Transp plane=13 > 41041areal_synth_sou_pad.H
	< 41041areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=41041areal_synth_sou_pad.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 std=2 dt=0.008 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ 

41042areal: # correct velocity areal shot random encode in z
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=41042areal_synth_rec.H sou_out=41042areal_synth_sou.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 combine=1> /dev/null 

41502arealmig_vel.H: # correct velocity areal shot random encode in z
	< 41042areal_synth_rec.H Transp plane=13 > 41042areal_synth_rec_pad.H
	< 41042areal_synth_sou.H Transp plane=13 > 41042areal_synth_sou_pad.H
	< 41042areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=41042areal_synth_sou_pad.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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ 

4105areal: # correct velocity areal shot random encode in z
#	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=4105areal_synth_rec.H sou_out=4105areal_synth_sou.H ngap=41 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 t0=10. combine=1 > /dev/null 

4160arealmig_vel.H: # correct velocity areal shot random encode in w
	< 4105areal_synth_rec.H Transp plane=13 > 4105areal_synth_rec_pad.H
	< 4105areal_synth_sou.H Transp plane=13 > 4105areal_synth_sou_pad.H
	< 4105areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4105areal_synth_sou_pad.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 std=2 dt=0.008 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ 

#!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
#
# Example of how to diminish cross-talk between reflectors by, similarly to what Biondo did 
# when time windowing the crosscorrelation in the imaging condition for 2Way, summing all
# frequencies of source and receiver wavefields prior to imaging condition. Imaging condition
# applied after propagating all frequencies.
#
4300areal:
	< 4008_mig_brad.H Window3d min1=0 max1=4 | Transp plane=13 > 4008_mig_brad_transp.H
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source.H vel=4000vel.H \
	rec_out=4301areal_synth_rec.H sou_out=4301areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \
	traces_in_taper=0 encode=0 > /dev/null 

4351arealmig_vel.H: # for receiver, upward = true
#	< 4301areal_synth_rec.H Pad beg1=200 end1=200 > 4301areal_synth_rec_pad.H
#	< 4301areal_synth_sou.H Pad beg1=200 end1=200 > 4301areal_synth_sou_pad.H
	< 4301areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4301areal_synth_sou_pad.H \
	vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 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 std=2 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ 

4351arealmig_vel_ang.H:
	Window3d < 4351arealmig_vel.H min3=0 max3=4 n4=1 f4=40| Transp plane=34 | Scale scale_to=1 | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=1 > $@

#
# Now, with wrong velocity
# 
4006areal_synth_sou.H: # For source wvfld use upgoing=0
#	< 4006_rot_sig_sou.H Transp plane=23 > 4006_rot_sig_sou_transp.H
	< 4006_rot_sig_sou_transp.H ./bin/OW_arealsht1.x sou=4004data_recip_source.H \
	vel=4000velow.H GFinput=4006_GGG.H upgoing=0 \
	ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6000 traces_in_taper=0 encode=0 encfreq=0 gfopt=0 > $@

4006areal_synth_rec.H: # For receiver wvfld use upgoing=1
#	< 4006_rot_sig_rec.H Transp plane=23 > 4006_rot_sig_rec_transp.H
	< 4006_rot_sig_rec_transp.H ./bin/OW_arealsht1.x sou=4004data_recip_source.H \
	vel=4000velow.H GFinput=4006_GGG.H upgoing=1 \
	ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6000 traces_in_taper=0 encode=0 encfreq=0 gfopt=0 > $@

4356arealmig_velow.H:
	< 4006areal_synth_rec.H Pad beg1=200 end1=200 | Transp plane=24 > 4006areal_synth_rec_pad.H
	< 4006areal_synth_sou.H Pad beg1=200 end1=200 | Transp plane=24 > 4006areal_synth_sou_pad.H
	< 4006areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4006areal_synth_sou_pad.H \
	vel=4000velow.H nrefvel=1 xmin=-2.4 xmax=7.25 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 std=2 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ 

43561arealmig_velow.H:
	< 4006areal_synth_rec.H Pad beg1=200 end1=200 | Transp plane=24 > 4006areal_synth_rec_pad.H
	< 4006areal_synth_sou.H Pad beg1=200 end1=200 | Transp plane=24 > 4006areal_synth_sou_pad.H
	< 4006areal_synth_rec_pad.H ./bin/OWmig2d_areal.x sou=4006areal_synth_sou_pad.H \
	vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 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 std=2 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ 

#4352arealmig_vel.H:
#	< 4301areal_synth_rec.H Pad beg1=200 end1=200 > 4301areal_synth_rec_pad.H
#	< 4301areal_synth_sou.H Pad beg1=200 end1=200 > 4301areal_synth_sou_pad.H
#	< 4301areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4301areal_synth_sou_pad.H \
#	vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
#	traces_in_taper=0 every=1 > $@ 


#####################################
#
# Sigsbee
#
sigareal:
	< ../2Wpar/04sig_mig_brad.H Window3d min1=12000 max1=32000 max5=20000 | Transp plane=13 > 5000temp.H
	./bin/OW_arealsht1.x rec_cig=5000temp.H sou_cig=5000temp.H sou=../2Wpar/02sig_sou1.H \
	vel=1000velmig_my_areal_ldtz.H rec_out=5001sig_rec.H sou_out=5001sig_sou.H GreenFct=5001sig_GF.H \ # GFinput=GGG.H \
	ngap=33 nrefvel=4 xmin=12000 xmax=62000 ymin=0 ymax=0 vsal=14000 traces_in_taper=100 encode=0 encfreq=0 gfopt=1 > /dev/null
	./bin/OW_arealsht1.x rec_cig=5000temp.H sou_cig=5000temp.H sou=../2Wpar/02sig_sou1.H \
	vel=1000velmig_my_areal_ldtz.H rec_out=5001sig_rec.H sou_out=5001sig_sou.H GFinput=50011sig_GF.H \
	ngap=33 nrefvel=4 xmin=12000 xmax=62000 ymin=0 ymax=0 vsal=14000 traces_in_taper=100 encode=0 encfreq=0 gfopt=0 > /dev/null

sigareal1:
	Window3d <../2Wpar/1000velmig_transp_z16000.H | Transp > 1000velmig_z16000.H
	< ../2Wpar/04sig_mig_brad_cut.H Window3d min1=10925 max1=30000 | Pad beg3=220 | Transp plane=13 > 5000temp.H
	./bin/OW_arealsht1.x rec_cig=5000temp.H sou_cig=5000temp.H sou=../2Wpar/02sig_sou1_cut.H \
	vel=1000velmig_z16000.H rec_out=5001sig_rec.H sou_out=5001sig_sou.H GreenFct=5001sig_GF.H \
	ngap=33 nrefvel=4 xmin=10925 xmax=30000 ymin=0 ymax=0 vsal=14700 traces_in_taper=100 encode=0 encfreq=0 gfopt=1 > /dev/null
	./bin/OW_arealsht1.x rec_cig=5000temp.H sou_cig=5000temp.H sou=../2Wpar/02sig_sou1_cut.H \
	vel=1000velmig_z16000.H rec_out=5001sig_rec.H sou_out=5001sig_sou.H GFinput=5001sig_GF.H \
	ngap=33 nrefvel=4 xmin=10925 xmax=30000 ymin=0 ymax=0 vsal=14700 traces_in_taper=100 encode=0 encfreq=0 gfopt=0 > /dev/null
#####################################
#
# Marmousi
#
marmousi1:
	< ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3000temp.H
	<3000temp.H ./bin/OW_arealsht1.x sou=../2Wpar/8102marm_source_brad.H \
	vel=../2Wpar/8000marmvel_24m.H rec_out=3001marm_rec.H sou_out=3001marm_sou.H GreenFct=3001marm_GF.H \
	ngap=33 nrefvel=4 xmin=1000 xmax=8000 ymin=0 ymax=0 vsal=6000 traces_in_taper=100 encode=0 encfreq=0 gfopt=1 > /dev/null

3001marm_rec.H: # For receiver wvfld use upgoing=1
	< ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3000temp.H
	< 3000temp.H ./bin/OW_arealsht1.x sou=../2Wpar/8102marm_source_brad.H \
	vel=../2Wpar/8000marmvel_24m.H GFinput=3001marm_GF.H upgoing=1 \
	ngap=33 nrefvel=4 xmin=1000 xmax=8000 ymin=0 ymax=0 vsal=6000 traces_in_taper=100 encode=0 encfreq=0 gfopt=0 > $@

3001marm_sou.H: # For source wvfld use upgoing=0
	< 3000temp.H ./bin/OW_arealsht1.x sou=../2Wpar/8102marm_source_brad.H \
	vel=../2Wpar/8000marmvel_24m.H GFinput=3001marm_GF.H upgoing=0 \
	ngap=33 nrefvel=4 xmin=1000 xmax=8000 ymin=0 ymax=0 vsal=6000 traces_in_taper=100 encode=0 encfreq=0 gfopt=0 > $@

3002marm_arealmig%.H:
	< 3001marm_rec.H Window3d n2=11 f2=$* squeeze=n | Pad beg1=100 end1=100 | Transp plane=24 > 3001marm_rectemp$*.H
	< 3001marm_sou.H Window3d n2=11 f2=$* squeeze=n | Pad beg1=100 end1=100 | Transp plane=24 > 3001marm_soutemp$*.H
	< 3001marm_rectemp$*.H ./bin/OWmig2d_areal.x sou=3001marm_soutemp$*.H \
	vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-1800 xmax=11664 ymin=0 ymax=0 zmin=0 zmax=3000. nh=1 81 oh=0 -.4 dh=24. wbott=0 \
	vsal=6000 traces_in_taper=100 every=1 deltat=0.02 std=2 dt=0.008 ws=3002marm_$*_GS.H wr=3002marm_$*_GR.H velout= velout.H > $@ 

3003marm_arealmig_full.H:
	Window3d < 3002marm_arealmig00.H n4=1 f4=10 min3=0 max3=9000 > 3003marm_z00.H
	Window3d < 3002marm_arealmig11.H n4=1 f4=10 min3=0 max3=9000 > 3003marm_z11.H
	Window3d < 3002marm_arealmig22.H n4=1 f4=10 min3=0 max3=9000 > 3003marm_z22.H
	Math file1=3003marm_z00.H file2=3003marm_z11.H file3=3003marm_z22.H \
	exp="file1+file2+file3" > $@

3004marm_arealmig%.H:
	< 3001marm_rectemp$*.H ./bin/OWmig2d_areal.x sou=3001marm_soutemp$*.H \
	vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-1800 xmax=11664 ymin=0 ymax=0 zmin=0 zmax=3000. nh=1 81 oh=0 -.4 dh=24. wbott=0 \
	vsal=6000 traces_in_taper=100 every=1 deltat=0.004 std=2 dt=0.004 ws= 3002marm_$*_GS.H wr= 3002marm_$*_GR.H velout= velout.H > $@ 

3005marm_arealmig_full.H:
	Window3d < 3004marm_arealmig00.H n4=1 f4=10 min3=0 max3=9000 > 3005marm_z00.H
	Window3d < 3004marm_arealmig11.H n4=1 f4=10 min3=0 max3=9000 > 3005marm_z11.H
	Window3d < 3004marm_arealmig22.H n4=1 f4=10 min3=0 max3=9000 > 3005marm_z22.H
	Math file1=3005marm_z00.H file2=3005marm_z11.H file3=3005marm_z22.H \
	exp="file1+file2+file3" > $@
		
3006marm_arealmig%.H:
	< 3001marm_rectemp$*.H ./bin/OWmig2d_areal.x sou=3001marm_soutemp$*.H \
	vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-1800 xmax=11664 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-400 dh=24. wbott=0 \
	vsal=6000 traces_in_taper=100 every=1 deltat=0.02 std=2 dt=0.004 ws=3002marm_$*_GS.H wr=3002marm_$*_GR.H velout= velout.H > $@ 

3007marm_arealmig_full.H:
	Window3d < 3006marm_arealmig00.H n4=1 f4=10 min3=0 max3=9000 > 3007marm_z00.H
	Window3d < 3006marm_arealmig11.H n4=1 f4=10 min3=0 max3=9000 > 3007marm_z11.H
	Window3d < 3006marm_arealmig22.H n4=1 f4=10 min3=0 max3=9000 > 3007marm_z22.H
	Math file1=3007marm_z00.H file2=3007marm_z11.H file3=3007marm_z22.H \
	exp="file1+file2+file3" > $@