@essrt/physics 23.3.3 → 24.0.0

package/types/model/surface.d.ts

@@ -0,0 +1,281 @@
+/** @param {import('./cardinal-instants.js').Ctx} ctx @returns the published surface (frozen) */
+export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Readonly<{
+    time: Readonly<{
+        yearFromJD: (jd: number) => number;
+        jdFromYear: (year: number) => number;
+        siTropicalYearDays: number;
+    }>;
+    massChain: Readonly<{
+        moonDistanceKm: any;
+        moonSiderealMonthDays: any;
+        siderealYearDays: any;
+        dayLengthSeconds: any;
+        massRatioEarthMoon: any;
+        lunarSmallParameter: number;
+        moonMassFraction: number;
+        solarTidalShiftKm: number;
+        moonDistanceCorrectedKm: number;
+        gmEarthMoonSystemKm3S2: number;
+        gmEarthKm3S2: number;
+        gmMoonKm3S2: number;
+        gmSunPlusEarthKm3S2: number;
+        gmSunKm3S2: number;
+        sunMassKg: number;
+        earthMassKg: number;
+        gravitationalConstantKm3KgS2: any;
+    }>;
+    sun: Readonly<{
+        /** Mean longitude of date, degrees, unwrapped — L0 + 360·∫dt/T_trop (the one-family tropical year). @param {number} year @returns {number} */
+        meanLongitudeDegAt: (year: number) => number;
+        /** The Sun's mean perihelion longitude of date (Earth's ϖ, mean elements), degrees, wrapped to [0, 360). @param {number} year @returns {number} */
+        perihelionLongitudeDegAt: (year: number) => number;
+    }>;
+    epoch: Readonly<{
+        yearToTMa: (year: number) => number;
+        lodSecondsAtYear: (year: number) => number | null;
+        lodRealSecondsAtYear: (year: number) => number;
+        alphaAtYear: (year: number) => number;
+        moonDistanceKmAtYear: (year: number) => number;
+        siderealYearSecondsAtYear: (year: number) => number;
+        deltaTSecondsAtYear: (year: number) => number;
+        axialPrecessionYearsAtYear: (year: number) => number;
+        tropicalYearSecondsAtYear: (year: number) => number;
+    }>;
+    lunisolar: Readonly<{
+        /** The mean lunisolar precession period T_p(t), years — 1,296,000/the precession RATE of date (25,771.6 at J2000): THE published deep-time clock. @param {number} year @returns {number} */
+        meanPeriodYearsAtYear: (year: number) => number;
+        /** THE PRECESSION RATE OF DATE ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ″/yr — the constant's scaling on the spin and the lunar torque times the cited mean obliquity's cosine (plan 12 §8 option (a); p₀ = 50.2880, the IAU 2006 input). @param {number} year @returns {number} */
+        meanRateArcsecPerYrAtYear: (year: number) => number;
+        /** THE PRECESSION CONSTANT OF DATE as a rate on the J2000 obliquity, ″/yr — the spin and torque scaling alone (the hybrid obliquity integration's integrand: its own ŝ·n̂ carries the cosine). @param {number} year @returns {number} */
+        torqueRateArcsecPerYrAtYear: (year: number) => number;
+        /** 1,296,000/the torque rate, years. @param {number} year @returns {number} */
+        torquePeriodYearsAtYear: (year: number) => number;
+        /** cos ε̄(t)/cos ε̄(0) — the rate of date's obliquity factor on the constant (1 at J2000; 1.0075 at the Devonian, 1.0245 at 2.46 Ga). @param {number} year @returns {number} */
+        precessionObliquityFactorAtYear: (year: number) => number;
+        /** THE CITED MEAN OBLIQUITY OF DATE ε̄(t), degrees — Farhat et al. 2022's own column of the Earth–Moon history (a Myr-scale MEAN: 23.25° at age 0 against the J2000 instant 23.44°; 22.23° Devonian, 19.73° at 2.46 Ga); read as the shape cos ε̄(t)/cos ε̄(0) on the precession rate. NaN beyond the input's span (4,429 Ma); forward of J2000 the J2000 slope is held. @param {number} year @returns {number} */
+        meanObliquityCitedDegAtYear: (year: number) => number;
+        /** The of-date period — the one-family beat inside ±2 Myr, the composed mean beyond. @param {number} year @returns {number} */
+        ofDatePeriodYearsAtYear: (year: number) => number;
+        /** f_S, the solar fraction of the J2000 precession torque. */
+        solarShareJ2000: number;
+        /** THE PRECESSION OF DATE — the solar-torque factor 1 + f_S·([(1 − e²)/(1 − e₀²)]^(−3/2) − 1) at the eccentricity of date (1 at J2000; +3.2·10⁻⁴ on the million-year mean). The spin integration multiplies it into the secular rate. @param {number} year @returns {number} */
+        solarTorqueFactorOfDateAtYear: (year: number) => number;
+        /** THE PRECESSION OF DATE — J₂(t)/J₂₀, the oblateness of date from the GIA channel (ψ̇ ∝ J₂; 1 at J2000, bounded ±5·10⁻⁴; dJ₂/dt at J2000 = the observed satellite rate). @param {number} year @returns {number} */
+        ellipticityRatioOfDateAtYear: (year: number) => number;
+        /** (a₀/a_M(t))³ — the lunar torque's growth on the recession history. @param {number} year @returns {number} */
+        lunarTorqueFactorAtYear: (year: number) => number;
+        /** f_S + (1 − f_S)(a₀/a_M)³ — the torque term of the composed precession rate. @param {number} year @returns {number} */
+        torqueTermAtYear: (year: number) => number;
+        /** The hybrid's precession constant α = p₀ / cos ε₀, ″/yr (p₀ the derived J2000 rate, ε₀ the J2000 obliquity input; 54.81). */
+        torqueConstantJ2000ArcsecPerYr: number;
+        /** |s₃|, the dominant nodal mode of Earth's orbit — the largest-amplitude ζ mode of the banked deep secular modes, ″/yr (18.85; the obliquity beat's partner). ONE home for the registry, the browser and the API. */
+        nodalModeS3ArcsecPerYr: number;
+        /** The ecliptic (nodal) precession period 1,296,000/|s₃|, years (68,751) — an orbital quantity, fixed at every epoch at the two-body level. */
+        nodalPeriodYears: number;
+        /** THE SPIN AXIS'S INERTIAL PRECESSION RATE OF DATE, ″/yr — the composed rate scaled by the
+         *  hybrid's self-anchor factor (the inertial-to-general ratio at J2000, the multiplicative form
+         *  the obliquity integration itself applies): 50.385 at J2000 (series-driven). The obliquity beat and the
+         *  obliquity lines ride THIS rate; the general precession (meanRateArcsecPerYrAtYear) is the
+         *  equinox along the ecliptic of date. @param {number} year @returns {number} */
+        spinAxisInertialRateArcsecPerYrAtYear: (year: number) => number;
+        /** The obliquity beat 2π/(ψ̇ᵢₙ(t) − |s₃|), years — the SHIPPED deep-time obliquity period at an INSTANT
+         *  (falsification leg 1; 41,098 at J2000). ψ̇ᵢₙ is the spin axis's INERTIAL rate (above): the equator
+         *  against the ecliptic's nodal mode in the inertial frame, the rate the hybrid obliquity integration
+         *  self-anchors on — not the general precession, which reads the beat 0.3 % long (41,224). v23.0: a
+         *  dated record sees the beat of the pole rate's mean over its own window, with the pole's of-date
+         *  factors (obliquityPeriodRealizedYears below: 40,927 over the LR04 record, the climate formula's p+s₀
+         *  line; La2004's own series reads 40,900 on the same window — compare realized with realized, never
+         *  this instant with Laskar's long-window precession frequency 50.476). @param {number} year @returns {number} */
+        obliquityBeatYearsAtYear: (year: number) => number;
+        /** THE REALIZED POLE RATE OVER A WINDOW, ″/yr — see the one home above (the mean of the inertial rate of
+         *  date with its three of-date factors); the rate the obliquity lines of a record spanning the window carry. */
+        spinAxisInertialRateMeanArcsecPerYr: (yearFrom: number, yearTo: number, stepYears?: number, opts?: {
+            ellipticity?: boolean;
+        }) => number;
+        /** THE REALIZED OBLIQUITY PERIOD OVER A WINDOW, years — 1,296,000/(⟨ψ̇ᵢₙ⟩_window − |s₃|): the period a record
+         *  spanning [yearFrom, yearTo] sees (40,930 over the LR04 record 0–5.32 Ma, 41,018 over the last Myr); the
+         *  J2000 beat (obliquityBeatYearsAtYear) is its instantaneous limit.
+         *  @param {number} yearFrom @param {number} yearTo @returns {number} */
+        obliquityPeriodRealizedYears: (yearFrom: number, yearTo: number) => number;
+        /** 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), 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;
+        /** The published of-date window, years from 2000. */
+        publishedWindowYears: 2000000;
+        /** THE SPIN AXIS'S INERTIAL PRECESSION PERIOD at J2000, years (25,722 in the shipped,
+         *  series-driven configuration) — the rate the hybrid obliquity integration self-anchors
+         *  ŝ on (earth/deep-orbital-history, D4e): the general-precession anchor p₀ PLUS the
+         *  planetary mean the N-body ecliptic's own motion carries (+0.097 ″/yr ≡ IAU's planetary
+         *  precession ψ_A − p_A; the bare 18-mode table anchors 0.18 — the g₂-gap trap: pass the
+         *  secular series). Two rates, two uses — the obliquity line frequencies p + s_i (the
+         *  spin axis against the ecliptic's nodal modes, in the inertial frame) carry THIS
+         *  rate; the climatic-precession lines p + g_i (the perihelion from the equinox along
+         *  the ecliptic of date) carry p₀. */
+        spinAxisInertialPeriodYearsJ2000: number;
+        /** The same as a rate, ″/yr (50.385 at J2000 with the series; p₀ = 50.288 is the general precession). */
+        spinAxisInertialRateArcsecPerYrJ2000: number;
+    }>;
+    earth: Readonly<{
+        perihelionLongitudeDeg: (year: number) => number;
+        obliquityDeg: (year: number) => number;
+        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (plan 12 §8.5: the published surfaces serve `obliquityDeg` inside it only; null without a series, i.e. a bare createModel()). */
+        obliquityHybridDomainYears: [number, number] | null;
+        /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column; plan 12 §8 option (a)), degrees — a Myr-scale mean; NaN beyond 4,429 Ma. @param {number} year @returns {number} */
+        meanObliquityCitedDeg: (year: number) => number;
+        eccentricity: (year: number) => number;
+        inclinationDeg: (year: number) => number;
+        ascendingNodeDeg: (year: number) => number;
+    }>;
+    lengths: Readonly<{
+        tropicalYearDays: (year: number) => number;
+        siderealYearDays: (year: number) => number;
+        anomalisticYearDays: (year: number) => number;
+        dayLengthSeconds: (year: number) => number;
+        lodRealSeconds: (year: number) => number;
+        siderealDaySeconds: (year: number) => number;
+        stellarDaySeconds: (year: number) => number;
+    }>;
+    cardinal: Readonly<{
+        jd: (year: number, type: string) => number;
+        jdNearUT: (jdUT: number, type: string) => number;
+        raDeg: (year: number, type: string) => number;
+        yearLengthDays: (year: number, type: string) => number;
+    }>;
+    cardinalStructure: Readonly<{
+        yearLengthSeconds: (year: number, type: "VE" | "SS" | "AE" | "WS") => number;
+        eocOffsetSeconds: (year: number, type: "VE" | "SS" | "AE" | "WS") => number;
+        spreadSeconds: (year: number) => {
+            VE: number;
+            SS: number;
+            AE: number;
+            WS: number;
+            meanSeconds: number;
+        };
+        anomalisticYearSeconds: (year: number) => number;
+    }>;
+    yearLengths: Readonly<{
+        tropicalYearSecondsAtYear: (year: number) => number;
+        siderealYearSecondsAtYear: (year: number) => number;
+        anomalisticYearSecondsAtYear: (year: number) => number;
+        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;
+            spreadSeconds: (y: number) => {
+                VE: number;
+                SS: number;
+                AE: number;
+                WS: number;
+                meanSeconds: number;
+            };
+        }>;
+        planetaryRelAtYear: (year: number) => number;
+    }>;
+    moon: Readonly<{
+        distanceKmAtYear: (year: number) => number;
+        siderealMonthDaysAtYear: (year: number) => number;
+        synodicMonthDays: number;
+        lonDegAtJD: (jd: number) => number;
+        betaDegAtJD: (jd: number) => number;
+        distanceKmAtJD: (jd: number) => number;
+    }>;
+    eclipse: Readonly<{
+        sunLonDegAtJD: (jd: number) => number;
+        sunLonCompletedDegAtJD: (jd: number) => number;
+        /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings (plan 06 R1/I1). @param {number} jdUT @returns {number} */
+        sunApparentLonDegAtJD: (jdUT: number) => number;
+        findLunarInRange: (jdStart: number, jdEnd: number) => {
+            jd: number;
+            beta: number;
+            moonDistance_km: number;
+            type: string;
+            magnitudeUmbral: number;
+            magnitudePenumbral: number;
+        }[];
+        findSolarInRange: (jdStart: number, jdEnd: number) => {
+            jd: number;
+            beta: number;
+            moonDistance_km: number;
+            type: string;
+            moonAppR_topo: number;
+            sunAppR: number;
+            moonSunRatio: number;
+        }[];
+        deltaTSecondsAtJD: (jd: number) => number;
+        /** The ΔT trend anchor at J2000 (seconds): the TT − UT bridge every finder axis adds to a
+         *  true-UT JD. The joint ΔT fit's output, shipped in its own artifact (FIT.DT_ANCHOR —
+         *  plan 14: never a constants-file input), exposed here so no consumer reaches into the
+         *  constants for it. */
+        ttBridgeSeconds: any;
+        frameworkSunDeps: Readonly<{
+            sunMeanLongitudeJ2000Deg: any;
+            tropicalRateDegPerCy: number;
+            eccentricityAt: (year: number) => any;
+            perihelionLongitudeDegAt: (year: number) => number;
+            meanLongitudeDegAt: (year: number) => number;
+        }>;
+        umbraGroundAtJD: (jd: number) => {
+            latDeg: number;
+            lonDeg: number;
+        } | null;
+        /** THE ONE GMST — mean sidereal time at Greenwich, degrees, at a UT JD (eclipse/besselian.cjs computeGmstDeg; the observer tier's hour angle). @param {number} jdUT @returns {number} */
+        gmstDegAtJD: (jdUT: number) => number;
+        solarLocalCircumstances: (jdGreatest: number, latDeg: number, lonDeg: number) => {
+            kind: "none" | "partial" | "annular" | "total";
+            magnitude: number;
+            maxJd: number;
+            contacts: {
+                c1: number | null;
+                c2: number | null;
+                c3: number | null;
+                c4: number | null;
+            };
+            centralDurationSeconds: number | null;
+        };
+    }>;
+    climate: Readonly<{
+        l1OrbitalPermil: (year: number) => number;
+    }>;
+    planets: Readonly<{
+        keys: readonly string[];
+        record: (k: string) => Record<string, any> | undefined;
+        perihelionLongitudeDeg: (k: string, year: number) => number;
+        ascendingNodeInvPlaneDeg: (k: string, year: number) => number;
+        invPlaneInclinationDeg: (k: string, year: number) => number;
+        eccentricity: (k: string, year: number) => number;
+        spin: (k: string) => Readonly<{
+            key: string;
+            cassiniLocked: boolean;
+            momentOfInertiaFactor: number;
+            momentOfInertiaFactorClass: string;
+            alphaArcsecPerYr: number;
+            satelliteQuadrupole: number;
+            satelliteAngularMomentum: number;
+            spinRetrograde: boolean;
+            meanMotionRadPerYr: number;
+            obliquityJ2000Deg: number;
+            spinPrecessionRateArcsecPerYrJ2000: number | null;
+            axialPrecessionPeriodYearsJ2000: number | null;
+            obliquityDegAtYear: (year: number) => number | null;
+            spinPrecessionRateArcsecPerYrAtYear: (year: number) => number | null;
+            obliquityEnvelopeDeg: (year: number, spanYr: number) => {
+                minDeg: number;
+                meanDeg: number;
+                maxDeg: number;
+            } | null;
+            spinAxisEclipticJ2000AtYear: (year: number) => number[] | null;
+            orbitNormalEclipticJ2000AtYear: (year: number) => number[];
+            stepYr: number;
+            maxSpanYr: 10000000;
+        }>;
+    }>;
+}>;

