Newer
Older
#!/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/bin/marchenko
#zero angle
if [[ "$1" == "plane0" ]];
then
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
exit
fi
#5 degrees dip without (set src_angle=0)
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
138
139
140
141
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
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