iterations.scr

#!/bin/bash

#second reflector at time:
# 800/1800+600/2300
# .70531400966183574878 => sample 176
#third reflector at time model.scr:
# 800/1800+600/2300+800/2000
# 1.10531400966183574878 sample 276
#third reflector at time modepm.scr
#800/1800+600/3200+520/2000
#.96531400966183574878 sample 241

export ROOT=`grep ROOT= ../../../Make_include | head -1 | sed 's/ROOT=//' `
#adjust the PATH to where the code is installed
export PATH=$ROOT/bin:$PATH:

R=shotsdx5_rp.su
makewave fp=20 dt=0.004 file_out=wave.su nt=1024 t0=0.0 scale=0 scfft=1
select=451

if [[ "$1" == "Figure34910" ]];
then
for istart in 276
do
(( iend = istart + 1 ))

marchenko_primaries file_shot=$R ishot=$select file_src=wave.su \
	nshots=901 verbose=2 istart=$istart iend=$iend fmax=90 \
	pad=512 niter=33 smooth=8 niterskip=600 shift=16 file_rr=pred_rr_${istart}.su T=0 file_iter=iter.su 
done
fi

if [[ "$1" == "Figure6" ]];
then
istart=200
(( iend = istart + 1 ))
marchenko_primaries file_shot=$R ishot=$select file_src=wave.su \
	nshots=901 verbose=2 istart=$istart iend=$iend fmax=90 \
	pad=512 niter=30 smooth=8 niterskip=600 shift=16 file_rr=pred_rr_${istart}.su T=0 file_iter=iter.su 
fi

if [[ "$1" == "Figure11" ]];
then
for istart in 246 256 266 276 286 296 306 316
do
(( iend = istart + 1 ))

marchenko_primaries file_shot=$R ishot=$select file_src=wave.su \
	nshots=901 verbose=2 istart=$istart iend=$iend fmax=90 \
	pad=512 niter=33 smooth=8 niterskip=600 shift=16 file_rr=pred_rr_${istart}.su T=0 file_iter=iter.su 
done
fi

if [[ "$1" == "Figure13" ]];
then
istart=276
(( iend = istart + 1 ))
marchenko_primaries file_shot=$R ishot=$select file_src=wave.su \
	nshots=901 verbose=2 istart=$istart iend=$iend fmax=90 \
	pad=512 niter=33 smooth=8 niterskip=600 shift=16 file_rr=pred_rr_${istart}.su T=1 file_iter=iterT.su 
fi