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
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
#!/bin/bash
#PBS -l nodes=1
#PBS -N InterfModeling
#PBS -q fourweeks
#PBS -V
#
# receivers and source placed on model with topography, 161 hours
echo " !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
echo " This Figure can not be reproduced completely"
echo " The progrom corrsmp used is still in test phase"
echo " Mail j.w.thorbecke@tudelft.nl if you want to reproduce this Figure."
echo " Only the reference response will be modeled"
export PATH=../../bin:$PATH
dt=0.0004
ntap=120
fmax=45
makemod sizex=10000 sizez=4100 dx=5 dz=5 cp0=0 ro0=1000 file_base=real2.su \
orig=0,-800 gradunit=0 \
intt=def poly=2 cp=2450 ro=1000 gradcp=14 grad=0 \
x=0,1000,1700,1800,2000,3000,4000,4500,6000,6800,7000,7500,8100,8800,10000 \
z=-100,-200,-250,-200,-200,-120,-300,-600,-650,-500,-350,-200,-200,-150,-200 \
intt=rough var=200,3.2,1 poly=2 x=0,3000,8000,10000 \
z=400,250,300,500 cp=4500,4200,4800,4500 ro=1400 gradcp=5 grad=0 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=1100,1100,1100,1600,1100,1100,1100 cp=4000 ro=2000 gradcp=8 grad=0 \
intt=def poly=0 x=0,10000 z=1750,2050 cp=4500,5100 ro=1500 gradcp=13 grad=0 \
intt=def poly=0 x=0,10000 z=1850,2150 cp=6000,4200 ro=1500 gradcp=14 grad=0 \
intt=def poly=0 x=0,10000 z=1950,2250 cp=4800,4800 ro=1500 gradcp=5 grad=0 \
intt=def poly=0 x=0,10000 z=2000,2300 cp=6100,5000 ro=1500 gradcp=13 grad=0 \
intt=def poly=0 x=0,10000 z=2100,2400 cp=3800,5000 ro=1500 gradcp=20 grad=0 \
intt=def poly=0 x=0,10000 z=2150,2450 cp=5000 ro=1500 gradcp=14 grad=0 \
intt=def poly=0 x=0,10000 z=2350,2650 cp=5800 ro=1500 gradcp=5 grad=0 \
intt=def poly=0 x=0,10000 z=2600,2600 cp=5500 ro=2200 gradcp=5 grad=0
sushw key=f1 a=0 < real2_cp.su | \
sushw key=f1 a=0 | \
supsimage hbox=6 wbox=8 labelsize=12 f2num=-5000 d2num=1000 \
wrgb=0,0,1.0 grgb=0,1.0,0 brgb=1.0,0,0 legend=1 lnice=1 lstyle=vertright \
bclip=7053.02 wclip=0 label1="depth [m]" label2="lateral position [m]" \
> model2_cp.eps
makewave w=g2 fmax=45 t0=0.10 dt=$dt nt=4096 db=-40 file_out=G2.su verbose=1
extendModel file_in=real2_ro.su nafter=$ntap nbefore=$ntap nabove=0 nbelow=$ntap > vel2_edge_ro.su
extendModel file_in=real2_cp.su nafter=$ntap nbefore=$ntap nabove=0 nbelow=$ntap > vel2_edge_cp.su
#reference
fdelmodc \
file_cp=vel2_edge_cp.su ischeme=1 \
file_den=vel2_edge_ro.su \
file_rcv=shot_real2_x5000_topo.su \
file_src=G2.su \
dtrcv=0.004 \
verbose=4 \
tmod=3.004 \
dxrcv=20.0 \
zrcv1=-800 \
zrcv2=-800 \
xrcv1=0 \
xrcv2=10000 \
sinkdepth=1 \
src_random=0 \
wav_random=0 \
xsrc1=5000 \
xsrc2=5000 \
zsrc1=-800 \
tsrc1=0.0 \
dipsrc=1 \
ntaper=$ntap \
left=4 right=4 top=1 bottom=4
sushw key=f1,delrt a=0.0,0.0 < shot_real2_x5000_topo_rvz.su | \
basop choice=1 shift=-0.1 | \
supsimage clip=2e-12 f1=0 f2=-5000 x1end=3.004 hbox=8 wbox=6 \
label1="time (s)" label2="lateral position (m)" \
labelsize=10 f2num=-5000 d2num=1000 d1num=0.5 > shot_real2_x5000_topo.eps
#clip minclip=0 < SrcRecPositions.su > nep.su
#addmul file_in1=nep.su file_in2=vel2_edge_cp.su a=8000 | \
exit;
tlength=50
nsrc=1500
xsrc1=100
xsrc2=9900
zsrc1=1750
zsrc2=2600
fmax=30
tmod=3600
tsrc1=0.1
tsrc2=3600
file_shot=shot_real2_T${tmod}_S${nsrc}_Dt${tlength}_F${fmax}.su
fdelmodc \
file_cp=vel2_edge_cp.su ischeme=1 \
file_den=vel2_edge_ro.su \
file_rcv=$file_shot \
dtrcv=0.004 \
dt=$dt \
verbose=4 \
rec_ntsam=15000 \
tmod=$tmod \
dxrcv=20.0 \
zrcv1=-800 \
zrcv2=-800 \
xrcv1=0 \
xrcv2=10000 \
sinkdepth=1 \
src_random=1 \
wav_random=1 \
fmax=$fmax \
xsrc1=$xsrc1 \
xsrc2=$xsrc2 \
zsrc1=$zsrc1 \
zsrc2=$zsrc2 \
tsrc1=0.0 \
tsrc2=$tsrc2 \
tlength=$tlength \
nsrc=$nsrc \
dipsrc=0 \
ntaper=$ntap \
left=4 right=4 top=1 bottom=4
suop2 SrcRecPositions.su vel2_edge_cp.su w1=8000 w2=1.0 op=sum | \
sugain nclip=0 | \
supsimage hbox=6 wbox=8 labelsize=10 \
f2=-5600 f2num=-5000 d2num=1000 f1num=-500 d1num=500 \
x2beg=-5000 x2end=5000 legend=1 \
bclip=8200 wclip=0 label1="depth (m)" label2="lateral position (m)" \
> real_sources_cp.eps
supsimage < vel2_edge_cp.su hbox=6 wbox=8 labelsize=10 \
f2=-5600 f2num=-5000 d2num=1000 f1num=-500 d1num=500 \
x2beg=-5000 x2end=5000 legend=1 bps=24 lstyle=vertright \
bclip=6500 label1="depth (m)" label2="lateral position (m)" \
brgb=0.1,0.1,0.1 grgb=0.9,0.9,0.9 wrgb=1.0,1.0,1.0 \
curve=SrcPositions1500.txt npair=1500 curvecolor=black curvedash=-1 curvewidth=3 \
> real_sources_gray_cp.eps
supsimage hbox=6 wbox=8 labelsize=12 < SrcRecPositions.su \
f2=-5600 f2num=50000 d2num=1000 f1num=5000 d1num=500 \
legend=1 lnice=1 \
label1=" " label2=" " bclip=1 wclip=0 \
> rec_pos_real2.eps
export OMP_NUM_THREADS=4
export MKL_SERIAL=yes
export DFTI_NUMBER_OF_USER_THREADS=1
$HOME/src/CorrSMP_FD_Files/corrsmp file_base=shot_real2_T3600_S1500_Dt50_F30_001_rvz.su ntc=2048 fmax=60 \
verbose=3 nstation=501 \
file_out=corrsmp_out.su nc=1 fullcorr=1 dt=0.004 nt=15000
suwind key=fldr min=251 max=251 < corrsmp_out_cc1.su | \
supsimage clip=5e-13 f1=0 f2=-5000 d2=20 x1end=3.004 hbox=8 wbox=6 \
label1="time (s)" label2="lateral position (m)" \
labelsize=10 f2num=-5000 d2num=1000 d1num=0.5 > corr_real2_x5000_topo.eps