prepare_solution

Calculates helper columns from an astronomical solution input.

Easily calculate precession and obliquity from an orbital solution (defaults to ZB18a from Zeebe and Lourens (2019) <doi:10.1126/science.aax0612>) and assumed or reconstructed values for tidal dissipation (Td) and dynamical ellipticity (Ed). This is a translation and adaptation of the 'C'-code in the supplementary material to Zeebe and Lourens (2022) <doi:10.1029/2021PA004349>, with further details on the methodology described in Zeebe (2022) <doi:10.3847/1538-3881/ac80f8>. The name of the 'C'-routine is 'snvec', which refers to the key units of computation: spin vector s and orbit normal vector n.

Ilja Kocken

snvecR

Calculate Earth’s Obliquity and Precession in the Past

Richard Zeebe

prepare_solution function

<dl><dt>data</dt>
<dd>A data frame with the following columns:<ul>
<li><code>t</code> Time \(t\) (days).</li>
<li><code>ee</code> Eccentricity \(e\) (unitless).</li>
<li><code>lph</code> Longitude of perihelion \(\varpi\) (degrees).</li>
<li><code>lan</code> Longitude of the ascending node \(\Omega\) (degrees).</li>
<li><code>inc</code> Inclination \(I\) (degrees).</li>
</ul></dd>The easiest way to get this is with <code>get_solution()</code>.
<dt>quiet</dt>
<dd>Be quiet?<ul>
<li>If <code>TRUE</code>, hide info messages.</li>
<li>If <code>FALSE</code> (the default) print info messages and timing.</li>
</ul></dd></dl>

Arguments

Prepare Astronomical Solution — prepare_solution

<dl>

<dt>data</dt>
<dd>A data frame with the following columns:<ul>
<li><code>t</code> Time \(t\) (days).</li>
<li><code>ee</code> Eccentricity \(e\) (unitless).</li>
<li><code>lph</code> Longitude of perihelion \(\varpi\) (degrees).</li>
<li><code>lan</code> Longitude of the ascending node \(\Omega\) (degrees).</li>
<li><code>inc</code> Inclination \(I\) (degrees).</li>
</ul>

The easiest way to get this is with <code>get_solution()</code>.</dd>


<dt>quiet</dt>
<dd>Be quiet?<ul>
<li>If <code>TRUE</code>, hide info messages.</li>
<li>If <code>FALSE</code> (the default) print info messages and timing.</li>
</ul></dd>

</dl>

prepare_solution: Prepare Astronomical Solution

Description

Usage

Value

Arguments

Details

See Also