package/types/model/time-axis.d.ts

@@ -0,0 +1,12 @@
+/**
+ * The time axis — exact JD ↔ model-year conversion on the SI axis.
+ *
+ * One slice of `assembleModel` (../model.js): receives the context the earlier slices built, returns the bindings the later ones read.
+ */
+/** @param {import('./inv-plane-elements.js').Ctx} ctx @returns the bindings the later slices read */
+export function assembleTimeAxis(ctx: import("./inv-plane-elements.js").Ctx): {
+    jdFromYear: (year: number) => number;
+    siTropicalYearDays: number;
+    yearFromJD: (jd: number) => number;
+};
+export type Ctx = import("./inv-plane-elements.js").Ctx & ReturnType<typeof assembleTimeAxis>;

package/types/model/year-day-lengths.d.ts

@@ -0,0 +1,8 @@
+/** @param {import('./earth-scalars.js').Ctx} ctx @returns the bindings the later slices read */
+export function assembleYearDayLengths(ctx: import("./earth-scalars.js").Ctx): {
+    anomalisticYearDays: (year: number) => number;
+    certifiedAxialPrecessionJ2000Years: () => number;
+    dayLengthSeconds: (year: number) => number;
+    siderealYearDays: (year: number) => number;
+};
+export type Ctx = import("./earth-scalars.js").Ctx & ReturnType<typeof assembleYearDayLengths>;

