@essrt/physics 20.0.0 → 21.0.0

package/src/deltat/deep-time.cjs

@@ -12,11 +12,6 @@
  *    factory wiring).
  *  - moiFactorAtAge — the L1-orbital-coupled GIA α(t). The lattice-α pin
  *    machinery (_withLatticeAlpha) is ENGINE STATE and stays there.
- *  - siderealYearDaysFourierAt — the IAU-base integrated-phase Fourier
- *    evaluator (Phase D matched pair). 8.4-1 S-D2: the browser's old Actual
- *    divided by its epoch-aware f(Y) evaluator instead, double-counting the
- *    deep-time drift already in mean(t) — 12% wrong at −400 Ma. Both engines
- *    now inject the ripple-form evaluator.
  *  - cycleLodSumAt / swingLodAt / swingLodRateAt — the engine-GATED
  *    cycle-stack sums (flags are engine dispatch).
  *
@@ -50,11 +45,9 @@ const { createComposedPrecession } = require('../earth/precession-composed.cjs')
  * @property {number} aMoonNowMetres - Farhat polynomial anchor
  * @property {number} alpha1PerMa @property {number} alpha3PerMa3
  * @property {number} alpha4PerMa4
- * @property {number} holisticYearJ2000
- * @property {number} lodNowH13Seconds - the H/13-identity mean LOD at J2000
+ * @property {number} lodNowH13Seconds - ω₀, the J2000 day length (the SI day since plan 09; identifier kept)
  * @property {number} meanSiderealYearJ2000Seconds
  * @property {number} solarMassLossFracPerYear
- * @property {number} siderealYearDaysKinematicJ2000 - Actual-LOD numerator
  * @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
  */
@@ -64,7 +57,6 @@ const { createComposedPrecession } = require('../earth/precession-composed.cjs')
  * @property {DeepTimeLodConstants} constants
  * @property {(tMa: number) => number} moonDistanceMetresAtAge
  * @property {(tMa: number) => number} moiFactorAtAge - GIA α(t)
- * @property {(year: number) => number} siderealYearDaysFourierAt
  * @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
@@ -117,26 +109,18 @@ function createDeepTimeLod(deps) {
     return (s === null) ? null : s / 3600;
   }
 
-  /** The FROZEN ERA CLOCK's phase convention: H_era(t) = H_J2000 · LOD(t)/LOD_J2000
-   * (pure spin scaling — the pre-Phase-3 "H/13 identity"). The frozen devices
-   * (the cardinal era clock, the year-length comb family, the ∫dt/H phase
-   * table) were fitted against THIS counter; it ships with their coefficients
-   * as a named device constant (plan 06 D8: two named counters). It is NOT
-   * the unit H(t) below.
-   * @param {number} t_Ma @returns {number|null} */
-  function eraClockHAtAge(t_Ma) {
-    return eraClockHFromLod(lodSecondsAtAge(t_Ma));
-  }
+  // (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).
-  /** @param {number|null} LOD_s @returns {number|null} */
-  function eraClockHFromLod(LOD_s) {
-    if (LOD_s === null) return null;
-    return K.holisticYearJ2000 * LOD_s / K.lodNowH13Seconds;
-  }
 
   // The composed lunisolar precession rate (plan 06 D6 → Phase 3 → S5):
   // ψ̇(t) = [ω(t)/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M(t))³], p₀ the model's
@@ -155,23 +139,6 @@ function createDeepTimeLod(deps) {
     yearToTMa: /** @param {number} year */ (year) => (2000 - year) / 1e6,
   });
 
