Skip to content
Snippets Groups Projects
Commit 64edae23 authored by joeri.brackenhoff's avatar joeri.brackenhoff
Browse files

3D

parent 7d0ae15a
No related branches found
No related tags found
No related merge requests found
Showing
with 3105 additions and 0 deletions
This diff is collapsed.
/* http://en.wikipedia.org/wiki/Multiply-with-carry */
#include<stdlib.h>
#include<limits.h>
/* random number generator which can be used as an alternative for drand48() */
/* http://school.anhb.uwa.edu.au/personalpages/kwessen/shared/Marsaglia03.html*/
static unsigned long Q[4096],c=362436; /* choose random initial c<809430660 and */
/* 4096 random 32-bit integers for Q[] */
void seedCMWC4096(void)
{
int i;
for (i=0; i<4096; i++) {
Q[i] = lrand48();
}
return;
}
unsigned long CMWC4096(void)
{
unsigned long long t, a=18782LL;
static unsigned long i=4095;
unsigned long x,r=0xfffffffe;
i=(i+1)&4095;
t=a*Q[i]+c;
c=(t>>32);
x=t+c;
if(x<c){x++;c++;}
return(Q[i]=r-x);
}
/* replace defaults with five random seed values in calling program */
static unsigned long x=123456789,y=362436069,z=521288629,w=88675123,v=886756453;
unsigned long xorshift(void)
{unsigned long t;
t=(x^(x>>7)); x=y; y=z; z=w; w=v;
v=(v^(v<<6))^(t^(t<<13)); return (y+y+1)*v;}
double dcmwc4096(void)
{
double rd;
// rd = ((double)xorshift())/((double)ULONG_MAX);
rd = ((double)CMWC4096())/((double)ULONG_MAX);
return rd;
}
#!/bin/bash
#PBS -l nodes=1
#PBS -N InterfModeling
#PBS -q fourweeks
#PBS -V
#
#reference and 2 SI results for visco-acoustic media, 2x200 s. + 2x1 hours
export PATH=../../bin:$PATH
makewave w=g1 fmax=30 t0=0.10 dt=0.0008 nt=4096 db=-40 file_out=G1.su verbose=1
makemod sizex=10000 sizez=4100 dx=10 dz=10 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=400,400 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 z=1100,1100,1100,1600,1100,1100,1100 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=2100,2100 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=2600,2600 cp=5500 ro=2200
xsrc1=1000
xsrc2=9000
zsrc1=3300
zsrc2=3900
dxsrc=10
base=fw5000_Q
../fdelmodc \
file_cp=simple_cp.su ischeme=2 \
Qp=15 \
file_den=simple_ro.su \
file_rcv=${base}.su \
file_src=G1.su \
rec_delay=0.1 \
dtrcv=0.008 \
verbose=3 \
tmod=4.108 \
dxrcv=50.0 \
src_random=0 \
wav_random=0 \
fmax=30 \
xsrc=5000 \
zsrc=0.0 \
dipsrc=0 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
makemod sizex=10000 sizez=5000 dx=10 dz=10 cp0=1500 ro0=1000 file_base=hom.su
../fdelmodc \
file_cp=hom_cp.su ischeme=2 \
Qp=15 \
file_den=hom_ro.su \
file_rcv=${base}_D.su \
file_src=G1.su \
rec_delay=0.1 \
dtrcv=0.008 \
verbose=3 \
tmod=4.108 \
dxrcv=50.0 \
src_random=0 \
wav_random=0 \
fmax=30 \
xsrc=5000 \
zsrc=0.0 \
dipsrc=0 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
suop2 ${base}_rvz.su ${base}_D_rvz.su > Reffw5000_Q_rvz.su
suwind s=1 j=1 tmax=4 f1=0.0 < Reffw5000_Q_rvz.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > shotRefQ_5000_0.eps
tmod=120
tsrc1=0.1
tsrc2=120
tlength=120
nsrc=8000
fmax=30
file_shot=shotRQ_volume_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
echo $file_shot
#volume
zsrc1=500
zsrc2=4090
../fdelmodc \
file_cp=simple_cp.su ischeme=2 \
Qp=15 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=1 \
tmod=$tmod \
dxrcv=10.0 \
plane_wave=0 \
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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 10
#deep
zsrc1=2700
zsrc2=4090
file_shot=shotRQ_deep_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
../fdelmodc \
file_cp=simple_cp.su ischeme=2 \
file_den=simple_ro.su \
Qp=15 \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=1 \
tmod=$tmod \
dxrcv=10.0 \
plane_wave=0 \
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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 10
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -q long
#PBS -V
#
export PATH=../../bin:$PATH
makewave w=g2 fp=10 t0=0.15 dt=0.0010 nt=4096 file_out=G2.su verbose=1
makemod sizex=10000 sizez=4100 dx=10 dz=10 cp0=1500 ro0=1000 file_base=hom.su
xsrc1=100
xsrc2=9900
dxsrc=10
tmod=120
tsrc1=0.1
tsrc2=120
tlength=120
nsrc=8000
fmax=30
for wav_random in 0 1;
do
file_shot=shotHRS${wav_random}_volume_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
echo $file_shot
#volume
zsrc1=500
zsrc2=4090
../fdelmodc \
file_cp=hom_cp.su ischeme=1 \
file_den=hom_ro.su \
file_src=G2.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=1 \
tmod=$tmod \
dxrcv=10.0 \
plane_wave=0 \
src_random=1 \
wav_random=${wav_random} \
fmax=$fmax \
xsrc1=$xsrc1 \
xsrc2=$xsrc2 \
zsrc1=$zsrc1 \
zsrc2=$zsrc2 \
tsrc1=0.0 \
tsrc2=$tsrc2 \
tlength=$tlength \
nsrc=$nsrc \
dipsrc=0 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 10
#deep
zsrc1=2700
zsrc2=4090
file_shot=shotHRS${wav_random}_deep_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=hom_cp.su ischeme=1 \
file_den=hom_ro.su \
file_src=G2.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=1 \
tmod=$tmod \
dxrcv=10.0 \
plane_wave=0 \
src_random=1 \
wav_random=${wav_random} \
fmax=$fmax \
xsrc1=$xsrc1 \
xsrc2=$xsrc2 \
zsrc1=$zsrc1 \
zsrc2=$zsrc2 \
tsrc1=0.0 \
tsrc2=$tsrc2 \
tlength=$tlength \
nsrc=$nsrc \
dipsrc=0 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 10
#plane
zsrc1=2700
zsrc2=2700
file_shot=shotHRS${wav_random}_plane_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=hom_cp.su ischeme=1 \
file_den=hom_ro.su \
file_src=G2.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=1 \
tmod=$tmod \
dxrcv=10.0 \
plane_wave=1 \
xsrc=5000 zsrc=2700 \
src_random=0 \
wav_random=${wav_random} \
fmax=$fmax \
tsrc1=0.0 \
tsrc2=$tsrc2 \
tlength=$tlength \
nsrc=$nsrc \
dipsrc=0 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 10
done
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -q long
#PBS -V
#
# calls fdelmodc_long.scr, can not be reproduced; software in test phase
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."
exit;
./fdelmodc_long.scr
export OMP_NUM_THREADS=1
~/src/CorrSMP_FD_Files/corrsmp \
file_base=long/shotAcoustic_T3600_S1500_Dt500_F30_001_rvz.su \
nc=1 nstation=201 fullcorr=1 dt=0.008 \
nt=16384 verbose=1 ntc=4096 file_out=outI.su causal=1
#results long recordings and using corrsmp
suwind s=1 j=1 tmax=4 f1=0.0 key=fldr min=101 max=101 < outI_cc1.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > long_Corr_I.eps
~/src/CorrSMP_FD_Files/corrsmp \
file_base=long/shotAcoustic_T3600_S1500_Dt500_F30_001_rvz.su \
nc=1 nstation=201 fullcorr=1 dt=0.008 \
nt=16384 verbose=1 ntc=4096 file_out=outII causal=2
suwind s=1 j=1 tmax=4 f1=0.0 key=fldr min=101 max=101 < outII_cc1.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > long_Corr_II.eps
~/src/CorrSMP_FD_Files/corrsmp \
file_base=long/shotAcoustic_T3600_S1500_Dt500_F30_001_rvz.su \
nc=1 nstation=201 fullcorr=1 dt=0.008 \
nt=16384 verbose=1 ntc=4096 file_out=outIII causal=4
suwind s=1 j=1 tmax=4 f1=0.0 key=fldr min=101 max=101 < outIII_cc1.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > long_Corr_III.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -q long
#PBS -V
#
# calls fdelmodc_amplitude.scr, can not be reproduced; software in test phase
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."
exit;
./fdelmodc_amplitude.scr
#amplitude distribution of sources
supsgraph wbox=3 hbox=4 style=seimic < src_ampl.su \
f1=-250 d1=5.05 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitudeDistr.eps
export OMP_NUM_THREADS=1
~/src/CorrSMP_FD_Files/corrsmp \
file_base=long/shotAcousticA500_T3600_S1500_Dt500_F30_001_rvz.su \
nc=1 nstation=201 fullcorr=1 dt=0.008 \
nt=16384 verbose=1 ntc=4096 file_out=outA.su causal=1
suwind key=fldr min=101 max=101 < outA.su_cc1.su > outA_Station101_Comp1.su
suwind s=1 j=1 tmax=4 f1=0.0 < outA_Station101_Comp1.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > longA_Corr_I.eps
suwind s=1 j=1 tmax=4 f1=0.0 < outA_Station101_Comp1.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
sunormalize norm=rms | \
supsimage perc=99 f1=0 f2=-5000 hbox=4 wbox=3 \
label1='time (s)' label2='lateral position (m)' \
labelsize=10 f2num=-5000 d2num=2500 > longA_norm_Corr_I.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -q long
#PBS -V
#
# amplitude variations on source strength, 3x1500 s.
export PATH=../../bin:$PATH
makemod sizex=10000 sizez=4100 dx=10 dz=10 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=400,400 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 z=1100,1100,1100,1600,1100,1100,1100 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=2100,2100 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=2600,2600 cp=5500 ro=2200
xsrc1=500
xsrc2=9500
zsrc1=500
zsrc2=4090
tmod=120
tsrc1=0.1
tsrc2=120
tlength=120
nsrc=150
fmax=30
#Gaussian amplitude distribution
A=10000
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=$A \
distribution=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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
#use clip=500
SIrand.scr ${base}_rvz.su 50
#amplitude distribution of sources
f1=`surange < src_ampl.su | grep f1 | awk '{print $2 }'`
d1=`surange < src_ampl.su | grep d1 | awk '{print $2 }'`
supsgraph wbox=3 hbox=1 style=normal < src_ampl.su \
f1=$f1 d1=$d1 d1num=10000 x1end=50000 d2num=1 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitudeDistrGauss${A}.eps
#Flat amplitude distribution
A=50000
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=$A \
distribution=0 \
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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
#use clip=900
SIrand.scr ${base}_rvz.su 50
#amplitude distribution of sources
#amplitude distribution of sources
f1=$(echo "scale=3; -5*$A" | bc -l)
d1=`surange < src_ampl.su | grep d1 | awk '{print $2 }'`
supsgraph wbox=3 hbox=1 style=normal < src_ampl.su \
f1=$f1 d1=$d1 d1num=10000 x1end=50000 d2num=1 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitudeDistrFlat${A}.eps
#No amplitude distribution
A=0
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=3 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
#use clip=4e-6
SIrand.scr ${base}_rvz.su 50
suspike nt=1000 ntr=1 nspk=1 ix1=1 it1=500 | sugain scale=$nsrc | \
supsgraph wbox=3 hbox=1 style=normal \
f1=-50000 d1=100.1 d1num=10000 d2num=50 x1end=50000 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitude${A}.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -q long
#PBS -V
#
# computes only the amplitude distributions pictures, 5 s.
export PATH=../../bin:$PATH
makemod sizex=10000 sizez=4100 dx=10 dz=10 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=400,400 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 z=1100,1100,1100,1600,1100,1100,1100 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=2100,2100 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=2600,2600 cp=5500 ro=2200
xsrc1=500
xsrc2=9500
zsrc1=500
zsrc2=4090
tmod=0.2
tsrc1=0.1
tsrc2=0.2
tlength=0.2
nsrc=150
fmax=30
#Gaussian amplitude distribution
A=10000
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
../fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=$A \
distribution=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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
#amplitude distribution of sources
f1=`surange < src_ampl.su | grep f1 | awk '{print $2 }'`
d1=`surange < src_ampl.su | grep d1 | awk '{print $2 }'`
supsgraph wbox=3 hbox=1 style=normal < src_ampl.su \
f1=$f1 d1=$d1 d1num=10000 x1end=50000 d2num=1 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitudeDistrGauss${A}.eps
#Flat amplitude distribution
A=50000
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0010 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=$A \
distribution=0 \
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=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
#amplitude distribution of sources
#amplitude distribution of sources
f1=$(echo "scale=3; -5*$A" | bc -l)
d1=`surange < src_ampl.su | grep d1 | awk '{print $2 }'`
supsgraph wbox=3 hbox=1 style=normal < src_ampl.su \
f1=$f1 d1=$d1 d1num=10000 x1end=50000 d2num=1 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitudeDistrFlat${A}.eps
#No amplitude distribution
A=0
file_shot=shotRS_A${A}_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
suspike nt=1000 ntr=1 nspk=1 ix1=1 it1=500 | sugain scale=$nsrc | \
supsgraph wbox=3 hbox=1 style=normal \
f1=-50000 d1=100.1 d1num=10000 d2num=50 x1end=50000 linecolor=black labelsize=10 titlesize=10 \
label1=amplitude label2=occurence > amplitude${A}.eps
#!/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
#!/bin/bash
#PBS -l nodes=1
#PBS -N InterfModeling
#PBS -q fourweeks
#PBS -V
#
export PATH=../../bin:$PATH
which fdelmodc
makemod sizex=10000 sizez=4100 dx=10 dz=10 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=400,400 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=1100,1100,1100,1600,1100,1100,1100 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=2100,2100 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=2600,2600 cp=5500 ro=2200
makewave w=g2 fmax=45 t0=0.10 dt=0.001 nt=4096 db=-40 file_out=G2.su verbose=1
xsrc1=100
xsrc2=9900
zsrc1=2100
zsrc2=4000
file_shot=shotRandomPos${xsrc}_${zsrc1}.su
fdelmodc \
file_cp=simple_cp.su ischeme=1 \
file_den=simple_ro.su \
file_rcv=$file_shot \
dtrcv=0.008 \
dt=0.0010 \
verbose=4 \
tmod=10.000 \
dxrcv=20.0 \
zrcv1=10 \
zrcv2=10 \
xrcv1=0 \
xrcv2=10000 \
src_random=1 \
wav_random=1 \
fmax=30 \
xsrc1=$xsrc1 \
xsrc2=$xsrc2 \
zsrc1=$zsrc1 \
zsrc2=$zsrc2 \
tsrc1=0.0 \
tsrc2=5.0 \
nsrc=20 \
dipsrc=0 \
ntaper=45 \
tsnap1=0.1 tsnap2=6.0 dtsnap=0.2 \
left=4 right=4 top=1 bottom=4 \
nxmax=2500 nzmax=1400 ntmax=10000
psgraph < srcTimeLengthN\=10001.bin n1=10001 \
labelsize=12 d1=0.001 style=normal linecolor=blue \
label1="start time (s)" label2="source duration (s)" \
d1num=1 d2num=1 wbox=8 hbox=4 x1end=5 > srcTimeLength.eps
supswigp < src_nwav.su \
labelsize=10 label1='time (s)' label2='source number' x1end=6 \
d2=1 d2num=1 hbox=4 wbox=6 fill=0 \
titlesize=-1 > src_nwav.eps
suwind < src_nwav.su key=tracl min=11 max=11 | \
supsgraph hbox=2 wbox=4 style=normal \
labelsize=10 label2='amplitude' label1='time (s)' \
titlesize=-1 d1num=1.0 > src11_wiggle.eps
suwind < src_nwav.su key=tracl min=11 max=11 | \
supsgraph hbox=2 wbox=4 style=normal \
labelsize=10 label2='amplitude' label1='time (s)' \
titlesize=-1 x1end=0.05 > src11_wiggle_zbeg.eps
suwind < src_nwav.su key=tracl min=11 max=11 | \
supsgraph hbox=2 wbox=4 style=normal \
labelsize=10 label2='amplitude' label1='time (s)' \
titlesize=-1 x1beg=3.60 x1end=3.65 > src11_wiggle_zend.eps
suwind < src_nwav.su key=tracl min=11 max=11 | \
sufft | suamp| supsgraph hbox=2 wbox=4 style=normal \
labelsize=10 label2='amplitude' label1='frequency (Hz)' \
titlesize=-1 x1end=100 d1num=10 > src11_ampl.eps
fconv file_in1=src_nwav.su auto=1 shift=1 mode=cor1 | \
sugain qbal=1 | \
supswigp x1beg=-1 x1end=1 d2num=1 hbox=4 wbox=6 \
labelsize=10 label2='source number' label1='time (s)' \
titlesize=-1 fill=0 > src_nwav_autoCorr_Norm.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# starts fdelmodc only to compute the source positions, 1 s.
export PATH=../../bin:$PATH
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=0 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=900,900 cp=1500 ro=1000 \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc1=100
xsrc2=9900
dxsrc=10
#volume
zsrc1=500
zsrc2=4090
tmod=0.01
tsrc2=0.1
tlength=0.01
nsrc=1000
fmax=30
#Figure 2b,c,d,e,f
file_shot=shotR_T${tmod}_S${nsrc}_Dt${tlength}_F${fmax}.su
echo $file_shot
#dummy modeling just to generate the source positions to be used in the Figure
../fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0005 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=5000 zsrc=2700 \
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 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
#Figure 2
suop2 SrcRecPositions.su vel_cp.su w1=7000 w2=1.0 op=sum | \
sugain nclip=1500 | \
sushw key=f1,f2 a=0,-5000 | \
supsimage hbox=6 wbox=8 labelsize=10 \
f2num=-5000 d2num=1000 f1num=-500 d1num=500 \
x2beg=-5000 x2end=5000 x1beg=-900 x1end=4000 \
bclip=7000 wclip=0 label1="depth (m)" label2="lateral position (m)" \
x2beg=-5000 x2end=5000 \
> simple_sources_cp.eps
# laternative color scheme for above picture
# wrgb=0,0,1.0 grgb=0,1.0,0 brgb=1.0,0,0 x2beg=-5000 x2end=5000 \
supsimage hbox=6 wbox=8 labelsize=18 < vel_cp.su \
f2=-5000 f2num=-5000 d2num=1000 f1num=-500 d1num=500 \
x2beg=-5000 x2end=5000 x1beg=-900 x1end=4000 \
label1="depth (m)" label2="lateral position (m)" \
wrgb=0,0,1.0 grgb=0,1.0,0 brgb=1.0,0,0 x2beg=-5000 x2end=5000 \
> simple_cp.eps
# use adapted psimage from su to plot points using curve function
supsimage hbox=6 wbox=8 labelsize=10 < simple_cp.su \
f1=-900 f2=-5000 f2num=-5000 d2num=1000 f1num=-500 d1num=500 \
x2beg=-5000 x2end=5000 x1beg=-900 x1end=4000 \
label1="depth (m)" label2="lateral position (m)" \
wrgb=0,0,1.0 grgb=0,1.0,0 brgb=1.0,0,0 x2beg=-5000 x2end=5000 \
curve=SrcPositions1000.txt npair=1000 curvecolor=black curvedash=-1 curvewidth=1 \
> simple_srcpos_cp.eps
supsimage hbox=6 wbox=8 labelsize=18 < SrcRecPositions.su\
f2num=50000 d2num=1000 f1num=50000 d1num=500 \
bclip=1 wclip=0 label1=" " label2=" "\
> sources.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
cp=2000
rho=1000
dx=2.5
dt=0.0001
halfdt=`perl -e "print 0.5*$dt;"`
export PATH=../../bin:$PATH
makemod sizex=2000 sizez=2000 dx=$dx dz=$dx cp0=$cp ro0=$rho orig=-1000,0 file_base=simple.su
makewave fp=15 dt=$dt file_out=wave.su nt=4096 t0=0.1
######### MONOPOLE #######
../fdelmodc \
file_cp=simple_cp.su ischeme=1 iorder=4 \
file_den=simple_ro.su \
file_src=wave.su \
file_rcv=shot_fd.su \
src_type=1 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0005 \
verbose=2 \
tmod=0.5115 \
dxrcv=5.0 \
xrcv1=-500 xrcv2=500 \
zrcv1=500 zrcv2=500 \
xsrc=0 zsrc=1000 \
ntaper=80 \
left=4 right=4 top=4 bottom=4
makewave fp=15 dt=0.0005 file_out=wave.su nt=4096 t0=0.1
green c=$cp rho=$rho file_src=wave.su zsrc1=500 xrcv=-500,500 dxrcv=5 nt=4096 dip=0 | suwind nt=1024 > shot_green_rp.su
green c=$cp rho=$rho file_src=wave.su zsrc1=500 xrcv=-500,500 dxrcv=5 nt=4096 dip=0 p_vz=1 | suwind nt=1024 > shot_green_rvz.su
# rp
(suwind key=tracl min=101 max=101 < shot_fd_rp.su | basop choice=shift shift=$halfdt; suwind key=tracl min=101 max=101 < shot_green_rp.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=4 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 d2num=0.04 x2end=0.080 f2num=-0.4 > mon_rp.eps
(suwind key=tracl min=101 max=101 < shot_fd_rp.su | basop choice=shift shift=$halfdt; suwind key=tracl min=101 max=101 < shot_green_rp.su;) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=2 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green x1beg=0.255 x1end=0.258 x2beg=0.0785 x2end=0.0797 d2num=0.0002 f2num=0.0786 > mon_zoom_rp.eps
suwind key=tracl min=101 max=101 < shot_fd_rp.su | basop choice=shift shift=$halfdt > trace_fd.su
suwind key=tracl min=101 max=101 < shot_green_rp.su > trace_green.su
sumax < trace_green.su outpar=nep
Mmax=`cat nep | awk '{print $1}'`
a=`perl -e "print 100.0/$Mmax;"`
echo $a
sudiff trace_green.su trace_fd.su | basop choice=shift shift=-0.1 | sugain scale=$a | suop op=abs | supsgraph style=normal labelsize=12 wbox=4 hbox=2 label1='time in seconds' label2="Relative error in percentage of peak" linewidth=0.1 linecolor=red x2beg=0.0 x2end=1.0 f2num=0.0 d2num=0.5 > mon_diff_dx${dx}_rp.eps
(suwind key=tracl min=101 max=101 < shot_fd_rp.su | basop choice=shift shift=$halfdt; suwind key=tracl min=101 max=101 < shot_green_rp.su ) | basop choice=shift shift=-0.1 | suxgraph
# rvz
(suwind key=tracl min=101 max=101 < shot_fd_rvz.su ; suwind key=tracl min=101 max=101 < shot_green_rvz.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=4 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 d2num=2e-8 f2num=-4e-8 x2end=4e-8 x2beg=-2.7e-8 > mon_rvz.eps
(suwind key=tracl min=101 max=101 < shot_fd_rvz.su ; suwind key=tracl min=101 max=101 < shot_green_rvz.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=2 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green d2num=2e-10 f2num=3.86e-8 x2end=3.96e-8 x2beg=3.86e-8 x1beg=0.255 x1end=0.258 > mon_zoom_rvz.eps
suwind key=tracl min=101 max=101 < shot_fd_rvz.su > trace_fd.su
suwind key=tracl min=101 max=101 < shot_green_rvz.su > trace_green.su
sumax < trace_green.su outpar=nep
Mmax=`cat nep | awk '{print $1}'`
a=`perl -e "print 100.0/$Mmax;"`
echo $a
sudiff trace_green.su trace_fd.su | basop choice=shift shift=-0.1 | sugain scale=$a | suop op=abs | supsgraph style=normal labelsize=12 wbox=4 hbox=2 label1='time in seconds' label2="Relative error in percentage of peak" linewidth=0.1 linecolor=red x2beg=0.0 x2end=1.0 f2num=0.0 d2num=0.5 > mon_diff_dx${dx}_rvz.eps
(suwind key=tracl min=101 max=101 < shot_fd_rvz.su ; suwind key=tracl min=101 max=101 < shot_green_rvz.su) | basop choice=shift shift=-0.1 | suxgraph
exit;
######### DIPOLE #######
makewave fp=15 dt=$dt file_out=wave.su nt=4096 t0=0.1
fdelmodc \
file_cp=simple_cp.su ischeme=1 iorder=4 \
file_den=simple_ro.su \
file_src=wave.su \
file_rcv=shot_fd_dip.su \
src_type=1 \
src_orient=2 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0005 \
verbose=2 \
tmod=0.5115 \
dxrcv=5.0 \
xrcv1=-500 xrcv2=500 \
zrcv1=500 zrcv2=500 \
xsrc=0 zsrc=1000 \
ntaper=80 \
left=4 right=4 top=4 bottom=4
makewave fp=15 dt=0.0005 file_out=wave.su nt=4096 t0=0.1
green c=$cp rho=$rho file_src=wave.su zsrc1=500 xrcv=-500,500 dxrcv=5 nt=4096 dip=1 | suwind nt=1024 > shot_green_dip_rp.su
green c=$cp rho=$rho file_src=wave.su zsrc1=500 xrcv=-500,500 dxrcv=5 nt=4096 p_vz=1 dip=1 | suwind nt=1024 > shot_green_dip_rvz.su
shift=`perl -e "print 0.5*$dx/$cp;"`
# rp
(suwind key=tracl min=101 max=101 < shot_fd_dip_rp.su | basop choice=shift shift=-$shift; suwind key=tracl min=101 max=101 < shot_green_dip_rp.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal titlesize=-1 labelsize=10 wbox=4 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green d2num=0.001 x2end=0.004001 f2num=-0.03 > dip_rp.eps
(suwind key=tracl min=101 max=101 < shot_fd_dip_rp.su | basop choice=shift shift=-$shift; suwind key=tracl min=101 max=101 < shot_green_dip_rp.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=2 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green x2beg=0.00379 x2end=0.00388 d2num=0.00002 f2num=0.00380 x1beg=0.2430 x1end=0.24601 f1num=0.243 d1num=0.002 > dip_zoom_rp.eps
suwind key=tracl min=101 max=101 < shot_fd_dip_rp.su | basop choice=shift shift=-$shift > trace_fd.su
suwind key=tracl min=101 max=101 < shot_green_dip_rp.su > trace_green.su
sumax < trace_green.su outpar=nep
Mmax=`cat nep | awk '{print $1}'`
a=`perl -e "print 100.0/$Mmax;"`
echo $a
sudiff trace_green.su trace_fd.su | basop choice=shift shift=-0.1 | sugain scale=$a | suop op=abs | supsgraph style=normal labelsize=12 wbox=4 hbox=2 label1='time in seconds' label2="Relative error in percentage of peak" linewidth=0.1 linecolor=red x2beg=0.0 x2end=1.0 f2num=0.0 d2num=0.5 > dip_diff_dx${dx}_rp.eps
(suwind key=tracl min=101 max=101 < shot_fd_dip_rp.su | basop choice=shift shift=-$shift ; suwind key=tracl min=101 max=101 < shot_green_dip_rp.su ) | basop choice=shift shift=-0.1 | suxgraph
# rvz
(suwind key=tracl min=101 max=101 < shot_fd_dip_rvz.su | basop choice=shift shift=-$shift; suwind key=tracl min=101 max=101 < shot_green_dip_rvz.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=4 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green d2num=6e-10 f2num=-1.2e-9 x2end=2e-9 x2beg=-1.3e-9 > dip_rvz.eps
(suwind key=tracl min=101 max=101 < shot_fd_dip_rvz.su | basop choice=shift shift=-$shift; suwind key=tracl min=101 max=101 < shot_green_dip_rvz.su ) | basop choice=shift shift=-0.1 | supsgraph style=normal labelsize=10 wbox=2 hbox=2 label1='time in seconds' label2="Amplitude" linewidth=0.1 linecolor=red,green d2num=2e-11 f2num=1.85e-9 x2end=1.944e-9 x2beg=1.85e-9 x1beg=0.2430 x1end=0.24601 f1num=0.243 d1num=0.002 > dip_zoom_rvz.eps
suwind key=tracl min=101 max=101 < shot_fd_dip_rvz.su | basop choice=shift shift=-$shift > trace_fd.su
suwind key=tracl min=101 max=101 < shot_green_dip_rvz.su > trace_green.su
sumax < trace_green.su outpar=nep
Mmax=`cat nep | awk '{print $1}'`
a=`perl -e "print 100.0/$Mmax;"`
echo $a
sudiff trace_green.su trace_fd.su | basop choice=shift shift=-0.1 | sugain scale=$a | suop op=abs | supsgraph style=normal labelsize=12 wbox=4 hbox=2 label1='time in seconds' label2="Relative error in percentage of peak" linewidth=0.1 linecolor=red x2beg=0.0 x2end=1.0 f2num=0.0 d2num=0.5 > dip_diff_dx${dx}_rvz.eps
(suwind key=tracl min=101 max=101 < shot_fd_dip_rvz.su | basop choice=shift shift=-$shift; suwind key=tracl min=101 max=101 < shot_green_dip_rvz.su ) | basop choice=shift shift=-0.1 | suxgraph
#!/bin/bash
#
#Forward model all source positions one by one (this takes a very long time)
#
# calls Simple_model_base, and Simple_model_sides.scr, 122 hours!
export PATH=../../bin:$PATH
./Simple_model_base.scr
./Simple_model_sides.scr
# Correlation horizontal layer
export dxsrc=10
export xsrc1=1000
export xsrc2=9000
export xrcv1=0
export xrcv2=10000
export file_out=/tmp/corr.su
export file_in=T_simple_dxrcv50_dx10.su
export dxrcv=50
(( nshots = ($xsrc2-$xsrc1)/$dxsrc + 1 ))
(( nrecv = ($xrcv2-$xrcv1)/$dxrcv + 1 ))
(( middle = $xsrc1 + (($nshots-1)/2)*$dxsrc ))
(( xa = ($nrecv+1)/2 ))
(( xsrc_sim = $xrcv1 + ($xa-1)*$dxrcv ))
echo xsrc_sim=$xsrc_sim
echo nshot=$nshots
echo nrec=$nrecv
echo xa=$xa
echo middle=$middle
# first do the correlation with the middle trace of each shot
rm -rf $file_out
for (( xsrc = $xsrc1; xsrc<=$xsrc2; xsrc+=$dxsrc )) do
echo $xsrc $middle;
sxkey=${xsrc}000
suwind key=sx min=$sxkey max=$sxkey < $file_in > shot.su
suwind key=tracl min=$xa max=$xa < shot.su | \
fconv mode=cor2 file_in1=shot.su shift=1 >> $file_out
done;
#then sort on constant gx values and stack all traces belonging to the same gx (=sum over sx)
(( f2 = -($middle-$xrcv1) ))
susort < $file_out gx > /tmp/corr_sort.su
sustack key=gx < /tmp/corr_sort.su | sushw key=f2 a=$f2 > corr_stack${dxrcv}_${xsrc_sim}_${dxsrc}.su
#Correlation sides
export file_in=T_simple_sides_dxrcv50_dx10.su
nshots=720
export fldr1=902
export fldr2=1621
echo $xsrc_sim
echo $nshots
echo $nrecv
echo $xa
echo $middle
rm -rf $file_out
for (( fldr = $fldr1; fldr<=$fldr2; fldr+=1 )) do
echo $fldr;
suwind key=fldr min=$fldr max=$fldr < $file_in > shot.su
suwind key=tracl min=$xa max=$xa < shot.su | \
fconv mode=cor2 file_in1=shot.su shift=1 >> $file_out
done;
(( f2 = -($middle-$xrcv1) ))
susort < $file_out gx > corr_sort.su
sustack key=gx < corr_sort.su | sushw key=f2 a=$f2 > corr_stack${dxrcv}_sides_$xsrc_sim.su
rm $file_out /tmp/corr_sort.su shot.su
# sources at left and right sides of the model receivers on top gray scales dxrc
cat corr_stack50_sides_5000.su | \
suwind tmax=4 tmin=0 f1=-8.008 | \
sushw key=f1,delrt a=0.0,0.0 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage clip=2e-15 f1=0 f2=-5000 hbox=4 wbox=3 \
labelsize=10 f2num=-5000 d2num=2500 \
label1='time (s)' label2='lateral position (m)' > SimCorr50_sides_sx5000_simple.eps
cat corr_stack50_5000_10.su | \
suwind tmax=4 tmin=0 f1=-8.008 | \
sushw key=f1,delrt a=0.0,0.0 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage clip=2e-15 f1=0 f2=-5000 hbox=4 wbox=3 \
labelsize=10 f2num=-5000 d2num=2500 \
label1='time (s)' label2='lateral position (m)'> SimCorr50_sx5000_simple.eps
susum corr_stack50_5000_10.su corr_stack50_sides_5000.su | \
suwind tmax=4 tmin=0 f1=-8.008 | \
sushw key=f1,delrt a=0.0,0.0 | \
sufilter amps=0,0.5,1,1,0 f=0,2,3,50,60 | \
supsimage clip=2e-15 f1=0 f2=-5000 hbox=4 wbox=3 \
labelsize=10 f2num=-5000 d2num=2500 \
label1='time (s)' label2='lateral position (m)'> SimCorr50_add_sides_sx5000_simple.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# direct modeled reference result Figure 3d, 500 s.
export PATH=../../bin:$PATH
makewave w=g1 fmax=45 t0=0.15 dt=0.0005 nt=4096 db=-40 file_out=G1.su verbose=1
cp G1.su G1_c.su
fconv file_in1=G1.su file_in2=G1_c.su mode=cor1 verbose=1 > corrw.su
basop choice=shift shift=0.2 dx=1 file_in=corrw.su > G1corr.su
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
sushw key=f1 a=-900 < simple_cp.su > vel_cp.su
sushw key=f1 a=-900 < simple_ro.su > vel_ro.su
xsrc=5000
file_shot=shotRefNoFree_x${xsrc}.su
../fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_src=G1corr.su \
fmax=45 \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
verbose=4 \
rec_delay=0.2 \
tmod=4.2 \
dxrcv=50.0 \
zrcv1=0 \
zrcv2=0 \
xrcv1=0 xrcv2=10000 \
plane_wave=0 \
amplitude=0 \
xsrc=${xsrc} zsrc=0 \
ntaper=45 \
left=4 right=4 top=4 bottom=4
suwind s=1 j=1 tmax=4.008 f1=0.0 < shotRefNofree_x${xsrc}_rvz.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sugain scale=-1 | \
supsimage f1=0 f2=-5000 hbox=4 wbox=3 x1end=4.0 \
labelsize=10 f2num=-5000 d2num=2500 verbose=1 clip=2e-7 \
label1='time (s)' label2='lateral position (m)' > shotRefNofree_5000_0.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# direct modeled reference result Figure 3d, 500 s.
export PATH=../../bin:$PATH
makewave w=g1 fmax=45 t0=0.15 dt=0.0005 nt=4096 db=-40 file_out=G1.su verbose=1
cp G1.su G1_c.su
fconv file_in1=G1.su file_in2=G1_c.su mode=cor1 verbose=1 > corrw.su
basop choice=shift shift=0.2 dx=1 file_in=corrw.su > G1corr.su
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc=5000
file_shot=shotRef_x${xsrc}.su
../fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_src=G1corr.su \
fmax=45 \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
verbose=4 \
rec_delay=0.2 \
tmod=0.5 \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=${xsrc} zsrc=5 \
ntaper=45 \
left=4 right=4 top=1 bottom=4 verbose=1
suwind s=1 j=1 tmax=4.008 f1=0.0 < shotRef_x${xsrc}_rvz.su | \
sushw key=f1,delrt,d2 a=0.0,0.0,50 | \
sugain scale=-1 | \
supsimage f1=0 f2=-5000 hbox=4 wbox=3 x1end=4.0 \
labelsize=10 f2num=-5000 d2num=2500 verbose=1 clip=2e-7 \
label1='Time (s)' label2='Lateral position (m)' > shotRef_5000_0.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# inside the length loop the script SIrand.scr is called to compute the retrieved response from the modeled data
# 5 different source signature lengths, 5x3.5 hours
export PATH=../../bin:$PATH
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc1=100
xsrc2=9900
dxsrc=10
#volume
zsrc1=500
zsrc2=4090
tmod=120
tsrc2=120
tlength=120
nsrc=1000
fmax=30
#Figure 7b,c,d,e,f
for tlength in 120 60 30 10 5;
do
file_shot=shotR_T${tmod}_S${nsrc}_Dt${tlength}_F${fmax}.su
echo $file_shot
fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0005 \
verbose=2 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=5000 zsrc=2700 \
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 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su 50
done
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# simulates 8000 short (2.5 s) sources, 3.5 hours
export PATH=../../bin:$PATH
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc1=100
xsrc2=9900
#volume
tmod=120
tsrc2=120
fmax=30
zsrc1=500
zsrc2=4090
tlength=5
nsrc=8000
file_shot=shotR_T${tmod}_S${nsrc}_Dt${tlength}_F${fmax}.su
fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0005 \
verbose=2 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=5000 zsrc=2700 \
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 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# 5 different number of random sources, 5x3.5 hours
export PATH=../../bin:$PATH
#makewave w=g1 fp=10 t0=0.15 dt=0.0010 nt=4096 file_out=G2.su verbose=1
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc1=100
xsrc2=9900
#volume
zsrc1=500
zsrc2=4090
tmod=120
tsrc2=120
tlength=120
fmax=30
#Figure 6a,b,c,d,e
rm Trace.su
for nsrc in 8000 1000 500 100 50;
do
file_shot=shotR_T${tmod}_S${nsrc}_Dt${tsrc2}_F${fmax}.su
echo $file_shot
fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0005 \
verbose=2 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=5000 zsrc=2700 \
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 \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su
ntraces=`surange < ${base}_rvz.su | grep traces| awk '{print $1 }'`
echo $ntraces
middle=$(echo "scale=0; ($ntraces+1)/2"| bc -l)
echo $middle
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,8000 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
done
#add reference trace to trace comparison
suwind key=tracf min=101 max=101 < shotRef_x5000_rvz.su | sugain scale=-1 >> Trace.su
suwind tmin=1.5 tmax=3.0 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supsgraph \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
linecolor=red,green,blue,emerald,black f1=1.5 d1num=0.5 \
label1='time (s)' > shotTraces_T120_S_Dt120_F30.eps
suwind tmin=1.5 tmax=3.1 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supswigp \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
f1=1.5 d1num=0.5 linewidth=1 fill=0 f2=1 d2=1 d2num=1 \
label1='time (s)' label2='number of sources' > shotWiggles_T120_S_Dt120_F30.eps
#!/bin/bash
#
# make postscript file of middle trace after Figure6.scr, 1 s.
shot=shotR_T120_S8000_Dt120_F30_rvz.su
base=`echo $shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
ntraces=`surange < $shot | grep traces| awk '{print $1 }'`
echo $ntraces
middle=$(echo "scale=0; ($ntraces+1)/2"| bc -l)
echo $middle
./SIrand.scr $shot
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,8000 | \
suwind key=tracl min=$middle max=$middle > Trace.su
shot=shotR_T120_S1000_Dt120_F30_rvz.su
./SIrand.scr $shot
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,1000 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
shot=shotR_T120_S500_Dt120_F30_rvz.su
./SIrand.scr $shot
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,500 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
shot=shotR_T120_S100_Dt120_F30_rvz.su
./SIrand.scr $shot
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,100 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
suwind s=1 j=1 tmax=4 f1=0.0 < shotRef_5000_0_rvz.su | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,100 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
suwind tmin=1.5 tmax=3.0 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supsgraph \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
linecolor=red,green,blue,emerald,black f1=1.5 d1num=0.5 \
> shotTraces_T120_S_Dt120_F30.eps
suwind tmin=1.5 tmax=3.1 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supswigp \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
f1=1.5 d1num=0.5 linewidth=1 fill=0 f2=1 d2=1 \
label1='time (s)' label2='number of sources' > shotWiggles_T120_S_Dt120_F30.eps
#!/bin/bash
#PBS -l nodes=1:ppn=2
#PBS -N InterfModeling
#PBS -V
#
# alternative not used in paper, fixed source signature length, 5x3.5 hours
export PATH=../../bin:$PATH
#makewave w=g1 fp=10 t0=0.15 dt=0.0010 nt=4096 file_out=G2.su verbose=1
makemod sizex=10000 sizez=5000 dx=5 dz=5 cp0=1500 ro0=1000 file_base=simple.su \
intt=def poly=0 x=0,10000 z=1300,1300 cp=2000 ro=1400 \
intt=def poly=2 x=0,2000,3000,5000,7000,8000,10000 \
z=2000,2000,2000,2500,2000,2000,2000 cp=4000 ro=2000 \
intt=def poly=0 x=0,10000 z=3000,3000 cp=3000 ro=1500 \
intt=def poly=0 x=0,10000 z=3500,3500 cp=5500 ro=2200
suwind itmin=181 < simple_cp.su | sushw key=f1 a=0 > vel_cp.su
suwind itmin=181 < simple_ro.su | sushw key=f1 a=0 > vel_ro.su
xsrc1=100
xsrc2=9900
#volume
zsrc1=500
zsrc2=4090
tmod=120
tsrc2=120
tlength=120
fmax=30
#Figure 6a,b,c,d,e
rm Trace.su
for nsrc in 8000 1000 500 100 50;
do
file_shot=shotRVt_T${tmod}_S${nsrc}_Dt${tlength}_F${fmax}.su
echo $file_shot
fdelmodc \
file_cp=vel_cp.su ischeme=1 \
file_den=vel_ro.su \
file_rcv=$file_shot \
rec_type_p=0 \
dtrcv=0.008 \
rec_ntsam=16384 \
dt=0.0005 \
verbose=4 \
tmod=$tmod \
dxrcv=50.0 \
plane_wave=0 \
amplitude=0 \
xsrc=5000 zsrc=2700 \
src_random=1 \
wav_random=1 \
length_random=0 \
fmax=$fmax \
xsrc1=$xsrc1 \
xsrc2=$xsrc2 \
zsrc1=$zsrc1 \
zsrc2=$zsrc2 \
tsrc1=0.0 \
tsrc2=$tsrc2 \
tlength=$tlength \
nsrc=$nsrc \
ntaper=45 \
left=4 right=4 top=1 bottom=4
base=`echo $file_shot | awk 'BEGIN { FS = "." } ; { print $1 }'`
echo $base
SIrand.scr ${base}_rvz.su
ntraces=`surange < ${base}_rvz.su | grep traces| awk '{print $1 }'`
echo $ntraces
middle=$(echo "scale=0; ($ntraces+1)/2"| bc -l)
echo $middle
susum causal.su noncausal.su | \
suwind s=1 j=1 tmax=4 f1=0.0 | \
sushw key=f1,delrt,d2,fldr a=0.0,0.0,50,8000 | \
suwind key=tracl min=$middle max=$middle >> Trace.su
done
suwind tmin=1.5 tmax=3.0 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supsgraph \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
linecolor=red,green,blue,emerald,black f1=1.5 d1num=0.5 \
label1='time (s)' > shotTracesVt_T120_S_Dt120_F30.eps
suwind tmin=1.5 tmax=3.1 f1=0.0 < Trace.su | \
sunormalize norm=max | \
supswigp \
hbox=4 wbox=3 labelsize=10 x1beg=1.5 x1end=3.0 \
f1=1.5 d1num=0.5 linewidth=1 fill=0 f2=1 d2=1 \
label1='time (s)' label2='number of sources' > shotWigglesVt_T120_S_Dt120_F30.eps
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment