@essrt/physics 10.0.0 → 11.0.1

package/src/earth/sidereal-year-channel.cjs

@@ -18,10 +18,20 @@
 // J2000 node lamDotRel ≡ 1 by construction, so the caller's anchor is
 // preserved exactly.
 //
-// Scope: this channel feeds the model/chart surface for the sidereal year
-// OF DATE. The scene's sidereal frame stays on the certified H/13
-// identity — a 3e-8 relative frame change is invisible in the render and
-// would disturb the certified frozen-clock pair.
+// Scope: this channel is the sidereal year OF DATE everywhere — the
+// model/chart surfaces AND the scene: since plan 06 R1 the Sun's mean
+// longitude of date integrates the one-family tropical year (this channel
+// minus the movement's own general precession), so the scene's measured
+// sidereal angle reproduces this law (the scene-sampled gate,
+// tools/verify/scene-year-lengths.js). The frozen era clock's sidereal
+// year (mass-loss only, no planetary drift) is the certified era device,
+// not a member of this family.
+//
+// ESTIMATOR (2026-09): each banked node is the least-squares slope of a
+// quadratic fitted to the dump's unwrapped L over ±5 kyr — the former
+// 2-kyr boxcar of per-step rates telescoped to the difference of the
+// OSCULATING longitude at the window's two ends and left 0.1–0.3 s of
+// node noise (−266 ms at +4000 against DE441; +18…28″ of Sun by 6500).
 
 const { SIDEREAL_CHANNEL_ARTIFACT } = require('./sidereal-channel-artifact.cjs');
 

package/src/earth/year-lengths.cjs

@@ -10,81 +10,112 @@
 // pair's members is amplified ~20×/s by P = T_s/(T_s − T_t).
 //
 // The rules this factory owns:
-//   • ONE FAMILY: tropical = massLoss·(1 − p_geom/360) with p_geom the
-//     movement's own equinox rate (wobble included); sidereal = the D6
-//     λ̇ channel; anomalistic = the SECULAR mean-element construction on
-//     the chain's apsidal tangent (the injected fn — the same rate family
-//     the panel's Prec. cell shows; structure fallback when absent); every
-//     value λ̇-corrected coherently (rate form — beats invariant to
-//     second order).
-//   • ONE BASIS: SI seconds (days = seconds/86400 at the caller).
+//   • ONE FAMILY: sidereal = the D6 λ̇ channel; tropical = the sidereal
+//     year minus the movement's own GENERAL-PRECESSION rate of date (the
+//     broken-angle equinox longitude, wobble included), in RATE form with
+//     the sampler's Julian-year unit named; anomalistic = the tropical
+//     year minus the movement's own equinox-referenced apsidal rate of
+//     date (the same year-over-year reading — the scene's perihelion
+//     passages ARE this movement); every value λ̇-corrected coherently
+//     (rate form — beats invariant to second order).
+//   • ONE BASIS: SI seconds (days = seconds/86400 at the caller); periods
+//     in JULIAN years (the beat's duration ÷ 365.25 d — the unit the
+//     frame's precession rate is quoted in, so 360°/rate and the pair's
+//     beat are the SAME number; measured 2026-09: the dimensionless
+//     T_sid/(T_sid − T_trop) counts TROPICAL years, 2.1e-5 off — 0.55 yr
+//     on the axial period, 26 ms on the tropical year against the frame).
 //   • BEATS FROM THE PAIR: each precession period is computed from the
 //     two years THIS factory returns — recomputing a beat from the
 //     displayed/served values reproduces it exactly.
-//   • Scene MEASUREMENTS validate this family (sub-0.5 s, the report's
-//     measured columns); they are instruments, never mixed into beats.
+//   • THE SCENE IS THE REFERENCE (owner, 2026-09-30: "all calculations
+//     should match what we measure in the scene"): the scene's Sun
+//     integrates the tropical year below and turns with this movement's
+//     frame and apsidal line, so its measured sidereal, tropical and
+//     anomalistic years reproduce this family to the millisecond once the
+//     planetary completion terms are removed per event (the scene-sampled
+//     gate, tools/verify/scene-year-lengths.js). Before this restatement
+//     the scene read the law's sidereal year 26 ms short at J2000 and up
+//     to 1 s off at depth (the projected equinox rate), and the published
+//     anomalistic year rode the planet chain's secular tangent — seconds
+//     from the scene's own perihelion passages within a few millennia,
+//     395 s at +25 kyr where e passes 0.0034.
 //   • The frozen era clock stays the certified era device (opt-out /
 //     era-certification record), not a member of this family.
 
 const { createCardinalStructure } = require('../cardinal/one-source-structure.cjs');
 const { createSiderealYearChannel } = require('./sidereal-year-channel.cjs');
 