-  /** The unit H(t) = H_era / [f_S + (1 − f_S)(a₀/a_M)³] — the internal unit
-   * scales WITH the composed lunisolar precession period (H(t)/T_p(t) =
-   * H₀/axial0 = 13.011 at every epoch) but is NOT 13 of them: H₀ was fitted
-   * on the perihelion-of-date beat, and the "H = 13·T_p" claim is retired
-   * (plan 06 S5; docs/retired-record.md). Spelled in the SAME operations as
-   * layer0's holisticHCore so the layer0 gate holds the twins bit-identical.
-   * ≡ H_era wherever (a₀/a_M)³ ≈ 1. Identifier kept (plan 06 P4).
-   * @param {number} t_Ma @returns {number|null} */
-  function hAtAge(t_Ma) {
-    return hFromLod(t_Ma, lodSecondsAtAge(t_Ma));
-  }
-  /** @param {number} t_Ma @param {number|null} LOD_s @returns {number|null} */
-  function hFromLod(t_Ma, LOD_s) {
-    const hEra = eraClockHFromLod(LOD_s);
-    const term = composed.torqueTermAtAge(t_Ma);
-    return hEra === null || term === null ? null : hEra / term;
-  }
 
   /** The composed lunisolar precession rate, ″/yr — p₀·[ω/ω₀]·[f_S + (1 − f_S)
    * (a₀/a_M)³] on the derived J2000 anchor (S5). ≡ 1,296,000·(H₀/axial0)/H(t).
@@ -196,32 +163,32 @@ 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, and the entire source of the
-    // day-count invariant's 0.17 ppm residual: with the product form,
-    // H·(days/yr)·(AU₀/AU)² = TOTAL_DAYS_IN_H holds EXACTLY
-    // (tools/explore/deep-time-sensitivity.js §5; docs/99 §near-invariant).
+    // 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).
     return K.meanSiderealYearJ2000Seconds * (1 - mass_loss_fraction) * (1 - mass_loss_fraction);
   }
 
-  /** @param {number} t_Ma @returns {number} tropical = sidereal · (1 − 13/H(t)).
-   * The 13/H(t) here is the unit's CALENDAR convention (one turn per H₀/13
-   * of the unit — the kinematic day/year identities on H/(H − 13) and the
-   * deep JD↔year calendar ride it), NOT a precession claim: the published
-   * precession period is the composed T_p (S5), 0.086 % apart. Kept
-   * unchanged in S5 so nothing certified moves (calendar, goldens, CSV);
-   * plan-06 D2 DECIDED: relabel, not re-base — the 0.086 % is the fit
-   * anchor's convention (H₀/T_p = 13.011, constant at every epoch beyond
-   * ±2 Myr, measured), and re-basing would move every comb coefficient
-   * for a rounding-level gain. This tier stays a named device. */
+  /** 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).
+   * 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} */
   function tropicalYearSecondsAtAge(t_Ma) {
     return tropicalYearSecondsFromLod(t_Ma, lodSecondsAtAge(t_Ma));
   }
   /** @param {number} t_Ma @param {number|null} LOD_s @returns {number} */
   function tropicalYearSecondsFromLod(t_Ma, LOD_s) {
     const sidSec = siderealYearSecondsAtAge(t_Ma);
-    const Ht = hFromLod(t_Ma, LOD_s);
-    if (Ht === null) return sidSec * (1 - 13 / K.holisticYearJ2000);
-    return sidSec * (1 - 13 / Ht);
+    const term = composed.torqueTermAtAge(t_Ma);
+    const p0 = composed.p0ArcsecPerYr;
+    const rate = (LOD_s === null || !(LOD_s > 0) || term === null) ? p0 : (K.lodNowH13Seconds / LOD_s) * p0 * term;
+    return 1 / (1 / sidSec + rate / (1296000 * 365.25 * 86400));
   }
 
   /** @param {number} t_Ma @returns {number|null} */
@@ -237,25 +204,9 @@ function createDeepTimeLod(deps) {
     return tropicalYearSecondsAtAge(t_Ma) / LOD_s;
   }
 
