@essrt/physics 23.3.3 → 24.0.0

package/src/deltat/deep-time.cjs

@@ -1,11 +1,10 @@
 /**
  * Deep-time LOD / ΔT core — THE shared implementation (Phase 8.4, slice 3).
  *
- * The angular-momentum-conservation LOD (Layer 1/2), the H(t) identity, the
- * mass-loss year chain, the Simpson ΔT integrator with the H/5 ecliptic
- * "missing motion" integrand, the corrected-LOD composite, the Actual-LOD
- * Fourier ripple, and the dLOD/dt driver decomposition (tidal / GIA /
- * stack / resonator channels).
+ * The angular-momentum-conservation LOD (Layer 1/2), the composed precession
+ * rate, the mass-loss year chain, the Simpson ΔT integrator with the ecliptic
+ * "missing motion" integrand, the corrected-LOD composite, and the dLOD/dt
+ * driver decomposition (tidal / GIA / stack / resonator channels).
  *
  * INJECTED, by design:
  *  - moonDistanceMetresAtAge — the 8.2 moon month-chain (each engine's own
@@ -44,11 +43,11 @@ const { createComposedPrecession } = require('../earth/precession-composed.cjs')
  * @property {number} aLockMetres - tidal-lock distance (LOD denominator → 0)
  * @property {number} aMoonNowMetres - Farhat polynomial anchor
  * @property {number} alpha1PerMa - the LLR recession rate (the J2000 slope of the cited history)
- * @property {number} lodNowH13Seconds - ω₀, the J2000 day length (the SI day since plan 09; identifier kept)
+ * @property {number} lodNowH13Seconds - ω₀, the J2000 day length (the SI day; identifier kept)
  * @property {number} meanSiderealYearJ2000Seconds
  * @property {number} solarMassLossFracPerYear
  * @property {number} precessionSolarShareJ2000 - f_S, the solar fraction of the J2000
- *   precession torque (derive-params) — the unit H(t) = 13·T_p,composed needs it
+ *   precession torque (derive-params) — the composed precession rate needs it
  */
 
 /**
@@ -59,23 +58,20 @@ const { createComposedPrecession } = require('../earth/precession-composed.cjs')
  * @property {(year: number) => number} cycleLodSumAt - gated δLOD sum (incl. swing)
  * @property {(year: number) => number} swingLodAt - gated swing δLOD alone
  * @property {(year: number) => number} swingLodRateAt - gated analytic swing rate
- * @property {() => number} precessionPeriodJ2000YearsFn - the model's DERIVED
- *   J2000 axial precession period: the certified year-length laws' beat at
- *   2000, sid/(sid − trop) ≈ 25,771.4 yr (plan 06 S5 — one J2000 precession
- *   reading; the same anchor the hybrid obliquity self-anchors on). Read
+ * @property {() => number} precessionPeriodJ2000YearsFn - the model's J2000 axial
+ *   precession period, 1,296,000/p_A = 25,771.575 yr (the IAU 2006 general
+ *   precession; the same anchor the hybrid obliquity self-anchors on). Read
  *   LAZILY on the first composed-rate use, never at construction: every
  *   runtime's year laws read THIS factory's bases.
  * @property {() => number} nodalPeriodYearsFn - the period of Earth's orbit
  *   plane's nodal regression on the invariable plane, 1,296,000/|s₃| with s₃
  *   the dominant ζ mode of the banked deep secular modes (≈ 68,751 yr). An
  *   ORBITAL quantity (μ-tier: it does not follow Earth's spin), the reference
- *   the solar day's "ecliptic missing motion" term rides (plan 06 T2 item —
- *   formerly the spin unit's H/5, 67,063 yr, 2.5 % off and scaling with the
- *   spin at deep time). Lazy, like the precession anchor.
+ *   the solar day's "ecliptic missing motion" term rides. Lazy, like the
+ *   precession anchor.
  * @property {(tMa: number) => number} [lEmAtAgeKgm2S] - OPTIONAL time-dependent
  *   Earth-Moon angular momentum. Absent → the J2000 constant: conservation
- *   (plan 12 option A — the Driver-1½ solar-channel budget is deleted;
- *   measured: conservation on the cited history reproduces Farhat's own LOD
+ *   (measured: conservation on the cited history reproduces Farhat's own LOD
  *   to 0.01 h at 100 Ma, 0.14 h at 2.46 Ga).
  * @property {(tMa: number) => number|null} [moonDistanceSlopeMetresPerMaAtAge] -
  *   OPTIONAL analytic slope of the distance history (metres per Ma, positive
@@ -84,7 +80,7 @@ const { createComposedPrecession } = require('../earth/precession-composed.cjs')
  * @property {(tMa: number) => number|null} [meanObliquityCosRatioAtAge] -
  *   OPTIONAL cos ε̄(t)/cos ε̄(0) of the cited mean obliquity (earth/mean-
  *   obliquity-history) — the factor the precession RATE of date carries on the
- *   constant (plan 12 §8 option (a)); absent → 1, the constant alone, bit-identical.
+ *   constant; absent → 1, the constant alone, bit-identical.
  */
 
 /** @param {DeepTimeLodDeps} deps */
