A Fortran Code for Null Geodesic Solutions in the Lemaitre-Tolman-Bondi Spacetime
Abstract
This paper describes the Fortran 77 code SIMU, version 1.1, designed for numerical simulations of observational relations along the past null geodesic in the Lemaitre-Tolman-Bondi (LTB) spacetime. SIMU aims at finding scale invariant solutions of the average density, but due to its full modularity it can be easily adapted to any application which requires LTB's null geodesic solutions. In version 1.1 the numerical output can be read by the GNUPLOT plotting package to produce a fully graphical output, although other plotting routines can be easily adapted. Details of the code's subroutines are discussed, and an example of its output is shown.
Cite
@article{arxiv.gr-qc/0205095,
title = {A Fortran Code for Null Geodesic Solutions in the Lemaitre-Tolman-Bondi Spacetime},
author = {Marcelo B. Ribeiro},
journal= {arXiv preprint arXiv:gr-qc/0205095},
year = {2015}
}
Comments
13 pages, 10 figures, LaTeX. Fortran code included with the LaTeX source code (also available at http://www.if.ufrj.br/~mbr/codes). Accepted for publication in "Computer Physics Communications"