-  // ── The frozen era clock's twins (device tier) ──────────────────────────
-  // The same forms on H_era — the bases the frozen cardinal clock (its drift
-  // integrand and real-LOD convention) and the year-length comb family were
-  // fitted against. Bit-identical to the pre-Phase-3 tropicalYearSecondsAtAge /
-  // yearInDaysAtAge; the physical (unit) forms above are what the model
-  // publishes.
-  /** @param {number} t_Ma @returns {number} */
-  function eraClockTropicalYearSecondsAtAge(t_Ma) {
-    const sidSec = siderealYearSecondsAtAge(t_Ma);
-    const Ht = eraClockHAtAge(t_Ma);
-    if (Ht === null) return sidSec * (1 - 13 / K.holisticYearJ2000);
-    return sidSec * (1 - 13 / Ht);
-  }
-  /** @param {number} t_Ma @returns {number|null} */
-  function eraClockYearInDaysAtAge(t_Ma) {
-    const LOD_s = lodSecondsAtAge(t_Ma);
-    if (LOD_s === null) return null;
-    return eraClockTropicalYearSecondsAtAge(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
@@ -323,17 +274,9 @@ function createDeepTimeLod(deps) {
     return tidal + deps.cycleLodSumAt(year);
   }
 
-  /** ACTUAL LOD: tidal mean × Fourier ripple —
-   *  actual(t) = mean(t) × Y_days_kinematic_J2000 / Y_days_fourier(year).
-   *  Fixed bases (the ripple isolates year-to-year variation; the deep-time
-   *  drift is already in mean(t)). @param {number} t_Ma @returns {number|null} */
-  function lodSecondsActualAtAge(t_Ma) {
-    const mean_t = lodSecondsAtAge(t_Ma);
-    if (mean_t === null) return null;
-    const year_at_t = 2000 - t_Ma * 1e6;
-    const Y_days_fourier = deps.siderealYearDaysFourierAt(year_at_t);
-    return mean_t * K.siderealYearDaysKinematicJ2000 / Y_days_fourier;
-  }
+  // (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.
@@ -401,18 +344,19 @@ function createDeepTimeLod(deps) {
   }
 
   return {
-    lodSecondsAtAge, lodSecondsAtAgeWithAlpha, lodHoursAtAge, hAtAge,
+    lodSecondsAtAge, lodSecondsAtAgeWithAlpha, lodHoursAtAge,
     lunisolarPrecessionRateArcsecPerYrAtAge, lunisolarPrecessionPeriodYearsAtAge,
     // the composition's terms (plan 06 Phase 3 S3, 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)
+    tropicalYearSecondsFromLod,
     // the solar day's ecliptic missing-motion term on the nodal period (plan 06 T2 item; one home)
     eclipticLodCorrectionSecondsAtAge,
     yearInDaysAtAge, deltaTRawSecondsAtAge, lodSecondsWithCorrectionsAtAge,
-    lodSecondsActualAtAge, dLodDtDecompositionAtAge,
-    // the frozen era clock's named counter and bases (device tier, plan 06 D8)
-    eraClockHAtAge, eraClockTropicalYearSecondsAtAge, eraClockYearInDaysAtAge,
+    dLodDtDecompositionAtAge,
+    // (plan 09 step 7a: the frozen era clock's counter and bases left with the last comb)
   };
 }
 

package/src/deltat/recession-history.cjs

@@ -4,7 +4,8 @@
  *
  * Background (docs/99 §"The ungated mid-Precambrian window", docs/106 §D,
  * and the research record in tools/explore/{farhat-divergence-probe,
- * regime-aware-recession-study,fit-regime-recession}.js): Farhat 2022's
+ * fit-regime-recession}.js — the hand-shaped regime-aware-recession-study.js
+ * was deleted at plan 09 step 7e): Farhat 2022's
  * central result is that lunar recession followed a resonant STAIRCASE —
  * a single smooth quartic anchored at the modern LLR rate, the gated
  * 0–650 Ma era and the Roche crossing cannot follow it through the

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

@@ -16,10 +16,10 @@
  * shim tools/lib/deep-orbital-history.js (the deep insolation extractor),
  * and the browser's published ε surfaces (VFP chart, Earth-panel row) via
  * the generated deep-modes embed. The ERA machinery (frame conversions,
- * lunar D5, besselian, the scene's A-solve and Step 6b) stays on the
- * certified fitted law — measured: the hybrid differs from it by 5″ rms in
- * 1900–2100 / 19.5″ in 1600–2400, mostly hybrid-side there — the same
- * certification split as e and ϖ.
+ * lunar D5, besselian) rode the comb-era fitted ε law until plan 09 deleted
+ * it (measured at the time: the hybrid differed from that law by 5″ rms in
+ * 1900–2100 / 19.5″ in 1600–2400, mostly hybrid-side); since then every
+ * surface reads this hybrid — the one obliquity of date.
  *
  * The climatic precession is the PHYSICAL of-date angle from the moving
  * equinox γ = unit(ŝ×n̂) to the perihelion direction, measured in the orbit

package/src/earth/precession-composed.cjs

@@ -26,9 +26,9 @@
  * composed period — H(t)/T_p(t) = H₀/axial0 = 13.011 at every epoch — but
  * is not 13 of them; the 13 is a fit-era label that leaves the physics.
  * The pre-Phase-3 structural clock 13·1,296,000/H_era(t) (pure ω-scaling,
- * H_era = H₀·LOD/LOD₀) survives only as the FROZEN ERA CLOCK's named phase
- * convention (eraClockHAtAge) and in docs/retired-record.md, with its
- * record: it read 70.9 ″/yr against Lantink et al. 2022's 108.6 ± 8.5 at
+ * H_era = H₀·LOD/LOD₀) survived as the FROZEN ERA CLOCK's named phase
+ * convention (eraClockHAtAge) until plan 09 step 7a deleted it with the last
+ * comb; its record (doc 110 chain 2, the paleo-anchors gate rows): it read 70.9 ″/yr against Lantink et al. 2022's 108.6 ± 8.5 at
  * 2.46 Ga and 65.2 against Meyers & Malinverno 2018's 85.79 ± 2.72 at
  * 1.4 Ga, where the composed rate reads 104.6 and 86.6.
  *

