From d4d850495207094dabc667d59f0f2651afa8dfbc Mon Sep 17 00:00:00 2001
From: JanThorbecke <janth@xs4all.nl>
Date: Thu, 16 Jun 2016 09:26:02 +0200
Subject: [PATCH] demo Marchenko
---
marchenko/demo/README | 9 ++
marchenko/demo/direct.scr | 37 ++++++
marchenko/demo/first_arrival.scr | 92 ++++++++++++++
marchenko/demo/marchenko.scr | 39 ++++++
marchenko/demo/model.scr | 212 +++++++++++++++++++++++++++++++
marchenko/demo/referenceShot.scr | 44 +++++++
marchenko/demo/remove_direct.scr | 38 ++++++
marchenko/demo/shots.scr | 77 +++++++++++
8 files changed, 548 insertions(+)
create mode 100644 marchenko/demo/README
create mode 100755 marchenko/demo/direct.scr
create mode 100755 marchenko/demo/first_arrival.scr
create mode 100755 marchenko/demo/marchenko.scr
create mode 100755 marchenko/demo/model.scr
create mode 100755 marchenko/demo/referenceShot.scr
create mode 100755 marchenko/demo/remove_direct.scr
create mode 100755 marchenko/demo/shots.scr
diff --git a/marchenko/demo/README b/marchenko/demo/README
new file mode 100644
index 0000000..5fc5036
--- /dev/null
+++ b/marchenko/demo/README
@@ -0,0 +1,9 @@
+Description of files:
+1) shots.scr create the shots
+2) model.scr computes the model
+3) direct_wave.scr crate the direct wave to be removed from the shots
+4) remove_direct.scr remove the direct wave from the shots and scale them
+5) first_arrival.scr computes the first arrival
+6) marchenko.scr perform the Marchenko scheme
+7) referenceShot.scr creates the reference Green's function
+
diff --git a/marchenko/demo/direct.scr b/marchenko/demo/direct.scr
new file mode 100755
index 0000000..a27121e
--- /dev/null
+++ b/marchenko/demo/direct.scr
@@ -0,0 +1,37 @@
+#!/bin/bash
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
+
+dx=2.5
+dt=0.0005
+
+makemod sizex=12000 sizez=4000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
+ orig=-6000,-1000 file_base=noContrast.su
+
+export OMP_NUM_THREADS=8
+makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0
+
+$HOME/bin/fdelmodc \
+ file_cp=noContrast_cp.su ischeme=1 iorder=4 \
+ file_den=noContrast_ro.su \
+ file_src=wavefw.su \
+ file_rcv=direct.su \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=1 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ rec_delay=0.3 \
+ dtrcv=0.004 \
+ verbose=2 \
+ tmod=4.394 \
+ dxrcv=10.0 \
+ xrcv1=-6000 xrcv2=6000 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=0 \
+ ntaper=400 \
+ left=4 right=4 top=4 bottom=4
+
diff --git a/marchenko/demo/first_arrival.scr b/marchenko/demo/first_arrival.scr
new file mode 100755
index 0000000..721721d
--- /dev/null
+++ b/marchenko/demo/first_arrival.scr
@@ -0,0 +1,92 @@
+#!/bin/bash
+#PBS -N fdelmod
+#PBS -q verylong
+#PBS -l nodes=1
+#PBS -k eo
+#PBS -j eo
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
+
+dx=2.5
+dt=0.0005
+
+makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
+ orig=-3000,0 file_base=synclDown.su verbose=2 \
+ intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
+ intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
+ intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
+ intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
+ intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
+
+# intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
+# intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
+# intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
+# intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
+# intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
+
+makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
+
+#smooth file_in=synclDown_cp.su power=-1.0 ntsm=29 nxsm=29 niter=15 file_out=syncls_cp.su
+#smooth file_in=synclDown_ro.su power=-1.0 ntsm=29 nxsm=29 niter=15 file_out=syncls_ro.su
+
+dxshot=10
+ishot=300
+nshots=301
+
+export OMP_NUM_THREADS=1
+mkdir -p shots
+mkdir -p jobs
+
+while (( ishot < nshots ))
+do
+
+ (( xsrc = -3000 + ${ishot}*${dxshot} ))
+# (( xsrc = -1100 + ${ishot}*${dxshot} ))
+ echo xsrc=$xsrc
+ file_rcv=shots/shotsmonPz1100_${xsrc}.su
+
+ cat << EOF > jobs/pbs_$ishot.job
+#!/bin/bash
+#
+#PBS -q medium
+#PBS -N mod_${xsrc}
+#PBS -j eo
+#PBS -m n
+#PBS -l nodes=1
+#PBS -V
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
+export OMP_NUM_THREADS=1
+
+ $HOME/bin64/fdelmodc \
+ file_cp=synclDown_cp.su ischeme=1 iorder=4 \
+ file_den=synclDown_ro.su \
+ file_src=wave.su \
+ file_rcv=$file_rcv \
+ src_type=1 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ rec_delay=0.1 \
+ dtrcv=0.004 \
+ verbose=2 \
+ tmod=2.100 \
+ dxrcv=10.0 \
+ xrcv1=-3000 xrcv2=3000 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=$xsrc zsrc=1100 \
+ ntaper=300 \
+ left=4 right=4 top=4 bottom=4
+EOF
+
+ qsub jobs/pbs_$ishot.job
+
+ (( ishot = $ishot + 1))
+done
+
+
+
diff --git a/marchenko/demo/marchenko.scr b/marchenko/demo/marchenko.scr
new file mode 100755
index 0000000..bee3690
--- /dev/null
+++ b/marchenko/demo/marchenko.scr
@@ -0,0 +1,39 @@
+#!/bin/bash -x
+#PBS -N fdelmod
+#PBS -q verylong
+#PBS -l nodes=1
+#PBS -k eo
+#PBS -j eo
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+which makewave
+which makemod
+which fmute
+which syn2d
+export OMP_NUM_THREADS=8
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
+
+tmpdir=/tmp/shotI
+mkdir -p $tmpdir
+#for dt=0.004 with modeling at 0.0005
+scale=1.0
+w1=1
+smooth=3
+#smooth=0
+
+fmute file_shot=shots/shotsmonPz1100_0_rp.su file_out=p0plus.su above=-1 shift=-10 verbose=1 check=1 hw=4
+
+suwind itmax=1023 < p0plus.su | \
+suwind key=gx min=-3000000 max=3000000 | \
+sushw key=fldr a=1 > p0plussx.su
+
+#~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen_1.su verbose=1 tap=0 ntap=0 reci=0 niter=1 hw=8 shift=7 smooth=$smooth w=1 file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0_1.su file_f1min=f1min0_1.su file_pplus=Pplus0_1.su
+#
+#~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen_4.su verbose=1 tap=0 ntap=0 reci=0 niter=4 hw=8 shift=7 smooth=$smooth w=1 file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0_4.su file_f1min=f1min0_4.su file_pplus=Pplus0_4.su
+
+#for backpropagating pplus in marchenko scheme must be written to file
+~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen.su verbose=1 tap=0 ntap=10 niter=15 hw=8 shift=7 smooth=$smooth file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0.su file_f1min=f1min0.su file_pplus=Pplus0.su
+
+exit;
+
diff --git a/marchenko/demo/model.scr b/marchenko/demo/model.scr
new file mode 100755
index 0000000..777a556
--- /dev/null
+++ b/marchenko/demo/model.scr
@@ -0,0 +1,212 @@
+#!/bin/bash
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+
+dx=2.5
+dt=0.0005
+
+makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
+ orig=-3000,0 file_base=syncl.su verbose=2 \
+ intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
+ intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
+ intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
+ intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
+ intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
+ intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
+ intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
+ intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
+ intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
+ intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
+
+ #intt=diffr x=-2000,-1000,0,1000,2000 z=1800,1800,1800,1800,1800 cp=0,0,0,0,0 ro=5000,5000,5000,5000,5000
+
+
+makemod sizex=6000 sizez=2000 dx=10 dz=5 cp0=1900 ro0=1200 \
+ orig=-3000,0 file_base=syncl_migr.su verbose=2 \
+ intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
+ intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
+ intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
+ intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
+ intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
+ intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
+ intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
+ intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
+ intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
+ intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
+
+exit
+
+#example FD modeling with model defined above
+makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0
+
+export OMP_NUM_THREADS=1
+
+makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
+
+~/bin/fdelmodc \
+ file_cp=syncl_cp.su ischeme=1 iorder=4 \
+ file_den=syncl_ro.su \
+ file_src=wave.su \
+ file_rcv=shot_fd.su \
+ src_type=1 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=4.195 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=0 \
+ ntaper=400 \
+ tsnap1=3.1 tsnap2=2.5 dtsnap=0.1 \
+ left=4 right=4 top=4 bottom=4
+
+
+exit
+
+makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
+ orig=-3000,-1000 file_base=hom.su
+
+~/bin/fdelmodc \
+ file_cp=hom_cp.su ischeme=1 iorder=4 \
+ file_den=hom_ro.su \
+ file_src=wave.su \
+ file_rcv=shot_hom_fd.su \
+ src_type=1 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=4.195 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=0 \
+ ntaper=400 \
+ left=4 right=4 top=4 bottom=4
+
+sudiff shot_fd_rp.su shot_hom_fd_rp.su > shot_rp.su
+
+migr file_shot=shot_rp.su file_vel=syncl_migr_cp.su imc=0 file_image=image.su verbose=1
+
+exit
+
+makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=$cp cs0=$cs ro0=$rho \
+ orig=-3000,-1000 file_base=synclTop.su \
+ intt=def x=-3000,0,3000 z=200,200,200 poly=0 cp=1800 ro=5000 \
+ intt=def x=-3000,-2000,-1000,-400,0,200,900,1800,3000 z=520,520,560,670,950,790,600,520,500 poly=2 cp=2300 ro=1800 \
+
+~/bin/fdelmodc \
+ file_cp=synclTop_cp.su ischeme=1 iorder=4 \
+ file_den=synclTop_ro.su \
+ file_src=wave.su \
+ file_rcv=p0.su \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=2.100 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=1100 \
+ ntaper=300 \
+ left=4 right=4 top=4 bottom=4 &
+
+~/bin/fdelmodc \
+ file_cp=synclTop_cp.su ischeme=1 iorder=4 \
+ file_den=synclTop_ro.su \
+ file_src=wave.su \
+ file_rcv=shot_top.su \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=4.195 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=0 \
+ ntaper=300 \
+ left=4 right=4 top=4 bottom=4 &
+
+makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=2300 ro0=1800 \
+ orig=-3000,-1000 file_base=synclBot.su \
+ intt=def x=-3000,0,3000 z=1310,1310,1310 poly=0 cp=2450 ro=1950 \
+ intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2460 ro=1820 \
+ intt=def x=-3000,0,3000 z=1490,1490,1570 poly=0 cp=2470 ro=2100 \
+ intt=def x=-3000,3000 z=1580,1580 poly=0 cp=2480 ro=2000 \
+ intt=def x=-3000,3000 z=1680,1680 poly=0 cp=2490 ro=1850
+
+~/bin/fdelmodc \
+ file_cp=synclBot_cp.su ischeme=1 iorder=4 \
+ file_den=synclBot_ro.su \
+ file_src=wave.su \
+ file_rcv=pRef.su \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=2.100 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=1100 zrcv2=1100 \
+ xsrc=0 zsrc=1100 \
+ ntaper=300 \
+ left=4 right=4 top=4 bottom=4 &
+
+makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=2300 ro0=1800 \
+ orig=-3000,-1000 file_base=synclBotHom.su
+
+~/bin/fdelmodc \
+ file_cp=synclBotHom_cp.su ischeme=1 iorder=4 \
+ file_den=synclBotHom_ro.su \
+ file_src=wave.su \
+ file_rcv=pRefHom.su \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=0 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.0010 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=2.100 \
+ dxrcv=10.0 \
+ xrcv1=-2250 xrcv2=2250 \
+ zrcv1=1100 zrcv2=1100 \
+ xsrc=0 zsrc=1100 \
+ ntaper=300 \
+ left=4 right=4 top=4 bottom=4 &
+
+wait
+
+sudiff shot_fd_rp.su shot_hom_fd_rp.su > shot_rp.su
+sudiff shot_top_rp.su shot_hom_fd_rp.su > shotTop_rp.su
+sudiff pRef_rp.su pRefHom_rp.su > pref_rp.su
+
diff --git a/marchenko/demo/referenceShot.scr b/marchenko/demo/referenceShot.scr
new file mode 100755
index 0000000..b1ea70e
--- /dev/null
+++ b/marchenko/demo/referenceShot.scr
@@ -0,0 +1,44 @@
+#!/bin/bash
+#PBS -N fdelmod
+#PBS -q verylong
+#PBS -l nodes=1
+#PBS -k eo
+#PBS -j eo
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
+
+#makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3
+
+dx=2.5
+dt=0.0005
+
+makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
+
+export OMP_NUM_THREADS=8
+
+$HOME/bin/fdelmodc \
+ file_cp=../syncl_cp.su ischeme=1 iorder=4 \
+ file_den=../syncl_ro.su \
+ file_src=wave.su \
+ file_rcv=virtual_shot_fd_P_zsrc1100.su \
+ src_type=1 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_ud=1 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ dtrcv=0.004 \
+ rec_delay=0.1 \
+ verbose=2 \
+ tmod=2.144 \
+ dxrcv=10.0 \
+ xrcv1=-3000 xrcv2=3000 \
+ zrcv1=0 zrcv2=0 \
+ xsrc=0 zsrc=1100 \
+ file_snap=backpropref.su tsnap1=0.1 dtsnap=0.010 tsnap2=2.100 dxsnap=10 dzsnap=10 zsnap1=0 zsnap2=2000 xsnap1=-2250 xsnap2=2250 \
+ ntaper=400 \
+ left=4 right=4 top=4 bottom=4
+
+
diff --git a/marchenko/demo/remove_direct.scr b/marchenko/demo/remove_direct.scr
new file mode 100755
index 0000000..24a3f41
--- /dev/null
+++ b/marchenko/demo/remove_direct.scr
@@ -0,0 +1,38 @@
+#!/bin/bash
+#PBS -N fdelmod
+#PBS -q verylong
+#PBS -l nodes=1
+#PBS -k eo
+#PBS -j eo
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
+
+dxshot=10
+ishot=0
+nshots=601
+
+rm shots/refl_rp.su
+
+while (( ishot < nshots ))
+do
+ (( xsrc = -3000 + ${ishot}*${dxshot} ))
+ (( iishot = ${ishot}*${dxshot}/10 ))
+ (( tr1 = 601 - ${iishot} ))
+ (( tr2 = ${tr1} + 600 ))
+ echo xsrc=$xsrc tr1=$tr1 tr2=$tr2
+ suwind < direct_rp.su key=tracl min=$tr1 max=$tr2 > direct.su
+
+ file_rcv=shots/shots_${xsrc}_rp.su
+ suwind key=tracl min=1 max=601 < $file_rcv > shotz0.su
+
+ sudiff shotz0.su direct.su > refl.su
+
+ (( ishot = $ishot + 1))
+
+ sushw < refl.su key=fldr a=$ishot | \
+ suwind itmax=1023 >> shots/refl_rp.su
+
+done
+
diff --git a/marchenko/demo/shots.scr b/marchenko/demo/shots.scr
new file mode 100755
index 0000000..907744a
--- /dev/null
+++ b/marchenko/demo/shots.scr
@@ -0,0 +1,77 @@
+#!/bin/bash
+#PBS -N fdelmod
+#PBS -q long
+#PBS -l nodes=1
+#PBS -k eo
+#PBS -j eo
+
+export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
+which makewave
+which makemod
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
+
+dt=0.0005
+makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0 scfft=1
+
+./model.scr
+
+mkdir -p shots
+mkdir -p jobs
+
+dxshot=10
+ishot=0
+nshots=601
+zsrc=0
+
+while (( ishot < nshots ))
+do
+
+ (( xsrc = -3000 + ${ishot}*${dxshot} ))
+
+ echo ishot=$ishot xsrc=$xsrc zsrc=$zsrc
+
+ cat << EOF > jobs/pbs_$ishot.job
+#!/bin/bash
+#
+#PBS -q medium
+#PBS -N mod_${xsrc}
+#PBS -j eo
+#PBS -m n
+#PBS -l nodes=1
+#PBS -V
+
+cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
+
+ export OMP_NUM_THREADS=1
+ file_rcv=shots/shots_${xsrc}.su
+
+ $HOME/bin/fdelmodc \
+ file_cp=syncl_cp.su ischeme=1 iorder=4 \
+ file_den=syncl_ro.su \
+ file_src=wavefw.su \
+ file_rcv=\$file_rcv \
+ src_type=7 \
+ src_orient=1 \
+ src_injectionrate=1 \
+ rec_type_vz=1 \
+ rec_type_p=1 \
+ rec_int_vz=2 \
+ rec_delay=0.3 \
+ dtrcv=0.004 \
+ verbose=2 \
+ tmod=4.394 \
+ dxrcv=10.0 \
+ xrcv1=-3000,-3000,-3000 xrcv2=3000,3000,3000 \
+ zrcv1=0,1000,1600 zrcv2=0,1000,1600 \
+ xsrc=$xsrc zsrc=$zsrc \
+ ntaper=400 \
+ left=4 right=4 top=4 bottom=4
+EOF
+
+qsub jobs/pbs_$ishot.job
+
+ (( ishot = $ishot + 1))
+
+done
+
--
GitLab