@essrt/physics 8.0.0 → 9.1.1

package/README.md

@@ -1,7 +1,7 @@
 # @essrt/physics
 
 The physics core of **ESSRT** — the Expanding Solar System Resonance Theory —
-as implemented by the [Holistic Universe Model](https://3d.holisticuniverse.com)
+as implemented by the [ESSRT simulator](https://3d.holisticuniverse.com)
 ([source, AGPL-3.0](https://github.com/dvansonsbeek/3d)). Pure computation:
 no I/O, no globals, no DOM, constants injected. The package ships the
 **complete model**, including the fitted coefficients, so an installed copy

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "8.0.0",
+  "version": "9.1.1",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core — pure, injectable-constants, no I/O, no globals, no DOM. The complete model incl. fitted coefficients; every published version immutably ships one recorded model identity.",
   "essrt": {
     "modelVersion": "v14.0",
@@ -55,6 +55,8 @@
     "./planets/asc-node-integrator": { "types": "./types/planets/asc-node-integrator.d.cts", "default": "./src/planets/asc-node-integrator.cjs" },
     "./planets/orbit-chain": { "types": "./types/planets/orbit-chain.d.cts", "default": "./src/planets/orbit-chain.cjs" },
     "./planets/keplerian-chain": { "types": "./types/planets/keplerian-chain.d.cts", "default": "./src/planets/keplerian-chain.cjs" },
+    "./planets/spin-channel": { "types": "./types/planets/spin-channel.d.cts", "default": "./src/planets/spin-channel.cjs" },
+    "./planets/secular-shape": { "types": "./types/planets/secular-shape.d.cts", "default": "./src/planets/secular-shape.cjs" },
     "./planets/inv-plane-frame": { "types": "./types/planets/inv-plane-frame.d.cts", "default": "./src/planets/inv-plane-frame.cjs" },
     "./planets/secular-series": { "types": "./types/planets/secular-series.d.cts", "default": "./src/planets/secular-series.cjs" },
     "./planets/model": { "types": "./types/planets/model.d.cts", "default": "./src/planets/model.cjs" },

package/src/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: "2f0bb1d372d35ed9";
+export declare const CONSTANTS_HASH: "b9c4d78647577e67";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "2f0bb1d372d35ed9";
+  readonly hash: "b9c4d78647577e67";
   readonly additionalBodies: {
     "pluto": {
       "name": string;
@@ -372,6 +372,113 @@ export declare const DEFAULT_CONSTANTS: {
       "rotationPeriodDays": number;
     };
   };
+  readonly planetSpinPhysical: {
+    "mercury": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": unknown[];
+    };
+    "venus": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": unknown[];
+    };
+    "mars": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "jupiter": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "momentOfInertiaFactorAlternative": number;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "saturn": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "momentOfInertiaFactorAlternative": [number, number];
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "uranus": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "neptune": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+        "retrograde": boolean;
+      }>;
+    };
+  };
   readonly planets: {
     "mercury": {
       "name": string;
@@ -386,9 +493,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "venus": {
       "name": string;
@@ -403,9 +507,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": null;
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "mars": {
       "name": string;
@@ -420,9 +521,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "jupiter": {
       "name": string;
@@ -437,9 +535,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "saturn": {
       "name": string;
@@ -454,9 +549,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "uranus": {
       "name": string;
@@ -471,10 +563,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
-      "eccentricityCycleApprox": boolean;
     };
     "neptune": {
       "name": string;
@@ -489,10 +577,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": null;
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
-      "eccentricityCycleApprox": boolean;
     };
   };
   readonly timeReference: {

package/src/constants/generated.js

@@ -33,7 +33,7 @@
  * carry it so a counterfactual is reproducible (§2d).
  * @type {string}
  */
-export const CONSTANTS_HASH = "2f0bb1d372d35ed9";
+export const CONSTANTS_HASH = "b9c4d78647577e67";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
 export const MODEL_VERSION = "v14.0";
@@ -43,7 +43,7 @@ export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";
 
 /** @type {Readonly<Record<string, unknown>>} */
 export const DEFAULT_CONSTANTS = Object.freeze({
-  hash: "2f0bb1d372d35ed9",
+  hash: "b9c4d78647577e67",
   additionalBodies: {
     "pluto": {
       "name": "Pluto",
@@ -430,6 +430,216 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "rotationPeriodDays": 0.671300001591743
     }
   },
+  planetSpinPhysical: {
+    "mercury": {
+      "j2": 0.0000503,
+      "j2ReferenceRadiusKm": 2440,
+      "momentOfInertiaFactor": 0.346,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": 6.1385108,
+      "poleRaJ2000Deg": 281.0103,
+      "poleDecJ2000Deg": 61.4155,
+      "cassiniLocked": true,
+      "satellites": []
+    },
+    "venus": {
+      "j2": 0.0000044044,
+      "j2ReferenceRadiusKm": 6051,
+      "momentOfInertiaFactor": 0.337,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": -1.4813688,
+      "poleRaJ2000Deg": 272.76,
+      "poleDecJ2000Deg": 67.16,
+      "cassiniLocked": false,
+      "satellites": []
+    },
+    "mars": {
+      "j2": 0.001956608,
+      "j2ReferenceRadiusKm": 3396,
+      "momentOfInertiaFactor": 0.3644,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": 350.891982443297,
+      "poleRaJ2000Deg": 317.68143,
+      "poleDecJ2000Deg": 52.8865,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Phobos",
+          "gmKm3S2": 0.0007087,
+          "semiMajorAxisKm": 9375,
+          "orbitalPeriodDays": 0.3187
+        },
+        {
+          "name": "Deimos",
+          "gmKm3S2": 0.0000962,
+          "semiMajorAxisKm": 23457,
+          "orbitalPeriodDays": 1.2625
+        }
+      ]
+    },
+    "jupiter": {
+      "j2": 0.014696572,
+      "j2ReferenceRadiusKm": 71492,
+      "momentOfInertiaFactor": 0.26393,
+      "momentOfInertiaFactorClass": "gravity-constrained",
+      "momentOfInertiaFactorAlternative": 0.2756,
+      "rotationRateDegPerDay": 870.536,
+      "poleRaJ2000Deg": 268.056595,
+      "poleDecJ2000Deg": 64.495303,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Io",
+          "gmKm3S2": 5959.91547,
+          "semiMajorAxisKm": 421800,
+          "orbitalPeriodDays": 1.762732
+        },
+        {
+          "name": "Europa",
+          "gmKm3S2": 3202.7121,
+          "semiMajorAxisKm": 671100,
+          "orbitalPeriodDays": 3.525463
+        },
+        {
+          "name": "Ganymede",
+          "gmKm3S2": 9887.83275,
+          "semiMajorAxisKm": 1070400,
+          "orbitalPeriodDays": 7.155588
+        },
+        {
+          "name": "Callisto",
+          "gmKm3S2": 7179.2834,
+          "semiMajorAxisKm": 1882700,
+          "orbitalPeriodDays": 16.69044
+        }
+      ]
+    },
+    "saturn": {
+      "j2": 0.016290573,
+      "j2ReferenceRadiusKm": 60330,
+      "momentOfInertiaFactor": 0.2181,
+      "momentOfInertiaFactorClass": "gravity-constrained",
+      "momentOfInertiaFactorAlternative": [
+        0.2257,
+        0.2438
+      ],
+      "rotationRateDegPerDay": 810.7939024,
+      "poleRaJ2000Deg": 40.589,
+      "poleDecJ2000Deg": 83.537,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Mimas",
+          "gmKm3S2": 2.50349,
+          "semiMajorAxisKm": 186000,
+          "orbitalPeriodDays": 0.942422
+        },
+        {
+          "name": "Enceladus",
+          "gmKm3S2": 7.21037,
+          "semiMajorAxisKm": 238400,
+          "orbitalPeriodDays": 1.370218
+        },
+        {
+          "name": "Tethys",
+          "gmKm3S2": 41.21353,
+          "semiMajorAxisKm": 295000,
+          "orbitalPeriodDays": 1.887802
+        },
+        {
+          "name": "Dione",
+          "gmKm3S2": 73.11607,
+          "semiMajorAxisKm": 377700,
+          "orbitalPeriodDays": 2.736916
+        },
+        {
+          "name": "Rhea",
+          "gmKm3S2": 153.94175,
+          "semiMajorAxisKm": 527200,
+          "orbitalPeriodDays": 4.517503
+        },
+        {
+          "name": "Titan",
+          "gmKm3S2": 8978.1371,
+          "semiMajorAxisKm": 1221900,
+          "orbitalPeriodDays": 15.945448
+        },
+        {
+          "name": "Iapetus",
+          "gmKm3S2": 120.51511,
+          "semiMajorAxisKm": 3561700,
+          "orbitalPeriodDays": 79.331002
+        }
+      ]
+    },
+    "uranus": {
+      "j2": 0.00351068,
+      "j2ReferenceRadiusKm": 25559,
+      "momentOfInertiaFactor": 0.225,
+      "momentOfInertiaFactorClass": "interior-model",
+      "rotationRateDegPerDay": -501.1600928,
+      "poleRaJ2000Deg": 257.311,
+      "poleDecJ2000Deg": -15.175,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Miranda",
+          "gmKm3S2": 4.3,
+          "semiMajorAxisKm": 129846,
+          "orbitalPeriodDays": 1.413479
+        },
+        {
+          "name": "Ariel",
+          "gmKm3S2": 83.5,
+          "semiMajorAxisKm": 190929,
+          "orbitalPeriodDays": 2.520379
+        },
+        {
+          "name": "Umbriel",
+          "gmKm3S2": 85.1,
+          "semiMajorAxisKm": 265986,
+          "orbitalPeriodDays": 4.144177
+        },
+        {
+          "name": "Titania",
+          "gmKm3S2": 226.9,
+          "semiMajorAxisKm": 436298,
+          "orbitalPeriodDays": 8.705869
+        },
+        {
+          "name": "Oberon",
+          "gmKm3S2": 205.3,
+          "semiMajorAxisKm": 583511,
+          "orbitalPeriodDays": 13.463237
+        }
+      ]
+    },
+    "neptune": {
+      "j2": 0.00340843,
+      "j2ReferenceRadiusKm": 25225,
+      "momentOfInertiaFactor": 0.23,
+      "momentOfInertiaFactorClass": "interior-model",
+      "rotationRateDegPerDay": 541.1397757,
+      "poleRaJ2000Deg": 299.36,
+      "poleDecJ2000Deg": 43.46,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Triton",
+          "gmKm3S2": 1428.49546,
+          "semiMajorAxisKm": 354800,
+          "orbitalPeriodDays": 5.876994,
+          "retrograde": true
+        },
+        {
+          "name": "Proteus",
+          "gmKm3S2": 2.58342,
+          "semiMajorAxisKm": 117600,
+          "orbitalPeriodDays": 1.122315
+        }
+      ]
+    }
+  },
   planets: {
     "mercury": {
       "name": "Mercury",
@@ -446,19 +656,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 32.83,
       "inclinationCycleAnchor": 234.52,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 9,
-      "obliquityCycleFraction": [
-        8,
-        3
-      ],
-      "axialPrecessionFraction": [
-        -8,
-        9
-      ],
-      "eccentricityCycleFraction": [
-        2,
-        21
-      ]
+      "ascendingNodeCyclesIn8H": 9
     },
     "venus": {
       "name": "Venus",
@@ -475,16 +673,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 54.7,
       "inclinationCycleAnchor": 218.64,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 1,
-      "obliquityCycleFraction": null,
-      "axialPrecessionFraction": [
-        8,
-        91
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        19
-      ]
+      "ascendingNodeCyclesIn8H": 1
     },
     "mars": {
       "name": "Mars",
@@ -501,19 +690,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 354.87,
       "inclinationCycleAnchor": 236.07,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 64,
-      "obliquityCycleFraction": [
-        8,
-        21
-      ],
-      "axialPrecessionFraction": [
-        -1,
-        2
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        52
-      ]
+      "ascendingNodeCyclesIn8H": 64
     },
     "jupiter": {
       "name": "Jupiter",
@@ -530,19 +707,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 312.89,
       "inclinationCycleAnchor": 287.06,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 36,
-      "obliquityCycleFraction": [
-        1,
-        2
-      ],
-      "axialPrecessionFraction": [
-        -8,
-        21
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        44
-      ]
+      "ascendingNodeCyclesIn8H": 36
     },
     "saturn": {
       "name": "Saturn",
@@ -559,19 +724,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 118.81,
       "inclinationCycleAnchor": 116.26,
       "antiPhase": true,
-      "ascendingNodeCyclesIn8H": 36,
-      "obliquityCycleFraction": [
-        1,
-        3
-      ],
-      "axialPrecessionFraction": [
-        -4,
-        3
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        163
-      ]
+      "ascendingNodeCyclesIn8H": 36
     },
     "uranus": {
       "name": "Uranus",
@@ -588,20 +741,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 307.8,
       "inclinationCycleAnchor": 21.33,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 11,
-      "obliquityCycleFraction": [
-        1,
-        2
-      ],
-      "axialPrecessionFraction": [
-        610,
-        1
-      ],
-      "eccentricityCycleFraction": [
-        1,
-        10
-      ],
-      "eccentricityCycleApprox": true
+      "ascendingNodeCyclesIn8H": 11
     },
     "neptune": {
       "name": "Neptune",
@@ -618,17 +758,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 192.04,
       "inclinationCycleAnchor": 174.04,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 3,
-      "obliquityCycleFraction": null,
-      "axialPrecessionFraction": [
-        -68,
-        1
-      ],
-      "eccentricityCycleFraction": [
-        2,
-        25
-      ],
-      "eccentricityCycleApprox": true
+      "ascendingNodeCyclesIn8H": 3
     }
   },
   timeReference: {

package/src/deltat/deep-time.cjs

@@ -125,7 +125,15 @@ function createDeepTimeLod(deps) {
    * the unit H(t) below.
    * @param {number} t_Ma @returns {number|null} */
   function eraClockHAtAge(t_Ma) {
-    const LOD_s = lodSecondsAtAge(t_Ma);
+    return eraClockHFromLod(lodSecondsAtAge(t_Ma));
+  }
+  // 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;
   }
@@ -156,7 +164,11 @@ function createDeepTimeLod(deps) {
    * ≡ H_era wherever (a₀/a_M)³ ≈ 1. Identifier kept (plan 06 P4).
    * @param {number} t_Ma @returns {number|null} */
   function hAtAge(t_Ma) {
-    const hEra = eraClockHAtAge(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;
   }
@@ -202,8 +214,12 @@ function createDeepTimeLod(deps) {
    * ±2 Myr, measured), and re-basing would move every comb coefficient
    * for a rounding-level gain. This tier stays a named device. */
   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 = hAtAge(t_Ma);
+    const Ht = hFromLod(t_Ma, LOD_s);
     if (Ht === null) return sidSec * (1 - 13 / K.holisticYearJ2000);
     return sidSec * (1 - 13 / Ht);
   }
@@ -281,15 +297,16 @@ function createDeepTimeLod(deps) {
       const tau = i * h;
       const lodMean = lodSecondsAtAge(tau);
       if (lodMean === null) return NaN;
-      const yearS = tropicalYearSecondsAtAge(tau);
+      const yearS = tropicalYearSecondsFromLod(tau, lodMean);
       // Ecliptic "missing motion" — the solar day is measured against the
       // Sun on the ECLIPTIC, whose plane turns on the invariable plane at the
       // nodal rate s₃ (period ≈ 68,751 yr, an orbital quantity), not in the
       // inclination frame. Adds ~3.4 ms at J2000; the fitted cycle stack
       // closes Layer-4 LOD_real onto the USNO anchor. Non-null: the helper
       // can only be null when lodSecondsAtAge(tau) is null, which already
-      // returned NaN above — the checker can't see the chain.
-      const lodH5Raw = lodMean + /** @type {number} */ (eclipticLodCorrectionSecondsAtAge(tau));
+      // returned NaN above. Spelled inline from the node's own lodMean and
+      // yearS — the same operations as eclipticLodCorrectionSecondsAtAge(tau).
+      const lodH5Raw = lodMean + lodMean / (deps.nodalPeriodYearsFn() * (yearS / 86400));
       const integrand = (86400 - lodH5Raw) * yearS * 1e6 / 86400;
       const w = (i === 0 || i === n) ? 1 : (i % 2 === 1 ? 4 : 2);
       sum += w * integrand;

package/src/index.js

@@ -127,9 +127,13 @@ export { integrateAscendingNode } from './planets/asc-node-integrator.cjs';
 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';
 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).
+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).