@essrt/physics 22.1.0 → 23.2.0

package/src/layer0/index.js

@@ -32,9 +32,7 @@
 /**
  * @typedef {Object} EpochParams
  * @property {number} epochYear                  year at t_Ma = 0 (2000, NOT startmodelYear)
- * @property {number} alpha1PerMa                Moon recession, LLR-anchored
- * @property {number} alpha3PerMa3
- * @property {number} alpha4PerMa4
+ * @property {number} alpha1PerMa                Moon recession, LLR-anchored (the J2000 slope of the cited history)
  * @property {number} moonDistanceNowM           a_Moon at J2000, metres
  * @property {number} moonLockDistanceM          tidal-lock asymptote, metres
  * @property {number} totalAngularMomentumKgM2S  L_tot of the Earth-Moon system
@@ -53,17 +51,24 @@
 /**
  * @typedef {Object} EpochPrimitives
  * @property {(year: number) => number} tMa
- * @property {(year: number) => number} moonDistanceMetres
+ * @property {(year: number) => number|null} moonDistanceMetres - null beyond the cited history's span
  * @property {(year: number) => number|null} lodSeconds
  * @property {(year: number) => number} siderealYearSeconds
- * @property {(year: number) => number} tropicalYearSeconds      sidereal less one equinox turn per T_p(t) (the composed precession period)
+ * @property {(year: number) => number} tropicalYearSeconds      sidereal less one equinox turn per T_p(t) (the composed precession period on the RATE of date)
  */
 
 /**
- * @param {{params: EpochParams, alphaAtAgeMa: (tMa: number) => number}} deps
+ * @param {{params: EpochParams, alphaAtAgeMa: (tMa: number) => number,
+ *   moonDistanceMetresAtAge: (tMa: number) => number|null,
+ *   meanObliquityCosRatioAtAge?: (tMa: number) => number|null}} deps —
+ *   moonDistanceMetresAtAge is the engine's cited distance history
+ *   (deltat/moon-distance-history: Farhat et al. 2022 as a shape on the
+ *   J2000 anchor; plan 12 option A — the layer0 quartic twin is deleted);
+ *   meanObliquityCosRatioAtAge the cited mean obliquity's cos ε̄(t)/cos ε̄(0)
+ *   (earth/mean-obliquity-history, plan 12 §8 option (a)) — OPTIONAL, absent → 1.
  * @returns {EpochPrimitives}
  */