+/** SI seconds in one Julian year — the sampler's time unit (t = years from J2000 TT, 365.25 d). */
+const JULIAN_YEAR_SECONDS = 365.25 * 86400;
+
 /**
  * @param {{
- *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number},
+ *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number, generalPrecessionLonDeg: number},
  *   massLossSiderealSecondsAtYearFn: (year: number) => number,
- *   apsidalSecularDegPerYrFn?: (year: number) => number,
  * }} opts
- *   sampleAt: the one-source movement sampler, calendar-year keyed
- *   (createDeepOrbitalHistory build().at wrapped by the caller — series
- *   tier where the artifact is available, mode tail beyond).
+ *   sampleAt: the one-source movement sampler, keyed in the sampler's own
+ *   year (Julian years from J2000 TT, 2000 = J2000; createDeepOrbitalHistory
+ *   build().at wrapped by the caller — series tier where the artifact is
+ *   available, mode tail beyond).
  *   massLossSiderealSecondsAtYearFn: the caller's secular mass-loss
  *   sidereal-year law, SI seconds, IAU-anchored at J2000.
- *   apsidalSecularDegPerYrFn: the chain's SECULAR apsidal tangent of date
- *   (keplerian-chain computeApsidalSecularDegPerYr for Earth — deg/yr,
- *   J2000-ecliptic frame). When given, the anomalistic year is the
- *   mean-element construction on this smooth rate; when absent, it falls
- *   back to the cardinal structure's measured year-over-year interval
- *   (which carries the banked series' century-scale ϖ̇ wobble).
  */