package/types/model.d.ts

@@ -1,9 +1,3 @@
-/**
- * RA-day-offset Fourier amplitudes (ms). KNOWN EXCEPTION carried over from
- * the reference adapters: these are fit results that live as literals in the
- * engine too — packaging them into FITTED_COEFFICIENTS is the remaining
- * §7a-step-1b move; until then this is their single packaged home.
- */
 /**
  * Assemble the model surfaces from a resolved constants context + fitted
  * coefficients. Internal: `createModel` in index.js composes this with the
@@ -11,28 +5,15 @@
  *
  * @param {Readonly<Record<string, any>>} C  the frozen constants context
  * @param {Readonly<Record<string, any>>} F  the fitted coefficients
- * @param {{ eccentricityAt?: (year: number) => number, eccentricityRateAt?: (year: number) => number, perihelionLongitudeDegAt?: (year: number) => number }} [laws]
- *   RESEARCH OVERRIDES for Earth's orbit laws (doc 109 §7): an alternative e(t),
- *   de/dt(t) and ϖ_of-date(t) flow through every consumer (the eclipse Sun, the
- *   Moon's E-factor, the cardinal points) exactly as the shipped laws do. Default
- *   {} = the shipped laws, bit-identical to before this parameter existed. Not a
- *   counterfactual in the §2d sense (no hash change) — callers must say when
- *   they used it; the generators refuse to --write under an override.
- * @param {{ giaEllipticityOfDate?: boolean }} [options]
- *   `giaEllipticityOfDate: false` builds the one-source history WITHOUT the J₂(t)/J₂₀ factor of date
- *   (the ice-driven GIA channel, which rides the climate formula's own post-MPT layer): the orbital
- *   FORCING alone. The T1 exporter derives the climate lines' obliquity rate on it, so the lines
- *   cannot depend on their own fit (v23.0; measured 3e-5 of the rate, 2–3 yr of period). Default true =
- *   the shipped physics (v15.0), bit-identical to before this parameter existed.
- * @returns the assembled surfaces (epoch, earth, lengths, cardinal, moon) — type inferred so ReturnType stays precise
+ * @param {import('./model/context.js').Laws} [laws]  research overrides for Earth's orbit laws
+ *   (model/context.js); default {} = the shipped laws, bit-identical to before this parameter existed
+ * @param {any} [secularSeriesArtifact]  the governed N-body secular series (the API/Node path; the
+ *   npm package alone has no 8 MB series)
+ * @param {import('./model/context.js').Options} [options]  `giaEllipticityOfDate: false` = the orbital
+ *   forcing alone (model/context.js); default true = the shipped physics, bit-identical to before
+ * @returns the assembled surfaces (epoch, earth, lengths, cardinal, moon, …) — type inferred so ReturnType stays precise
  */
