From 2f6835f02bd3c56dae7344723d161b884de16823 Mon Sep 17 00:00:00 2001
From: JBrackenhoff <J.A.Brackenhoff@tudelft.nl>
Date: Wed, 29 May 2019 15:29:56 +0200
Subject: [PATCH] WS15 Joeri
---
marchenko/demo/WS15/README.5 | 19 ++++++++++++
marchenko/demo/twoD/homg_reference.scr | 43 ++++++++++++++++++++++++++
marchenko/demo/twoD/homgpng.scr | 6 ++++
marchenko/demo/twoD/homgview.scr | 4 +++
marchenko/demo/twoD/marchenko_ray.scr | 3 ++
5 files changed, 75 insertions(+)
create mode 100644 marchenko/demo/WS15/README.5
create mode 100755 marchenko/demo/twoD/homg_reference.scr
create mode 100755 marchenko/demo/twoD/homgpng.scr
create mode 100755 marchenko/demo/twoD/homgview.scr
diff --git a/marchenko/demo/WS15/README.5 b/marchenko/demo/WS15/README.5
new file mode 100644
index 0000000..cec5287
--- /dev/null
+++ b/marchenko/demo/WS15/README.5
@@ -0,0 +1,19 @@
+To generate the raytimes and the estimated amplitudes for the virtual source and receivers, run:
+
+rayvsp.scr
+
+After the raytimes have been determined you can create the source function of the homogeneous Green's function by running:
+
+marchenko_ray.scr
+
+A reference Green's function can be modeled by running
+
+homg_reference.scr
+
+To view the retrieved and reference Green's function, run:
+
+homgview.scr
+
+Or to convert the data to eps and png, run:
+
+homgpng.scr
diff --git a/marchenko/demo/twoD/homg_reference.scr b/marchenko/demo/twoD/homg_reference.scr
new file mode 100755
index 0000000..30d5678
--- /dev/null
+++ b/marchenko/demo/twoD/homg_reference.scr
@@ -0,0 +1,43 @@
+#!/bin/bash
+
+export PATH=$HOME/src/OpenSource/bin:$PATH:
+
+dx=2.5
+dt=0.0005
+
+makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
+
+
+export KMP_AFFINITY=disabled
+export OMP_NUM_THREADS=16
+
+
+
+fdelmodc \
+ file_cp=syncl_cp.su ischeme=1 iorder=4 \
+ file_den=syncl_ro.su \
+ file_src=wave.su \
+ file_rcv=homg_ref.su \
+ 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=0 zsrc=1000 \
+ file_snap=homg_ref.su \
+ xsnap1=0 xsnap2=0 dxsnap=10 \
+ zsnap1=100 zsnap2=1900 dzsnap=10 \
+ tsnap=0.1 tsnap2=2.044 dtsnap=0.004 \
+ ntaper=200 npml=100 \
+ left=2 right=2 top=2 bottom=2
+
+sushw < homg_ref_sp.su key=ntr a=0 > tmp.su
+mv tmp.su homg_ref_sp.su
diff --git a/marchenko/demo/twoD/homgpng.scr b/marchenko/demo/twoD/homgpng.scr
new file mode 100755
index 0000000..dfd55d0
--- /dev/null
+++ b/marchenko/demo/twoD/homgpng.scr
@@ -0,0 +1,6 @@
+#!/bin/bash -x
+
+suwind < homg_rayz.su key=fldr min=513 | supsimage title='Retrieved' perc=99 > homg_retrieved.eps
+/usr/bin/convert -quality 90 -antialias homg_retrieved.eps homg_retrieved.png
+supsimage < homg_ref_sp.su title='Reference' perc=99 > homg_reference.eps
+/usr/bin/convert -quality 90 -antialias homg_reference.eps homg_reference.png
diff --git a/marchenko/demo/twoD/homgview.scr b/marchenko/demo/twoD/homgview.scr
new file mode 100755
index 0000000..d5efb59
--- /dev/null
+++ b/marchenko/demo/twoD/homgview.scr
@@ -0,0 +1,4 @@
+#!/bin/bash -x
+
+suwind < homg_rayz.su key=fldr min=513 | suximage title='Retrieved' perc=99 &
+suximage < homg_ref_sp.su title='Reference' perc=99 &
diff --git a/marchenko/demo/twoD/marchenko_ray.scr b/marchenko/demo/twoD/marchenko_ray.scr
index f6728fd..ad12049 100755
--- a/marchenko/demo/twoD/marchenko_ray.scr
+++ b/marchenko/demo/twoD/marchenko_ray.scr
@@ -17,3 +17,6 @@ marchenko3D file_shot=shots/refl_rp.su verbose=1 file_inp=greensrc.su \
tap=0 niter=6 hw=8 shift=7 smooth=3 geomspread=1 ampest=1 \
file_ray=rayz_time.su file_amp=rayz_amp.su file_wav=wavefp.su \
file_homg=homg_rayz.su file_green=green_rayz.su
+
+sushw < homg_rayz.su key=ntr a=0 > tmp.su
+mv tmp.su homg_rayz.su
--
GitLab