-function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, apsidalSecularDegPerYrFn }) {
+function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }) {
   if (typeof sampleAt !== 'function' || typeof massLossSiderealSecondsAtYearFn !== 'function') {
     throw new Error('createYearLengths: sampleAt and massLossSiderealSecondsAtYearFn are required');
   }
   const w180 = (/** @type {number} */ d) => ((d + 540) % 360) - 180;
   const chan = createSiderealYearChannel({ massLossSiderealSecondsAtYearFn });
+  /** A year length from a longer year and a retrograde angular rate of the
+   *  reference point, RATE form: 1/T = 1/T_long + p/(360°·JY) — exact (the
+   *  linear T_long·(1 − p/360) drops (p/360)², 47 ms on the tropical year),
+   *  with p in degrees per JULIAN year (the sampler's unit; read as "per
+   *  tropical year" it was 26 ms short). @param {number} tLong @param {number} pDegPerJulianYr */
+  const shorterYear = (tLong, pDegPerJulianYr) => 1 / (1 / tLong + pDegPerJulianYr / (360 * JULIAN_YEAR_SECONDS));
+  /** The movement's own general precession of date, degrees per Julian
+   *  year, retrograde positive — the year-over-year advance of the
+   *  broken-angle equinox longitude (wobble included). @param {number} year */
+  const generalPrecessionDegPerJulianYr = (year) =>
+    w180(sampleAt(year - 0.5).generalPrecessionLonDeg - sampleAt(year + 0.5).generalPrecessionLonDeg);
+  /** The movement's own apsidal rate of date RELATIVE TO THE EQUINOX (the
+   *  climatic-precession rate, ~0.017°/yr at J2000), degrees per Julian
+   *  year, prograde positive. @param {number} year */
+  const apsidalVsEquinoxDegPerJulianYr = (year) =>
+    w180(sampleAt(year + 0.5).periOfDateDeg - sampleAt(year - 0.5).periOfDateDeg);
   // The mean tropical year OF DATE, raw (pre-λ̇): the movement's own
-  // equinox-rate form — carries the real precession-rate wobble, unlike
-  // the smooth secular mean (a wobble-free mean under a wobbling family
-  // was one of the duplicated implementations this factory retires).
-  const tropicalRawSeconds = (/** @type {number} */ year) => {
-    const p = w180(sampleAt(year - 0.5).equinoxLonJ2000Deg
-      - sampleAt(year + 0.5).equinoxLonJ2000Deg);          // retrograde → positive
-    return massLossSiderealSecondsAtYearFn(year) * (1 - p / 360);
-  };
+  // general-precession rate — carries the real precession-rate wobble,
+  // unlike the smooth secular mean (a wobble-free mean under a wobbling
+  // family was one of the duplicated implementations this factory retires).
+  const tropicalRawSeconds = (/** @type {number} */ year) =>
+    shorterYear(massLossSiderealSecondsAtYearFn(year), generalPrecessionDegPerJulianYr(year));
   const structure = createCardinalStructure({
     sampleAt,
     tropicalYearSecondsAtYearFn: tropicalRawSeconds,
   });
   const trop = (/** @type {number} */ year) => chan.correctedYearSeconds(year, tropicalRawSeconds(year));
   const sid = (/** @type {number} */ year) => chan.siderealYearSecondsAtYear(year);
-  // The anomalistic year OF DATE — the SECULAR (mean-element) construction
-  // when the chain tangent is injected: T_anom_raw = T_sid_raw·360/(360 − w)
-  // with w = dϖ/dt (deg/yr, J2000-ecliptic, the SAME rate family the panel's
-  // Prec. cell shows), then λ̇-corrected coherently. The displayed beat
-  // anom/(anom − sid) then equals the chain period 360/w EXACTLY at J2000
-  // (the rate-form correction leaves beats invariant), and the J2000 value
-  // lands +0.07 s from the published mean-element 365.259636 d. The former
-  // structure form (the year-over-year interval on the banked series) read
-  // the series' century-scale ϖ̇ wobble — ±0.015″/yr ≈ ±0.35 s, an alias at
-  // the series' 250-yr sampling — and sat +0.20 s above the smooth tangent
-  // at J2000 (owner-found: the panel beat read 111,491 where the Prec. cell
-  // read 111,570). The measured-interval form stays on the cardinal
-  // structure for callers that want the wobble.
-  const anom = apsidalSecularDegPerYrFn
-    ? (/** @type {number} */ year) => chan.correctedYearSeconds(year,
-        massLossSiderealSecondsAtYearFn(year) * 360 / (360 - apsidalSecularDegPerYrFn(year)))
-    : (/** @type {number} */ year) => chan.correctedYearSeconds(year, structure.anomalisticYearSeconds(year));
+  // The anomalistic year OF DATE: the tropical year lengthened by the
+  // movement's own apsidal advance against the equinox (a perihelion
+  // passage is a pure mean-anomaly event — no equation-of-centre term), in
+  // rate form, then λ̇-corrected coherently. This IS the scene: its apsidal
+  // wheel turns with periOfDateDeg, and the scene's perihelion-to-perihelion
+  // intervals read this form to 0.03–0.13 s at every sampled epoch. The
+  // former construction on the planet chain's secular apsidal tangent (the
+  // same family as the Perihelion Longitudes Prec. cell) agreed at J2000
+  // (+0.07 s) and parted from the scene by +3.3 s at 0 AD, −6.6 s at 6000,
+  // −14 s at 10,000 and −395 s at +25,000 (e = 0.0034, the apsidal line
+  // swings fast) — the chain tangent's dϖ̇/dt runs ~40 % steep against the
+  // secular mean elements, the series' reads them to ~1 s over ±4 kyr. The
+  // rate the beats and the Prec. cell quote is apsidalRateFixedDegPerJulianYr
+  // below — ONE family, so 360°/rate and anom/(anom − sid) are the same
+  // number. The year-over-year reading carries the banked series' own
+  // century-scale ripple (±0.015″/yr ≈ ±0.35 s; the cadence is named in the
+  // artifact) — it is the movement's, not smoothed away.
+  const anom = (/** @type {number} */ year) =>
+    chan.correctedYearSeconds(year, shorterYear(tropicalRawSeconds(year), -apsidalVsEquinoxDegPerJulianYr(year)));
+  /** A beat period in JULIAN years from two year lengths (SI s): the time for
+   *  the shorter year to gain one turn on the longer, ÷ 365.25 d. @param {number} tLong @param {number} tShort */
+  const beatJulianYears = (tLong, tShort) => tLong * tShort / ((tLong - tShort) * JULIAN_YEAR_SECONDS);
 
   return Object.freeze({
     /** Mean tropical year of date, SI seconds. @param {number} year */
@@ -93,12 +124,16 @@ function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, apsidalS
     siderealYearSecondsAtYear: sid,
     /** Anomalistic year of date, SI seconds. @param {number} year */
     anomalisticYearSecondsAtYear: anom,
-    /** Axial precession period of date (equinox regression beat), years. @param {number} y */
-    axialPrecessionYearsAtYear: (y) => { const s = sid(y), t = trop(y); return s / (s - t); },
-    /** Perihelion (apsidal-vs-equinox) precession period of date, years. @param {number} y */
-    perihelionPrecessionYearsAtYear: (y) => { const a = anom(y), t = trop(y); return a / (a - t); },
-    /** Apsidal precession period of date vs the fixed stars, years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
-    inclinationPrecessionYearsAtYear: (y) => { const a = anom(y), s = sid(y); return a / (a - s); },
+    /** Axial precession period of date (equinox regression beat), JULIAN years. @param {number} y */
+    axialPrecessionYearsAtYear: (y) => beatJulianYears(sid(y), trop(y)),
+    /** Perihelion (apsidal-vs-equinox) precession period of date, JULIAN years. @param {number} y */
+    perihelionPrecessionYearsAtYear: (y) => beatJulianYears(anom(y), trop(y)),
+    /** Apsidal precession period of date vs the fixed stars, JULIAN years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
+    inclinationPrecessionYearsAtYear: (y) => beatJulianYears(anom(y), sid(y)),
+    /** The movement's own general precession of date, ″ per Julian year (retrograde positive) — 360°/this ≡ axialPrecessionYearsAtYear to second order. @param {number} y */
+    generalPrecessionArcsecPerJulianYr: (y) => generalPrecessionDegPerJulianYr(y) * 3600,
+    /** The movement's own apsidal rate of date against the FIXED frame, ″ per Julian year (prograde positive) — the rate the Perihelion Longitudes Prec. cell shows for Earth; 360°/this ≡ inclinationPrecessionYearsAtYear exactly. @param {number} y */
+    apsidalRateFixedArcsecPerJulianYr: (y) => { const a = anom(y), s = sid(y); return 360 * JULIAN_YEAR_SECONDS * (1 / s - 1 / a) * 3600; },
     /** The per-cardinal layer (D6 rules: year lengths λ̇-corrected; offsets/spreads raw — μs-class/second-order there). */
     cardinal: Object.freeze({
       /** @param {number} y @param {'VE'|'SS'|'AE'|'WS'} type */

package/src/model.js

@@ -24,7 +24,6 @@ import { createPhaseMachinery } from './phase/index.cjs';
 import { createYearLengths, ONE_FAMILY_WINDOW_YEARS } from './earth/year-lengths.cjs';
 import { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 import { CHAIN_ARTIFACT } from './planets/chain-artifact.js';
-import { buildPlanetChainsFromArtifactData, computeApsidalSecularDegPerYr } from './planets/keplerian-chain.cjs';
 import { computeSecularShape } from './planets/secular-shape.cjs';
 import { createPlanetSpinChannelFromArtifacts, computeObliquityJ2000Deg } from './planets/spin-channel.cjs';
 import { createDeltaTCycles } from './deltat/cycles.cjs';
@@ -533,9 +532,6 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // path; the npm package alone has no 8 MB series), the movement runs the
   // SERIES tier exactly like the browser — killing the mode-vs-series
   // value split (measured: 6 s on the anomalistic year at J2000).
-  // The engine-D planet chains (ONE build; the one-family route's apsidal
-  // tangent and the lunisolar surface's n_aps read the same instance).
-  const kcChainsM = buildPlanetChainsFromArtifactData(CHAIN_ARTIFACT);
   // The ONE-SOURCE movement: the hybrid's ε(t) (plan 06 Phase 3 S3b — THE
   // published obliquity) and the one-family year lengths, on the same
   // per-tier sampler.
@@ -611,13 +607,12 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       }
       return samplers.get(key).at(t);
     };
-    // The anomalistic rides the chain's SECULAR apsidal tangent (the same
-    // rate family the panel's Prec. cell shows) — ONE helper, keplerian-chain.
-    const kcChains = kcChainsM;
+    // The anomalistic year rides THIS movement's own apsidal line (the
+    // scene's perihelion) — the factory owns the construction; the planet
+    // chain's secular tangent is no longer a member (year-lengths.cjs).
     const yearLengths = createYearLengths({
       sampleAt,
       massLossSiderealSecondsAtYearFn: (year) => deepLod.siderealYearSecondsAtAge(yearToTMa(year)),
-      apsidalSecularDegPerYrFn: (year) => computeApsidalSecularDegPerYr(year, kcChains.earth, kcChains),
     });
     return {
       yearLengths,
@@ -1505,23 +1500,23 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
         const psiDot = r === null ? 1296000 / certifiedAxialPrecessionJ2000Years() : r;
         return 1296000 / (psiDot - s3ArcsecPerYr);
       },
-      /** The apsidal (perihelion vs the stars) period from the engine-D chain's secular tangent, years — inside the published window only (the tangent is an extrapolation beyond the banked series: it turns negative at −5 Myr); null beyond. @param {number} year @returns {number|null} */
-      apsidalPeriodYearsAtYear: (year) => (Math.abs(year - 2000) <= ONE_FAMILY_WINDOW_YEARS ? 360 / computeApsidalSecularDegPerYr(year, kcChainsM.earth, kcChainsM) : null),
+      /** The apsidal (perihelion vs the stars) period of date, JULIAN years — the ONE movement's own apsidal rate (the scene's perihelion; anom/(anom − sid) of the one family, ≡ 360°/apsidalRateFixedArcsecPerJulianYr), inside the published window only; null beyond. Before 2026-09 this rode the planet chain's secular tangent, which parted from the scene's perihelion passages by seconds of anomalistic year within a few millennia (year-lengths.cjs). @param {number} year @returns {number|null} */
+      apsidalPeriodYearsAtYear: (year) => (Math.abs(year - 2000) <= ONE_FAMILY_WINDOW_YEARS ? yearLengthsM.inclinationPrecessionYearsAtYear(year) : null),
       /** T_aps(t) / T_p(t) — the apsidal period in of-date precession periods (4.33 at J2000, a reading; 0.84 … 9.9 across ±26 kyr, measured); null beyond the published window. @param {number} year @returns {number|null} */
       apsidalPerPrecessionAtYear: (year) => {
         if (Math.abs(year - 2000) > ONE_FAMILY_WINDOW_YEARS) return null;
-        return (360 / computeApsidalSecularDegPerYr(year, kcChainsM.earth, kcChainsM)) / publishedAxialPrecessionYearsAtYear(year);
+        return yearLengthsM.inclinationPrecessionYearsAtYear(year) / publishedAxialPrecessionYearsAtYear(year);
       },
-      /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), years, on the of-date T_p; null beyond the published window. @param {number} year @returns {number|null} */
+      /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), JULIAN years, on the of-date T_p — ≡ the one family's anom/(anom − trop); null beyond the published window. @param {number} year @returns {number|null} */
       periOfDatePeriodYearsAtYear: (year) => {
         if (Math.abs(year - 2000) > ONE_FAMILY_WINDOW_YEARS) return null;
-        const tp = publishedAxialPrecessionYearsAtYear(year), taps = 360 / computeApsidalSecularDegPerYr(year, kcChainsM.earth, kcChainsM);
+        const tp = publishedAxialPrecessionYearsAtYear(year), taps = yearLengthsM.inclinationPrecessionYearsAtYear(year);
         return 1 / (1 / tp + 1 / taps);
       },
       /** T_peri(t) / T_p(t) (0.812 at J2000 — the J2000 reading); null beyond the published window. @param {number} year @returns {number|null} */
       periOfDatePerPrecessionAtYear: (year) => {
         if (Math.abs(year - 2000) > ONE_FAMILY_WINDOW_YEARS) return null;
-        const tp = publishedAxialPrecessionYearsAtYear(year), taps = 360 / computeApsidalSecularDegPerYr(year, kcChainsM.earth, kcChainsM);
+        const tp = publishedAxialPrecessionYearsAtYear(year), taps = yearLengthsM.inclinationPrecessionYearsAtYear(year);
         return (1 / (1 / tp + 1 / taps)) / tp;
       },
       /** The published of-date window, years from 2000. */

package/types/earth/deep-orbital-history.d.cts

@@ -71,6 +71,7 @@ export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSerie
             epsDeg: number;
             e: number;
             periOfDateDeg: number;
+            periOfDateRateDegPerYr: number;
             eSinPeri: number;
             eCosPeri: number;
             inclEclDeg: number;
@@ -78,6 +79,8 @@ export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSerie
             orbitNormalY: number;
             equinoxLonJ2000Deg: number;
             equinoxLonRateDegPerYr: number;
+            generalPrecessionLonDeg: number;
+            generalPrecessionLonRateDegPerYr: number;
         };
     };
     alphaArcsecPerYr: number;

