######################################## # # Sigsbee # 2000areal: # correct velocity areal shot random encode in z < ../2Wpar/04sig_mig_brad.H Window3d max1=30000 min1=15000 max5=16000| Transp plane=13 > 2000temp.H ./bin/OW_arealsht.x rec_cig=2000temp.H sou_cig=2000temp.H sou=../2Wpar/02sig_sou1.H \ vel=../2Wpar/1000velmig_my_areal.H rec_out=2000sig_rec_comb1_enc0_vel.H sou_out=2000sig_sou_comb1_enc0_vel.H \ ngap=15 nrefvel=4 vsal=14500 padtraces=100 wbott=6000 encode=0 funct=0 eps=0.1 encw=0 encz=1 combine=1> /dev/null Rm 2000temp.H 2000sig_mig_comb1_enc01_vel.H: < 2000sig_rec_comb1_enc0_vel.H Transp plane=13 > 2000rec_temp.H < 2000sig_sou_comb1_enc0_vel.H Transp plane=13 > 2000sou_temp.H < 2000rec_temp.H ./bin/OWmig2d_areal.x sou=2000sou_temp.H \ vel=../2Wpar/1000velmig_my_areal.H nrefvel=4 xmin=600 xmax=4400 ymin=0 ymax=0 zmin=0 zmax=16000. nh=33 oh=-1200 dh=75 wbott=0 vsal=14500 \ traces_in_taper=0 every=1 encode=1 deltat=0.04 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 2000*temp.H 2000sig_mig_comb1_enc00_vel.H: < 2000sig_rec_comb1_enc0_vel.H Transp plane=13 > 2000rec_temp.H < 2000sig_sou_comb1_enc0_vel.H Transp plane=13 > 2000sou_temp.H < 2000rec_temp.H ./bin/OWmig2d_areal.x sou=2000sou_temp.H \ vel=../2Wpar/1000velmig_my_areal.H nrefvel=4 xmin=600 xmax=4400 ymin=0 ymax=0 zmin=0 zmax=16000. nh=33 oh=-1200 dh=75 wbott=0 vsal=14500 \ traces_in_taper=0 every=1 encode=0 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 2000*temp.H 2001areal%: # correct velocity areal shot random encode in z # < ../2Wpar/04sig_mig_brad.H Window3d min1=15000 max1=30000 max5=16000| Transp plane=13 > 2001temp.H ./bin/OW_arealsht.x rec_cig=2001temp.H sou_cig=2001temp.H sou=../2Wpar/02sig_sou1.H \ vel=../2Wpar/1000velmig_my_areal.H rec_out=2001sig_rec_comb1_enc1z$*_vel.H sou_out=2001sig_sou_comb1_enc1z$*_vel.H \ ngap=15 nrefvel=4 vsal=14500 padtraces=100 wbott=6000 encode=1 funct=0 eps=$* encw=0 encz=1 combine=1> /dev/null Rm 3001temp.H 2001sig_mig_comb1_enc1z_vel%.H: < 2001sig_rec_comb1_enc1z$*_vel.H Transp plane=13 > 2001rec_temp$*.H < 2001sig_sou_comb1_enc1z$*_vel.H Transp plane=13 > 2001sou_temp$*.H < 2001rec_temp$*.H ./bin/OWmig2d_areal.x sou=2001sou_temp$*.H \ vel=../2Wpar/1000velmig_my_areal.H nrefvel=4 xmin=600 xmax=4400 ymin=0 ymax=0 zmin=0 zmax=16000. nh=33 oh=-1200 dh=75 wbott=0 vsal=14500 \ traces_in_taper=0 every=1 encode=1 wbott=6000 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 2001*temp.H 2002areal: # correct velocity areal shot random encode in z < ../2Wpar/04sig_mig_brad.H Window3d min1=15000 max1=30000 max5=16000| Transp plane=13 > 2002temp.H ./bin/OW_arealsht.x rec_cig=2002temp.H sou_cig=2002temp.H sou=../2Wpar/02sig_sou1.H \ vel=../2Wpar/1000velmig_my_areal.H rec_out=2002sig_rec_comb1_enc1w_vel.H sou_out=2002sig_sou_comb1_enc1w_vel.H \ ngap=15 nrefvel=4 vsal=14500 wbott=6000 padtraces=100 encode=1 funct=0 eps=1. encw=1 encz=0 combine=1> /dev/null Rm 2002temp.H 2002sig_mig_comb1_enc1w_vel.H: < 2002sig_rec_comb1_enc1w_vel.H Transp plane=13 > 2002rec_temp$*.H < 2002sig_sou_comb1_enc1w_vel.H Transp plane=13 > 2002sou_temp$*.H < 2002rec_temp$*.H ./bin/OWmig2d_areal.x sou=2002sou_temp$*.H \ vel=../2Wpar/1000velmig_my_areal.H nrefvel=4 xmin=600 xmax=4400 ymin=0 ymax=0 zmin=0 zmax=16000. nh=33 oh=-1200 dh=75 wbott=0 vsal=14500 \ traces_in_taper=0 every=1 encode=1 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 2002*temp.H ##################################### # # Marmousi # 3001areal%: # correct velocity areal shot random encode in z < ../2Wpar/80042marm_mig_brad_fdmod.H Window3d min1=2000 max1=6000 max5= 1000| Transp plane=13 > 3001temp.H ./bin/OW_arealsht.x rec_cig=3001temp.H sou_cig=3001temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3001marm_rec_comb1_enc1z_vel$*.H sou_out=3001marm_sou_comb1_enc1z_vel$*.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=0 encode=1 funct=0 eps=0.1 encw=0 encz=1 combine=1> /dev/null Rm 3001temp.H 3001marm_mig_comb1_enc1z_vel%.H: < 3001marm_rec_comb1_enc1z_vel$*.H Transp plane=13 > 3001rec_temp$*.H < 3001marm_sou_comb1_enc1z_vel$*.H Transp plane=13 > 3001sou_temp$*.H < 3001rec_temp$*.H ./bin/OWmig2d_areal.x sou=3001sou_temp$*.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=600 xmax=4400 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=1 encode=1 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3001*temp.H 3002areal: # correct velocity areal shot no encode < ../2Wpar/80042marm_mig_brad_fdmod.H Window3d min1=2000 max1=6000 max5= 1000| Transp plane=13 > 3001temp.H ./bin/OW_arealsht.x rec_cig=3001temp.H sou_cig=3001temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3002marm_rec_comb1_enc0_vel.H sou_out=3002marm_sou_comb1_enc0_vel.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=0 encode=0 funct=0 eps=0.1 encw=0 encz=1 combine=1> /dev/null Rm 3001temp.H 3002marm_mig_comb1_enc0_vel.H: < 3002marm_rec_comb1_enc0_vel.H Transp plane=13 > 3002rec_temp.H < 3002marm_sou_comb1_enc0_vel.H Transp plane=13 > 3002sou_temp.H < 3002rec_temp.H ./bin/OWmig2d_areal.x sou=3002sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=1624 xmax=6376 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=1 encode=1 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3001*temp.H 3003marm_mig_comb1_enc0_vel.H: < 3002marm_rec_comb1_enc0_vel.H Transp plane=13 > 3002rec_temp.H < 3002marm_sou_comb1_enc0_vel.H Transp plane=13 > 3002sou_temp.H < 3002rec_temp.H ./bin/OWmig2d_areal.x sou=3002sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=1624 xmax=6376 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=1 encode=0 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3001*temp.H 3004marm_mig_comb1_enc0_vel.H: < 3002marm_rec_comb1_enc0_vel.H Transp plane=13 > 3002rec_temp.H < 3002marm_sou_comb1_enc0_vel.H Transp plane=13 > 3002sou_temp.H < 3002rec_temp.H ./bin/OWmig2d_areal.x sou=3002sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=1624 xmax=6376 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=1 encode=0 deltat=0.04 dt=0.004 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3001*temp.H 3001marm_mig_comb1_enc1z_vel.H: # correct velocity areal shot random encode in z < ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3001temp.H ./bin/OW_arealsht.x rec_cig=3001temp.H sou_cig=3001temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3001marm_rec_comb1_enc1z1_vel.H sou_out=3001marm_sou_comb1_enc1z1_vel.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=61 encode=1 funct=1 encw=0 encz=1 combine=1> /dev/null Rm 3001temp.H < 3001marm_rec_comb1_enc1z1_vel.H Transp plane=13 > 3001rec_temp.H < 3001marm_sou_comb1_enc1z1_vel.H Transp plane=13 > 3001sou_temp.H < 3001rec_temp.H ./bin/OWmig2d_areal.x sou=3001sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-864 xmax=9816 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3001*temp.H 3002marm_mig_comb1_enc1z_vel.H: # correct velocity areal shot random encode in z < ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3002temp.H ./bin/OW_arealsht.x rec_cig=3002temp.H sou_cig=3002temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3002marm_rec_comb1_enc1z2_vel.H sou_out=3002marm_sou_comb1_enc1z2_vel.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=61 encode=1 funct=2 encw=0 encz=1 combine=1> /dev/null Rm 3002temp.H < 3002marm_rec_comb1_enc1z1_vel.H Transp plane=13 > 3002rec_temp.H < 3002marm_sou_comb1_enc1z1_vel.H Transp plane=13 > 3002sou_temp.H < 3002rec_temp.H ./bin/OWmig2d_areal.x sou=3002sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-864 xmax=9816 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3002*temp.H 3003marm_mig_comb1_enc1z_vel.H: # correct velocity areal shot random encode in z < ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3003temp.H ./bin/OW_arealsht.x rec_cig=3003temp.H sou_cig=3003temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3003marm_rec_comb1_enc1z3_vel.H sou_out=3003marm_sou_comb1_enc1z3_vel.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=61 encode=1 funct=3 encw=0 encz=1 combine=1> /dev/null Rm 3003temp.H < 3003marm_rec_comb1_enc1z1_vel.H Transp plane=13 > 3003rec_temp.H < 3003marm_sou_comb1_enc1z1_vel.H Transp plane=13 > 3003sou_temp.H < 3003rec_temp.H ./bin/OWmig2d_areal.x sou=3003sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-864 xmax=9816 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3002*temp.H 3004marm_mig_comb1_enc1z_vel.H: # correct velocity areal shot random encode in z < ../2Wpar/8004marm_mig_brad_fdmod.H Window3d min1=1000 max1=8000 | Transp plane=13 > 3004temp.H ./bin/OW_arealsht.x rec_cig=3004temp.H sou_cig=3004temp.H sou=../2Wpar/8102marm_source_brad.H \ vel=../2Wpar/8000marmvel_24m.H rec_out=3004marm_rec_comb1_enc1z4_vel.H sou_out=3004marm_sou_comb1_enc1z4_vel.H \ ngap=21 nrefvel=4 vsal=6000 padtraces=61 encode=1 funct=4 encw=0 encz=1 combine=1> /dev/null Rm 3004temp.H < 3004marm_rec_comb1_enc1z1_vel.H Transp plane=13 > 3004rec_temp.H < 3004marm_sou_comb1_enc1z1_vel.H Transp plane=13 > 3004sou_temp.H < 3004rec_temp.H ./bin/OWmig2d_areal.x sou=3004sou_temp.H \ vel=../2Wpar/8000marmvel_24m.H nrefvel=4 xmin=-864 xmax=9816 ymin=0 ymax=0 zmin=0 zmax=3000. nh=33 oh=-384 dh=24 wbott=0 vsal=6000 \ traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3004*temp.H 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 > $@ ################################################## # # Plane reflector 2km/s - Modeling and preparing # 4000data_hor.H: susynlv ref="1:0,.5.;10.,.5" nt=1501 dt=0.002 dxs=0.01 nxs=401 fpeak=20 fxs=0 dxo=0.01 nxo=301 fxo=-3 v00=2 > junk1.su < junk1.su Su2sep > junk1.H susynlv ref="1:0,0.;10.,2.5" nt=1501 dt=0.002 dxs=0.01 nxs=401 fpeak=20 fxs=0 dxo=0.01 nxo=301 fxo=-3 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=301 o2=-3 d2=0.01 n3=401 o3=0 d3=0.01 >> $@ 4000data_prep.H: < 4000data_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 > $@ 4000data_source.H: Window3d < 4000data_hor.H n3=1 f3=380 max1=.75 n2=1 min2=-.4| Pad n1out=376 > 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 > 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 > $@ 4000_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=4000data_prep.H D_d=4000data_source.H V=4000vel_brad_pad10.H R=$@ \ pverb=2 nx_shots=401 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-4.5 maxx=5.5 & #restart=y 4000_mig_brad_velow.H: # With 401 shots separated 0.01km correct velocity @date nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \ D=4000data_prep.H D_d=4000data_source.H V=4000vtemp_0.01.H R=$@ \ pverb=2 nx_shots=401 ox_shots=0 \ par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-4.5 maxx=5.5 & #restart=y 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 >> $@ 4001data_hor2.H: susynlv ref="1:0,.5;10.,.5" nt=1501 dt=0.002 dxs=0.02 nxs=401 fpeak=20 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su < junk1.su Su2sep > junk1.H susynlv ref="1:0,1.;10.,1." nt=1501 dt=0.002 dxs=0.02 nxs=401 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=2 > $@; echo n2=201 o2=-2 d2=0.02 n3=401 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 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 > $@ 4003data_recip_prep2.H: < 4001data_hor2.H Window3d min3=2 max3=6 | /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_prep2.H: < 4001data_hor.H Pad beg1=0 end1=1502 > 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 > 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 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 > 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 > $@ 4004data_recip_source2.H: Window3d < 4001data_hor.H n3=1 f3=200 | Mute vmute=2 tmute=0.7 | Window3d n2=1 min2=-1 | Pad beg1=0 end1=1502 >x.H y.H < y.H Window3d | Cabs > 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=200 ox_shots=0 \ par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y 4008_mig_brad2.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_prep2.H D_d=4004data_recip_source2.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 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 ######################## # # On 04/04/08, back to # exploding reflector # propagation # ######################## # # Show crosstalk problem ngap=41 # 3890areal: < 4008_mig_brad.H Window3d min1=1.5 n1=1 max5=1 | Transp > junk1.H ./bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk1.H sou=4004data_recip_source.H vel=4000vel.H \ rec_out=3890areal_rec_comb1_enc0_vel.H sou_out=3890areal_sou_comb1_enc0_vel.H ngap=1 nrefvel=1 vsal=6 \ padtraces=81 encode=0 combine=1 > /dev/null 3890mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 < 3890areal_rec_comb1_enc0_vel.H Transp plane=13 > 3890areal_rec_temp.H < 3890areal_sou_comb1_enc0_vel.H Transp plane=13 > 3890areal_sou_temp.H < 3890areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3890areal_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= xxx_GS.H wr= xxx_GR.H > $@ 3891mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 < 3890areal_rec_comb1_enc0_vel.H Transp plane=13 > 3890areal_rec_temp.H < 3890areal_sou_comb1_enc0_vel.H Transp plane=13 > 3890areal_sou_temp.H < 3890areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3890areal_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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ 3892areal: < 4008_mig_brad.H Window3d min1=1.5 n1=1 max5=1 | Transp > junk1.H ./bin/OW_arealsht.x rec_cig=junk1.H sou_cig=junk1.H sou=4004data_recip_source.H vel=4000vel.H \ rec_out=3890areal_rec_comb1_enc0_vel.H sou_out=3890areal_sou_comb1_enc0_vel.H ngap=1 nrefvel=1 vsal=6 \ padtraces=81 encode=0 combine=1 > /dev/null 3892mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 < 3890areal_rec_comb1_enc0_vel.H Transp plane=13 > 3890areal_rec_temp.H < 3890areal_sou_comb1_enc0_vel.H Transp plane=13 > 3890areal_sou_temp.H < 3890areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3890areal_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= xxx_GS.H wr= xxx_GR.H > $@ 3990areal: < 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=3990areal_rec_comb1_enc0_vel.H sou_out=3990areal_sou_comb1_enc0_vel.H ngap=21 nrefvel=1 vsal=6 \ padtraces=0 encode=0 combine=1 > /dev/null 3990mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 < 3990areal_rec_comb1_enc0_vel.H Transp plane=13 > 3990areal_rec_temp.H < 3990areal_sou_comb1_enc0_vel.H Transp plane=13 > 3990areal_sou_temp.H < 3990areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3990areal_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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ # Rm 3990areal*temp.H 3990zang_comb1_enc0_vel.H: Window3d < 3990mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 3991mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15 < 3990areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3990areal_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= xxx_GS.H wr= xxx_GR.H > $@ Rm 3990areal*temp.H 3992ang_comb1_enc0_vel.H: Window3d < 3991mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 4100areal: # correct velocity areal shots # < 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=4100areal_rec_comb1_enc0_vel.H sou_out=4100areal_sou_comb1_enc0_vel.H ngap=51 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=1 > /dev/null 4100mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 5:39.28 < 4100areal_rec_comb1_enc0_vel.H Transp plane=13 > 4100areal_rec_temp.H < 4100areal_sou_comb1_enc0_vel.H Transp plane=13 > 4100areal_sou_temp.H < 4100areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4100areal_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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4100areal*temp.H 4100zang_comb1_enc0_vel.H: Window3d < 4100mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 4101mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took 22:26.15 < 4100areal_rec_comb1_enc0_vel.H Transp plane=13 > 4100areal_rec_temp.H < 4100areal_sou_comb1_enc0_vel.H Transp plane=13 > 4100areal_sou_temp.H < 4100areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4100areal_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= xxx_GS.H wr= xxx_GR.H > $@ Rm 4100areal*temp.H 4102ang_comb1_enc0_vel.H: Window3d < 4101mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 4105areal: # correct velocity one synth shot separately ./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_rec_comb0_enc0_vel.H sou_out=4105areal_sou_comb0_enc0_vel.H ngap=41 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=0 > /dev/null 4105mig_comb0_enc0_vel.H: # one synth shot separately migration using window (encode=0) < 4105areal_rec_comb0_enc0_vel.H Transp plane=13 > 4105areal_rec_temp.H < 4105areal_sou_comb0_enc0_vel.H Transp plane=13 > 4105areal_sou_temp.H < 4105areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4105areal_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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4105areal_*_temp.H 4106ang_comb0_enc0_vel.H: Window3d < 4105mig_comb0_enc0_vel.H min3=0 max3=4 | 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; $@ # # Finished showing crosstalk problem ngap=41 # ############################################################ # # Show crosstalk problem ngap=11 # Make the adcigs # # 4200areal: # correct velocity areal shots < 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=4200areal_rec_comb1_enc0_vel.H sou_out=4200areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=1 > /dev/null ##### # # Stereographic i.cond. 4200mig_comb1_enc0_vel_stereo.H: # correct velocity areal shots ### It took 1:31:12 < 4200areal_rec_comb1_enc0_vel.H Transp plane=13 > 4200areal_rec_temp.H < 4200areal_sou_comb1_enc0_vel.H Transp plane=13 > 4200areal_sou_temp.H < 4200areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4200areal_sou_temp.H \ vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 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 deltat=0.02 dt=0.008 encode=0 stereo=1 geodip=30. azn=31 axn=31 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4200areal*temp.H 4200mig_comb1_enc0_vel.H: # correct velocity areal shots ### It took 1:31:12 < 4200areal_rec_comb1_enc0_vel.H Transp plane=13 > 4200areal_rec_temp.H < 4200areal_sou_comb1_enc0_vel.H Transp plane=13 > 4200areal_sou_temp.H < 4200areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4200areal_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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4200areal*temp.H 4200zang_comb1_enc0_vel.H: Window3d < 4200mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 4201mig_comb1_enc0_vel.H: # areal shots with encode=1 to let migration generate crosstalk between reflectors ### It took < 4200areal_rec_comb1_enc0_vel.H Transp plane=13 > 4200areal_rec_temp1.H < 4200areal_sou_comb1_enc0_vel.H Transp plane=13 > 4200areal_sou_temp1.H < 4200areal_rec_temp1.H ./bin/OWmig2d_areal.x sou=4200areal_sou_temp1.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= xxx_GS.H wr= xxx_GR.H > $@ Rm 4200areal*temp1.H 4202ang_comb1_enc0_vel.H: Window3d < 4201mig_comb1_enc0_vel.H min3=0 max3=4 | 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; $@ 4205areal: # correct velocity one synth shot separately ./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=4205areal_rec_comb0_enc0_vel.H sou_out=4205areal_sou_comb0_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=0 > /dev/null 4205mig_comb0_enc0_vel.H: # one synth shot separately migration using window (encode=0) < 4205areal_rec_comb0_enc0_vel.H Transp plane=13 > 4205areal_rec_temp2.H < 4205areal_sou_comb0_enc0_vel.H Transp plane=13 > 4205areal_sou_temp2.H < 4205areal_rec_temp2.H ./bin/OWmig2d_areal.x sou=4205areal_sou_temp2.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=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4205areal_*_temp2.H 4206ang_comb0_enc0_vel.H: Window3d < 4205mig_comb0_enc0_vel.H min3=0 max3=4 | 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; $@ 4190areal: # correct velocity areal shots < 4008_mig_brad2.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > 4008_mig_brad2_transp.H ./bin/OW_arealsht.x rec_cig=4008_mig_brad2_transp.H sou_cig=4008_mig_brad2_transp.H sou=4004data_recip_source2.H vel=4000vel.H \ rec_out=4190areal_rec_comb1_enc0_vel.H sou_out=4190areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=1 > /dev/null 4190mig_comb1_enc0_vel.H: # correct velocity areal shots ### encode 1 to force fre IC < 4190areal_rec_comb1_enc0_vel.H Transp plane=13 > 4190areal_rec_temp.H < 4190areal_sou_comb1_enc0_vel.H Transp plane=13 > 4190areal_sou_temp.H < 4190areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4190areal_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= xxx_GS.H wr= xxx_GR.H > $@ Rm 4190areal*temp.H # # Finished showing crosstalk problem ngap=11 # ############################################################ # # Encoding from the model space for ngap=11 # # Make this: # 43005mig_sum_comb1_encz1_vel_ang.H: # 4302mig_sum_comb1_encw1_vel_ang.H: # # 4333areal: # correct velocity areal shot linear phase encode < 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_source2.H \ vel=4000vel.H rec_out=4300areal_rec_comb1_enclin_vel.H sou_out=4300areal_sou_comb1_enclin_vel.H ngap=41 nrefvel=1 vsal=6 \ padtraces=41 encode=1 t0=25. enclin=1 encw=0 encz=0 combine=1> /dev/null 4333mig_comb1_enclin_vel.H: # correct velocity areal shot random encode in z < 4300areal_rec_comb1_enclin_vel.H Transp plane=13 > 4300areal_rec_temp.H < 4300areal_sou_comb1_enclin_vel.H Transp plane=13 > 4300areal_sou_temp.H < 4300areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4300areal_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=1 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4300*temp.H 4334mig_comb1_enclin_vel.H: # correct velocity areal shot random encode in z < 4300areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp.H \ vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 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 encode=0 deltat=0.02 dt=0.008 ws=xxx_GS.H wr= xxx_GR.H > $@ # Rm 4300*temp.H 4300areal%: # correct velocity areal shot random encode in z # < 4008_mig_brad.H Window3d min1=1 max1=3 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=4300areal_rec_comb1_enc1z_vel_rnd$*.H sou_out=4300areal_sou_comb1_enc1z_vel_rnd$*.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=$* combine=1> /dev/null 4300mig_comb1_enc1z_vel_rnd%.H: # correct velocity areal shot random encode in z < 4300areal_rec_comb1_enc1z_vel_rnd$*.H Transp plane=13 > 4300areal_rec_temp_rnd$*.H < 4300areal_sou_comb1_enc1z_vel_rnd$*.H Transp plane=13 > 4300areal_sou_temp_rnd$*.H < 4300areal_rec_temp_rnd$*.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp_rnd$*.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 > $@ Rm 4300*temp_rnd$*.H 4300mig_comb1_enc1z_vel_rndall.H: Math file1=4300mig_comb1_enc1z_vel_rnd127.H file2=4300mig_comb1_enc1z_vel_rnd331.H file3=4300mig_comb1_enc1z_vel_rnd449.H \ file4=4300mig_comb1_enc1z_vel_rnd237.H file5=4300mig_comb1_enc1z_vel_rnd551.H file6=4300mig_comb1_enc1z_vel_rnd657.H file7=4300mig_comb1_enc1z_vel_rnd779.H \ file8=4300mig_comb1_enc1z_vel_rnd857.H file9=4300mig_comb1_enc1z_vel_rnd991.H file10=4300mig_comb1_enc1z_vel_rnd000.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4300ang_comb1_enc1z_vel_rndall.H: Window3d < 4300mig_comb1_enc1z_vel_rndall.H min3=1 max3=3 | 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; $@ 4301areal%: # 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=4301areal_rec_comb1_enc1w_vel_$*.H sou_out=4301areal_sou_comb1_enc1w_vel_$*.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4301mig_comb1_enc1w_vel_%.H: # correct velocity areal shot random encode in w < 4301areal_rec_comb1_enc1w_vel_$*.H Transp plane=13 > 4301areal_rec_temp_$*.H < 4301areal_sou_comb1_enc1w_vel_$*.H Transp plane=13 > 4301areal_sou_temp_$*.H < 4301areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4301areal_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4301*temp_$*.H 4301mig_comb1_enc1w_vel_rndall.H: Math file1=4301mig_comb1_enc1w_vel_01.H file2=4301mig_comb1_enc1w_vel_03.H file3=4301mig_comb1_enc1w_vel_04.H \ file4=4301mig_comb1_enc1w_vel_02.H file5=4301mig_comb1_enc1w_vel_05.H file6=4301mig_comb1_enc1w_vel_06.H file7=4301mig_comb1_enc1w_vel_07.H \ file8=4301mig_comb1_enc1w_vel_08.H file9=4301mig_comb1_enc1w_vel_09.H file10=4301mig_comb1_enc1w_vel_10.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4302mig_comb1_enc1w_vel_%.H: # correct velocity areal shot random encode in w < 4301areal_rec_comb1_enc1w_vel_$*.H Transp plane=13 > 4302areal_rec_temp_$*.H < 4301areal_sou_comb1_enc1w_vel_$*.H Transp plane=13 > 4302areal_sou_temp_$*.H < 4302areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4302areal_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 encode=0 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4302*temp_$*.H 4302mig_comb1_enc1w_vel_rndall.H: Math file1=4302mig_comb1_enc1w_vel_01.H file2=4302mig_comb1_enc1w_vel_03.H file3=4302mig_comb1_enc1w_vel_04.H \ file4=4302mig_comb1_enc1w_vel_02.H file5=4302mig_comb1_enc1w_vel_05.H file6=4302mig_comb1_enc1w_vel_06.H file7=4302mig_comb1_enc1w_vel_07.H \ file8=4302mig_comb1_enc1w_vel_08.H file9=4302mig_comb1_enc1w_vel_09.H file10=4302mig_comb1_enc1w_vel_10.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4301ang_comb1_encw1_vel.H: Window3d < 4301mig_comb1_enc1w_vel.H min3=0 max3=4 | 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; $@ 4305areal%: # 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=4305areal_rec_comb1_enc1w_vel_$*.H sou_out=4305areal_sou_comb1_enc1w_vel_$*.H ngap=1 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4305mig_comb1_enc1w_vel_%.H: # correct velocity areal shot random encode in w < 4305areal_rec_comb1_enc1w_vel_$*.H Transp plane=13 > 4305areal_rec_temp_$*.H < 4305areal_sou_comb1_enc1w_vel_$*.H Transp plane=13 > 4305areal_sou_temp_$*.H < 4305areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4305areal_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4305*temp_$*.H 4305mig_comb1_enc1w_vel_rndall.H: Math file1=4305mig_comb1_enc1w_vel_01.H file2=4305mig_comb1_enc1w_vel_03.H file3=4305mig_comb1_enc1w_vel_04.H \ file4=4305mig_comb1_enc1w_vel_02.H file5=4305mig_comb1_enc1w_vel_05.H file6=4305mig_comb1_enc1w_vel_06.H file7=4305mig_comb1_enc1w_vel_07.H \ file8=4305mig_comb1_enc1w_vel_08.H file9=4305mig_comb1_enc1w_vel_09.H file10=4305mig_comb1_enc1w_vel_10.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ # # Encoding in z reduces cross-talks but destroys the signal. # So, I used an epsilon to control the amount of randomness. # The results are better. I will try later linear encoding. # Edit the flows. # ################################### # # Slower velocity # # Make: # 4490areal: #slower velocity Window3d <4006_rot_sig_sou.H max1=1 min2=3 n2=1 > junks.H Window3d <4006_rot_sig_rec.H max1=1 min2=3 n2=1 > junkr.H ./bin/OW_arealsht.x rec_cig=junkr.H sou_cig=junks.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4490areal_rec_comb1_enc0_velow.H sou_out=4490areal_sou_comb1_enc0_velow.H ngap=1 nrefvel=1 vsal=6 \ padtraces=81 encode=0 combine=1 > /dev/null 4490mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4490areal_rec_comb1_enc0_velow.H Transp plane=13 > 4490areal_rec_temp.H < 4490areal_sou_comb1_enc0_velow.H Transp plane=13 > 4490areal_sou_temp.H < 4490areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4490areal_sou_temp.H \ vel=4000velow.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 vsal=5 \ traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4490areal*temp.H 44900mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4490areal_rec_comb1_enc0_velow.H Transp plane=13 > 4490areal_rec_temp.H < 4490areal_sou_comb1_enc0_velow.H Transp plane=13 > 4490areal_sou_temp.H < 4490areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4490areal_sou_temp.H \ vel=4000velow.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 vsal=5 \ traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4490areal*temp.H 4490zang_comb1_enc0_velow.H: Window3d < 4490mig_comb1_enc0_velow.H | 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; $@ 4500areal: #slower velocity Window3d <4006_rot_sig_sou.H max1=1 min2=1 max2=3 |Transp plane=23 > 4500_rot_sig_sou.H Window3d <4006_rot_sig_rec.H max1=1 min2=1 max2=3 |Transp plane=23 > 4500_rot_sig_rec.H ./bin/OW_arealsht.x rec_cig=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4500areal_rec_comb1_enc0_velow.H sou_out=4500areal_sou_comb1_enc0_velow.H ngap=51 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=1 > /dev/null 4500mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4500areal_rec_comb1_enc0_velow.H Transp plane=13 > 4500areal_rec_temp.H < 4500areal_sou_comb1_enc0_velow.H Transp plane=13 > 4500areal_sou_temp.H < 4500areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4500areal_sou_temp.H \ vel=4000velow.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 vsal=5 \ traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4500areal*temp.H 4500zang_comb1_enc0_velow.H: Window3d < 4500mig_comb1_enc0_velow.H | 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; $@ 4501mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4500areal_rec_comb1_enc0_velow.H Transp plane=13 > 4500areal_rec_temp.H < 4500areal_sou_comb1_enc0_velow.H Transp plane=13 > 4500areal_sou_temp.H < 4500areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4500areal_sou_temp.H \ vel=4000velow.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= xxx_GS.H wr= xxx_GR.H > $@ Rm 4500areal*temp.H 4501zang_comb1_enc0_velow.H: Window3d < 4501mig_comb1_enc0_velow.H | 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; $@ 4550areal: #slower velocity # Window3d <4006_rot_sig_sou.H max1=1 min2=1 max2=3 |Transp plane=23 > 4500_rot_sig_sou.H # Window3d <4006_rot_sig_rec.H max1=1 min2=1 max2=3 |Transp plane=23 > 4500_rot_sig_rec.H ./bin/OW_arealsht.x rec_cig=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4550areal_rec_comb1_enc0_velow.H sou_out=4550areal_sou_comb1_enc0_velow.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=0 combine=1 > /dev/null 4550mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4550areal_rec_comb1_enc0_velow.H Transp plane=13 > 4550areal_rec_temp.H < 4550areal_sou_comb1_enc0_velow.H Transp plane=13 > 4550areal_sou_temp.H < 4550areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4550areal_sou_temp.H \ vel=4000velow.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 vsal=5 \ traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4550areal*temp.H 4550zang_comb1_enc0_velow.H: Window3d < 4550mig_comb1_enc0_velow.H | 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; $@ 4551mig_comb1_enc0_velow.H: # slower velocity areal shots ### It took 5:39.28 < 4550areal_rec_comb1_enc0_velow.H Transp plane=13 > 4550areal_rec_temp.H < 4550areal_sou_comb1_enc0_velow.H Transp plane=13 > 4550areal_sou_temp.H < 4550areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4550areal_sou_temp.H \ vel=4000velow.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= xxx_GS.H wr= xxx_GR.H > $@ Rm 4500areal*temp.H 4551zang_comb1_enc0_velow.H: Window3d < 4551mig_comb1_enc0_velow.H | 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; $@ ############################################################ # # Encoding from the model space for ngap=11 # # Make: # # 46005mig_sum_comb1_encz1_velow_ang.H: # 4602ang_comb1_encw1_velow.H: # 4600areal%: # correct velocity areal shot random encode in z # < 4008_mig_brad.H Window3d min1=1 max1=3 max5=1 | Transp plane=13 > 4008_mig_brad_transp.H ./bin/OW_arealsht.x rec_cig=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4600areal_rec_comb1_enc1z_velow_rnd$*.H sou_out=4600areal_sou_comb1_enc1z_velow_rnd$*.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=$* combine=1> /dev/null 4600mig_comb1_enc1z_velow_rnd%.H: # correct velocity areal shot random encode in z < 4600areal_rec_comb1_enc1z_velow_rnd$*.H Transp plane=13 > 4600areal_rec_temp_rnd$*.H < 4600areal_sou_comb1_enc1z_velow_rnd$*.H Transp plane=13 > 4600areal_sou_temp_rnd$*.H < 4600areal_rec_temp_rnd$*.H ./bin/OWmig2d_areal.x sou=4600areal_sou_temp_rnd$*.H \ vel=4000velow.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 > $@ Rm 4600*temp_rnd$*.H 4600mig_comb1_enc1z_velow_rndall.H: Math file1=4600mig_comb1_enc1z_velow_rnd127.H file2=4600mig_comb1_enc1z_velow_rnd331.H file3=4600mig_comb1_enc1z_velow_rnd449.H \ file4=4600mig_comb1_enc1z_velow_rnd237.H file5=4600mig_comb1_enc1z_velow_rnd551.H file6=4600mig_comb1_enc1z_velow_rnd657.H \ file7=4600mig_comb1_enc1z_velow_rnd779.H file8=4600mig_comb1_enc1z_velow_rnd857.H file9=4600mig_comb1_enc1z_velow_rnd991.H \ file10=4600mig_comb1_enc1z_velow_rnd000.H exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4600ang_comb1_enc1z_velow_rndall.H: Window3d < 4600mig_comb1_enc1z_velow_rndall.H min3=1 max3=3 | 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; $@ 4600mig_comb1_enc1z_velow_rndall_15.H: Math file1=4600mig_comb1_enc1z_velow_rnd127.H file2=4600mig_comb1_enc1z_velow_rnd331.H file3=4600mig_comb1_enc1z_velow_rnd449.H \ file4=4600mig_comb1_enc1z_velow_rnd237.H file5=4600mig_comb1_enc1z_velow_rnd551.H file6=4600mig_comb1_enc1z_velow_rnd657.H \ file7=4600mig_comb1_enc1z_velow_rnd779.H file8=4600mig_comb1_enc1z_velow_rnd857.H file9=4600mig_comb1_enc1z_velow_rnd991.H \ file10=4600mig_comb1_enc1z_velow_rnd000.H file11=4600mig_comb1_enc1z_velow_rnd1271.H file12=4600mig_comb1_enc1z_velow_rnd3311.H \ file13=4600mig_comb1_enc1z_velow_rnd4491.H file14=4600mig_comb1_enc1z_velow_rnd2371.H file15=4600mig_comb1_enc1z_velow_rnd5511.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15 > $@ 4600ang_comb1_enc1z_velow_rndall_15.H: Window3d < 4600mig_comb1_enc1z_velow_rndall_15.H min3=1 max3=3 | 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; $@ 4600mig_comb1_enc1z_velow_rndall_20.H: Math file1=4600mig_comb1_enc1z_velow_rnd127.H file2=4600mig_comb1_enc1z_velow_rnd331.H file3=4600mig_comb1_enc1z_velow_rnd449.H \ file4=4600mig_comb1_enc1z_velow_rnd237.H file5=4600mig_comb1_enc1z_velow_rnd551.H file6=4600mig_comb1_enc1z_velow_rnd657.H \ file7=4600mig_comb1_enc1z_velow_rnd779.H file8=4600mig_comb1_enc1z_velow_rnd857.H file9=4600mig_comb1_enc1z_velow_rnd991.H \ file10=4600mig_comb1_enc1z_velow_rnd000.H file11=4600mig_comb1_enc1z_velow_rnd1271.H file12=4600mig_comb1_enc1z_velow_rnd3311.H \ file13=4600mig_comb1_enc1z_velow_rnd4491.H file14=4600mig_comb1_enc1z_velow_rnd2371.H file15=4600mig_comb1_enc1z_velow_rnd5511.H \ file16=4600mig_comb1_enc1z_velow_rnd6571.H file17=4600mig_comb1_enc1z_velow_rnd7791.H file18=4600mig_comb1_enc1z_velow_rnd8571.H \ file19=4600mig_comb1_enc1z_velow_rnd9911.H file20=4600mig_comb1_enc1z_velow_rnd0001.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15+file16+file17+file18+file19+file20 > $@ 4600ang_comb1_enc1z_velow_rndall_20.H: Window3d < 4600mig_comb1_enc1z_velow_rndall_20.H min3=1 max3=3 | 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; $@ 4601areal%: # 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=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4601areal_rec_comb1_enc1w_velow_$*.H sou_out=4601areal_sou_comb1_enc1w_velow_$*.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4601mig_comb1_enc1w_velow_%.H: # correct velocity areal shot random encode in w < 4601areal_rec_comb1_enc1w_velow_$*.H Transp plane=13 > 4601areal_rec_temp_$*.H < 4601areal_sou_comb1_enc1w_velow_$*.H Transp plane=13 > 4601areal_sou_temp_$*.H < 4601areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4601areal_sou_temp_$*.H \ vel=4000velow.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 > $@ Rm 4601*temp_$*.H 4601ang_comb1_encw1_velow.H: Window3d < 4601mig_comb1_enc1w_velow.H min3=0 max3=4 | 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; $@ 4602areal%: # 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=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4602areal_rec_comb1_enc1w_velow_$*.H sou_out=4602areal_sou_comb1_enc1w_velow_$*.H ngap=1 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4602mig_comb1_enc1w_velow_%.H: # correct velocity areal shot random encode in w < 4602areal_rec_comb1_enc1w_velow_$*.H Transp plane=13 > 4602areal_rec_temp_$*.H < 4602areal_sou_comb1_enc1w_velow_$*.H Transp plane=13 > 4602areal_sou_temp_$*.H < 4602areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4602areal_sou_temp_$*.H \ vel=4000velow.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 > $@ Rm 4601*temp_$*.H 46020mig_comb1_enc1w_vel_%.H: # correct velocity areal shot random encode in w < 4602areal_rec_comb1_enc1w_velow_$*.H Transp plane=13 > 4602areal_rec_temp_$*.H < 4602areal_sou_comb1_enc1w_velow_$*.H Transp plane=13 > 4602areal_sou_temp_$*.H < 4602areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4602areal_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 encode=1 deltat=0.02 dt=0.008 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4602*temp_$*.H 4603areal%: # 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=4500_rot_sig_rec.H sou_cig=4500_rot_sig_sou.H sou=4004data_recip_source.H vel=4000velow.H \ rec_out=4603areal_rec_comb1_enc1w_velow_$*.H sou_out=4603areal_sou_comb1_enc1w_velow_$*.H ngap=5 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4603mig_comb1_enc1w_velow_%.H: # correct velocity areal shot random encode in w < 4603areal_rec_comb1_enc1w_velow_$*.H Transp plane=13 > 4603areal_rec_temp_$*.H < 4603areal_sou_comb1_enc1w_velow_$*.H Transp plane=13 > 4603areal_sou_temp_$*.H < 4603areal_rec_temp_$*.H ./bin/OWmig2d_areal.x sou=4603areal_sou_temp_$*.H \ vel=4000velow.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 > $@ Rm 4601*temp_$*.H 4603ang_comb1_encw1_velow.H: Window3d < 4603mig_comb1_enc1w_velow.H min3=0 max3=4 | 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; $@ # # # # ################################### # Finished showing crosstalk problem ngap=11 # ############################################################ # # Encoding from the model space for ngap=11 # # Make this: # 43005mig_sum_comb1_encz1_vel_ang.H: # 4302mig_sum_comb1_encw1_vel_ang.H: # # 4700areal%: # correct velocity areal shot random encode in z # < 4008_mig_brad.H Window3d min1=1 max1=3 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=4700areal_rec_comb1_enc1z_vel_rnd$*.H sou_out=4700areal_sou_comb1_enc1z_vel_rnd$*.H ngap=1 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=$* combine=1> /dev/null 4700mig_comb1_enc1z_vel_rnd%.H: # correct velocity areal shot random encode in z < 4700areal_rec_comb1_enc1z_vel_rnd$*.H Transp plane=13 > 4700areal_rec_temp_rnd$*.H < 4700areal_sou_comb1_enc1z_vel_rnd$*.H Transp plane=13 > 4700areal_sou_temp_rnd$*.H < 4700areal_rec_temp_rnd$*.H ./bin/OWmig2d_areal.x sou=4700areal_sou_temp_rnd$*.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 > $@ Rm 4700*temp_rnd$*.H 4700mig_comb1_enc1z_vel_rndall.H: Math file1=4700mig_comb1_enc1z_vel_rnd127.H file2=4700mig_comb1_enc1z_vel_rnd331.H file3=4700mig_comb1_enc1z_vel_rnd449.H \ file4=4700mig_comb1_enc1z_vel_rnd237.H file5=4700mig_comb1_enc1z_vel_rnd551.H file6=4700mig_comb1_enc1z_vel_rnd657.H file7=4700mig_comb1_enc1z_vel_rnd779.H \ file8=4700mig_comb1_enc1z_vel_rnd857.H file9=4700mig_comb1_enc1z_vel_rnd991.H file10=4700mig_comb1_enc1z_vel_rnd000.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4700ang_comb1_enc1z_vel_rndall.H: Window3d < 4700mig_comb1_enc1z_vel_rndall.H min3=1 max3=3 | 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; $@ 47000areal%: # correct velocity areal shot random encode in z # < 4008_mig_brad.H Window3d min1=1 max1=3 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=4700areal_rec_comb1_enc1z_vel_rnd$*_eps25.H sou_out=4700areal_sou_comb1_enc1z_vel_rnd$*_eps25.H ngap=1 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=0 encz=1 eps=0.25 rndcte=$* combine=1> /dev/null 4700mig_comb1_enc1z_vel_rnd%_eps25.H: # correct velocity areal shot random encode in z < 4700areal_rec_comb1_enc1z_vel_rnd$*_eps25.H Transp plane=13 > 4700areal_rec_temp_rnd$*.H < 4700areal_sou_comb1_enc1z_vel_rnd$*_eps25.H Transp plane=13 > 4700areal_sou_temp_rnd$*.H < 4700areal_rec_temp_rnd$*.H ./bin/OWmig2d_areal.x sou=4700areal_sou_temp_rnd$*.H \ vel=4000vel.H nrefvel=1 xmin=0 xmax=4 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 > $@ Rm 4700*temp_rnd$*.H 47001areal%: # correct velocity areal shot random encode in z # < 4008_mig_brad.H Window3d min1=1 max1=3 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=4700areal_rec_comb1_enc1z_vel_rnd$*_eps1.H sou_out=4700areal_sou_comb1_enc1z_vel_rnd$*_eps1.H ngap=1 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=0 encz=1 eps=0.1 rndcte=$* combine=1> /dev/null 4700mig_comb1_enc1z_vel_rnd%_eps1.H: # correct velocity areal shot random encode in z < 4700areal_rec_comb1_enc1z_vel_rnd$*_eps1.H Transp plane=13 > 4700areal_rec_temp_rnd$*.H < 4700areal_sou_comb1_enc1z_vel_rnd$*_eps1.H Transp plane=13 > 4700areal_sou_temp_rnd$*.H < 4700areal_rec_temp_rnd$*.H ./bin/OWmig2d_areal.x sou=4700areal_sou_temp_rnd$*.H \ vel=4000vel.H nrefvel=1 xmin=0 xmax=4 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 > $@ Rm 4700*temp_rnd$*.H 4700mig_comb1_enc1z_vel_rndall00.H: Math file1=4700mig_comb1_enc1z_vel_rnd01.H file2=4700mig_comb1_enc1z_vel_rnd03.H file3=4700mig_comb1_enc1z_vel_rnd04.H \ file4=4700mig_comb1_enc1z_vel_rnd02.H file5=4700mig_comb1_enc1z_vel_rnd05.H file6=4700mig_comb1_enc1z_vel_rnd06.H file7=4700mig_comb1_enc1z_vel_rnd07.H \ file8=4700mig_comb1_enc1z_vel_rnd08.H file9=4700mig_comb1_enc1z_vel_rnd09.H file10=4700mig_comb1_enc1z_vel_rnd10.H \ exp=file1+file2+file3+file4+file5+file6+file7+file8+file9+file10 > $@ 4701areal%: # 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=4301areal_rec_comb1_enc1w_vel_$*.H sou_out=4301areal_sou_comb1_enc1w_vel_$*.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 4701mig_comb1_enc1w_vel.H: # correct velocity areal shot random encode in w < 4301areal_rec_comb1_enc1w_vel.H Transp plane=13 > 4301areal_rec_temp.H < 4301areal_sou_comb1_enc1w_vel.H Transp plane=13 > 4301areal_sou_temp.H < 4301areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4301areal_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4301*temp.H 4701ang_comb1_encw1_vel.H: Window3d < 4301mig_comb1_enc1w_vel.H min3=0 max3=4 | 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; $@ # # Encoding in z reduces cross-talks but destroys the signal. # So, I used an epsilon to control the amount of randomness. # The results are better. I will try later linear encoding. # Edit the flows. # ################################### xxxareal1: # 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=1r.H sou_out=1s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 1mig.H: # correct velocity areal shot random encode in z < 1r.H Transp plane=13 > 1r_temp.H < 1s.H Transp plane=13 > 1s_temp.H < 1r_temp.H ./bin/OWmig2d_areal.x sou=1s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 1*_temp.H xxxareal2: # 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=2r.H sou_out=2s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 2mig.H: # correct velocity areal shot random encode in z < 2r.H Transp plane=13 > 2r_temp.H < 2s.H Transp plane=13 > 2s_temp.H < 2r_temp.H ./bin/OWmig2d_areal.x sou=2s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 2*_temp.H xxxareal3: # 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=3r.H sou_out=3s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 3mig.H: # correct velocity areal shot random encode in z < 3r.H Transp plane=13 > 3r_temp.H < 3s.H Transp plane=13 > 3s_temp.H < 3r_temp.H ./bin/OWmig2d_areal.x sou=3s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 3*_temp.H xxxareal4: # 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=4r.H sou_out=4s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 4mig.H: # correct velocity areal shot random encode in z < 4r.H Transp plane=13 > 4r_temp.H < 4s.H Transp plane=13 > 4s_temp.H < 4r_temp.H ./bin/OWmig2d_areal.x sou=4s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 4*_temp.H xxxareal5: # 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=5r.H sou_out=5s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 5mig.H: # correct velocity areal shot random encode in z < 5r.H Transp plane=13 > 5r_temp.H < 5s.H Transp plane=13 > 5s_temp.H < 5r_temp.H ./bin/OWmig2d_areal.x sou=5s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 5*_temp.H xxxareal6: # 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=6r.H sou_out=6s.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=0 encz=1 eps=0.05 funct=1 combine=1> /dev/null 6mig.H: # correct velocity areal shot random encode in z < 6r.H Transp plane=13 > 6r_temp.H < 6s.H Transp plane=13 > 6s_temp.H < 6r_temp.H ./bin/OWmig2d_areal.x sou=6s_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 encode=1 ws= xxx_GS.H wr= xxx_GR.H > $@ Rm 6*_temp.H wwwareal1: # 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=1rw.H sou_out=1sw.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=1 encz=0 eps=0.05 combine=1> /dev/null 1migw.H: # correct velocity areal shot random encode in z < 1rw.H Transp plane=13 > 1r_tempw.H < 1sw.H Transp plane=13 > 1s_tempw.H < 1r_tempw.H ./bin/OWmig2d_areal.x sou=1s_tempw.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 > $@ Rm 1*_tempw.H wwwareal2: # 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=2rw.H sou_out=2sw.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=1 encz=0 eps=1 0.5 combine=1> /dev/null 2migw.H: # correct velocity areal shot random encode in z < 2rw.H Transp plane=13 > 2r_tempw.H < 2sw.H Transp plane=13 > 2s_tempw.H < 2r_tempw.H ./bin/OWmig2d_areal.x sou=2s_tempw.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 > $@ Rm 2*_tempw.H wwwareal3: # 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=3rw.H sou_out=3sw.H ngap=11 nrefvel=1 vsal=6 \ padtraces=41 encode=1 encodset=4 encw=1 encz=0 eps=1 0.5 combine=1> /dev/null 3migw.H: # correct velocity areal shot random encode in z < 3rw.H Transp plane=13 > 3r_tempw.H < 3sw.H Transp plane=13 > 3s_tempw.H < 3r_tempw.H ./bin/OWmig2d_areal.x sou=3s_tempw.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 > $@ Rm 3*_tempw.H