@@ -117,26 +113,12 @@ function createDeepTimeLod(deps) {
     return (s === null) ? null : s / 3600;
   }
 
-  // (Plan 09 step 7a deleted the FROZEN ERA CLOCK's public counter
-  // `eraClockHAtAge` — H_era(t) = H₀·LOD(t)/LOD₀, plan 06 D8's second named
-  // counter — with its year twins; step 7b deleted the unit itself, `hAtAge`
-  // = H_era/[f_S + (1 − f_S)(a₀/a_M)³] — H(t) was 13.011 composed precession
-  // periods, a fit constant's multiple of T_p(t), never a period of anything
-  // the model computes. The physics it scaled is the composed precession
-  // below: ψ̇(t) and T_p(t) = 1,296,000/ψ̇ are the published quantities.)
-  // The `…FromLod` twins take an ALREADY-COMPUTED LOD(t): the Simpson ΔT
-  // integrand below evaluated the LOD four times and the tropical year twice
-  // per node through these helpers (measured 8.8 → 4.8 ms per cardinal panel
-  // year with the model's ΔT memo). Same operations, same inputs — every
-  // value bit-identical (11,958-value sweep, 0 differences).
-
-  // The composed lunisolar precession (plan 06 D6 → Phase 3 → S5; plan 12 §8):
-  // the CONSTANT's scaling [ω(t)/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M(t))³] and the
-  // RATE of date, that × cos ε̄(t)/cos ε̄(0) on the cited mean obliquity; p₀ the
-  // model's J2000 rate 1,296,000/axial0 (the IAU 2006 general precession since
-  // plan 09). ONE formula home: earth/precession-composed, built here on this
-  // factory's own LOD, recession and obliquity histories; the anchor resolves
-  // lazily on the first rate use.
+  // The composed lunisolar precession: the CONSTANT's scaling
+  // [ω(t)/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M(t))³] and the RATE of date, that ×
+  // cos ε̄(t)/cos ε̄(0) on the cited mean obliquity; p₀ the model's J2000 rate
+  // 1,296,000/axial0 (the IAU 2006 general precession). ONE formula home:
+  // earth/precession-composed, built here on this factory's own LOD, recession
+  // and obliquity histories; the anchor resolves lazily on the first rate use.
   const composed = createComposedPrecession({
     p0ArcsecPerYr: () => 1296000 / deps.precessionPeriodJ2000YearsFn(),
     solarShare: K.precessionSolarShareJ2000,
@@ -148,11 +130,10 @@ function createDeepTimeLod(deps) {
     meanObliquityCosRatioAtAge: deps.meanObliquityCosRatioAtAge,
   });
 
-
   /** THE PRECESSION RATE OF DATE ψ̇(t), ″/yr — the constant's scaling
    * p₀·[ω/ω₀]·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ε̄ the cited mean
-   * obliquity (plan 12 §8 option (a); the chain held cos ε̄ at J2000 before —
-   * 2.5 % of rate at 2.46 Ga, 2e-7 at 10 kyr). p₀ the J2000 anchor (S5).
+   * obliquity (its factor is 2.5 % of the rate at 2.46 Ga, 2e-7 at 10 kyr);
+   * p₀ the J2000 anchor.
    * @param {number} t_Ma @returns {number|null} */
   function lunisolarPrecessionRateArcsecPerYrAtAge(t_Ma) {
     return composed.composedRateArcsecPerYrAtAge(t_Ma);
@@ -192,23 +173,18 @@ function createDeepTimeLod(deps) {
     // a(t) = a₀·(1−Δm)): T ∝ M⁻² ⇒ T(t) = T₀·(1−Δm)² — the SAME Driver-2
     // law the planet chains use (orbit-chain driver2PeriodSecondsAtAge).
     // The previous (1 − 2Δm) was this law's first-order Taylor — the one
-    // form in the model inconsistent with it (measured on the retired day-count
-    // identity of the era counter: with the product form the counter's
-    // H_era·(days/yr)·(AU₀/AU)² held exactly — a tautology of spin scaling ∘
-    // Kepler, deleted with the counter and the unit at plan 09 steps 7a/7b).
+    // form in the model inconsistent with it.
     return K.meanSiderealYearJ2000Seconds * (1 - mass_loss_fraction) * (1 - mass_loss_fraction);
   }
 
   /** Tropical year of date, seconds — the sidereal year less the equinox's
    * own motion over one year: one turn per the composed lunisolar precession
    * period T_p(t), 1/T_trop = 1/T_sid + ψ̇(t)/(1,296,000″ · Julian-year s),
-   * ψ̇(t) the RATE of date (with the cited mean obliquity's cosine, plan 12 §8).
-   * Plan 09 step 5d: it was the unit's calendar convention sid·(1 − 13/H(t))
-   * — plan 06 D2 had kept it as a named device — 0.086 % slow on the equinox
-   * term (1.0 s of year at J2000), and the deep JD↔year calendar integrates
-   * THIS. The SAME operations as layer0's tropicalYearSecondsCore (the
-   * layer0 gate holds the twins bit-identical). Past the tidal chain's domain
-   * the J2000 rate stands in. @param {number} t_Ma @returns {number} */
+   * ψ̇(t) the RATE of date (with the cited mean obliquity's cosine). The deep
+   * JD↔year calendar integrates THIS. The SAME operations as layer0's
+   * tropicalYearSecondsCore (the layer0 gate holds the twins bit-identical).
+   * Past the tidal chain's domain the J2000 rate stands in.
+   * @param {number} t_Ma @returns {number} */
   function tropicalYearSecondsAtAge(t_Ma) {
     return tropicalYearSecondsFromLod(t_Ma, lodSecondsAtAge(t_Ma));
   }
@@ -235,10 +211,6 @@ function createDeepTimeLod(deps) {
     return tropicalYearSecondsAtAge(t_Ma) / LOD_s;
   }
 
-  // (Plan 09 step 7a: the frozen era clock's year twins on H_era —
-  // T_sid·(1 − 13/H_era) and its day count, the bases the deleted cardinal
-  // record and year-length combs were fitted against — are deleted.)
-
   /**
    * The solar day's "ecliptic missing motion" correction, seconds: the mean
    * solar day is counted against the Sun on the ECLIPTIC, whose plane turns
@@ -247,8 +219,7 @@ function createDeepTimeLod(deps) {
    * (≈ 3.44 ms at J2000). T_s₃ is an orbital quantity: constant at every
    * epoch (μ-tier), so the term rides only LOD_mean and the year length at
    * deep time. ONE home for the three runtimes' `h5Correction` twins (the
-   * identifier keeps its name; the divisor is the nodal period since the
-   * plan 06 T2 restatement, formerly the spin unit's H/5).
+   * identifier keeps its historical name; the divisor is the nodal period).
    * @param {number} t_Ma @returns {number|null} seconds (null past tidal lock)
    */
   function eclipticLodCorrectionSecondsAtAge(t_Ma) {
@@ -305,10 +276,6 @@ function createDeepTimeLod(deps) {
     return tidal + deps.cycleLodSumAt(year);
   }
 
-  // (Plan 09: the "Actual LOD" — the tidal mean × the retired sidereal-year
-  // comb's ripple — is deleted with the comb; the ripple was the sidereal
-  // year's λ̇ variation misattributed to the day length.)
-
   /**
    * dLOD/dt driver decomposition, ms/century per channel.
    * tidal: Moon-recession (Farhat da/dt, LLR-anchored) · gia: α(t) numeric
@@ -331,12 +298,11 @@ function createDeepTimeLod(deps) {
     const I_E = alpha * K.mEarthAloneKg * K.rEarthMetres * K.rEarthMetres;
     const year = 2000 - t_Ma * 1e6;
 
-    // The cited history's slope (plan 12 option A): the engine's
-    // moonDistanceSlopeMetresPerMaAtAge when it provides one (the Hermite's
-    // analytic derivative — C1, smooth), else a ±1-kyr central difference on
-    // the distance itself. t_Ma > 0 = past; forward-time sign flip. At J2000
-    // the slope IS the LLR observation (the history's boundary condition):
-    // +3.82 cm/yr ✓
+    // The cited history's slope: the engine's moonDistanceSlopeMetresPerMaAtAge
+    // when it provides one (the Hermite's analytic derivative — C1, smooth),
+    // else a ±1-kyr central difference on the distance itself. t_Ma > 0 =
+    // past; forward-time sign flip. At J2000 the slope IS the LLR observation
+    // (the history's boundary condition): +3.82 cm/yr ✓
     let slopePerMa = deps.moonDistanceSlopeMetresPerMaAtAge ? deps.moonDistanceSlopeMetresPerMaAtAge(t_Ma) : null;
     if (slopePerMa === null || slopePerMa === undefined) {
       const H_MA = 1e-3;
@@ -387,19 +353,18 @@ function createDeepTimeLod(deps) {
   return {
     lodSecondsAtAge, lodSecondsAtAgeWithAlpha, lodHoursAtAge,
     lunisolarPrecessionRateArcsecPerYrAtAge, lunisolarPrecessionPeriodYearsAtAge,
-    // the constant's scaling alone (the hybrid's integrand) and the rate's obliquity factor (plan 12 §8)
+    // the constant's scaling alone (the hybrid's integrand) and the rate's obliquity factor
     lunisolarTorqueRateArcsecPerYrAtAge, lunisolarTorquePeriodYearsAtAge, precessionObliquityFactorAtAge,
-    // the composition's terms (plan 06 Phase 3 S3, the lunisolar surface): (a₀/a_M)³ and f_S + (1 − f_S)(a₀/a_M)³
+    // the composition's terms (the lunisolar surface): (a₀/a_M)³ and f_S + (1 − f_S)(a₀/a_M)³
     lunarTorqueFactorAtAge: composed.lunarTorqueFactorAtAge,
     precessionTorqueTermAtAge: composed.torqueTermAtAge,
     siderealYearSecondsAtAge, tropicalYearSecondsAtAge, tropicalYearDaysAtAge,
-    // the same law on a caller-supplied LOD (the browser's L1-α sensitivity sweep; plan 09 step 7c)
+    // the same law on a caller-supplied LOD (the browser's L1-α sensitivity sweep)
     tropicalYearSecondsFromLod,
-    // the solar day's ecliptic missing-motion term on the nodal period (plan 06 T2 item; one home)
+    // the solar day's ecliptic missing-motion term on the nodal period (one home)
     eclipticLodCorrectionSecondsAtAge,
     yearInDaysAtAge, deltaTRawSecondsAtAge, lodSecondsWithCorrectionsAtAge,
     dLodDtDecompositionAtAge,
-    // (plan 09 step 7a: the frozen era clock's counter and bases left with the last comb)
   };
 }
 

package/src/earth/deep-orbital-history.cjs

@@ -30,6 +30,49 @@
 
 'use strict';
 
+/**
+ * The CO-ROTATING anchored mode sum — the z-side evaluator of the history
+ * (the ϖ̇ fix, measured): the plain DC anchor residual is a spurious
+ * zero-frequency mode that preserves the J2000 VALUE but dilutes the local
+ * argument rate by |sum(0)|/|anchor| (z: 11.47 → 7.93″/yr, a 92 s-class
+ * anomalistic-year bias in the mode tier). Here the residual ROTATES at the
+ * raw sum's own local arg-rate g* = Im(ż·z̄)/|z|² (analytic), so at t = 0
+ * both the value AND the rate are exact, and at depth the residual is one
+ * |R|-amplitude term at a physical frequency instead of a DC offset (deep
+ * envelope measured unchanged: [0.0001, 0.0672] → [0.0003, 0.0678] over
+ * 10–50 Myr). Exported: e(t) = hypot(sum(t)) needs no spin integration,
+ * so a fixture can read a deep eccentricity node here (bit-identical at the
+ * grid nodes; the published surface is the grid's linear interpolation and
+ * stays on the history).
+ * @param {ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>} modes
+ * @param {[number, number]} anchor the J2000 value (q, p) the sum is pinned to
+ * @returns {(t: number) => [number, number]} t in years from J2000
+ */
+function createCoRotatingAnchoredModeSum(modes, anchor) {
+  const sum = (/** @type {number} */ t) => {
+    let re = 0, im = 0;
+    for (const m of modes) {
+      const c = Math.cos(m.omegaRadPerYr * t), s = Math.sin(m.omegaRadPerYr * t);
+      re += m.re * c - m.im * s;
+      im += m.re * s + m.im * c;
+    }
+    return [re, im];
+  };
+  let q0 = 0, p0 = 0, dq0 = 0, dp0 = 0;
+  for (const m of modes) {
+    q0 += m.re; p0 += m.im;
+    dq0 += -m.omegaRadPerYr * m.im;
+    dp0 += m.omegaRadPerYr * m.re;
+  }
+  const gStar = (dp0 * q0 - dq0 * p0) / (q0 * q0 + p0 * p0);
+  const R = [anchor[0] - q0, anchor[1] - p0];
+  return (/** @type {number} */ t) => {
+    const [x, y] = sum(t);
+    const c = Math.cos(gStar * t), s = Math.sin(gStar * t);
+    return [x + R[0] * c - R[1] * s, y + R[0] * s + R[1] * c];
+  };
+}
+
 /**
  * @param {{
  *   zModes: ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>,
@@ -97,43 +140,14 @@ function createDeepOrbitalHistory({
   // CO-ROTATING anchor (the ϖ̇ fix, measured): the DC residual above is a
   // spurious zero-frequency mode — it preserves the J2000 VALUE but dilutes
   // the local argument rate by |sum(0)|/|anchor| (z: 11.47 → 7.93″/yr, a
-  // 92 s-class anomalistic-year bias in the mode tier). Here the residual
-  // ROTATES at the raw sum's own local arg-rate g* = Im(ż·z̄)/|z|²
-  // (analytic), so at t=0 both the value AND the rate are exact, and at
-  // depth the residual is one |R|-amplitude term at a physical frequency
-  // instead of a DC offset (deep envelope measured unchanged:
-  // [0.0001, 0.0672] → [0.0003, 0.0678] over 10–50 Myr).
+  // 92 s-class anomalistic-year bias in the mode tier). The co-rotating form
+  // is the module-level createCoRotatingAnchoredModeSum (exported — the
+  // plan-06 fixture reads the deep eccentricity from it directly).
   // z ONLY: ζ keeps the DC anchor DELIBERATELY — its anchor is ~0 by
   // definition (Earth's inclination to the J2000 ecliptic at J2000), i.e.
   // the coordinate origin; a rotating residual there would inject a fake
   // wobble into n̂(t).
-  const mkAnchoredCoRot = (
-    /** @type {ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>} */ modes,
-    /** @type {[number, number]} */ anchor,
-  ) => {
-    const sum = (/** @type {number} */ t) => {
-      let re = 0, im = 0;
-      for (const m of modes) {
-        const c = Math.cos(m.omegaRadPerYr * t), s = Math.sin(m.omegaRadPerYr * t);
-        re += m.re * c - m.im * s;
-        im += m.re * s + m.im * c;
-      }
-      return [re, im];
-    };
-    let q0 = 0, p0 = 0, dq0 = 0, dp0 = 0;
-    for (const m of modes) {
-      q0 += m.re; p0 += m.im;
-      dq0 += -m.omegaRadPerYr * m.im;
-      dp0 += m.omegaRadPerYr * m.re;
-    }
-    const gStar = (dp0 * q0 - dq0 * p0) / (q0 * q0 + p0 * p0);
-    const R = [anchor[0] - q0, anchor[1] - p0];
-    return (/** @type {number} */ t) => {
-      const [x, y] = sum(t);
-      const c = Math.cos(gStar * t), s = Math.sin(gStar * t);
-      return [x + R[0] * c - R[1] * s, y + R[0] * s + R[1] * c];
-    };
-  };
+  const mkAnchoredCoRot = createCoRotatingAnchoredModeSum;
   const s2h = Math.sin(anchorInclEclipticDeg / 2 * D2R);
   const zetaAnchor = /** @type {[number, number]} */ (
     [s2h * Math.cos(anchorAscNodeEclipticDeg * D2R), s2h * Math.sin(anchorAscNodeEclipticDeg * D2R)]);
@@ -676,4 +690,4 @@ function createDeepOrbitalHistory({
   };
 }
 
-module.exports = { createDeepOrbitalHistory };
+module.exports = { createDeepOrbitalHistory, createCoRotatingAnchoredModeSum };

package/src/earth/year-lengths.cjs

@@ -42,6 +42,7 @@
 //   • The frozen era clock stays the certified era device (opt-out /
 //     era-certification record), not a member of this family.
 
+const { wrap180 } = require('../math/angles.cjs');
 const { createCardinalStructure } = require('../cardinal/one-source-structure.cjs');
 const { createSiderealYearChannel } = require('./sidereal-year-channel.cjs');
 
@@ -64,7 +65,7 @@ function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }) {
   if (typeof sampleAt !== 'function' || typeof massLossSiderealSecondsAtYearFn !== 'function') {
     throw new Error('createYearLengths: sampleAt and massLossSiderealSecondsAtYearFn are required');
   }
-  const w180 = (/** @type {number} */ d) => ((d + 540) % 360) - 180;
+  const w180 = wrap180;   // math/angles — bit-identical to the former inline form for the sub-degree differences wrapped here
   const chan = createSiderealYearChannel({ massLossSiderealSecondsAtYearFn });
   /** A year length from a longer year and a retrograde angular rate of the
    *  reference point, RATE form: 1/T = 1/T_long + p/(360°·JY) — exact (the

package/src/eclipse/besselian.cjs

@@ -69,6 +69,8 @@
 
 'use strict';
 
+const { wrap180 } = require('../math/angles.cjs');
+
 /**
  * @typedef {Object} BesselianDeps
  * @property {(dDaysTT: number) => {lonDeg: number, latDeg: number, distKm: number}} moonFullAtDaysTT
@@ -88,6 +90,23 @@
  *   gmstMeanSiderealT2Deg: number, speedOfLightKmS: number }} constants
  */
 
+/** Mean sidereal time at Greenwich, degrees (mean equinox — pairs with the
+ *  mean-of-date body longitudes; coefficients single-sourced from
+ *  astro-reference.json physicalConstants). THE ONE GMST of the model: the
+ *  eclipse location tier, the API's observer tier, the eclipse and lunar
+ *  audits and the browser's diagnostics all read it. Meeus 12.4's
+ *  cubic term (−T³/38,710,000) is deliberately absent: it is 0.9″ in the
+ *  eclipse canon and at most 93″ (6 s of time) at the edge of the observer
+ *  era, and every recorded eclipse and lunar value was measured without it.
+ *  @param {number} jdUT
+ *  @param {{ j2000JD: number, julianCenturyDays: number, gmstMeanSiderealT0Deg: number, gmstMeanSiderealRateDegPerDay: number, gmstMeanSiderealT2Deg: number }} K
+ *  @returns {number} */
+function computeGmstDeg(jdUT, K) {
+  const T = (jdUT - K.j2000JD) / K.julianCenturyDays;
+  return ((K.gmstMeanSiderealT0Deg + K.gmstMeanSiderealRateDegPerDay * (jdUT - K.j2000JD)
+    + K.gmstMeanSiderealT2Deg * T * T) % 360 + 360) % 360;
+}
+
 /** @param {BesselianDeps} deps */
 function createBesselian(deps) {
   const K = deps.constants;
@@ -97,14 +116,8 @@ function createBesselian(deps) {
   const R_MOON_KM = K.moonDiameterKm / 2;
   const F = 1 / K.earthFlatteningInverse;
 
-  /** Mean sidereal time at Greenwich, degrees (mean equinox — pairs with the
-   *  mean-of-date body longitudes; coefficients single-sourced from
-   *  astro-reference.json physicalConstants). @param {number} jdUT @returns {number} */
-  function gmstDeg(jdUT) {
-    const T = (jdUT - K.j2000JD) / K.julianCenturyDays;
-    return ((K.gmstMeanSiderealT0Deg + K.gmstMeanSiderealRateDegPerDay * (jdUT - K.j2000JD)
-      + K.gmstMeanSiderealT2Deg * T * T) % 360 + 360) % 360;
-  }
+  /** @param {number} jdUT @returns {number} */
+  const gmstDeg = (jdUT) => computeGmstDeg(jdUT, K);
 
   /** Ecliptic-of-date → equatorial-of-date vector (km).
    *  @param {number} lonDeg @param {number} betDeg @param {number} distKm
@@ -150,7 +163,7 @@ function createBesselian(deps) {
     const moon = deps.moonFullAtDaysTT(dTT);
     const H_DAYS = 0.02;
     const dLon = deps.moonFullAtDaysTT(dTT + H_DAYS).lonDeg - deps.moonFullAtDaysTT(dTT - H_DAYS).lonDeg;
-    const moonRateDegPerDay = (((dLon + 540) % 360) - 180) / (2 * H_DAYS);
+    const moonRateDegPerDay = wrap180(dLon) / (2 * H_DAYS);
     const tauDays = moon.distKm / K.speedOfLightKmS / 86400;
     return {
       S: eclToEq(sunLon, 0, sunDistanceKm(year, sunLon), eps),
@@ -310,7 +323,7 @@ function createBesselian(deps) {
     return { kind, magnitude: Math.max(0, atMax.magnitude), maxJd, contacts, centralDurationSeconds };
   }
 
-  return { umbraGroundAt, shadowStateAt, localCircumstances };
+  return { umbraGroundAt, shadowStateAt, localCircumstances, gmstDeg };
 }
 
-module.exports = { createBesselian };
+module.exports = { createBesselian, computeGmstDeg };

package/src/index.js

@@ -61,7 +61,7 @@ const NEVER_INJECTABLE = Object.keys(REFERENCE_DATA);
 
 /**
  * Generated at build time from `public/input/astro-reference.json` (plus the
- * deep-time block the generator assembles; model-parameters.json is deleted, plan 11)
+ * deep-time block the generator assembles)
  * (§2g) by `tools/constants/generate.mjs`. 324 values across 15 blocks.
  *
  * Free parameters and measured anchors only. Validation targets are excluded by
@@ -95,29 +95,21 @@ export { FITTED_COEFFICIENTS, COEFFICIENTS_HASH } from './constants/index.js';
 export { createEpochPrimitives } from './layer0/index.js';
 export { deriveEpochParams } from './layer0/derive-params.js';
 export { createDerivedViews } from './layer1/index.js';
-// (Plan 09 step 7a: the Phase 7 integrated-phase machinery `createPhaseMachinery`
-// — the ∫dt/H_era table and the cycle counters every comb read — and the
-// frozen cardinal record evaluator `createCardinalModel` (the §10 derived
-// form on the deleted CARDINAL_POINT_* coefficients) are DELETED; the
-// cardinal instants are createModel().cardinal, the one Sun's crossings.)
 // The deep-time calendar (year ↔ JD under the drifting year length) — ONE
-// implementation for both engines since the year/day collapse (formerly the
-// PHASE-B-DUPLICATE family). CJS on purpose: tools/lib requires the same
+// implementation for both engines. CJS on purpose: tools/lib requires the same
 // file via the exports-map subpath; re-exported here so bundled ESM
 // consumers need only the package root.
 export { createDeepCalendar } from './phase/deep-calendar.cjs';
-// Plan 13 B1: the civil calendar ↔ JD (Meeus ch. 7 with the 1582-10-15
-// reform) and the angle constants/wraps — ONE home; the browser's eighteen
-// calendar copies and 27 local D2R constants import these.
+// The civil calendar ↔ JD (Meeus ch. 7 with the 1582-10-15 reform) and the
+// angle constants/wraps — ONE home for every runtime.
 export { GREGORIAN_START, GREGORIAN_START_JD, calendarToJD, jdToCalendar, jdToDecimalYear, dayFractionToHMS, julianCalendarToJD, jdToJulianCalendar } from './time/calendar.cjs';
 export { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wrap180Closed, wrap2pi, wrapPi } from './math/angles.cjs';
-// Plan 13 B2: THE deep-time chain assembly — the GIA α channel, the Layer-0
-// bundle, the cited Moon-distance history, the ΔT cycle stack, |s₃|, the
-// LOD/ΔT core, the lunar month chain and the mean*AtAge family, assembled
-// ONCE from the generated constants; the Node twin and the browser destructure it.
+// THE deep-time chain assembly — the GIA α channel, the Layer-0 bundle, the
+// cited Moon-distance history, the ΔT cycle stack, |s₃|, the LOD/ΔT core, the
+// lunar month chain and the mean*AtAge family, assembled ONCE from the
+// generated constants; the Node twin and the browser destructure it.
 export { createDeepTimeChain, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE } from './deltat/deep-time-chain.js';
-// Plan 13 B3: the J2000 mass chain and the Moon's J2000 periods (one home for
-// three copies each) and the unit constants.
+// The J2000 mass chain, the Moon's J2000 periods and the unit constants — ONE home.
 export { computeMassChainJ2000, computeMoonPeriodsJ2000 } from './earth/mass-chain.cjs';
 export { SECONDS_PER_DAY, DAYS_PER_JULIAN_YEAR, JULIAN_CENTURY_DAYS, SECONDS_PER_JULIAN_YEAR } from './constants/units.cjs';
 // D4b: the one-source cardinal structure — the EoC layer (year lengths,
@@ -127,7 +119,6 @@ export { createCardinalStructure } from './cardinal/one-source-structure.cjs';
 // Phase 8.2 — the lunar machinery, extracted layer by layer (survey order:
 // eccentricity channel → month chain → cycle tables → arguments → series →
 // apparent). Same CJS + root re-export convention.
-// (plan 09 step 6: createMoonEccChannel — the H/3 eccentricity law — is deleted; Earth's e of date is the N-body series on every path)
 export { createMoonMonthChain } from './moon/month-chain.cjs';
 export { createChainCycleIntegrator } from './chain-cycles/index.cjs';
 export { createMoonArguments } from './moon/arguments.cjs';
@@ -137,26 +128,26 @@ export { createMoonApparent } from './moon/apparent.cjs';
 // sample) and the offset-circle wheel-angle solver: the scene twins' one placement.
 export { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG } from './earth/frame-of-date.cjs';
 // Phase 8.3 — the planet machinery, extracted by LAW over body records
-// (survey order: geometry → channels → chains). The ψ/K amplitude laws
-// (planets/fibonacci-laws), the record-based orientation forms
-// (planets/orientation) and the ascending-node integrator
-// (planets/asc-node-integrator) are DELETED at plan 07 R5/R6: the planets'
-// elements of date have one home, the N-body chain.
+// (survey order: geometry → channels → chains). The planets' elements of date
+// have one home, the N-body chain.
 export { derivePlanetGeometry, scaffoldOrbitInputs } from './planets/geometry.cjs';
 // (plan 09 step 6: eccentricityFromCycles / computeEccentricityIntegrated — the H/16 law-of-cosines eccentricity — are deleted; no runtime read them)
 export * as planetOrbitChain from './planets/orbit-chain.cjs';
 // P5/K4.6c — the engine-D Keplerian chain (pure evaluator + the embedded
 // governed artifact; the browser flag path consumes exactly these).
 export { computeSecularShape } from './planets/secular-shape.cjs';
-export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, computeApsidalSecularDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from './planets/keplerian-chain.cjs';
+// `computeApsidalTangentDegPerYr` is the chain's ±150-yr central tangent of the perihelion
+// longitude (the panels' "Prec." cell — osculating, NOT the ±1-kyr secular stencil of the
+// "Secular" cell). Its former name `computeApsidalSecularDegPerYr` left the API at 24.0.0
+// (no reader in the simulator, the website or the paper).
+export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, computeApsidalTangentDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from './planets/keplerian-chain.cjs';
 export { computeEquatorNodeOriginSFrameDeg, convertNodeSFrameToEquatorOriginDeg } from './planets/inv-plane-frame.cjs';
 export { createSecularSeriesOverride } from './planets/secular-series.cjs';
-// Plan 06 Phase 7 — the planets' spin channel (the precession constant from
-// each planet's own torques, the spin integrated on its own ζ plane history).
+// The planets' spin channel (the precession constant from each planet's own
+// torques, the spin integrated on its own ζ plane history).
 export { computePlanetPrecessionConstant, createPlanetSpinChannel, createPlanetSpinChannelFromArtifacts, computeObliquityJ2000Deg } from './planets/spin-channel.cjs';
 export { CHAIN_ARTIFACT, CHAIN_ARTIFACT_HASH } from './planets/chain-artifact.js';
-// Engine-switch Stage B, T5d-(d): the deep-time Earth-z channel (the deep
-// lunar-chain modulation rides it; era consumers stay on the H/3 channel).
+// The deep-time Earth-z channel (the deep lunar-chain modulation rides it).
 export { createDeepEccChannel } from './moon/deep-ecc-channel.cjs';
 export { DEEP_MODES_ARTIFACT, DEEP_MODES_ARTIFACT_HASH } from './moon/deep-modes-artifact.cjs';
 // Stage C-3: the ONE home of the obliquity-hybrid / deep-orbital-history
@@ -166,28 +157,25 @@ export { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 export { createSiderealYearChannel } from './earth/sidereal-year-channel.cjs';
 // S2: THE ONE of-date year-length family (years + precession beats, SI seconds)
 export { createYearLengths, computeMeanTropicalYearJ2000Days, ONE_FAMILY_WINDOW_YEARS } from './earth/year-lengths.cjs';
-// Plan 06 D6: THE ONE home of the composed lunisolar precession rate ψ̇(t)
-// (spin ω(t) × [solar torque + lunar torque on the recession history]) —
-// leg-1's physical rate; the structural H(t)/13 clock is its named diagnostic.
+// THE ONE home of the composed lunisolar precession rate ψ̇(t) (spin ω(t) ×
+// [solar torque + lunar torque on the recession history]) — falsification
+// leg 1's physical rate.
 export { computeSolarTorqueShare, createComposedPrecession } from './earth/precession-composed.cjs';
-// Plan 12 §8 option (a): the cited mean obliquity of date — the second column read from the
-// Earth–Moon history; the precession RATE of date carries its cos ε̄(t)/cos ε̄(0) on the constant.
+// The cited mean obliquity of date — the second column read from the Earth–Moon
+// history; the precession RATE of date carries its cos ε̄(t)/cos ε̄(0) on the constant.
 export { createMeanObliquityHistory } from './earth/mean-obliquity-history.cjs';
 // (planets/predict.cjs — the planet predictive-precession feature basis — was
 // retired at plan 06 R8 with PREDICT_COEFFS_PHYSICAL.)
 // L10 — the composition front door: one law set, N body records. Thin by
 // design; engines keep their direct call sites (see planets/model.cjs).
 export { createPlanetModel } from './planets/model.cjs';
-// plan 11 step 6 — the planets' presentation names, the last leaf of the
-// deleted model-parameters.json (code, not an input: not injectable, not hashed).
+// the planets' presentation names (code, not an input: not injectable, not hashed).
 export { PLANET_DISPLAY_NAMES } from './planets/display-names.cjs';
 // Phase 8.4 — the climate/ΔT machinery, extracted layer by layer.
 export { createDeltaTCycles } from './deltat/cycles.cjs';
 export { createDeepTimeLod } from './deltat/deep-time.cjs';
-// plan 12 (option A): the Earth–Moon distance history is a CITED PHYSICAL
-// INPUT (Farhat et al. 2022, the AstroGeo22 ephemeris) read as a shape on the
-// model's J2000 anchor; the fitted quartic, the rock-fitted regime spline and
-// the solar-channel budget are deleted.
+// The Earth–Moon distance history is a CITED PHYSICAL INPUT (Farhat et al.
+// 2022, the AstroGeo22 ephemeris) read as a shape on the model's J2000 anchor.
 export { createMoonDistanceHistory } from './deltat/moon-distance-history.cjs';
 export { MOON_DISTANCE_ARTIFACT, MOON_DISTANCE_ARTIFACT_HASH } from './deltat/moon-distance-artifact.cjs';
 export { deltaTEspenakMeeusCanonSeconds } from './deltat/historical.cjs';
@@ -199,9 +187,6 @@ export { createEclipseFinders } from './eclipse/finders.cjs';
 // (reference/published-curves MIGRATED to @essrt/reference at the 4.0.0
 // major — the published model package carries no comparison references;
 // the K2 one-way wall is package-structural.)
-// Phase 9 — S-P8: the fitted sun-longitude harmonic stack (three copies → one).
-// (plan 09 step 5a: createSunLongitudeCorrection — the Step-0 Sun comb's evaluator,
-//  retired from the Sun node at plan 06 R4 — is deleted with the comb and its fitter.)
 
 /**
  * Build a model bound to a set of constants.
@@ -270,15 +255,12 @@ export const createModel = (constants = GENERATED, opts = {}) => {
     constants: ctx,
     hash,
 
-    // Plan 09: the fit-era lattice identities H/13, H/3, H/16
-    // (`computeLatticePeriodsYears`, DEVICE tier) left this surface — the
-    // published periods are `epoch.axialPrecessionYearsAtYear` and the
+    // The published periods are `epoch.axialPrecessionYearsAtYear` and the
     // lunisolar surface's apsidal and perihelion-of-date periods; the
-    // counterfactual gate now reads the obliquity anchor through the assembly.
+    // counterfactual gate reads the obliquity anchor through the assembly.
 
     /**
-     * Earth eccentricity at a year — the Phase-6 promise, fulfilled by the
-     * §7a assembly (law-of-cosines on the integrated H/16 phase).
+     * Earth eccentricity at a year — the one-source series through the §7a assembly.
      * @param {number} year
      * @returns {number}
      */