@essrt/physics 23.3.3 → 23.6.0

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;

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
@@ -172,3 +173,4 @@ export function buildPlanetChainsFromArtifactData(art: any, opts?: {
 } | undefined): {
     [x: string]: KcPlanetChain;
 };
+export { ANCHOR_EPOCH_JD };