package/src/earth/year-lengths.cjs

@@ -159,4 +159,19 @@ function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }) {
  */
 const ONE_FAMILY_WINDOW_YEARS = 2000000;
 
-module.exports = { createYearLengths, ONE_FAMILY_WINDOW_YEARS };
+/** The J2000 mean tropical year, DAYS — derived in ONE home from two IAU
+ *  inputs, the sidereal year and the general precession in longitude p_A
+ *  (″ per JULIAN year), in the exact rate form
+ *      1/T_trop = 1/T_sid + p_A / (1,296,000″ · 365.25 d)
+ *  — the same beat `shorterYear` evaluates of date. 365.242190469 d on the
+ *  IAU 2006 inputs. Plan 09: replaces the Meeus & Savoie 1992 input
+ *  365.2421897 (66 ms shorter, whose beat with the sidereal year read
+ *  25,770.7 yr against p_A's 25,771.6 — three inputs over-determining a
+ *  two-parameter identity). Counting p_A per tropical or sidereal year
+ *  instead moves this by 26 ms; IAU states p_A per Julian year.
+ *  @param {number} siderealYearDays @param {number} generalPrecessionArcsecPerJulianYr @returns {number} */
+function computeMeanTropicalYearJ2000Days(siderealYearDays, generalPrecessionArcsecPerJulianYr) {
+  return 1 / (1 / siderealYearDays + generalPrecessionArcsecPerJulianYr / (1296000 * 365.25));
+}
+
+module.exports = { createYearLengths, computeMeanTropicalYearJ2000Days, ONE_FAMILY_WINDOW_YEARS };

package/src/eclipse/sun-planetary-completion.cjs

@@ -67,8 +67,9 @@
  *       besselian.cjs), so a Sun-only nutation term is a frame
  *       inconsistency for elongation. Removed: the nutation-free syzygy
  *       fleet improves 3.877 → 3.756″ and the fuller-bridged Sun
- *       comparison 2.344 → 2.123″ (fq7s-jpl-preview.mjs --no-2le;
- *       fq7s-annual-channel.mjs). This file now carries ZERO fitted
+ *       comparison 2.344 → 2.123″ (measured by the fq7s preview and
+ *       annual-channel instruments, deleted at plan 09 step 7e — docs/109
+ *       §22). This file now carries ZERO fitted
  *       constants. The Earth-around-EMB wobble is DERIVED: amplitude
  *       a_M·μ/AU with μ = 1/(1+M_E/M_M), injected by the model wiring
  *       from live package constants (6.4399″ at current constants), sign
@@ -393,15 +394,13 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy, e
   return { sunPlanetaryCompletionDeg };
 }
 
-/** sha256/16 of JSON.stringify(FITTED_COEFFICIENTS.SUN_LONGITUDE_HARMONICS)
- *  at derivation time — the matched-pair fingerprint asserted by test:model.
- *  Unchanged from v1/v2: SUN_HARMONICS did not move in the D2 or N3
- *  landings. Plan 07 R10 moved the hash WITHOUT a refit: the three (sin, cos)
- *  pairs were rotated exactly from the t₀ phase origin to J2000
- *  (tools/fit/reorigin-combs-j2000.mjs), the finder Sun is bit-identical, so
- *  the residual this table was derived from is unchanged — the table stands
- *  (fidelity 0.616″, 70 terms), only the fingerprint is re-recorded. */
-const PAIRED_SUN_HARMONICS_SHA256 = '53530a28a60b3c4a';   // was cbc189cea1c20292 (eccentricity unification: Step-0 refit → N2/N3 re-derived); R10 rotation, same Sun
+/* (Plan 09 step 5a: PAIRED_SUN_HARMONICS_SHA256 — the sha256/16 of the Step-0
+ *  SUN_LONGITUDE_HARMONICS this table was derived beside, last 53530a28a60b3c4a
+ *  (the R10 rotation of a comb never refit since the eccentricity unification)
+ *  — is retired with the comb. The comb left every Sun path at plan 06 R4, so
+ *  the residual the table was derived from is the certified Sun's; the table
+ *  stands (fidelity 0.616″, 70 terms) and the carrier-rate pair below remains
+ *  its enforced matched pair.) */
 
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
  *  full-precision carrier rates (deg/cy TT) the tables pair with: since plan