-export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Record<string, any>>, laws?: {
-    eccentricityAt?: (year: number) => number;
-    eccentricityRateAt?: (year: number) => number;
-    perihelionLongitudeDegAt?: (year: number) => number;
-}, secularSeriesArtifact?: any, options?: {
-    giaEllipticityOfDate?: boolean;
-}): Readonly<{
+export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Record<string, any>>, laws?: import("./model/context.js").Laws, secularSeriesArtifact?: any, options?: import("./model/context.js").Options): Readonly<{
     time: Readonly<{
         yearFromJD: (jd: number) => number;
         jdFromYear: (year: number) => number;
@@ -58,9 +39,7 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         gravitationalConstantKm3KgS2: any;
     }>;
     sun: Readonly<{
-        /** Mean longitude of date, degrees, unwrapped — L0 + 360·∫dt/T_trop (the one-family tropical year). @param {number} year @returns {number} */
         meanLongitudeDegAt: (year: number) => number;
-        /** The Sun's mean perihelion longitude of date (Earth's ϖ, mean elements), degrees, wrapped to [0, 360). @param {number} year @returns {number} */
         perihelionLongitudeDegAt: (year: number) => number;
     }>;
     epoch: Readonly<{
@@ -75,90 +54,39 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         tropicalYearSecondsAtYear: (year: number) => number;
     }>;
     lunisolar: Readonly<{
-        /** The mean lunisolar precession period T_p(t), years — 1,296,000/the precession RATE of date (25,771.6 at J2000): THE published deep-time clock. @param {number} year @returns {number} */
         meanPeriodYearsAtYear: (year: number) => number;
-        /** THE PRECESSION RATE OF DATE ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ″/yr — the constant's scaling on the spin and the lunar torque times the cited mean obliquity's cosine (plan 12 §8 option (a); p₀ = 50.2880, the IAU 2006 input). @param {number} year @returns {number} */
         meanRateArcsecPerYrAtYear: (year: number) => number;
-        /** THE PRECESSION CONSTANT OF DATE as a rate on the J2000 obliquity, ″/yr — the spin and torque scaling alone (the hybrid obliquity integration's integrand: its own ŝ·n̂ carries the cosine). @param {number} year @returns {number} */
         torqueRateArcsecPerYrAtYear: (year: number) => number;
-        /** 1,296,000/the torque rate, years. @param {number} year @returns {number} */
         torquePeriodYearsAtYear: (year: number) => number;
-        /** cos ε̄(t)/cos ε̄(0) — the rate of date's obliquity factor on the constant (1 at J2000; 1.0075 at the Devonian, 1.0245 at 2.46 Ga). @param {number} year @returns {number} */
         precessionObliquityFactorAtYear: (year: number) => number;
-        /** THE CITED MEAN OBLIQUITY OF DATE ε̄(t), degrees — Farhat et al. 2022's own column of the Earth–Moon history (a Myr-scale MEAN: 23.25° at age 0 against the J2000 instant 23.44°; 22.23° Devonian, 19.73° at 2.46 Ga); read as the shape cos ε̄(t)/cos ε̄(0) on the precession rate. NaN beyond the input's span (4,429 Ma); forward of J2000 the J2000 slope is held. @param {number} year @returns {number} */
         meanObliquityCitedDegAtYear: (year: number) => number;
-        /** The of-date period — the one-family beat inside ±2 Myr, the composed mean beyond. @param {number} year @returns {number} */
         ofDatePeriodYearsAtYear: (year: number) => number;
-        /** f_S, the solar fraction of the J2000 precession torque. */
         solarShareJ2000: number;
-        /** THE PRECESSION OF DATE — the solar-torque factor 1 + f_S·([(1 − e²)/(1 − e₀²)]^(−3/2) − 1) at the eccentricity of date (1 at J2000; +3.2·10⁻⁴ on the million-year mean). The spin integration multiplies it into the secular rate. @param {number} year @returns {number} */
         solarTorqueFactorOfDateAtYear: (year: number) => number;
-        /** THE PRECESSION OF DATE — J₂(t)/J₂₀, the oblateness of date from the GIA channel (ψ̇ ∝ J₂; 1 at J2000, bounded ±5·10⁻⁴; dJ₂/dt at J2000 = the observed satellite rate). @param {number} year @returns {number} */
         ellipticityRatioOfDateAtYear: (year: number) => number;
-        /** (a₀/a_M(t))³ — the lunar torque's growth on the recession history. @param {number} year @returns {number} */
         lunarTorqueFactorAtYear: (year: number) => number;
-        /** f_S + (1 − f_S)(a₀/a_M)³ — the torque term the unit divides H_era by. @param {number} year @returns {number} */
         torqueTermAtYear: (year: number) => number;
-        /** The hybrid's precession constant α = p₀ / cos ε₀, ″/yr (p₀ the derived J2000 rate, ε₀ the J2000 obliquity input; 54.81). */
         torqueConstantJ2000ArcsecPerYr: number;
-        /** |s₃|, the dominant nodal mode of Earth's orbit — the largest-amplitude ζ mode of the banked deep secular modes, ″/yr (18.85; the obliquity beat's partner). ONE home for the registry, the browser and the API. */
         nodalModeS3ArcsecPerYr: number;
-        /** The ecliptic (nodal) precession period 1,296,000/|s₃|, years (68,751) — an orbital quantity, fixed at every epoch at the two-body level. */
         nodalPeriodYears: number;
-        /** THE SPIN AXIS'S INERTIAL PRECESSION RATE OF DATE, ″/yr — the composed rate scaled by the
-         *  hybrid's self-anchor factor (the inertial-to-general ratio at J2000, the multiplicative form
-         *  the obliquity integration itself applies): 50.385 at J2000 (series-driven). The obliquity beat and the
-         *  obliquity lines ride THIS rate; the general precession (meanRateArcsecPerYrAtYear) is the
-         *  equinox along the ecliptic of date. @param {number} year @returns {number} */
         spinAxisInertialRateArcsecPerYrAtYear: (year: number) => number;
-        /** The obliquity beat 2π/(ψ̇ᵢₙ(t) − |s₃|), years — the SHIPPED deep-time obliquity period at an INSTANT
-         *  (falsification leg 1; 41,098 at J2000). ψ̇ᵢₙ is the spin axis's INERTIAL rate (above): the equator
-         *  against the ecliptic's nodal mode in the inertial frame, the rate the hybrid obliquity integration
-         *  self-anchors on — not the general precession, which reads the beat 0.3 % long (41,224). v23.0: a
-         *  dated record sees the beat of the pole rate's mean over its own window, with the pole's of-date
-         *  factors (obliquityPeriodRealizedYears below: 40,927 over the LR04 record, the climate formula's p+s₀
-         *  line; La2004's own series reads 40,900 on the same window — compare realized with realized, never
-         *  this instant with Laskar's long-window precession frequency 50.476). @param {number} year @returns {number} */
         obliquityBeatYearsAtYear: (year: number) => number;
-        /** THE REALIZED POLE RATE OVER A WINDOW, ″/yr — see the one home above (the mean of the inertial rate of
-         *  date with its three of-date factors); the rate the obliquity lines of a record spanning the window carry. */
         spinAxisInertialRateMeanArcsecPerYr: (yearFrom: number, yearTo: number, stepYears?: number, opts?: {
             ellipticity?: boolean;
         }) => number;
-        /** THE REALIZED OBLIQUITY PERIOD OVER A WINDOW, years — 1,296,000/(⟨ψ̇ᵢₙ⟩_window − |s₃|): the period a record
-         *  spanning [yearFrom, yearTo] sees (40,930 over the LR04 record 0–5.32 Ma, 41,018 over the last Myr); the
-         *  J2000 beat (obliquityBeatYearsAtYear) is its instantaneous limit.
-         *  @param {number} yearFrom @param {number} yearTo @returns {number} */
         obliquityPeriodRealizedYears: (yearFrom: number, yearTo: number) => number;
-        /** 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), 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;
-        /** The published of-date window, years from 2000. */
         publishedWindowYears: 2000000;
-        /** THE SPIN AXIS'S INERTIAL PRECESSION PERIOD at J2000, years (25,722 in the shipped,
-         *  series-driven configuration) — the rate the hybrid obliquity integration self-anchors
-         *  ŝ on (earth/deep-orbital-history, D4e): the general-precession anchor p₀ PLUS the
-         *  planetary mean the N-body ecliptic's own motion carries (+0.097 ″/yr ≡ IAU's planetary
-         *  precession ψ_A − p_A; the bare 18-mode table anchors 0.18 — the g₂-gap trap: pass the
-         *  secular series). Two rates, two uses — the obliquity line frequencies p + s_i (the
-         *  spin axis against the ecliptic's nodal modes, in the inertial frame) carry THIS
-         *  rate; the climatic-precession lines p + g_i (the perihelion from the equinox along
-         *  the ecliptic of date) carry p₀. */
         spinAxisInertialPeriodYearsJ2000: number;
-        /** The same as a rate, ″/yr (50.385 at J2000 with the series; p₀ = 50.288 is the general precession). */
         spinAxisInertialRateArcsecPerYrJ2000: number;
     }>;
     earth: Readonly<{
         perihelionLongitudeDeg: (year: number) => number;
         obliquityDeg: (year: number) => number;
-        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (plan 12 §8.5: the published surfaces serve `obliquityDeg` inside it only; null without a series, i.e. a bare createModel()). */
         obliquityHybridDomainYears: [number, number] | null;
-        /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column; plan 12 §8 option (a)), degrees — a Myr-scale mean; NaN beyond 4,429 Ma. @param {number} year @returns {number} */
         meanObliquityCitedDeg: (year: number) => number;
         eccentricity: (year: number) => number;
         inclinationDeg: (year: number) => number;
@@ -224,7 +152,6 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
     eclipse: Readonly<{
         sunLonDegAtJD: (jd: number) => number;
         sunLonCompletedDegAtJD: (jd: number) => number;
-        /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings (plan 06 R1/I1). @param {number} jdUT @returns {number} */
         sunApparentLonDegAtJD: (jdUT: number) => number;
         findLunarInRange: (jdStart: number, jdEnd: number) => {
             jd: number;
@@ -244,6 +171,7 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
             moonSunRatio: number;
         }[];
         deltaTSecondsAtJD: (jd: number) => number;
+        ttBridgeSeconds: any;
         frameworkSunDeps: Readonly<{
             sunMeanLongitudeJ2000Deg: any;
             tropicalRateDegPerCy: number;
@@ -255,6 +183,7 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
             latDeg: number;
             lonDeg: number;
         } | null;
+        gmstDegAtJD: (jdUT: number) => number;
         solarLocalCircumstances: (jdGreatest: number, latDeg: number, lonDeg: number) => {
             kind: "none" | "partial" | "annular" | "total";
             magnitude: number;

package/types/moon/arguments.d.cts

@@ -38,17 +38,13 @@ export type MoonArgsDeg = {
  *   pair powers the (d′) of-date rate completion: pDyn = the DYNAMICAL
  *   axial precession (day-form beat of the engine's sidereal/solar year
  *   evaluators — the tweakpane identity, real at J2000, epoch-valid);
- *   pKin = the KINEMATIC pair's beat (the H/13-family mean the chains
- *   embed). Absent ⇒ the completion is disabled (pre-(d′) behaviour).
- *   Plan 09: the Sun of the arguments is the engine's ONE certified Sun —
+ *   pKin = the KINEMATIC pair's beat (the mean the chains embed). Absent ⇒
+ *   the completion is disabled (pre-(d′) behaviour).
+ *   The Sun of the arguments is the engine's ONE certified Sun —
  *   sunMeanLonDegAtJdTT (the mean longitude that integrates the one-family
  *   tropical year, unwrapped) and sunPerihelionLonDegAtJdTT (the mean
  *   perihelion of date, wrapped — the arguments wrap M, and the rate anchors
- *   difference it over ±25 days wrapped to ±180°). The former "model timeline" Sun — 360° per mean solar
- *   year plus the closed-form integrals of the retired year-length combs,
- *   the perihelion on the H/16 chain count — sat 0.2° behind the IAU mean
- *   Sun at −2000 (the combs could not carry the tropical year's secular
- *   drift); measured against Meeus L0 before the move.
+ *   difference it over ±25 days wrapped to ±180°).
  */
 export function createMoonArguments({ constants, fns }: {
     constants: {
@@ -128,8 +124,6 @@ import { jdToDecimalYear } from "../time/calendar.cjs";
 /** The secular-completion taper half-width, Julian years: the cos² envelope
  *  inside which the chain's Meeus-measured rate anchors and the Lp secular
  *  completion act (full at the anchor, zero beyond). A recorded window of the
- *  completion, not a law — the anchors were measured with it (plan 09 step
- *  5c: it was written as the era unit's twelfth, H/12, "lattice-derived";
- *  the value is kept bit-identical and named for what it is). */
+ *  completion, not a law — the anchors were measured with it. */
 export const SECULAR_COMPLETION_TAPER_YEARS: 27943.083333333332;
 export { jdToDecimalYear };

package/types/planets/keplerian-chain.d.cts

@@ -77,7 +77,8 @@ export type KcPlanetChain = {
     invariablePlane?: (KcInvariablePlane | undefined);
 };
 export const ANCHOR_EPOCH_YEAR: 2000;
-export const ANCHOR_EPOCH_JD: 2451545;
+import ANCHOR_EPOCH_JD_1 = require("../constants/units.cjs");
+import ANCHOR_EPOCH_JD = ANCHOR_EPOCH_JD_1.J2000_JD;
 /** Full osculating elements from a heliocentric ecliptic-J2000 state vector.
  *  @param {number[]} rKm  heliocentric position, km
  *  @param {number[]} vKmS heliocentric velocity, km/s
@@ -154,7 +155,7 @@ export function computePlanetElementsAtYear(year: number, planetChain: KcPlanetC
  * @param {Object<string,KcPlanetChain>=} allChains the full chain map
  * @returns {number} deg/yr, prograde positive
  */
-export function computeApsidalSecularDegPerYr(year: number, planetChain: KcPlanetChain, allChains?: {
+export function computeApsidalTangentDegPerYr(year: number, planetChain: KcPlanetChain, allChains?: {
     [x: string]: KcPlanetChain;
 } | undefined): number;
 /** The chain's canonical constructor from a PARSED governed artifact (K2.1 →
@@ -172,3 +173,4 @@ export function buildPlanetChainsFromArtifactData(art: any, opts?: {
 } | undefined): {
     [x: string]: KcPlanetChain;
 };
+export { ANCHOR_EPOCH_JD };

package/types/planets/model.d.cts

@@ -1,8 +1,8 @@
 export type PlanetModelBody = {
     /**
-     * - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
+     * - the body's orbital-period input, SI days — for the planets the period OF DATE (measured in sun-and-eclipses.js: their rates sit 1.4–1.6°/cy above the VSOP sidereal mean motions, the precession plus input rounding), not the sidereal period the field was once documented as; the scaffold period of the no-chain bodies — the chain planets ride the artifact's elements
      */
-    solarYearInput: number;
+    orbitalPeriodInputDays: number;
     invPlaneInclinationJ2000?: number | undefined;
     longitudePerihelion?: number | undefined;
     /**
@@ -43,7 +43,7 @@ export type PlanetModelRecord = {
 };
 /**
  * @typedef {Object} PlanetModelBody
- * @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} orbitalPeriodInputDays - the body's orbital-period input, SI days — for the planets the period OF DATE (measured in sun-and-eclipses.js: their rates sit 1.4–1.6°/cy above the VSOP sidereal mean motions, the precession plus input rounding), not the sidereal period the field was once documented as; 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)