This work received funding from the European Research Council (grant 742703) and the NWO Domain Applied and Engineering Sciences (grant 13939).
THE ACCOMPANYING PROGRAM IS PROVIDED UNDER THE TERMS OF THIS COMMON PUBLIC LICENSE ("AGREEMENT"). ANY USE, REPRODUCTION OR DISTRIBUTION OF THE PROGRAM CONSTITUTES RECIPIENT'S ACCEPTANCE OF THIS AGREEMENT.
A copy of this license can be found in the file 'Common_Public_License.txt' in the directory where you have found this README.
http://www.opensource.org/licenses/cpl1.0.php
Some routines are from Seismic Unix and include the SU LEGAL_STATEMENT in the source code.
Copyright (c) 2017 by the Society of Exploration Geophysicists. For more information, go to http://software.seg.org/2017/00XX . You must read and accept usage terms at: http://software.seg.org/disclaimer.txt before use.
-0- DOI reference of this software release https://zenodo.org/badge/latestdoi/23060862
-1- If the Finite Difference code has helped you in your research please refer to the following paper in your publications:
Finite-difference modeling experiments for seismic interferometry Jan Thorbecke and Deyan Draganov 2011, Geophysics, Vol. 76, no. 6 (November-December); p H1--H18, doi: 10.1190/GEO2010-0039.1 Download: https://janth.home.xs4all.nl/Publications/Articles/ThorbeckeDraganov2012.pdf
-2- If the Machenko code has helped you in your research please refer to this paper in your publications:
Implementation of the Marchenko method Jan Thorbecke, Evert Slob, Joeri Brackenhoff, Joost van der Neut, and Kees Wapenaar 2017, Geophysics, Vol. 82, no. 6 (November-December); p. WB29--WB45, doi: 10.1190/GEO2017-0108.1 Download: https://janth.home.xs4all.nl/Publications/Articles/ThorbeckeGPY2017.pdf
-3- When you are using the marchenko_primaries algorithm developed by Lele Zhang please refer to the following paper
Free-surface and internal multiple elimination in one step without adaptive subtraction Lele Zhang and Evert Slob 2019, Geophysics, Vol. 84, no. 1 (January-February); p. A7-A11, doi: 10.1190/GEO2018-0548.1 Download: http://homepage.tudelft.nl/t4n4v/BeyondInterferometry/geo_19h.pdf
-4- If you use the fdacrtmc code of Max Holicki please refer to the following paper:
Acoustic directional snapshot wavefield decomposition Holicki, M., Drijkoningen, G., and Wapenaar, K., 2019, Acoustic directional snapshot wavefield decomposition: Geophysical Prospecting, Vol. 67, 32-51 Download: http://homepage.tudelft.nl/t4n4v/4_Journals/Geophys.Prosp/GP_19a.pdf
-5-
- To compile and link the code you first have to set the ROOT variable in the Make_include file which can be found in the directory where you have found this README. You can use Make_include_template as a first start:
cp Make_include_template Make_include
-
Check the compiler and CFLAGS options in the file Make_include and adapt to the system you are using. The default options are set for a the GNU C-compiler on a Linux system. A Fortran or C++ compiler are not needed to compile the code. The Makefile has been tested with GNU make.
-
If the compiler options are set in the Make_include file you can type
make clean
make
and the Makefile will make:
- FFT library
- fdelmodc / 3D
- marchenko / 3D
- utils
The libraries will be placed in the lib/ directory and the executables in the bin/ directory.
To use the executables don't forget to include the pathname in your PATH:
bash/sh:
export PATH='path_to_this_directory'/bin:$PATH:
csh:
setenv PATH 'path_to_this_directory'/bin:$PATH:
If the compilation has finished without errors and produced an executable called bin/fdelmodc you can run one of the demo programs by running
./fdelmodc_plane.scr
in the directory fdelmodc/demo/
The demo directory contains many scripts which demonstrate the different possibilities of the modeling program.
To reproduce the Figures shown in the GEOPHYICS manuscript "Finite-difference modeling experiments for seismic interferometry" the scripts in FiguresPaper directory can be used. Please read the README in the FiguresPaper directory for more instructions and guidelines.
If the compilation has finished without errors and produced an executable called bin/marchenko you can run one of the demo programs by running a set of scripts that are explained in a README in one of the directories marchenko/demo/oneD or demo/twoD
To reproduce the Figures shown in the GEOPHYICS paper "Implementation of the Marchenko method" the scripts in marchenko/demo/oneD directory can be used. The README in this directory gives more instructions and guidelines.
To reproduce the Figures shown in the Scientific Reports paper "Virtual acoustics in inhomogeneous media with single-sided access" the scripts in marchenko/demo/ScientificReports directory can be used. The README in this directory gives more instructions and guidelines.
To reproduce the Figures shown in the GEOPHYICS paper "Implementation of the Marchenko Multiple Elimination algorithm" the scripts in marchenko/demo/oneD directory can be used. The README_PRIMARIES in this directory gives more instructions and guidelines.
The MDD kernels depend on BLAS and LAPACK calls. Free downloads of these libraries can be found on
https://www.netlib.org/blas/index.html https://www.netlib.org/lapack/index.html
If you are running on x64_86 processors you can download (for free) Intel's highly optimised MKL package:
https://software.intel.com/en-us/mkl/choose-download
These Libraries include highly optimised libraries of BLAS, LAPACK, and FFT(W). Usually the MKL libraries are installed in $MKLROOT. If that variable is not set, you can try to find the correct path by searching for one of the libraries:
find /opt/intel -name libmkl_gf_lp64.so
and adjust MKLROOT in Make_include accordingly. You can also completely disable the use of MKL by commenting out the MKL parts in Make_include.
If you want to use the .su files with SU from CWP: git clone https://github.com/JohnWStockwellJr/SeisUnix
==> Please make sure that SU is compiled without XDR (in $CWPROOT/Makefile.config make sure that XDRFLAG is NOT enabled). The SU output files of fdelmodc are all based on local IEEE data. To exclude the XDRFLAG in SU you have to use the following line in $CWPROOT/src/Makefile.config around line 35:
XDRFLAG =
so the XDRFLAG is empty. If you have made that change you have to remake SU with the commands: make remake make xtremake
The fdacrtmc and the 3D Marchenko code makes use of ZFP compression to store snaphots in CPU memory or do effient IO. This package is included in this repository for your convenience. The latest package and detailed explanation can be found on:
fdacrtmc uses FFTW and the wisdom computations are stored on disk for re-usage. This directory is defined in fdacrtmc.h
#ifndef WISDOMDIR #define WISDOMDIR "/tmp/fftw/" #endif
For the code to run properly the directory /tmp/fftw/ must exist. It is also possible to compile the fdacrtmc code that defines WISDOMDIR. In fdacrtmc/Makefile add:
CFLAGS += -DWISDOMDIR="/directory/that/exists"
or you can change the name of WISDOMDIR in fdacrtmc.h
Other make commands which can be useful:
make clean : removes all object files, but leaves libraries and executables make realclean: removes also object files, libraries and executables.
If you encounter missing trigometric / mathematical functions during compilation, for example;
defineSource.c:(.text+0x2356): undefined reference to sin getParameters.o: In function getParameters:
add '-lm -lc' around line 109 in Make_include:
LIBS= -L$L -lgenfft $(BLAS) -lm -lc
The latest version of the source code and manual can be found at:
http://www.xs4all.nl/~janth/Software/Software.html
or at github:
git clone https://github.com/JanThorbecke/OpenSource.git git clone git://github.com/JanThorbecke/OpenSource.git
The code is used by many different people and if there is a request for a new option in the code, then I will try to implement, test and make it available.