##################################### # # 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 ./bin/OW_arealsht1.x rec_cig=3000temp.H sou_cig=3000temp.H 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 rec_cig=3000temp.H sou_cig=3000temp.H 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" > $@ #################################### # # Plane reflector 1km/s # 4001data_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 # 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 > 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 >> $@ 4002data_recip.H: # Create3d < 4001data_hor.H keyname1=offset keyname2=sx verb=1 > junk.H # Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H # Sort3d < junk1.H nkeys=2 key2=cmp ng2=1011 og2=0 dg2=0.01 key1=offset ng1=101 og1=0 dg1=0.02 max_size=1000 > junk2.H # junk3.H; # junk30.H Interp d2out=0.01 o2out=0 n2out=201 maxsize=1000 junk31.H Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H Merge axis=2 space=n junk4.H junk31.H > junk5.H Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H Sort3d < junk7.H nkeys=2 key2=sx ng2=401 og2=0 dg2=0.01 key1=offset ng1=401 og1=-2 dg1=0.01 max_size=1000 > junk8.H junk9.H $@ 4003data_recip_prep.H: < 4001data_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=48 maxsize=1000 > j2.H < j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@ 4004data_recip_source.H: Window3d < 4001data_hor.H n3=1 f3=200 | Mute vmute=2 tmute=0.7 | Window3d min2=-1 | \ /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 | \ Window3d | Transp | Cabs | Stack > 4004_spec.H; Math file1=4004_spec.H exp="file1*0" > 4004_phse.H; Cmplx 4004_spec.H 4004_phse.H |\ Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@ 4006_mig_brad.H: # With 201 shots separated 0.02km slower velocity Window3d squeeze=n <4000velow_brad_pad.H j1=2 > 4000vtemp_0.02.H @date /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=4003data_recip_prep.H D_d=4004data_recip_source.H V=4000vtemp_0.02.H R=$@ \ pverb=2 nx_shots=201 ox_shots=0 \ par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y 4006_mig_brad_ang.H: Window3d < 4006_mig_brad.H min1=0 max1=4 max5=1 | Transp plane=13 | 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; $@ 4008_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=4003data_recip_prep.H D_d=4004data_recip_source.H V=4000vel_brad_pad.H R=$@ \ pverb=2 nx_shots=201 ox_shots=0 \ par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y 4006rot_vlow: < 4006_mig_brad.H Window3d min1=0 max1=4.02 | Transp plane=12 reshape=2,3 > 4006temp.H ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=4006temp.H Cig_tilt=4006_rot_sig_sou.H par=4000Cig-rot.P positive_rot=1 ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=4006temp.H Cig_tilt=4006_rot_sig_rec.H par=4000Cig-rot.P positive_rot=0 4008rot: < 4008_mig_brad.H Window3d min1=0 max1=4.02 | Transp plane=12 reshape=2,3 > 4008temp.H ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=4008temp.H Cig_tilt=4008_rot_sig_sou.H par=4000Cig-rot.P positive_rot=1 ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=4008temp.H Cig_tilt=4008_rot_sig_rec.H par=4000Cig-rot.P positive_rot=0 4006areal: # < 4006_rot_sig_sou.H Transp plane=23 > 4006_rot_sig_sou_transp.H # < 4006_rot_sig_rec.H Transp plane=23 > 4006_rot_sig_rec_transp.H # ./bin/OW_arealsht1.x rec_cig=4006_rot_sig_rec_transp.H sou_cig=4006_rot_sig_sou_transp.H sou=4004data_recip_source.H \ # vel=4000velow.H rec_out=4006areal_synth_rec.H sou_out=4006areal_synth_sou.H GreenFct=4006_GGG.H \ # ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=0 encfreq=1 gfopt=1 > /dev/null ./bin/OW_arealsht1.x rec_cig=4006_rot_sig_rec_transp.H sou_cig=4006_rot_sig_sou_transp.H sou=4004data_recip_source.H \ vel=4000velow.H rec_out=4006areal_synth_rec.H sou_out=4006areal_synth_sou.H GFinput=4006_GGG.H \ ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=0 gfopt=0 > /dev/null 4006mig01.H: < 4006areal_synth_rec.H Transp plane=24 > 4006areal_synth_rec_transp.H < 4006areal_synth_sou.H Transp plane=24 > 4006areal_synth_sou_transp.H < 4006areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4006areal_synth_sou_transp.H \ vel=4000velow.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 > $@ 4008areal: < 4008_rot_sig_sou.H Transp plane=23 > 4008_rot_sig_sou_transp.H < 4008_rot_sig_rec.H Transp plane=23 > 4008_rot_sig_rec_transp.H # ./bin/OW_arealsht1.x rec_cig=4008_rot_sig_rec_transp.H sou_cig=4008_rot_sig_sou_transp.H sou=4004data_recip_source.H \ # vel=4000vel.H rec_out=4008areal_synth_rec.H sou_out=4008areal_synth_sou.H GreenFct=4008_GGG.H \ # ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=0 encfreq=1 gfopt=1 > /dev/null ./bin/OW_arealsht1.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=4008areal_synth_rec.H sou_out=4008areal_synth_sou.H GFinput=4008_GGG.H \ ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=0 gfopt=0 > /dev/null 40081areal: # GF with diff rnd phase/freq ./bin/OW_arealsht1.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=4008areal_synth_rec_rndphase.H sou_out=4008areal_synth_sou_rndphase.H GFinput=4008_GGG.H \ ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=2 gfopt=0 > /dev/null 40082areal: # GF with equal rnd phase all freq ./bin/OW_arealsht1.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=4008areal_synth_rec_rndphase2.H sou_out=4008areal_synth_sou_rndphase2.H GFinput=4008_GGG.H \ ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=100 encode=2 gfopt=0 > /dev/null 4008mig01.H: < 4008areal_synth_rec.H Transp plane=24 > 4008areal_synth_rec_transp.H < 4008areal_synth_sou.H Transp plane=24 > 4008areal_synth_sou_transp.H < 4008areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4008areal_synth_sou_transp.H \ vel=4000vel.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 > $@ 4008mig02.H: < 4008areal_synth_rec_rndphase.H Transp plane=24 > 4008areal_synth_rec_rndphase_transp.H < 4008areal_synth_sou_rndphase.H Transp plane=24 > 4008areal_synth_sou_rndphase_transp.H < 4008areal_synth_rec_rndphase_transp.H ./bin/OWmig2d_areal.x sou=4008areal_synth_sou_rndphase_transp.H \ vel=4000vel.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 > $@ 4008mig03.H: < 4008areal_synth_rec_rndphase2.H Transp plane=24 > 4008areal_synth_rec_rndphase2_transp.H < 4008areal_synth_sou_rndphase2.H Transp plane=24 > 4008areal_synth_sou_rndphase2_transp.H < 4008areal_synth_rec_rndphase2_transp.H ./bin/OWmig2d_areal.x sou=4008areal_synth_sou_rndphase2_transp.H \ vel=4000vel.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 > $@ 4100areal_synth_sou.H: < 4008_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H cigset=4101cigset.H \ nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 encode=0 > $@ 4100areal_synth_rec.H: < 4008_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H \ nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 4101areal_synth_rec.H: < 4008_rot_sig_rec.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H \ nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 4150arealmig_velow.H: # < 4100areal_synth_rec.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4105areal_synth_rec_transp.H # < 4100areal_synth_sou.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4105areal_synth_sou_transp.H < 4105areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4105areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-0.4 xmax=5.25 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 4159arealmig_velow.H: # < 4100areal_synth_rec.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4105areal_synth_rec_transp.H # < 4100areal_synth_sou.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4105areal_synth_sou_transp.H < 4105areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4105areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=81 oh=-.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 41591arealmig_velow.H: < 4100areal_synth_rec.H Window3d min3=-0.4 max3=5.25 f4=29 squeeze=n| Transp plane=13 > 4105areal_synth_rec_transp.H < 4100areal_synth_sou.H Window3d min3=-0.4 max3=5.25 f4=29 squeeze=n| Transp plane=13 > 4105areal_synth_sou_transp.H < 4105areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4105areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-0.4 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=81 oh=-.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 4151arealmig_velow.H: # < 4101areal_synth_rec.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4106areal_synth_rec_transp.H # < 4100areal_synth_sou.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 > 4106areal_synth_sou_transp.H < 4106areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4106areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-0.4 xmax=5.25 ymin=0 ymax=0 zmin=0 zmax=1.35 nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 4200areal_synth_sou.H: < 4008_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H cigset=4101cigset.H \ ngap=82 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 encode=0 > $@ 4200areal_synth_rec.H: < 4008_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H \ ngap=82 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 4201areal_synth_rec.H: < 4008_rot_sig_rec.H Transp plane=23 | ./bin/OW_arealsht.x \ sou=4004data_recip_source.H vel=4000velow_pad1.H \ ngap=82 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 4250arealmig_velow.H: < 4200areal_synth_rec.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 | Pad beg1=200 end1=200 > 4205areal_synth_rec_transp.H < 4200areal_synth_sou.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 | Pad beg1=200 end1=200 > 4205areal_synth_sou_transp.H < 4205areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4205areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 4251arealmig_velow.H: < 4201areal_synth_rec.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 | Pad beg1=200 end1=200 > 4206areal_synth_rec_transp.H < 4200areal_synth_sou.H Window3d min3=-0.4 max3=5.25 squeeze=n| Transp plane=13 | Pad beg1=200 end1=200 > 4206areal_synth_sou_transp.H < 4206areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4206areal_synth_sou_transp.H \ vel=4000velow_pad1.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 4000velow_pad2.H: Math file1=4000velow_pad1.H exp="file1*0+2" > $@ 4250arealmig_vel.H: < 4205areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=4205areal_synth_sou_transp.H \ vel=4000velow_pad2.H nrefvel=1 xmin=-0.4 xmax=5.25 ymin=0 ymax=0 zmin=0 zmax=1.5 nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.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; $@ # # 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 > $@ #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 4302areal: ./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=4302areal_synth_rec.H sou_out=4302areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=1 > /dev/null #4352arealmig_vel.H: # < 4302areal_synth_rec.H Pad beg1=200 end1=200 > 4302areal_synth_rec_pad.H # < 4302areal_synth_sou.H Pad beg1=200 end1=200 > 4302areal_synth_sou_pad.H # < 4302areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4302areal_synth_sou_pad.H \ # vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ # traces_in_taper=0 every=1 > $@ 4303areal: ./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=4303areal_synth_rec.H sou_out=4303areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=3 t0=1. > /dev/null 4353arealmig_vel.H: < 4303areal_synth_rec.H Pad beg1=200 end1=200 > 4303areal_synth_rec_pad.H < 4303areal_synth_sou.H Pad beg1=200 end1=200 > 4303areal_synth_sou_pad.H < 4303areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4303areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4304areal: ./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=4304areal_synth_rec.H sou_out=4304areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=3 t0=2. > /dev/null 4354arealmig_vel.H: < 4304areal_synth_rec.H Pad beg1=200 end1=200 > 4304areal_synth_rec_pad.H < 4304areal_synth_sou.H Pad beg1=200 end1=200 > 4304areal_synth_sou_pad.H < 4304areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4304areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4306areal: ./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=4306areal_synth_rec.H sou_out=4306areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=3 t0=.1 > /dev/null 4356arealmig_vel.H: < 4306areal_synth_rec.H Pad beg1=200 end1=200 > 4306areal_synth_rec_pad.H < 4306areal_synth_sou.H Pad beg1=200 end1=200 > 4306areal_synth_sou_pad.H < 4306areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4306areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4307areal: ./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=4307areal_synth_rec.H sou_out=4307areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=3 t0=32. > /dev/null 4357arealmig_vel.H: < 4307areal_synth_rec.H Pad beg1=200 end1=200 > 4307areal_synth_rec_pad.H < 4307areal_synth_sou.H Pad beg1=200 end1=200 > 4307areal_synth_sou_pad.H < 4307areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4307areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4308areal: ./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=4308areal_synth_rec.H sou_out=4308areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=2 encfreq=1 t0=32. > /dev/null 4358arealmig_vel.H: < 4308areal_synth_rec.H Pad beg1=200 end1=200 > 4308areal_synth_rec_pad.H < 4308areal_synth_sou.H Pad beg1=200 end1=200 > 4308areal_synth_sou_pad.H < 4308areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4308areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4309areal: ./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=4309areal_synth_rec.H sou_out=4309areal_synth_sou.H ngap=41 nrefvel=1 xmin=0 xmax=4 ymin=0 ymax=0 vsal=6 \ traces_in_taper=0 encode=3 encfreq=1 t0=16. > /dev/null 4359arealmig_vel.H: < 4309areal_synth_rec.H Pad beg1=200 end1=200 > 4309areal_synth_rec_pad.H < 4309areal_synth_sou.H Pad beg1=200 end1=200 > 4309areal_synth_sou_pad.H < 4309areal_synth_rec_pad.H ../2Wpar/bin/OWmig2d_areal.x sou=4309areal_synth_sou_pad.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 every=1 > $@ 4921arealmig_vel.H: < 4920areal_synth_rec.H Transp plane=24 | Pad beg1=200 end1=200 > 4925areal_synth_rec_transp.H < 4920areal_synth_sou.H Transp plane=24 | Pad beg1=200 end1=200 > 4925areal_synth_sou_transp.H < 4925areal_synth_rec_transp.H ../2Wpar/bin/OWmig2d_areal.x sou=4925areal_synth_sou_transp.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 every=0 > $@ 4950arealmig_vel.H: < 4932areal_synth_rec.H Pad beg1=200 end1=200 > 4935areal_synth_rec_transp.H < 4932areal_synth_sou.H Pad beg1=200 end1=200 > 4935areal_synth_sou_transp.H < 4935areal_synth_rec_transp.H ../2Wpar/bin/OWmig2d_areal.x sou=4935areal_synth_sou_transp.H \ vel=4000vel.H nrefvel=1 xmin=-2.4 xmax=7.25 ymin=0 ymax=0 zmin=0 zmax=1. nh=1 81 oh=0 -.4 dh=0.01 wbott=0 vsal=5 \ traces_in_taper=0 every=0 > $@ #################################### # # Plane reflector # 6001data_hor.H: susynlv ref="1:0,1.5;10.,1.5" nt=251 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=201 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk1.H susynlv ref="1:0,-0.;10.,8." nt=251 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=201 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk2.H Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d > $@; echo n2=201 o2=0 d2=0.02 n3=201 o3=0 d3=0.02 >> $@ 6002data_recip.H: Create3d < 6001data_hor.H keyname1=offset keyname2=sx verb=1 > junk.H Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H Sort3d < junk1.H nkeys=2 key2=cmp ng2=611 og2=0 dg2=0.01 key1=offset ng1=101 og1=0 dg1=0.04 max_size=1000 > junk2.H junk3.H; junk30.H Interp d2out=0.02 o2out=0 n2out=201 maxsize=1000 junk31.H Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H Merge axis=2 space=n junk4.H junk31.H > junk5.H Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H Sort3d < junk7.H nkeys=2 key2=sx ng2=201 og2=0 dg2=0.02 key1=offset ng1=401 og1=-4 dg1=0.02 max_size=1000 > junk8.H junk9.H $@ 6003data_recip_prep.H: < 6002data_recip.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \ is2d=y wei=y f_min=6 f_min1=10 f_max1=50 f_max=60 maxsize=1000 > j2.H < j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@ 6004data_recip_source.H: < 6003data_recip_prep.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 6004_spec.H; Math file1=6004_spec.H exp="file1*0" > 6004_phse.H; Cmplx 6004_spec.H 6004_phse.H |\ Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@ 6100areal_synth_d%_a.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3 nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_a.H nphi=10 dphi=2. ophi=-70 encode=0 > $@ 6100areal_synth_d%_b.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_b.H nphi=10 dphi=2. ophi=-50 encode=0 > $@ 6100areal_synth_d%_c.H: ./bin/OW_arealsht.x < 6003data_recip_prep.H sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_c.H nphi=10 dphi=2. ophi=-30 encode=0 > $@ 6100areal_synth_d%_d.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_d.H nphi=10 dphi=2. ophi=-10 encode=0 > $@ 6100areal_synth_d%_e.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_e.H nphi=10 dphi=2. ophi=10 encode=0 > $@ 6100areal_synth_d%_f.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_f.H nphi=10 dphi=2. ophi=30 encode=0 > $@ 6100areal_synth_d%_g.H: < 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \ ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \ synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_g.H nphi=10 dphi=2. ophi=50 encode=0 > $@ 6120areal_mig_d1.5_%.H: < 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 squeeze=n n4=1 | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 squeeze=n n4=1 | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \ vel=7000vel.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -0.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@ 6130areal_mig_d0.500_ltd_full.H: Window3d < 6120areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_a.H Window3d < 6120areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_b.H Window3d < 6120areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_c.H Window3d < 6120areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_d.H Window3d < 6120areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_e.H Math file1=6130areal_mig_d0.500_ltd_a.H file2=6130areal_mig_d0.500_ltd_b.H file3=6130areal_mig_d0.500_ltd_c.H \ file4=6130areal_mig_d0.500_ltd_d.H file5=6130areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@ 6140areal_mig_d0500_ltd_angle.H: Transp plane=34 < 6130areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H < xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@ 7000velow.H: Math file1=7000vel.H exp="file1*0+2.7" > $@ 7000velhi.H: Math file1=7000vel.H exp="file1*0+3.3" > $@ 6220areal_mig_d1.5_%.H: < 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \ vel=7000velow.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 oh=0 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@ 6221areal_mig_d1.5_%.H: < 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H < 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \ vel=7000vel.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 oh=0 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@ 6230areal_mig_d0.500_ltd_full.H: Window3d < 6220areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_a.H Window3d < 6220areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_b.H Window3d < 6220areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_c.H Window3d < 6220areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_d.H Window3d < 6220areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_e.H Math file1=6230areal_mig_d0.500_ltd_a.H file2=6230areal_mig_d0.500_ltd_b.H file3=6230areal_mig_d0.500_ltd_c.H \ file4=6230areal_mig_d0.500_ltd_d.H file5=6230areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@ 6240areal_mig_d0500_angle.H: Transp plane=34 < 6230areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H < xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@ 6320areal_mig_d0.500_%.H: < 6101areal_synth_d0.500_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d0.500_$*_transp.H \ vel=7000velhi.H nrefvel=1 xmin=-6 xmax=10 ymin=0 ymax=0 zmin=0 zmax=1 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 6330areal_mig_d0.500_ltd_full.H: Window3d < 6320areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_a.H Window3d < 6320areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_b.H Window3d < 6320areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_c.H Window3d < 6320areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_d.H Window3d < 6320areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_e.H Math file1=6330areal_mig_d0.500_ltd_a.H file2=6330areal_mig_d0.500_ltd_b.H file3=6330areal_mig_d0.500_ltd_c.H \ file4=6330areal_mig_d0.500_ltd_d.H file5=6330areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@ 6340areal_mig_d0500_ltd_angle.H: Transp plane=34 < 6330areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H < xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@ 6420_mig_brad.H: @date /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=6003data_recip_prep.H D_d=6004data_recip_source.H V=7000vel_brad.H R=$@ \ pverb=2 nx_shots=201 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y 6425_mig_brad.H: @date /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=6003data_recip_prep.H D_d=6004data_recip_source.H V=7000velow_brad.H R=$@ \ pverb=2 nx_shots=201 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y 6501data_hor.H: susynlv ref="1:0,1.5;10.,1.5" nt=501 dt=0.008 dxs=0.02 nxs=401 fxs=0 dxo=0.02 nxo=401 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk1.H susynlv ref="1:0,-0.;10.,8." nt=501 dt=0.008 dxs=0.02 nxs=401 fxs=0 dxo=0.02 nxo=401 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk2.H Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d > $@; echo n2=401 o2=0 d2=0.02 n3=401 o3=0 d3=0.02 >> $@ 6502data_recip.H: Create3d < 6501data_hor.H keyname1=offset keyname2=sx verb=1 > junk.H Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H Sort3d < junk1.H nkeys=2 key2=cmp ng2=1211 og2=0 dg2=0.01 key1=offset ng1=201 og1=0 dg1=0.04 max_size=1000 > junk2.H junk3.H; junk30.H Interp d2out=0.02 o2out=0 n2out=401 maxsize=1000 junk31.H Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H Merge axis=2 space=n junk4.H junk31.H > junk5.H Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H Sort3d < junk7.H nkeys=2 key2=sx ng2=401 og2=0 dg2=0.02 key1=offset ng1=801 og1=-8 dg1=0.02 max_size=1000 > junk8.H junk9.H $@ 6503data_recip_prep.H: < 6502data_recip.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \ is2d=y wei=y f_min=6 f_min1=10 f_max1=50 f_max=60 maxsize=1000 > j2.H < j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@ 6504data_recip_source.H: < 6503data_recip_prep.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 6504_spec.H; Math file1=6504_spec.H exp="file1*0" > 6504_phse.H; Cmplx 6504_spec.H 6504_phse.H |\ Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@ 6525_mig_brad.H: @date /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=6503data_recip_prep.H D_d=6504data_recip_source.H V=7000velow_brad.H R=$@ \ pverb=2 nx_shots=401 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y 6530_mig_brad.H: @date /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=6503data_recip_prep.H D_d=6504data_recip_source.H V=7000vel_brad.H R=$@ \ pverb=2 nx_shots=401 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y #################################### # # Areal shots of rotated CIGs # rot_vel: < 6530_mig_brad.H Window3d min1=-1 max1=5 | Transp plane=12 reshape=2,3 > temp_brad.H ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=temp_brad.H Cig_tilt=6701_rot_sig_sou.H par=Cig-rot800.P positive_rot=1 ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=temp_brad.H Cig_tilt=6702_rot_sig_rec.H par=Cig-rot800.P positive_rot=0 7000areal_synth_sou.H: < 6701_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000vel.H cigset= 7101cigset.H \ nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 encode=0 > $@ 7000areal_synth_rec.H: < 6702_rot_sig_rec.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000vel.H cigset= 7101cigset.H \ nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 7210arealmig_vel.H: < 7000areal_synth_rec.H Transp plane=13 > 7050areal_synth_rec_transp.H < 7000areal_synth_sou.H Transp plane=13 > 7050areal_synth_sou_transp.H < 7050areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=7050areal_synth_sou_transp.H \ vel=7000vel.H nrefvel=1 xmin=-1.4 xmax=5.98 ymin=0 ymax=0 zmin=0 zmax=2 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ rot_vlow: < 6525_mig_brad.H Window3d min1=-1 max1=5 min3=-.4 max3=.41 | Transp plane=12 reshape=2,3 > temp_brad.H ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=temp_brad.H Cig_tilt=6601_rot_sig_sou.H par=Cig-rot800.P positive_rot=1 ~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=temp_brad.H Cig_tilt=6602_rot_sig_rec.H par=Cig-rot800.P positive_rot=0 7100areal_synth_sou.H: < 6601_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000velow.H cigset=7101cigset.H \ nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 cigset=cigset.H encode=0 > $@ 7100areal_synth_rec.H: < 6601_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000velow.H \ nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 cigset=cigset.H GGG=GGG_rec.H > $@ 7150arealmig_vel.H: < 7100areal_synth_rec.H Transp plane=13 > 7105areal_synth_rec_transp.H < 7100areal_synth_sou.H Transp plane=13 > 7105areal_synth_sou_transp.H < 7105areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=7105areal_synth_sou_transp.H \ vel=7000vel.H nrefvel=1 xmin=-1.4 xmax=5.98 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 7150arealmig_velow.H: # < 7100areal_synth_rec.H Transp plane=13 > 7105areal_synth_rec_transp.H # < 7100areal_synth_sou.H Transp plane=13 > 7105areal_synth_sou_transp.H < 7105areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=7105areal_synth_sou_transp.H \ vel=7000velow.H nrefvel=1 xmin=-1.4 xmax=5.98 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 71055areal_synth_sou.H: < 6601_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000velow.H \ ngap=82 nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 cigset=cigset.H encode=0 > $@ 71055areal_synth_rec.H: < 6601_rot_sig_sou.H Transp plane=23 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000velow.H \ ngap=82 nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 cigset=cigset.H GGG= GGG_rec.H > $@ 7155arealmig_velow.H: < 71055areal_synth_rec.H Transp plane=13 > 71055areal_synth_rec_transp.H < 71055areal_synth_sou.H Transp plane=13 > 71055areal_synth_sou_transp.H < 71055areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=71055areal_synth_sou_transp.H \ vel=7000velow.H nrefvel=1 xmin=-1.4 xmax=7.58 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@ 7300areal_synth_sou.H: < 6530_mig_brad.H Window3d min1=-1 max1=5 | Transp plane=13 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000vel.H \ cigset= 7101cigset.H nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=1 encode=0 > $@ 7300areal_synth_rec.H: < 6530_mig_brad.H Window3d min1=-1 max1=5 | Transp plane=13 | ./bin/OW_arealsht.x sou=6504data_recip_source.H vel=7000vel.H \ cigset= 7101cigset.H nrefvel=1 xmin=-1 xmax=4.98 ymin=0 ymax=0 vsal=6 traces_in_taper=0 souwfld=0 encode=0 > $@ 7350arealmig_vel.H: < 7300areal_synth_rec.H Transp plane=13 > 7301areal_synth_rec_transp.H < 7300areal_synth_sou.H Transp plane=13 > 7301areal_synth_sou_transp.H < 7301areal_synth_rec_transp.H ./bin/OWmig2d_areal.x sou=7301areal_synth_sou_transp.H \ vel=7000vel.H nrefvel=1 xmin=-1.4 xmax=5.98 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -.4 dh=0.02 wbott=0 vsal=5 \ traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@