@@ -421,4 +420,4 @@ const PAIRED_SUN_HARMONICS_SHA256 = '53530a28a60b3c4a';   // was cbc189cea1c2029
  *  was the I2 record-based carrier set. */
 const PAIRED_CARRIER_RATES_SHA256 = '26552b86ad69eaf2';   // running-mean series re-bank (mean=1 dump): carriers moved in the 9th digit, LONG_PERIOD_TERMS + the 8V−13E row re-derived on them; history: 8e11456fcf8db81c (I3), 2d066e92bae955e4 (I2)
 
-module.exports = { createSunPlanetaryCompletion, computeCarrierRatesDegPerCy, eccVectorOfDateEmbedded, PAIRED_SUN_HARMONICS_SHA256, PAIRED_CARRIER_RATES_SHA256 };
+module.exports = { createSunPlanetaryCompletion, computeCarrierRatesDegPerCy, eccVectorOfDateEmbedded, PAIRED_CARRIER_RATES_SHA256 };

package/src/index.js

@@ -54,7 +54,6 @@ const NEVER_INJECTABLE = Object.keys(REFERENCE_DATA);
  * @typedef {{
  *   constants: Constants,
  *   hash: string,
- *   computeLatticePeriodsYears: () => {axialPrecessionPeriodYears: number, inclinationPrecessionPeriodYears: number, perihelionPrecessionPeriodYears: number},
  *   eccentricity: (year: number) => number,
  *   identity: ModelIdentity,
  * } & ModelSurfaces} Model
@@ -95,24 +94,25 @@ 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';
-// Phase 7 — the shared integrated-phase and cardinal-point machinery (CJS on
-// purpose: tools/lib requires the same files via the exports-map subpaths;
-// re-exported here so bundled ESM consumers need only the package root).
-export { createPhaseMachinery } from './phase/index.cjs';
-// The deep-time calendar (year ↔ JD under the drifting year length, the
-// H-balanced event finder) — ONE implementation for both engines since the
-// year/day collapse (formerly the PHASE-B-DUPLICATE family).
+// (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
+// 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';
-export { createCardinalModel } from './cardinal/index.cjs';
 // D4b: the one-source cardinal structure — the EoC layer (year lengths,
 // crossing offsets, the e(t)-proportional spread) on the movement's own
-// e(t)/ϖ(t). Absolute dates deliberately excluded (per-renderer mean chain);
-// the frozen createCardinalModel stays the certified era device.
+// e(t)/ϖ(t). Absolute dates deliberately excluded (per-renderer mean chain).
 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.
-export { createMoonEccChannel } from './moon/ecc-channel.cjs';
+// (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';
@@ -127,8 +127,8 @@ export { computeEarthFrameOfDate, solveWheelAngleForLongitude } from './earth/fr
 // (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.
