Commit b23d519c authored by Roland Haas's avatar Roland Haas
Browse files

power: add some example parfile fragment

parent f13e117d
......@@ -365,6 +365,83 @@ explicitly or read from output files produced by
TwoPunctures~\cite{POWER-TwoPunctures:web} and
QuasiLocalMeasures~\cite{POWER-QuasiLocalMeasures:web}.
A typical parameter file fragment to produce all used output files looks like
this (taken from the GW150914 gallery
example~\cite{POWER-wardell_barry_2016_155394}:
\begin{verbatim}
ActiveThorns = "
Multipole
QuasiLocalMeasures
TwoPunctures
WeylScal4
"
################################################################
# Initial data
################################################################
ADMBase::initial_data = "twopunctures"
ADMBase::initial_lapse = "twopunctures-averaged"
################################################################
# Psi4 computation
################################################################
WeylScal4::fdOrder = 8
WeylScal4::calc_scalars = "psis"
WeylScal4::calc_invariants = "always"
################################################################
# Psi4 mode decomposition
################################################################
# Radii are chosen to be evenly spaced in 1/r as that is the
# variable extrapolation is performed in
Multipole::nradii = 7
Multipole::radius[0] = 100 # c
Multipole::radius[1] = 115
Multipole::radius[2] = 136
Multipole::radius[3] = 167
Multipole::radius[4] = 214
Multipole::radius[5] = 300
Multipole::radius[6] = 500
Multipole::ntheta = 120
Multipole::nphi = 240
Multipole::variables = "
WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}
"
# choose a number for $wave_extract_every here
Multipole::out_every = $wave_extract_every
Multipole::l_max = 8
Multipole::output_hdf5 = yes
# Disable ASCII output to avoid creating a large number of
# files
Multipole::output_ascii = no
################################################################
# Isolated Horizons (used only by Nakano method)
################################################################
QuasiLocalMeasures::verbose = no
QuasiLocalMeasures::veryverbose = no
QuasiLocalMeasures::interpolator =
"Lagrange polynomial interpolation"
QuasiLocalMeasures::interpolator_options = "order=4"
QuasiLocalMeasures::spatial_order = 4
QuasiLocalMeasures::num_surfaces = 3
QuasiLocalMeasures::surface_index [0] = 2
QuasiLocalMeasures::surface_index [1] = 3
QuasiLocalMeasures::surface_index [2] = 4
\end{verbatim}
Other parameter files and thorns can be used as long as required values for
ADM mass (both methods), and final mass and spin of the resulting object
(\texttt{Nakano} method only) are provided using command line options or
function call arguments.
\subsection{Support and Feedback}
Please use the Einstein Toolkit issue system
......@@ -454,6 +531,13 @@ Phys. Rev. Lett. \textbf{108}, 131101 (2012)
doi:10.1103/PhysRevLett.108.131101
[arXiv:1110.2938 [gr-qc]].
\bibitem{POWER-wardell_barry_2016_155394}
B.~Wardell, I.~Hinder and E.~Bentivegna,
``Simulation of GW150914 binary black hole merger using the Einstein Toolkit,''
doi:10.5281/zenodo.155394
\url{https://doi.org/10.5281/zenodo.155394}
\end{thebibliography}
% Do not delete next line
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment