Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
#!/bin/bash -x
ROOT=/vardim/home/thorbcke/R1/
ROOT=/vardim/home/thorbcke/src/
export PATH=$HOME/src/OpenSource/bin:$PATH:
export OMP_NUM_THREADS=4
exe=$ROOT/OpenSource/marchenko/marchenkowd
exe=$ROOT/OpenSource/marchenko/marchenko
#exe=$ROOT/OpenSource/bin/marchenko
#5 degrees dip without (set src_angle=0) tshift for Gmin in marchenko program (recompiled)
base=Source5a120_shiftedIter
base=Source5a120old
src_angle=5
#inifile=Source${src_angle}a120_rp.su
inifile=SourceNew${src_angle}a120_rp.su
base=${inifile%_rp.su}
fmute file_shot=$inifile file_out=fdplus${base}.su above=-1 shift=-6 verbose=1 check=0 hw=2
#for testing mute window in manual (B/W figures)
$exe file_shot=shots/refl_rp.su file_tinv=fdplus${base}.su nshots=601 verbose=2 \
tap=3 niter=16 hw=2 shift=7 smooth=3 plane_wave=1 src_angle=5 rotate=1 \
file_green=green${base}.su file_gplus=Gplus${base}.su file_gmin=Gmin${base}.su \
file_f1plus=f1plus${base}.su file_f1min=f1min${base}.su file_f2=f2${base}.su file_iter=iter${base}.su
#exit
#zero angle
src_angle=0
inifile=SourceNew${src_angle}a120_rp.su
base=${inifile%_rp.su}
fmute file_shot=$inifile file_out=fdplus${base}.su above=-1 shift=-6 verbose=1 check=0 hw=2
$exe file_shot=shots/refl_rp.su file_tinv=fdplus${base}.su nshots=601 verbose=2 \
tap=3 niter=16 hw=2 shift=7 smooth=3 plane_wave=0 src_angle=0 rotate=0 \
file_green=green${base}.su file_gplus=Gplus${base}.su file_gmin=Gmin${base}.su \
file_f1plus=f1plus${base}.su file_f1min=f1min${base}.su file_f2=f2${base}.su
for angle in 5 #10
do
for dip in "" -
do
src_angle=$dip$angle
#inifile=Source${src_angle}a120_rp.su
inifile=SourceNew${src_angle}a120_rp.su
base=${inifile%_rp.su}
fmute file_shot=$inifile file_out=fdplus${base}.su above=-1 shift=-6 verbose=1 check=0 hw=2
$exe file_shot=shots/refl_rp.su file_tinv=fdplus${base}.su nshots=601 verbose=2 \
tap=3 niter=16 hw=2 shift=7 plane_wave=1 src_angle=$src_angle \
file_green=green${base}.su file_gplus=Gplus${base}.su file_gmin=Gmin${base}.su \
file_f1plus=f1plus${base}.su file_f1min=f1min${base}.su file_f2=f2${base}.su
done
done
exit
#source at dipping line in model all fired at t=0
inifile=SourceDipa_rp.su
base=${inifile%_rp.su}
fmute file_shot=$inifile file_out=fdplus${base}.su above=-1 shift=-6 verbose=1 check=0 hw=2
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=fdplus${base}.su nshots=601 verbose=2 \
tap=3 niter=16 hw=2 shift=6 plane_wave=1 src_angle=5 \
file_green=green${base}.su file_gplus=Gplus${base}.su file_gmin=Gmin${base}.su \
file_f1plus=f1plus${base}.su file_f1min=f1min${base}.su file_f2=f2${base}.su
exit
fmute file_shot=TCorrect_Source150_rp.su file_out=p0plus.su above=-1 shift=-6 verbose=1 check=0 hw=8
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=p0plus.su nshots=601 verbose=2 \
tap=3 niter=12 hw=8 shift=12 plane_wave=1 \
file_green=pgreenPl.su file_gplus=GplusPl.su file_gmin=TGmin150.su \
file_f1plus=f1plusPl.su file_f1min=f1minPl.su file_f2=f2Pl.su
fmute file_shot=TCorrect_Source180_rp.su file_out=p0plus.su above=-1 shift=-6 verbose=1 check=0 hw=8
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=p0plus.su nshots=601 verbose=2 \
tap=3 niter=12 hw=8 shift=12 plane_wave=1 \
file_green=pgreenPl.su file_gplus=GplusPl.su file_gmin=TGmin180.su \
file_f1plus=f1plusPl.su file_f1min=f1minPl.su file_f2=f2Pl.su
# HCorrect DATASETS: src_angle=10, src_velo=1500:
fmute file_shot=HCorrect_Source120_rp.su file_out=p0plus.su above=-1 shift=-6 verbose=1 check=0 hw=8
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=p0plus.su nshots=601 verbose=2 \
tap=3 niter=12 hw=8 shift=12 plane_wave=1 \
file_green=pgreenPl.su file_gplus=GplusPl.su file_gmin=HGmin120.su \
file_f1plus=f1plusPl.su file_f1min=f1minPl.su file_f2=f2Pl.su
fmute file_shot=HCorrect_Source150_rp.su file_out=p0plus.su above=-1 shift=-6 verbose=1 check=0 hw=8
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=p0plus.su nshots=601 verbose=2 \
tap=3 niter=12 hw=8 shift=12 plane_wave=1 \
file_green=pgreenPl.su file_gplus=GplusPl.su file_gmin=HGmin150.su \
file_f1plus=f1plusPl.su file_f1min=f1minPl.su file_f2=f2Pl.su
fmute file_shot=HCorrect_Source180_rp.su file_out=p0plus.su above=-1 shift=-6 verbose=1 check=0 hw=8
$ROOT/OpenSource/marchenko/marchenko file_shot=shots/refl_rp.su file_tinv=p0plus.su nshots=601 verbose=2 \
tap=3 niter=12 hw=8 shift=12 plane_wave=1 \
file_green=pgreenPl.su file_gplus=GplusPl.su file_gmin=HGmin180.su \
file_f1plus=f1plusPl.su file_f1min=f1minPl.su file_f2=f2Pl.su
#suximage< TGmin120.su x1end=1.4 perc=99& img dTCorrect_Source120_rp.su x1end=1.4 perc=99&
#suximage< TGmin150.su x1end=1.4 perc=99& img dTCorrect_Source150_rp.su x1end=1.4 perc=99&
#suximage< TGmin180.su x1end=1.4 perc=99& img dTCorrect_Source180_rp.su x1end=1.4 perc=99&
#suximage< TGmin120.su x1end=1.4 perc=99& img dHCorrect_Source120_rp.su x1end=1.4 perc=99&
#suximage< TGmin150.su x1end=1.4 perc=99& img dHCorrect_Source150_rp.su x1end=1.4 perc=99&
#suximage< TGmin180.su x1end=1.4 perc=99& img dHCorrect_Source180_rp.su x1end=1.4 perc=99&
supsimage < TGmin120.su title="Gmin: Depth=600, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> TGmin120.eps
supsimage < TGmin150.su title="Gmin: Depth=750, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> TGmin150.eps
supsimage < TGmin180.su title="Gmin: Depth=900, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> TGmin180.eps
supsimage < HGmin120.su title="Gmin: Depth=600, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> HGmin120.eps
supsimage < HGmin150.su title="Gmin: Depth=750, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> HGmin150.eps
supsimage < HGmin180.su title="Gmin: Depth=900, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5> HGmin180.eps
supsimage < dTCorrect_Source120_rp.su title="Input: Depth=600, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
T120.eps
supsimage < dTCorrect_Source150_rp.su title="Input: Depth=750, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
T150.eps
supsimage < dTCorrect_Source180_rp.su title="Input: Depth=900, angle=5" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
T180.eps
supsimage < dHCorrect_Source120_rp.su title="Input: Depth=600, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
H120.eps
supsimage < dHCorrect_Source150_rp.su title="Input: Depth=750, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
H150.eps
supsimage < dHCorrect_Source180_rp.su title="Input: Depth=900, angle=10" label1="t(s)" label2="x(m)" x1beg=0 x1end=1.4 clip=4e5>
H180.eps