-export { derivePlanetGeometry } from './planets/geometry.cjs';
-export { eccentricityFromCycles, computeEccentricityIntegrated } from './planets/ecc-channel.cjs';
+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).
@@ -150,7 +150,7 @@ export { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 // D6: the sidereal-year-of-date channel (banked λ̇ ratio × the caller's mass-loss law)
 export { createSiderealYearChannel } from './earth/sidereal-year-channel.cjs';
 // S2: THE ONE of-date year-length family (years + precession beats, SI seconds)
-export { createYearLengths, ONE_FAMILY_WINDOW_YEARS } from './earth/year-lengths.cjs';
+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.
@@ -174,7 +174,8 @@ export { createEclipseFinders } from './eclipse/finders.cjs';
 // 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).
-export { createSunLongitudeCorrection } from './sun/longitude-correction.cjs';
+// (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.
@@ -243,36 +244,11 @@ export const createModel = (constants = GENERATED, opts = {}) => {
     constants: ctx,
     hash,
 
-    /**
-     * The fit-era lattice identities H/13, H/3, H/16, in years — DEVICE tier.
-     *
-     * PURE ALGEBRA OVER THE CONTEXT. No epoch, no formula, no fit. They exist
-     * here for one reason: a hash-only counterfactual test would still pass if
-     * `createModel` ignored its argument entirely. Something has to READ the
-     * context and return a number for injection to be demonstrated end to end.
-     *
-     * These are NOT the published periods (plan 06 S5): the axial precession
-     * period is `epoch.axialPrecessionYearsAtYear` (25,771.4 yr at J2000, the
-     * of-date year laws' beat); H/13 = 25,793.6 is the fit anchor's reading,
-     * 0.086 % slow — H₀ was fitted on the perihelion-of-date beat. Likewise
-     * H/3 and H/16 against the chain's apsidal and perihelion-of-date periods.
-     * The property names keep their historical form (plan 06 P4).
-     *
-     * `divisor` and `period` are never interchangeable (CLAUDE.md): 13 is the
-     * divisor, H/13 years is the identity's value. The names say which.
-     *
-     * @returns {{axialPrecessionPeriodYears: number, inclinationPrecessionPeriodYears: number, perihelionPrecessionPeriodYears: number}}
-     */
-    computeLatticePeriodsYears: () => {
-      const H = /** @type {{holisticyearLength: number}} */ (
-        /** @type {Record<string, unknown>} */ (ctx).foundational
-      ).holisticyearLength;
-      return {
-        axialPrecessionPeriodYears: H / 13,
-        inclinationPrecessionPeriodYears: H / 3,
-        perihelionPrecessionPeriodYears: H / 16,
-      };
-    },
+    // 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
+    // lunisolar surface's apsidal and perihelion-of-date periods; the
+    // counterfactual gate now reads the obliquity anchor through the assembly.
 
     /**
      * Earth eccentricity at a year — the Phase-6 promise, fulfilled by the

package/src/layer0/derive-params.js

@@ -27,9 +27,8 @@ import { computeSolarTorqueShare } from '../earth/precession-composed.cjs';
  * @property {number} alpha1PerMa
  * @property {number} alpha3PerMa3
  * @property {number} alpha4PerMa4
- * @property {number} holisticYearJ2000         H
  * @property {number} meanSiderealYearSeconds
- * @property {number} meanSiderealYearDaysKinematic
+ * @property {number} lodJ2000Seconds              ω₀: the J2000 day length, SI seconds (the SI day since plan 09)
  * @property {number} sunMassKg                 M_SUN
  * @property {number} gmEarthAloneKm3S2
  * @property {number} gmMoonAloneKm3S2
@@ -43,6 +42,7 @@ import { computeSolarTorqueShare } from '../earth/precession-composed.cjs';
  * @property {number} astronomicalUnitKm         1 AU (km)
  * @property {number} earthOrbitalEccentricityJ2000
  * @property {number} moonEclipticInclinationJ2000Deg
+ * @property {number} axialPrecessionJ2000Years    the model's ONE J2000 precession reading (1,296,000/p_A; plan 09)
  */
 
 /**
@@ -54,8 +54,10 @@ export const deriveEpochParams = (raw) => {
   const cSiMPerS = raw.speedOfLightKmPerS ? raw.speedOfLightKmPerS * 1000 : 299792458;
   const dmDtTotalKgS = raw.solarLuminosityW / (cSiMPerS * cSiMPerS) + raw.solarWindKgPerS;
 
-  // J2000 anchors derived from framework constants
-  const lodNowH13Seconds = raw.meanSiderealYearSeconds / raw.meanSiderealYearDaysKinematic;
+  // J2000 anchors. Plan 09: ω₀ is the SI day (an input, 86,400 s); the
+  // identifier keeps its historical name — the H/13 kinematic-day
+  // construction it once named (86,399.99968 s) is deleted.
+  const lodNowH13Seconds = raw.lodJ2000Seconds;
   const solarMassLossFracPerYear = dmDtTotalKgS * raw.meanSiderealYearSeconds / raw.sunMassKg;
 
   // Earth mass, moments
@@ -93,10 +95,13 @@ export const deriveEpochParams = (raw) => {
     moonEccentricityFactor,
     earthMassKg,
     earthRadiusM,
-    holisticYearJ2000: raw.holisticYearJ2000,
     lodNowH13Seconds,
     siderealYearJ2000Seconds: raw.meanSiderealYearSeconds,
     solarMassLossFracPerYear,
     precessionSolarShareJ2000,
+    // p₀, the composed precession rate's J2000 anchor — the tropical year of
+    // date's equinox term (plan 09 step 5d; the same expression the deep-time
+    // factory's composed module resolves, so the twins stay bit-identical)
+    precessionRateJ2000ArcsecPerYr: 1296000 / raw.axialPrecessionJ2000Years,
   });
 };