Reproducible paper explaining kirchoff modeling #13

This SConstruct generates a reproducible paper automatically after
running the numerical experiment of kirchoff modeling available in
the ./scons directory.

It uses references available in mybib.bib file and paper.tex is the
Latex file template.

More information about how to generate reproducicle papers and dependencies
are available in the Madagascar website:

http://www.ahay.org/wiki/Reproducible_computational_experiments_using_SCons

Author: Dirack

usage_examples/kirchoff_modeling/SConstruct

@@ -0,0 +1,26 @@
+from rsf.tex import *
+import os
+import glob
+
+# Run the experiment with scons
+print('Running experiment with scons...')
+os.chdir('scons/')
+os.system('scons')
+
+# Build figures
+print('Building figures...')
+os.chdir('Fig/')
+listOfFigures = glob.glob('*.vpl')
+
+for i in listOfFigures:
+	i = i.split('.')[0]
+	print('Building '+str(i)+'.pdf')
+	string = 'jpegpen <'+str(i)+'.vpl bgcolor=w >'+str(i)+'.jp2'
+	os.system(string)
+	string = 'img2pdf <'+str(i)+'.jp2 >'+str(i)+'.pdf'
+	os.system(string)
+
+print('Building reproducible papper...')
+os.chdir('../..')
+
+End(options='manuscript',use=['[bookmarks=true]hyperref','amsmath','color','listings'])

usage_examples/kirchoff_modeling/mybib.bib

@@ -0,0 +1,8 @@
+@Article{km,
+  author = 	 {B. Duquet and K. Marfut and J. Dellinger},
+  title = 	 {Kirchhoff modeling, inversion forreflectivity, and subsurface illumination},
+  journal = 	 {Geophysics},
+  year = 	 2000,
+  volume = 	 65,
+  number = 	 4,
+  pages = 	 {1195-1209}}

usage_examples/kirchoff_modeling/paper.tex

@@ -0,0 +1,54 @@
+\title{Reproducible experiment: Kirchhoff modeling}
+
+\author{Rodolfo A C Neves (Dirack)}
+\begin{abstract}
+To show how to implement Kirchoff modeling in Madagascar
+package. This paper is part of the current directory 
+documentation followed by a SConstruct script
+to run the numerical experiment.
+\end{abstract}
+
+\section{Introduction}
+
+\indent In this directory we build a seismic data cube from a gaussian reflector in a linear velocity model.
+\footnote{This numerical experiment was originaly 
+design by Sergey Fomel and Roman Kazinnik in 
+the paper Non-hyperbolic common reflection surface
+(2013). 
+This paper is avaliable as a reproducible paper in
+the \href{http://www.ahay.org}{Madagascar website}.}
+Velocity increases with depth with a 0.5 
+velocity gradient and near surface velocity 
+$v_0=1.5Km/s$.
+
+\indent A seismic data cube means that each sample
+in the data space can be identifyed by it's CMP, 
+half-offset and time coordinates. So, the seismic 
+data can be represented by an amplitude function
+A(m,h,t), where m is the CMP, h the half-offset
+and t is the time.
+
+Kirchhoff modeling \cite[]{km}:
+
+\begin{quote}
+Kirchhoff modeling can be defined as the 
+mathematical transpose of Kirchhoff migration.
+The resulting Kirchhoff modeling algorithm 
+has the same low computational cost as 
+Kirchhoff migration and, unlike expensive full 
+acoustic or elastic wave‐equation methods, 
+only models the events that Kirchhoff migration
+can image.
+\end{quote}
+
+\inputdir{scons}
+\plot{gaussianReflectorVelocityModel}{width=\textwidth}{Gaussian reflector model.}
+
+\plot{dataCube}{width=\textwidth}{Seismic data cube.}
+
+\plot{cmpGather}{width=\textwidth}{CMP gather 
+extracted from the seismic data cube, 
+5Km is the CMP Gather coordinate.}
+
+\bibliographystyle{seg}
+\bibliography{mybib}