@essrt/physics 20.0.0 → 21.0.0

package/src/planets/geometry.cjs

@@ -8,10 +8,19 @@
  * type-branched function over a body record — the shape the model's claim
  * has always had (six laws, per-body parameters).
  *
- * EXPRESSION FORMS ARE THE BROWSER'S (the certified engine) — e.g.
- * `Math.round((H·mSY)/input)`, `((H/N)**2)**(1/3)` — so browser bit-identity
- * holds by construction; Node's last-ulp variants dissolve here (the
- * cardinal-extraction pattern).
+ * EXPRESSION FORMS ARE THE BROWSER'S (the certified engine) — so browser
+ * bit-identity holds by construction; Node's last-ulp variants dissolve here
+ * (the cardinal-extraction pattern).
+ *
+ * PLAN 09 STEP 7b: the device count is gone. The geometry used to quantize
+ * every body's period onto the unit — N = round(H·mSY/input) orbits per
+ * anchor interval, period H/N, a = (H/N)^(2/3) — a rounding of up to half a
+ * count (2.5e-4 of Neptune's period). The caller now hands in the period in
+ * SI days directly (the chain's J2000 window mean motion for the seven chain
+ * planets, the observed input for the no-chain bodies) and, where a chain
+ * exists, the chain's J2000 semi-major axis; Kepler's a from the period is the
+ * fallback for the no-chain bodies. Rotation periods are the inputs (they were
+ * rounded onto the count too), Mercury's the 3:2 lock on its orbital period.
  *
  * Type classification (the survey's measured taxonomy):
  *  - Type I  (mercury, venus): peri = a·e_real·100, elip = peri/2
@@ -40,7 +49,8 @@
  * @typedef {{
  *   key: string,
  *   type?: string,
- *   solarYearInput: number,
+ *   periodDays: number,
+ *   semiMajorAxisAu?: number,
  *   orbitalEccentricityBase?: number,
  *   longitudePerihelion?: number,
  *   ascendingNode?: number,
@@ -50,19 +60,19 @@
  */
 
 /**
- * @param {PlanetGeometryBody} body
+ * @param {PlanetGeometryBody} body — `periodDays` the body's sidereal orbital
+ *   period in SI days; `semiMajorAxisAu` the chain's J2000 a where a chain
+ *   exists (else Kepler's a from the period, in the scene's solar-year units)
  * @param {{
- *   holisticYears: number,
  *   meanSolarYearDays: number,
  *   currentAUDistanceKm: number,
  *   earthEccentricityJ2000: number,
  *   earthPerihelionLongitudeJ2000Deg: number,
- * }} env — the engine's J2000 values (browser passes its live globals at
- *   load time, preserving its original initial-derivation semantics; the
- *   epoch machinery that later mutates the browser's `let` aliases is
- *   untouched by this layer).
+ * }} env — the engine's J2000 values; `meanSolarYearDays` is the scene's
+ *   year unit (`pos` counts mean solar years), so `periodYears` is the period
+ *   in that unit, the scaffold's `speed` = 2π/periodYears.
  * @returns {{
- *   solarYearCount: number, orbitDistance: number, periodYears: number,
+ *   orbitDistance: number, periodYears: number, periodDays: number,
  *   realOrbitalEccentricity: (number | undefined),
  *   elipticOrbit: (number | undefined), perihelionDistance: (number | undefined),
  *   speedKmh: number, rotationPeriodHours: (number | undefined),
@@ -70,17 +80,18 @@
  * }}
  */
 function derivePlanetGeometry(body, env) {
-  const { key, type, solarYearInput, orbitalEccentricityBase: base,
+  const { key, type, periodDays, semiMajorAxisAu, orbitalEccentricityBase: base,
           longitudePerihelion, ascendingNode, rotationPeriodDays, orbitDistanceOverride } = body;
-  const { holisticYears: H, meanSolarYearDays: mSY, currentAUDistanceKm: AU,
+  const { meanSolarYearDays: mSY, currentAUDistanceKm: AU,
           earthEccentricityJ2000, earthPerihelionLongitudeJ2000Deg } = env;
 
-  const solarYearCount = Math.round((H * mSY) / solarYearInput);
+  const periodYears = periodDays / mSY;   // in the scene's year unit (mean solar years)
   const orbitDistance = orbitDistanceOverride !== undefined
     ? orbitDistanceOverride
-    : ((H / solarYearCount) ** 2) ** (1 / 3);
-  const periodYears = H / solarYearCount;
-  const speedKmh = (orbitDistance * AU * Math.PI * 2) / (mSY * (H / solarYearCount)) / 24;
+    : semiMajorAxisAu !== undefined
+      ? semiMajorAxisAu
+      : (periodYears ** 2) ** (1 / 3);
+  const speedKmh = (orbitDistance * AU * Math.PI * 2) / periodDays / 24;
 
   /** @type {number | undefined} */
   let realOrbitalEccentricity;
@@ -127,19 +138,37 @@ function derivePlanetGeometry(body, env) {
   /** @type {number | undefined} */
   let rotationPeriodHours;
   if (key === 'mercury') {
-    // 3:2 spin-orbit lock — period from the orbit count itself.
-    rotationPeriodHours = 24 * (mSY * H) / (solarYearCount * 3 / 2);
+    // 3:2 spin-orbit lock — three sidereal rotations per two orbits.
+    rotationPeriodHours = 24 * periodDays / (3 / 2);
   } else if (rotationPeriodDays !== undefined) {
-    rotationPeriodHours = 24 * (mSY * H) / Math.round((mSY * H) / rotationPeriodDays);
+    rotationPeriodHours = 24 * rotationPeriodDays;
   }
 
   const lowestPoint = ascendingNode !== undefined ? 180 - ascendingNode : undefined;
 
   return {
-    solarYearCount, orbitDistance, periodYears, realOrbitalEccentricity,
+    orbitDistance, periodYears, periodDays, realOrbitalEccentricity,
     elipticOrbit, perihelionDistance, speedKmh, rotationPeriodHours,
     eccentricityPerihelion, lowestPoint,
   };
 }
 
-module.exports = { derivePlanetGeometry };
+/**
+ * The scaffold's orbit inputs for one body — ONE home for the rule (plan 09
+ * step 7b): where the governed engine-D chain artifact has the body, the
+ * chain's J2000 semi-major axis and its fitted era-window mean motion (deg
+ * per Julian year → the period in SI days); otherwise the body's observed
+ * period input (the no-chain bodies Pluto, Halley, Eros, Ceres).
+ * @param {string} key
+ * @param {any} chainArtifact the CHAIN_ARTIFACT (or the governed JSON it embeds)
+ * @param {number} fallbackPeriodDays the observed period input, SI days
+ * @returns {{ periodDays: number, semiMajorAxisAu: (number | undefined) }}
+ */
+function scaffoldOrbitInputs(key, chainArtifact, fallbackPeriodDays) {
+  const w = chainArtifact && chainArtifact.windowElementRates && chainArtifact.windowElementRates[key];
+  const a = chainArtifact && chainArtifact.j2000AnchorElements && chainArtifact.j2000AnchorElements[key];
+  if (w && a) return { periodDays: 360 / w.meanMotionDegPerYr * 365.25, semiMajorAxisAu: a.aAU };
+  return { periodDays: fallbackPeriodDays, semiMajorAxisAu: undefined };
+}
+
+module.exports = { derivePlanetGeometry, scaffoldOrbitInputs };

package/src/planets/model.cjs

@@ -31,12 +31,12 @@
 
 'use strict';
 
-const { derivePlanetGeometry } = require('./geometry.cjs');
+const { derivePlanetGeometry, scaffoldOrbitInputs } = require('./geometry.cjs');
 
 /**
  * @typedef {Object} PlanetModelBody
  * @property {string} [type] - geometry type branch ('I' | 'II' | 'III')
- * @property {number} solarYearInput
+ * @property {number} solarYearInput - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]
  * @property {number} [wobblePeriodYears] - carriers: the chain's g-mode beat (input since Phase 7 commit 2)
@@ -52,7 +52,7 @@ const { derivePlanetGeometry } = require('./geometry.cjs');
 
 /**
  * @typedef {Object} PlanetModelEnv
- * @property {number} holisticYears
+ * @property {any} [chainArtifact] - the governed engine-D chain artifact (CHAIN_ARTIFACT): the seven chain planets' scaffold a and period ride its J2000 elements (plan 09 step 7b)
  * @property {number} meanSolarYearDays
  * @property {number} balancedYear
  * @property {number} systemResetN
@@ -84,7 +84,6 @@ function createPlanetModel(env, bodies) {
   // from env.calibration) went with the inclination law.
 
   const geomEnv = {
-    holisticYears: env.holisticYears,
     meanSolarYearDays: env.meanSolarYearDays,
     currentAUDistanceKm: env.currentAUDistanceKm,
     earthEccentricityJ2000: env.earthEccentricityJ2000,
@@ -114,7 +113,8 @@ function createPlanetModel(env, bodies) {
 
     rec.geometry = derivePlanetGeometry({
       key, type: b.type,
-      solarYearInput: b.solarYearInput,
+      // plan 09 step 7b: the chain's J2000 elements where the artifact has the body, else the observed period input
+      ...scaffoldOrbitInputs(key, env.chainArtifact, b.solarYearInput),
       orbitalEccentricityBase: b.orbitalEccentricityBase !== undefined
         ? b.orbitalEccentricityBase : b.orbitalEccentricityJ2000,
       longitudePerihelion: b.longitudePerihelion,

package/types/chain-cycles/index.d.cts

@@ -5,7 +5,6 @@ export type PeriodSecondsAtAge = (tMa: number) => (number | null);
 /**
  * @param {{
  *   ageAnchorYear: number,
- *   tropicalYearSecondsAtAge: (tMa: number) => (number | null),
  *   tropicalYearJ2000Seconds: number,
  *   isDeepTime: () => boolean,
  *   gridMinYear?: number,
@@ -15,12 +14,18 @@ export type PeriodSecondsAtAge = (tMa: number) => (number | null);
  *   maxCacheEntries?: number,
  * }} deps — ageAnchorYear is the year whose age is 0 in the period
  *   functions' t_Ma coordinate (startmodelYear = 2000.5, the certified
- *   convention); tropicalYearJ2000Seconds feeds the snapshot branch;
+ *   convention); tropicalYearJ2000Seconds is the length of ONE tick of the
+ *   year axis the integrator is fed — the linear SI-year label, each tick
+ *   exactly one J2000 tropical year of seconds — so it is the integrand's
+ *   weight in BOTH branches (plan 09 step 7c′: the deep branch had weighted
+ *   each tick by the tropical year OF DATE, adding that year's deep-time
+ *   drift to every lunar count — ~5° of Moon at ±100 kyr, ~1.3 orbits at
+ *   1 Ma — although the tick itself never changes length; the owner's call:
+ *   a plain correctness fix, the tick already is an SI tick);
  *   isDeepTime is read PER CALL so the engine's runtime toggle works.
  */
-export function createChainCycleIntegrator({ ageAnchorYear, tropicalYearSecondsAtAge, tropicalYearJ2000Seconds, isDeepTime, gridMinYear, gridMaxYear, gridStepYears, gridAnchorYear, maxCacheEntries, }: {
+export function createChainCycleIntegrator({ ageAnchorYear, tropicalYearJ2000Seconds, isDeepTime, gridMinYear, gridMaxYear, gridStepYears, gridAnchorYear, maxCacheEntries, }: {
     ageAnchorYear: number;
-    tropicalYearSecondsAtAge: (tMa: number) => (number | null);
     tropicalYearJ2000Seconds: number;
     isDeepTime: () => boolean;
     gridMinYear?: number;

package/types/constants/coefficients.d.ts

@@ -2,73 +2,10 @@
 // The fitted coefficients' shapes for TypeScript consumers (§2g); values live
 // in coefficients.js, emitted VERBATIM from fitted-coefficients.json.
 
-export declare const COEFFICIENTS_HASH: "49147a1d09bdaabc";
+export declare const COEFFICIENTS_HASH: "91e2d041c7b6e038";
 
 export declare const FITTED_COEFFICIENTS: {
-  readonly ANOMALISTIC_YEAR_HARMONICS: Array<[number, number, number]>;
   readonly BALANCE_PRESETS: Array<Array<number | string>>;
-  readonly CARDINAL_POINT_ANCHORS: {
-    "SS": number;
-    "WS": number;
-    "VE": number;
-    "AE": number;
-  };
-  readonly CARDINAL_POINT_ANCHORS_ADJUSTED: {
-    "SS": number;
-    "WS": number;
-    "VE": number;
-    "AE": number;
-  };
-  readonly CARDINAL_POINT_DERIVED: {
-    "lincoef": number;
-    "h0": number;
-    "h1": number;
-    "eccOrderDivisors": [number, number];
-    "note": string;
-  };
-  readonly CARDINAL_POINT_ECC_TERMS: {
-    "SS": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "WS": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "VE": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "AE": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-  };
-  readonly CARDINAL_POINT_HARMONICS: {
-    "SS": Array<[number, number, number]>;
-    "WS": Array<[number, number, number]>;
-    "VE": Array<[number, number, number]>;
-    "AE": Array<[number, number, number]>;
-  };
-  readonly CARDINAL_POINT_JOINT_TERMS: {
-    "quadratureDeg": {
-      "SS": number;
-      "AE": number;
-      "WS": number;
-      "VE": number;
-    };
-    "terms": Array<{
-      "order": number;
-      "div": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "note": string;
-  };
   readonly CLIMATE_FORMULA_COEFFS: {
     "config": {
       "L1_lines": Array<{
@@ -303,105 +240,105 @@ export declare const FITTED_COEFFICIENTS: {
           "delta_l3": number;
         };
         "carbon_amplification_ratios": {
-          "19.0014": {
+          "19.0015": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "19.1554": {
+          "19.1555": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "22.4452": {
+          "22.4453": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "23.1931": {
+          "23.1933": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "23.7599": {
+          "23.7601": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "27.2689": {
+          "27.2691": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "28.8424": {
+          "28.8426": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.0047": {
+          "29.0050": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.6172": {
+          "29.6174": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.8091": {
+          "29.8094": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "39.5560": {
+          "39.5564": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.2292": {
+          "40.2297": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.9239": {
+          "40.9244": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "41.6445": {
+          "41.6450": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "53.6174": {
+          "53.6182": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
@@ -634,17 +571,4 @@ export declare const FITTED_COEFFICIENTS: {
   };
   readonly MEEUS_LATITUDE_TERMS: Array<number[]>;
   readonly MEEUS_LONGITUDE_TERMS: Array<number[]>;
-  readonly PERI_HARMONICS_RAW: Array<[number, number, number]>;
-  readonly PERI_OFFSET: number;
-  readonly RA_DAY_OFFSET_MS: {
-    "mean": number;
-    "terms": Array<[number, number, number]>;
-    "note": string;
-  };
-  readonly SIDEREAL_YEAR_HARMONICS: Array<[number, number, number]>;
-  readonly SOLSTICE_OBLIQUITY_HARMONICS: Array<[number, number, number]>;
-  readonly SOLSTICE_OBLIQUITY_MEAN_FITTED: number;
-  readonly SUN_LONGITUDE_HARMONICS: Array<[number, number, number]>;
-  readonly SUN_LONGITUDE_MEAN: number;
-  readonly TROPICAL_YEAR_HARMONICS: Array<[number, number, number]>;
 };

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: "6c2beb70bb3d5000";
+export declare const CONSTANTS_HASH: "c161ec19e5984b9e";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "6c2beb70bb3d5000";
+  readonly hash: "c161ec19e5984b9e";
   readonly additionalBodies: {
     "pluto": {
       "name": string;
@@ -148,6 +148,14 @@ export declare const DEFAULT_CONSTANTS: {
     "earthInvPlaneInclinationAmplitude": number;
     "eccentricityBase": number;
     "eccentricityAmplitude": number;
+    "kDeviceInclinationPhaseJ2000Deg": number;
+    "kDeviceWheelPeriodYears": {
+      "3": number;
+      "5": number;
+      "8": number;
+      "13": number;
+      "16": number;
+    };
     "kDeviceWheelPhaseAtJ2000Cycles": number;
   };
   readonly earthOrbital: {
@@ -156,13 +164,13 @@ export declare const DEFAULT_CONSTANTS: {
     "juneSolstice2000_JD": number;
     "obliquityJ2000_deg": number;
     "obliquityRate_arcsecPerCentury": number;
+    "generalPrecessionJ2000_arcsecPerYr": number;
     "earthEccentricityJ2000": number;
     "earthEccentricityDotJ2000": number;
     "earthEccentricityDotDotJ2000": number;
     "earthPerihelionLongitudeJ2000": number;
     "sunMeanLongitudeJ2000_deg": number;
     "earthAscendingNodeInvPlane": number;
-    "earthInclinationCycleAnchor": number;
     "perihelionPassageJ2000_JD": number;
     "earthInclinationJ2000_deg": number;
     "earthInclinationRate_arcsecPerCentury": number;
@@ -170,9 +178,6 @@ export declare const DEFAULT_CONSTANTS: {
     "sunTilt": number;
   };
   readonly foundational: {
-    "holisticyearLength": number;
-    "stepYears": number;
-    "inputmeanlengthsolaryearindays": number;
     "startmodelJD": number;
     "startmodelYear": number;
     "correctionDays": number;
@@ -480,7 +485,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "venus": {
       "name": string;
@@ -490,7 +494,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "mars": {
       "name": string;
@@ -500,7 +503,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "jupiter": {
       "name": string;
@@ -510,7 +512,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "saturn": {
       "name": string;
@@ -520,7 +521,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "uranus": {
       "name": string;
@@ -530,7 +530,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
     "neptune": {
       "name": string;
@@ -540,7 +539,6 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "mirrorPair": string;
       "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
     };
   };
   readonly timeReference: {
@@ -556,10 +554,7 @@ export declare const DEFAULT_CONSTANTS: {
     "tropicalYearSS": number;
     "tropicalYearAE": number;
     "tropicalYearWS": number;
-    "tropicalYearMean": number;
-    "anomalisticYear": number;
     "siderealYear": number;
-    "iauPrecessionJ2000": number;
     "tropicalYearRateSecPerCentury": number;
     "solarDay": number;
     "siderealDay": number;
@@ -716,6 +711,7 @@ export declare const REFERENCE_DATA: {
     "meeusPeriLong2000AD": number;
     "meeusPeriLong2500AD": number;
     "meeusPeriLong3000AD": number;
+    "simon1994AnomalisticYearDays": number;
     "mainstreamAxialPrecKyr": number;
     "mainstreamPeriPrecKyr": number;
     "mainstreamObliqCycleKyr": number;

package/types/deltat/deep-time.d.cts

@@ -25,17 +25,12 @@ export type DeepTimeLodConstants = {
     alpha1PerMa: number;
     alpha3PerMa3: number;
     alpha4PerMa4: number;
-    holisticYearJ2000: number;
     /**
-     * - the H/13-identity mean LOD at J2000
+     * - ω₀, the J2000 day length (the SI day since plan 09; identifier kept)
      */
     lodNowH13Seconds: number;
     meanSiderealYearJ2000Seconds: number;
     solarMassLossFracPerYear: number;
-    /**
-     * - Actual-LOD numerator
-     */
-    siderealYearDaysKinematicJ2000: number;
     /**
      * - f_S, the solar fraction of the J2000
      * precession torque (derive-params) — the unit H(t) = 13·T_p,composed needs it
@@ -49,7 +44,6 @@ export type DeepTimeLodDeps = {
      * - GIA α(t)
      */
     moiFactorAtAge: (tMa: number) => number;
-    siderealYearDaysFourierAt: (year: number) => number;
     /**
      * - gated δLOD sum (incl. swing)
      */
@@ -101,11 +95,9 @@ export type DeepTimeLodDeps = {
  * @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
  */
@@ -114,7 +106,6 @@ export type DeepTimeLodDeps = {
  * @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
@@ -141,7 +132,6 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
     lodSecondsAtAge: (t_Ma: number) => number | null;
     lodSecondsAtAgeWithAlpha: (t_Ma: number, alpha: number) => number | null;
     lodHoursAtAge: (t_Ma: number) => number | null;
-    hAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
     lunarTorqueFactorAtAge: (tMa: number) => (number | null);
@@ -149,11 +139,11 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
     siderealYearSecondsAtAge: (t_Ma: number) => number;
     tropicalYearSecondsAtAge: (t_Ma: number) => number;
     tropicalYearDaysAtAge: (t_Ma: number) => number | null;
+    tropicalYearSecondsFromLod: (t_Ma: number, LOD_s: number | null) => number;
     eclipticLodCorrectionSecondsAtAge: (t_Ma: number) => number | null;
     yearInDaysAtAge: (t_Ma: number) => number | null;
     deltaTRawSecondsAtAge: (t_Ma: number) => number;
     lodSecondsWithCorrectionsAtAge: (t_Ma: number) => number | null;
-    lodSecondsActualAtAge: (t_Ma: number) => number | null;
     dLodDtDecompositionAtAge: (t_Ma: number) => {
         tidal: number | null;
         gia: number | null;
@@ -163,7 +153,4 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
         net_L3: number | null;
         net_L4: number | null;
     };
-    eraClockHAtAge: (t_Ma: number) => number | null;
-    eraClockTropicalYearSecondsAtAge: (t_Ma: number) => number;
-    eraClockYearInDaysAtAge: (t_Ma: number) => number | null;
 };

package/types/earth/year-lengths.d.cts

@@ -53,6 +53,18 @@ export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }:
     /** The planetary λ̇ ratio (diagnostic; ≡1 at J2000). */
     planetaryRelAtYear: (year: number) => number;
 }>;
+/** 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} */
+export function computeMeanTropicalYearJ2000Days(siderealYearDays: number, generalPrecessionArcsecPerJulianYr: number): number;
 /**
  * The published of-date window (plan 06 Phase 3 S2): inside |year − 2000| ≤
  * this span the one-family route is THE published year-length / precession

package/types/eclipse/sun-planetary-completion.d.cts

@@ -63,15 +63,6 @@ export function computeCarrierRatesDegPerCy({ meanSiderealYearDays, moonSidereal
  * @returns {[number, number]}
  */
 export function eccVectorOfDateEmbedded(bodyIndex: number, T: number): [number, number];
-/** 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. */
-export const PAIRED_SUN_HARMONICS_SHA256: "53530a28a60b3c4a";
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
  *  full-precision carrier rates (deg/cy TT) the tables pair with: since plan
  *  06 I3 the model's own banked J2000 sidereal mean motions for the planets