package/types/earth/sidereal-channel-artifact.d.cts

@@ -1,7 +1,7 @@
 export const SIDEREAL_CHANNEL_ARTIFACT: Readonly<{
     t0Yr: -10000000;
     stepYr: 2000;
-    windowYr: 2000;
+    windowYr: 10000;
     lamDotRel: number[];
     meta: {
         dumpSha256: string;
@@ -9,4 +9,4 @@ export const SIDEREAL_CHANNEL_ARTIFACT: Readonly<{
         chaprontGate: number;
     };
 }>;
-export const SIDEREAL_CHANNEL_ARTIFACT_HASH: "528101ca3410dc9b";
+export const SIDEREAL_CHANNEL_ARTIFACT_HASH: "a6f4729c2de722af";

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

@@ -1,29 +1,23 @@
 /**
  * @param {{
- *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number},
+ *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number, generalPrecessionLonDeg: number},
  *   massLossSiderealSecondsAtYearFn: (year: number) => number,
- *   apsidalSecularDegPerYrFn?: (year: number) => number,
  * }} opts
- *   sampleAt: the one-source movement sampler, calendar-year keyed
- *   (createDeepOrbitalHistory build().at wrapped by the caller — series
- *   tier where the artifact is available, mode tail beyond).
+ *   sampleAt: the one-source movement sampler, keyed in the sampler's own
+ *   year (Julian years from J2000 TT, 2000 = J2000; createDeepOrbitalHistory
+ *   build().at wrapped by the caller — series tier where the artifact is
+ *   available, mode tail beyond).
  *   massLossSiderealSecondsAtYearFn: the caller's secular mass-loss
  *   sidereal-year law, SI seconds, IAU-anchored at J2000.
- *   apsidalSecularDegPerYrFn: the chain's SECULAR apsidal tangent of date
- *   (keplerian-chain computeApsidalSecularDegPerYr for Earth — deg/yr,
- *   J2000-ecliptic frame). When given, the anomalistic year is the
- *   mean-element construction on this smooth rate; when absent, it falls
- *   back to the cardinal structure's measured year-over-year interval
- *   (which carries the banked series' century-scale ϖ̇ wobble).
  */
-export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, apsidalSecularDegPerYrFn }: {
+export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }: {
     sampleAt: (year: number) => {
         e: number;
         periOfDateDeg: number;
         equinoxLonJ2000Deg: number;
+        generalPrecessionLonDeg: number;
     };
     massLossSiderealSecondsAtYearFn: (year: number) => number;
-    apsidalSecularDegPerYrFn?: (year: number) => number;
 }): Readonly<{
     /** Mean tropical year of date, SI seconds. @param {number} year */
     tropicalYearSecondsAtYear: (year: number) => number;
@@ -31,12 +25,16 @@ export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, a
     siderealYearSecondsAtYear: (year: number) => number;
     /** Anomalistic year of date, SI seconds. @param {number} year */
     anomalisticYearSecondsAtYear: (year: number) => number;
-    /** Axial precession period of date (equinox regression beat), years. @param {number} y */
+    /** Axial precession period of date (equinox regression beat), JULIAN years. @param {number} y */
     axialPrecessionYearsAtYear: (y: number) => number;
-    /** Perihelion (apsidal-vs-equinox) precession period of date, years. @param {number} y */
+    /** Perihelion (apsidal-vs-equinox) precession period of date, JULIAN years. @param {number} y */
     perihelionPrecessionYearsAtYear: (y: number) => number;
-    /** Apsidal precession period of date vs the fixed stars, years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
+    /** Apsidal precession period of date vs the fixed stars, JULIAN years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
     inclinationPrecessionYearsAtYear: (y: number) => number;
+    /** The movement's own general precession of date, ″ per Julian year (retrograde positive) — 360°/this ≡ axialPrecessionYearsAtYear to second order. @param {number} y */
+    generalPrecessionArcsecPerJulianYr: (y: number) => number;
+    /** The movement's own apsidal rate of date against the FIXED frame, ″ per Julian year (prograde positive) — the rate the Perihelion Longitudes Prec. cell shows for Earth; 360°/this ≡ inclinationPrecessionYearsAtYear exactly. @param {number} y */
+    apsidalRateFixedArcsecPerJulianYr: (y: number) => number;
     /** The per-cardinal layer (D6 rules: year lengths λ̇-corrected; offsets/spreads raw — μs-class/second-order there). */
     cardinal: Readonly<{
         /** @param {number} y @param {'VE'|'SS'|'AE'|'WS'} type */

package/types/model.d.ts

@@ -61,11 +61,11 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         nodalPeriodYears: number;
         /** The obliquity beat 2π/(ψ̇(t) − |s₃|), years — the SHIPPED deep-time obliquity period (falsification leg 1; 41,224 at J2000, on the composed rate). @param {number} year @returns {number} */
         obliquityBeatYearsAtYear: (year: number) => number;
-        /** The apsidal (perihelion vs the stars) period from the engine-D chain's secular tangent, years — inside the published window only (the tangent is an extrapolation beyond the banked series: it turns negative at −5 Myr); null beyond. @param {number} year @returns {number|null} */
+        /** The apsidal (perihelion vs the stars) period of date, JULIAN years — the ONE movement's own apsidal rate (the scene's perihelion; anom/(anom − sid) of the one family, ≡ 360°/apsidalRateFixedArcsecPerJulianYr), inside the published window only; null beyond. Before 2026-09 this rode the planet chain's secular tangent, which parted from the scene's perihelion passages by seconds of anomalistic year within a few millennia (year-lengths.cjs). @param {number} year @returns {number|null} */
         apsidalPeriodYearsAtYear: (year: number) => number | null;
         /** T_aps(t) / T_p(t) — the apsidal period in of-date precession periods (4.33 at J2000, a reading; 0.84 … 9.9 across ±26 kyr, measured); null beyond the published window. @param {number} year @returns {number|null} */
         apsidalPerPrecessionAtYear: (year: number) => number | null;
-        /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), years, on the of-date T_p; null beyond the published window. @param {number} year @returns {number|null} */
+        /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), JULIAN years, on the of-date T_p — ≡ the one family's anom/(anom − trop); null beyond the published window. @param {number} year @returns {number|null} */
         periOfDatePeriodYearsAtYear: (year: number) => number | null;
         /** T_peri(t) / T_p(t) (0.812 at J2000 — the J2000 reading); null beyond the published window. @param {number} year @returns {number|null} */
         periOfDatePerPrecessionAtYear: (year: number) => number | null;
@@ -116,6 +116,8 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         axialPrecessionYearsAtYear: (y: number) => number;
         perihelionPrecessionYearsAtYear: (y: number) => number;
         inclinationPrecessionYearsAtYear: (y: number) => number;
+        generalPrecessionArcsecPerJulianYr: (y: number) => number;
+        apsidalRateFixedArcsecPerJulianYr: (y: number) => number;
         cardinal: Readonly<{
             yearLengthSeconds: (y: number, type: "VE" | "SS" | "AE" | "WS") => number;
             eocOffsetSeconds: (y: number, type: "VE" | "SS" | "AE" | "WS") => number;