@essrt/physics 13.0.0 → 15.0.0

package/src/planets/model.cjs

@@ -2,31 +2,29 @@
  * createPlanetModel — the planet composition front door (Phase 8.3, L10).
  *
  * ONE law set, N body records: this factory binds an environment once and
- * runs the certified derivation chain over every body record, in the order
- * the chain requires (each step feeds the next):
+ * runs the certified derivation chain over every body record:
  *
- *   ψ constant → inclination law (amplitude, mean)
- *             → wobble period (beat of |axial| and |ICRF|)
- *   K constant → obliquity mean (snapshot form) → eccentricity law
- *             (amplitude, base, J2000 phase)
- *   geometry   (the type-branched ellipse family — the law-derived
- *              eccentricity base feeds geometry for the seven carriers;
- *              minor bodies use their record's base)
+ *   inputs     the chain's g-mode beat (wobble) and the spin channel's
+ *              derived J2000 obliquity, carried on the record since plan 06
+ *              Phase 7 commit 2
+ *   geometry   (the type-branched ellipse family on the OBSERVED J2000
+ *              eccentricity; minor bodies use their record's base)
+ *
+ * The ψ and K amplitude laws that opened this chain (ψ constant →
+ * inclination amplitude/mean; K constant → eccentricity amplitude/base/
+ * phase) are DELETED — K at plan 07 R6, ψ at plan 07 R5. A planet's
+ * inclination and eccentricity of date have one home, the N-body chain
+ * (model.js planetChainElementsAt); docs/retired-record.md is the record.
  *
  * THIN BY DESIGN. This is the composition surface, not a rewiring: both
- * engines keep their existing direct call sites into the law modules, and
- * the runtime channels (eccentricity-at-year, orientation, the ascending-
- * node integrator, predictive precession) stay direct module calls because
- * they consume engine-owned state (scene JD, epoch machinery, fitted
- * tables). What this factory adds is the seam future bodies plug into:
- * a new body is a record + (at most) one new geometry branch — see the
- * minor-body placeholder note in geometry.cjs (Phase 18 perturbation
- * types land the same way).
+ * engines keep their existing direct call sites, and the runtime channels
+ * stay direct module calls because they consume engine-owned state (scene
+ * JD, epoch machinery, fitted tables). What this factory adds is the seam
+ * future bodies plug into: a new body is a record + (at most) one new
+ * geometry branch — see the minor-body placeholder note in geometry.cjs
+ * (Phase 18 perturbation types land the same way).
  *
- * Guard semantics mirror tools/lib/constants.js verbatim: the law steps
- * run only where the record carries the required fields (fibonacciD +
- * mass fraction for the ψ/K families, perihelion + axial periods for
- * wobble); geometry runs for every body. The identity gate
+ * Geometry runs for every body. The identity gate
  * (test/planet-model-identity.test.mjs) holds this factory bit-exact
  * against the shipped Node derivation.
  */
@@ -34,18 +32,13 @@
 'use strict';
 
 const { derivePlanetGeometry } = require('./geometry.cjs');
-const FL = require('./fibonacci-laws.cjs');
 
 /**
  * @typedef {Object} PlanetModelBody
  * @property {string} [type] - geometry type branch ('I' | 'II' | 'III')
  * @property {number} solarYearInput
- * @property {number} [fibonacciD]
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]
- * @property {number} [inclinationCycleAnchor]
- * @property {boolean} [antiPhase]
- * @property {number} [perihelionEclipticYears]
  * @property {number} [wobblePeriodYears] - carriers: the chain's g-mode beat (input since Phase 7 commit 2)
  * @property {number} [obliquityMeanDeg] - carriers: the derived J2000 obliquity (input since Phase 7 commit 2)
  * @property {number} [axialTiltJ2000]
@@ -66,15 +59,10 @@ const FL = require('./fibonacci-laws.cjs');
  * @property {number} currentAUDistanceKm
  * @property {number} earthEccentricityJ2000
  * @property {number} earthPerihelionLongitudeJ2000Deg
- * @property {{ earthInvPlaneInclinationAmplitude: number,
- *   massEarthAlone: number, massSun: number }} calibration
- * @property {Record<string, number>} massFractions
  */
 
 /**
  * @typedef {Object} PlanetModelRecord
- * @property {number} [invPlaneInclinationAmplitude]
- * @property {number} [invPlaneInclinationMean]
  * @property {number} [wobblePeriodYears]
  * @property {number} [obliquityMeanDeg]
  * @property {ReturnType<typeof derivePlanetGeometry>} geometry
@@ -87,18 +75,13 @@ const FL = require('./fibonacci-laws.cjs');
  * @param {Record<string, PlanetModelBody>} bodies - keyed by body name
  *   (the key selects the body-unique geometry branches: mercury, pluto,
  *   halleys, ceres)
- * @returns {{ psiConstant: number,
- *   bodies: Record<string, PlanetModelRecord> }}
+ * @returns {{ bodies: Record<string, PlanetModelRecord> }}
  */
 function createPlanetModel(env, bodies) {
-  const psiConstant = FL.computePsiConstant({
-    earthInvPlaneInclinationAmplitude: env.calibration.earthInvPlaneInclinationAmplitude,
-    massEarthAlone: env.calibration.massEarthAlone,
-    massSun: env.calibration.massSun,
-  });
   // Plan 07 R6: kConstant, the System-Reset eccentricityAnchor
   // (balancedYear − systemResetN·H) and the t2000 phase it fed went with the
-  // eccentricity law.
+  // eccentricity law; plan 07 R5: psiConstant (ψ = 3·A_earth·√(m_E/m_☉),
+  // from env.calibration) went with the inclination law.
 
   const geomEnv = {
     holisticYears: env.holisticYears,
@@ -111,22 +94,9 @@ function createPlanetModel(env, bodies) {
   /** @type {Record<string, PlanetModelRecord>} */
   const out = {};
   for (const [key, b] of Object.entries(bodies)) {
-    const massFrac = env.massFractions[key];
     /** @type {PlanetModelRecord} */
     const rec = /** @type {PlanetModelRecord} */ ({});
 
-    if (b.fibonacciD && massFrac && b.invPlaneInclinationJ2000 !== undefined) {
-      const il = FL.computeInclinationLaw({
-        fibonacciD: b.fibonacciD, massFrac,
-        invPlaneInclinationJ2000: b.invPlaneInclinationJ2000,
-        longitudePerihelion: /** @type {number} */ (b.longitudePerihelion),
-        inclinationCycleAnchor: /** @type {number} */ (b.inclinationCycleAnchor),
-        antiPhase: /** @type {boolean} */ (b.antiPhase),
-      }, psiConstant);
-      rec.invPlaneInclinationAmplitude = il.amplitude;
-      rec.invPlaneInclinationMean = il.mean;
-    }
-
     // Plan 06 Phase 7 commit 2: the K law's cycle period and obliquity input
     // are INPUTS of the record now — the chain's own g-mode beat
     // (keplerian-chain computeSecularShape) and the derived J2000 obliquity
@@ -156,7 +126,7 @@ function createPlanetModel(env, bodies) {
     out[key] = rec;
   }
 
-  return { psiConstant, bodies: out };
+  return { bodies: out };
 }
 
 module.exports = { createPlanetModel };

package/types/constants/generated.d.ts

@@ -2,14 +2,14 @@
 // Gives the TypeScript website full type safety at the boundary (§2g) while
 // packages/physics stays JavaScript.
 
-export declare const CONSTANTS_HASH: "726b10e23d8ace09";
+export declare const CONSTANTS_HASH: "5ce921f2efa2e845";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "726b10e23d8ace09";
+  readonly hash: "5ce921f2efa2e845";
   readonly additionalBodies: {
     "pluto": {
       "name": string;
@@ -17,9 +17,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -29,9 +27,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -41,9 +37,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -53,9 +47,7 @@ export declare const DEFAULT_CONSTANTS: {
       "orbitDistanceOverride": number;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -486,13 +478,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "venus": {
@@ -500,13 +488,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "mars": {
@@ -514,13 +498,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "jupiter": {
@@ -528,13 +508,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "saturn": {
@@ -542,13 +518,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "uranus": {
@@ -556,13 +528,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "neptune": {
@@ -570,13 +538,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
   };

package/types/index.d.ts

@@ -12,9 +12,6 @@ export { createMoonArguments } from "./moon/arguments.cjs";
 export { createMoonSeries } from "./moon/series.cjs";
 export { createMoonApparent } from "./moon/apparent.cjs";
 export { derivePlanetGeometry } from "./planets/geometry.cjs";
-export * as planetFibonacciLaws from "./planets/fibonacci-laws.cjs";
-export * as planetOrientation from "./planets/orientation.cjs";
-export { integrateAscendingNode } from "./planets/asc-node-integrator.cjs";
 export * as planetOrbitChain from "./planets/orbit-chain.cjs";
 export { computeSecularShape } from "./planets/secular-shape.cjs";
 export { createSecularSeriesOverride } from "./planets/secular-series.cjs";

package/types/planets/model.d.cts

@@ -4,12 +4,8 @@ export type PlanetModelBody = {
      */
     type?: string | undefined;
     solarYearInput: number;
-    fibonacciD?: number | undefined;
     invPlaneInclinationJ2000?: number | undefined;
     longitudePerihelion?: number | undefined;
-    inclinationCycleAnchor?: number | undefined;
-    antiPhase?: boolean | undefined;
-    perihelionEclipticYears?: number | undefined;
     /**
      * - carriers: the chain's g-mode beat (input since Phase 7 commit 2)
      */
@@ -37,16 +33,8 @@ export type PlanetModelEnv = {
     currentAUDistanceKm: number;
     earthEccentricityJ2000: number;
     earthPerihelionLongitudeJ2000Deg: number;
-    calibration: {
-        earthInvPlaneInclinationAmplitude: number;
-        massEarthAlone: number;
-        massSun: number;
-    };
-    massFractions: Record<string, number>;
 };
 export type PlanetModelRecord = {
-    invPlaneInclinationAmplitude?: number | undefined;
-    invPlaneInclinationMean?: number | undefined;
     wobblePeriodYears?: number | undefined;
     obliquityMeanDeg?: number | undefined;
     geometry: ReturnType<typeof derivePlanetGeometry>;
@@ -55,12 +43,8 @@ export type PlanetModelRecord = {
  * @typedef {Object} PlanetModelBody
  * @property {string} [type] - geometry type branch ('I' | 'II' | 'III')
  * @property {number} solarYearInput
- * @property {number} [fibonacciD]
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]
- * @property {number} [inclinationCycleAnchor]
- * @property {boolean} [antiPhase]
- * @property {number} [perihelionEclipticYears]
  * @property {number} [wobblePeriodYears] - carriers: the chain's g-mode beat (input since Phase 7 commit 2)
  * @property {number} [obliquityMeanDeg] - carriers: the derived J2000 obliquity (input since Phase 7 commit 2)
  * @property {number} [axialTiltJ2000]
@@ -80,14 +64,9 @@ export type PlanetModelRecord = {
  * @property {number} currentAUDistanceKm
  * @property {number} earthEccentricityJ2000
  * @property {number} earthPerihelionLongitudeJ2000Deg
- * @property {{ earthInvPlaneInclinationAmplitude: number,
- *   massEarthAlone: number, massSun: number }} calibration
- * @property {Record<string, number>} massFractions
  */
 /**
  * @typedef {Object} PlanetModelRecord
- * @property {number} [invPlaneInclinationAmplitude]
- * @property {number} [invPlaneInclinationMean]
  * @property {number} [wobblePeriodYears]
  * @property {number} [obliquityMeanDeg]
  * @property {ReturnType<typeof derivePlanetGeometry>} geometry
@@ -99,11 +78,9 @@ export type PlanetModelRecord = {
  * @param {Record<string, PlanetModelBody>} bodies - keyed by body name
  *   (the key selects the body-unique geometry branches: mercury, pluto,
  *   halleys, ceres)
- * @returns {{ psiConstant: number,
- *   bodies: Record<string, PlanetModelRecord> }}
+ * @returns {{ bodies: Record<string, PlanetModelRecord> }}
  */
 export function createPlanetModel(env: PlanetModelEnv, bodies: Record<string, PlanetModelBody>): {
-    psiConstant: number;
     bodies: Record<string, PlanetModelRecord>;
 };
 import { derivePlanetGeometry } from "./geometry.cjs";

package/src/planets/asc-node-integrator.cjs

@@ -1,120 +0,0 @@
-/**
- * Dynamic ascending node — rate-based segment integration (Phase 8.3, L5).
- *
- * dΩ/dε = −sin(Ω)/tan(i), integrated over obliquity change with segment
- * handling at obliquity extrema and Earth-inclination crossovers (the
- * inclination DIRECTION flips the sign of the effect per segment).
- *
- * Extracted from tools/lib/orbital-engine.js
- * calculateDynamicAscendingNodeFromTilts (itself the mirror of the
- * browser's 6-arg variant — S-P5). The BROWSER variant remains engine-side
- * this slice: neither copy is fixture-pinned, so the signature unification
- * is deferred to the factory pass, WITH probes first (recorded follow-up).
- *
- * §2h BOUNDARY: this module takes {ascendingNodeDeg, inclinationDeg} — the
- * Tychosium `orbitTilta/orbitTiltb` scheme names never enter the package;
- * the engines own the atan2/hypot decomposition at their edge.
- *
- * All time-dependent inputs are INJECTED per call: the engine's obliquity
- * and Earth-inclination evaluators, the precomputed obliquity-extrema
- * lookup, the inclination-crossing finder, and (optionally) the per-planet
- * dynamic ecliptic inclination.
- */
-
-'use strict';
-
-/**
- * @param {{ ascendingNodeDeg: number, inclinationDeg: number }} tilt —
- *   static node and inclination (engine-decomposed from its scene scheme)
- * @param {number} currentYear
- * @param {{
- *   obliquityAt: (year: number) => number,
- *   earthInclinationAt: (year: number) => number,
- *   obliquityExtremaInRange: (yearMin: number, yearMax: number) => number[],
- *   inclinationCrossingsInRange: (inclinationDeg: number, yearMin: number, yearMax: number) => number[],
- *   eclipticInclinationAt: ((year: number) => number) | null,
- *   earthInclinationMeanDeg: number,
- *   earthInclinationAmplitudeDeg: number,
- *   epochYear?: number,
- * }} deps — eclipticInclinationAt non-null enables the dynamic-inclination
- *   crossover search and per-segment rates (the planetName path).
- * @returns {number} dynamic ascending node longitude, degrees 0–360
- */
-function integrateAscendingNode(tilt, currentYear, deps) {
-  const DEG2RAD = Math.PI / 180;
-  const RAD2DEG = 180 / Math.PI;
-
-  const staticOmega = ((tilt.ascendingNodeDeg % 360) + 360) % 360;
-  const planetInclination = tilt.inclinationDeg;
-
-  if (planetInclination < 1e-6) return staticOmega;
-
-  const i = planetInclination * DEG2RAD;
-  const OmegaRad = staticOmega * DEG2RAD;
-  const tanI = Math.tan(i);
-  if (Math.abs(tanI) < 1e-10) return staticOmega;
-
-  const sinOmega = Math.sin(OmegaRad);
-  const EPOCH_YEAR = deps.epochYear !== undefined ? deps.epochYear : 2000;
-
-  /** @param {number} fromYear @param {number} toYear @returns {number} */
-  const integrateEffect = (fromYear, toYear) => {
-    if (Math.abs(toYear - fromYear) < 0.1) return 0;
-
-    const yearMin = Math.min(fromYear, toYear);
-    const yearMax = Math.max(fromYear, toYear);
-    const dir = toYear >= fromYear ? 1 : -1;   // integration direction — signs the total
-
-    /** @type {number[]} */
-    let criticalYears = [yearMin, yearMax];
-
-    criticalYears.push(...deps.obliquityExtremaInRange(yearMin, yearMax));
-
-    // Find ALL inclination crossings
-    const minEarthIncl = deps.earthInclinationMeanDeg - deps.earthInclinationAmplitudeDeg;
-    const maxEarthIncl = deps.earthInclinationMeanDeg + deps.earthInclinationAmplitudeDeg;
-
-    // With a dynamic ecliptic inclination the planet can enter Earth's range
-    // even if the static value is outside — always search for crossovers.
-    if (deps.eclipticInclinationAt || (planetInclination >= minEarthIncl && planetInclination <= maxEarthIncl)) {
-      const crossIncl = deps.eclipticInclinationAt
-        ? deps.eclipticInclinationAt((yearMin + yearMax) / 2)
-        : planetInclination;
-      criticalYears.push(...deps.inclinationCrossingsInRange(crossIncl, yearMin, yearMax));
-    }
-
-    criticalYears = [...new Set(criticalYears)].sort((a, b) => a - b);
-
-    // Integrate over segments
-    let effect = 0;
-    for (let idx = 0; idx < criticalYears.length - 1; idx++) {
-      const segStart = criticalYears[idx];
-      const segEnd = criticalYears[idx + 1];
-
-      const oblStart = deps.obliquityAt(segStart);
-      const oblEnd = deps.obliquityAt(segEnd);
-      const deltaObl = (oblEnd - oblStart) * DEG2RAD;
-
-      const midYear = (segStart + segEnd) / 2;
-      const earthInclAtMid = deps.earthInclinationAt(midYear);
-
-      const dynIncl = deps.eclipticInclinationAt
-        ? deps.eclipticInclinationAt(midYear)
-        : planetInclination;
-      const inclDirection = earthInclAtMid > dynIncl ? 1 : -1;
-      const dynTanI = Math.tan(dynIncl * DEG2RAD);
-      if (Math.abs(dynTanI) < 1e-10) continue; // skip near-zero inclination
-      const segRate = -sinOmega / dynTanI;
-
-      effect += segRate * inclDirection * deltaObl * RAD2DEG;
-    }
-    return effect * dir;
-  };
-
-  const effectFromEpoch = integrateEffect(EPOCH_YEAR, currentYear);
-
-  const newOmega = staticOmega + effectFromEpoch;
-  return ((newOmega % 360) + 360) % 360;
-}
-
-module.exports = { integrateAscendingNode };

package/src/planets/fibonacci-laws.cjs

@@ -1,95 +0,0 @@
-/**
- * The ψ/K amplitude relations + the retired law-loop constructions —
- * THE shared implementation (Phase 8.3, layer L2). Status per relation
- * (doc 10 "The Six Fibonacci Relations"; the restatement: the N-body
- * chain is the planet path):
- *  - ψ (computePsiConstant/computeInclinationLaw) and K (computeKConstant/
- *    computeEccentricityLaw's amplitude) are OPEN PREDICTIONS — untested
- *    empirical amplitude laws, derived from Earth alone, kept as stated.
- *  - The wobble beat (computeWobblePeriodYears) is kinematic frame
- *    arithmetic — the permanent class.
- *  - The obliquity-mean snapshot and the eccentricity base/phase
- *    (System-Reset construction) belong to the RETIRED law framework —
- *    today serving the legacy scene scaffolding, the no-chain bodies
- *    (Pluto/Halley/Eros) and the o.fib* diagnostics (doc 72 — the closed
- *    loop; docs/37 is archived, see docs/retired-record.md).
- *
- * Extracted VERBATIM in the BROWSER's expression forms from src/script.js
- * (the ψ loop, calcWobblePeriod, the obliquity-cycle aliases,
- * calcObliquityMean's load-time snapshot branch, the K loop), which
- * tools/lib/constants.js + constants/utils.js hand-mirrored with last-ulp
- * operation-order variants (measured: 32/42 law outputs bit-exact, 10 at
- * the 1e-11..1e-15 relative class) — those variants dissolve here.
- *
- * LOAD-BEARING conventions:
- *  - ψ = 3·A_earth·√(M_EARTH_ALONE/M_SUN): the ALONE/SYSTEM mass asymmetry
- *    is calibration convention (doc 24) — switching Earth to SYSTEM shifts
- *    ψ by 0.612% and would require re-calibrating A_earth.
- *  - Wobble = beat of |axial| and |ICRF| RATES (sign-free — Venus's
- *    prograde-axial/retrograde-ICRF case); |axial| > 8H ⇒ frozen ⇒
- *    wobble = |ICRF| exactly.
- *  - The obliquity-mean law here is the SNAPSHOT form (the browser's
- *    module-load TDZ fallback — the value both engines actually ship at
- *    load); the browser's runtime integrated path stays engine-side with
- *    its Phase-8 anchors. Note it uses 1/(H/13), where the runtime
- *    integrated form uses 13/H — historical operation orders, preserved.
- *  - K-law phase offset: 90° in-phase, 270° anti-phase (Saturn) — the n=7
- *    System Reset state (all planets at mean e, Saturn falling).
- *  - The Node mirror carried an `else` branch for bodies without a wobble
- *    period — DEAD code (the loop is fibonacciD-guarded and all seven
- *    carriers have wobble periods); dropped here, recorded in the commit.
- */
-
-'use strict';
-
-/**
- * ψ constant from Earth's calibration. @param {{
- *   earthInvPlaneInclinationAmplitude: number,
- *   massEarthAlone: number, massSun: number }} c @returns {number} */
-function computePsiConstant(c) {
-  return 3 * c.earthInvPlaneInclinationAmplitude * Math.sqrt(c.massEarthAlone / c.massSun);
-}
-
-/**
- * ψ law: invariable-plane inclination amplitude and mean.
- * @param {{ fibonacciD: number, massFrac: number,
- *   invPlaneInclinationJ2000: number, longitudePerihelion: number,
- *   inclinationCycleAnchor: number, antiPhase: boolean }} b
- * @param {number} psiConstant
- * @returns {{ amplitude: number, mean: number }} */
-function computeInclinationLaw(b, psiConstant) {
-  const amplitude = psiConstant / (b.fibonacciD * Math.sqrt(b.massFrac));
-  const antiPhase = b.antiPhase ? -1 : 1;
-  const mean = b.invPlaneInclinationJ2000
-    - antiPhase * amplitude * Math.cos((b.longitudePerihelion - b.inclinationCycleAnchor) * Math.PI / 180);
-  return { amplitude, mean };
-}
-
-// (The device's wobble beat of integer axial and obliquity fractions, the
-// Venus/Neptune |ICRF| obliquity-cycle fallback and the snapshot "mean
-// obliquity" left with plan 06 Phase 7 commit 2: the K law's cycle period is
-// the chain's own g-mode beat (keplerian-chain computeSecularShape) and its
-// obliquity input the derived J2000 obliquity (spin-channel
-// computeObliquityJ2000Deg). docs/retired-record.md carries the record.)
-
-// Plan 07 R6: computeKConstant and computeEccentricityLaw stood here — the
-// K relation (K = e_amp·√m·a^1.5 / (sin(tilt)·√d), inverted per planet) and
-// the System-Reset eccentricity construction it fed
-// (base = A·cosθ + √(e_J2000² − A²sin²θ) at the J2000 phase, θ measured on
-// the wobble beat with a 90°/270° anti-phase offset).
-//
-// They are DELETED, not reclassified. Measured first (plan 07 §9e): scaling
-// their outputs by 1.10 moved no scene position and no published record
-// field — both twins place planets from the N-body chain, which overwrites
-// the construction entirely — and the registry keys built on them had no
-// consumer in the docs, on the website or in the paper. The nine balance and
-// Fibonacci-law instruments that computed from them are retired with them;
-// docs/retired-record.md is the record.
-//
-// A planet's eccentricity of date now has one home: the chain
-// (`model.planets.eccentricity(k, year)`), and its long-term shape is the
-// chain's own g-mode structure.
-
-module.exports = {
-  computePsiConstant, computeInclinationLaw,
-};