-export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
+export const createEpochPrimitives = ({ params: p, alphaAtAgeMa, moonDistanceMetresAtAge, meanObliquityCosRatioAtAge }) => {
   /**
    * Years before the epoch, in Ma. The ONE place this conversion lives.
    * @param {number} year
@@ -72,24 +77,29 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
   const tMa = (year) => (p.epochYear - year) / 1e6;
 
   /**
-   * Moon semi-major axis. Cubic in t_Ma: alpha1 is LLR-anchored, alpha3/alpha4
-   * are the Farhat 2022 deep-time fit.
+   * Moon distance at age t — the cited history (null beyond its span).
    * @param {number} t age in Ma
-   * @returns {number} metres
+   * @returns {number|null} metres
    */
-  const moonDistanceMetresCore = (t) =>
-    p.moonDistanceNowM
-      * (1 + p.alpha1PerMa * t + p.alpha3PerMa3 * t * t * t + p.alpha4PerMa4 * t * t * t * t);
+  const moonDistanceMetresCore = (t) => moonDistanceMetresAtAge(t);
+
+  /**
+   * cos ε̄(t)/cos ε̄(0) of the cited mean obliquity — the precession RATE's
+   * factor on the constant (plan 12 §8); 1 when no history is injected.
+   * @param {number} t age in Ma
+   * @returns {number|null}
+   */
+  const meanObliquityCosRatioCore = (t) => (meanObliquityCosRatioAtAge ? meanObliquityCosRatioAtAge(t) : 1);
 
   /**
    * Length of day from angular-momentum conservation: as the Moon recedes it
    * takes orbital angular momentum and Earth's spin slows to compensate.
    * @param {number} t age in Ma
-   * @returns {number|null} seconds; null past the tidal-lock asymptote
+   * @returns {number|null} seconds; null past the tidal-lock asymptote or the history's span
    */
   const lodSecondsCore = (t) => {
     const a = moonDistanceMetresCore(t);
-    if (a <= 0 || a >= p.moonLockDistanceM) return null;
+    if (a === null || a <= 0 || a >= p.moonLockDistanceM) return null;
     const iEarth = alphaAtAgeMa(t) * p.earthMassKg * p.earthRadiusM * p.earthRadiusM;
     return (2 * Math.PI * iEarth)
       / (p.totalAngularMomentumKgM2S
@@ -123,7 +133,9 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
    * Tropical year of date — the sidereal year less the equinox's own motion
    * over one year: one turn per the composed lunisolar precession period
    * T_p(t), i.e. 1/T_trop = 1/T_sid + ψ̇(t)/(1,296,000″ · Julian-year s) with
-   * ψ̇(t) = [ω₀/ω(t)]⁻¹… spelled as (LOD₀/LOD(t))·p₀·[f_S + (1 − f_S)(a₀/a_M)³].
+   * ψ̇(t) the RATE of date, spelled as (LOD₀/LOD(t))·p₀·[f_S + (1 − f_S)(a₀/a_M)³]
+   * × cos ε̄(t)/cos ε̄(0) — the cited mean obliquity's factor (plan 12 §8 option (a);
+   * injected, 1 when absent).
    * Plan 09 step 5d: it was the unit's calendar convention sid·(1 − 13/H(t))
    * (plan 06 D2 kept it as a named device), 0.086 % slow on the equinox term
    * — 1.0 s of year at J2000, integrated by the deep JD↔year calendar. The
@@ -143,7 +155,8 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
     const lf = (a === null || !(a > 0)) ? null : Math.pow(p.moonDistanceNowM / a, 3);
     const term = lf === null ? null : p.precessionSolarShareJ2000 + (1 - p.precessionSolarShareJ2000) * lf;
     const p0 = p.precessionRateJ2000ArcsecPerYr;
-    const rate = (lod === null || !(lod > 0) || term === null) ? p0 : (p.lodNowH13Seconds / lod) * p0 * term;
+    const f = meanObliquityCosRatioCore(t);
+    const rate = (lod === null || !(lod > 0) || term === null || f === null) ? p0 : (p.lodNowH13Seconds / lod) * p0 * term * f;
     return 1 / (1 / sid + rate / (1296000 * 365.25 * 86400));
   };
 

package/src/math/angles.cjs

@@ -0,0 +1,38 @@
+/**
+ * Angles — ONE home for the degree ↔ radian constants and the wraps
+ * (holisticuniverse plan 13 B1). Before B1 the browser carried 27 local
+ * `D2R`-style constants and twelve local wrap helpers in three ±180 forms;
+ * the package itself a dozen private `D2R`s.
+ *
+ * The two ±180 forms in use differ only AT the boundary: `wrap180` maps
+ * +180 to −180 (the half-open [−180, 180) the longitude-difference sites
+ * use), `wrap180Closed` keeps +180 (the (−180, 180] of the residual
+ * printers). A site keeps the form it had — B1 moves the formula, never
+ * the convention.
+ *
+ * Pure, no model dependencies.
+ */
+
+'use strict';
+
+/** π/180 — degrees to radians. */
+const DEG2RAD = Math.PI / 180;
+/** 180/π — radians to degrees. */
+const RAD2DEG = 180 / Math.PI;
+/** Arcseconds in a full turn. */
+const ARCSEC_PER_TURN = 1296000;
+/** Arcseconds per radian, 206264.806…  */
+const ARCSEC_PER_RAD = (180 / Math.PI) * 3600;
+
+/** @param {number} x degrees @returns {number} the same angle in [0, 360) */
+const wrap360 = (x) => ((x % 360) + 360) % 360;
+/** @param {number} x degrees @returns {number} the same angle in [−180, 180) */
+const wrap180 = (x) => ((x + 540) % 360 + 360) % 360 - 180;
+/** @param {number} x degrees @returns {number} the same angle in (−180, 180] */
+const wrap180Closed = (x) => { const w = ((x % 360) + 360) % 360; return w > 180 ? w - 360 : w; };
+/** @param {number} x radians @returns {number} the same angle in [0, 2π) */
+const wrap2pi = (x) => { const tau = 2 * Math.PI; return ((x % tau) + tau) % tau; };
+/** @param {number} x radians @returns {number} the same angle in [−π, π) */
+const wrapPi = (x) => { const tau = 2 * Math.PI; return ((x + Math.PI) % tau + tau) % tau - Math.PI; };
+
+module.exports = { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wrap180Closed, wrap2pi, wrapPi };