@essrt/physics 23.0.0 → 23.2.1

package/types/deltat/deep-time-chain.d.ts

@@ -0,0 +1,266 @@
+/**
+ * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
+ * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin
+ * gates it on BOTH env vars, DT_CORRECTIONS and DT_RESONATOR), else by
+ * `resonator`. The four cycles' δLOD terms follow their flags only when
+ * `cycleLodFollowsFlags` is true (the browser's research toggles switch a
+ * cycle's ΔT and δLOD together); the Node twin sums them unconditionally.
+ * Their ΔT terms always follow their flags. Identical with every flag ON.
+ * @typedef {{bond: boolean, hallstatt: boolean, jose5: boolean, jose4: boolean, resonator: boolean, resonatorDeltaT?: boolean, cycleLodFollowsFlags?: boolean}} DeltaTFlags
+ */
+/**
+ * @param {{
+ *   constants: any,
+ *   fitted: any,
+ *   deepModes: any,
+ *   moonDistanceArtifact: any,
+ *   ageAnchorYear?: number,
+ *   flagsAt?: () => DeltaTFlags,
+ *   planetAuRatios?: Record<string, number>,
+ *   stellarDayRaProjection?: number,
+ * }} deps — `constants` the generated DEFAULT_CONSTANTS, `fitted` the
+ *   generated FITTED_COEFFICIENTS, `deepModes` the DEEP_MODES_ARTIFACT,
+ *   `moonDistanceArtifact` the MOON_DISTANCE_ARTIFACT — the cited distance
+ *   AND mean-obliquity columns, plan 12 option A / §8 (a) (injected so a
+ *   counterfactual context can replace any of them); `ageAnchorYear` the
+ *   scene epoch (default constants.timeReference.sceneEpochYear);
+ *   `flagsAt` the ΔT stack's switches (default all ON); `planetAuRatios`
+ *   the planets' J2000 semi-major axes in AU (for the per-planet family);
+ *   `stellarDayRaProjection` cos ε₀-class factor for the projected stellar day.
+ */
+export function createDeepTimeChain(deps: {
+    constants: any;
+    fitted: any;
+    deepModes: any;
+    moonDistanceArtifact: any;
+    ageAnchorYear?: number;
+    flagsAt?: () => DeltaTFlags;
+    planetAuRatios?: Record<string, number>;
+    stellarDayRaProjection?: number;
+}): Readonly<{
+    meanSolarYearDays: number;
+    meanSiderealYearDays: any;
+    meanSiderealYearSeconds: number;
+    meanLengthOfDay: any;
+    meanTropicalYearJ2000Seconds: number;
+    siTropicalYearDays: number;
+    axialPrecessionJ2000Years: number;
+    axialPrecessionJ2000Days: number;
+    G_CONSTANT: any;
+    MASS_RATIO_EARTH_MOON: any;
+    currentAUDistance: any;
+    earthMoiFactorJ2000: any;
+    moonDistanceKm: any;
+    moonSiderealMonthInput: any;
+    moonOrbitalShift: number;
+    moonDistanceCorrected: number;
+    GM_EARTH_MOON_SYSTEM: number;
+    GM_EARTH_ALONE: number;
+    GM_MOON_ALONE: number;
+    GM_SUN_PLUS_EARTH: number;
+    GM_SUN: number;
+    M_SUN: number;
+    moonApsidalJ2000Seconds: number;
+    moonNodalJ2000Seconds: number;
+    moonSiderealMonthJ2000Seconds: number;
+    moonPeriods: {
+        moonApsidalPrecessionDaysEarth: number;
+        moonNodalPrecessionDaysEarth: number;
+        moonApsidalJ2000Seconds: number;
+        moonNodalJ2000Seconds: number;
+        moonSiderealMonthJ2000Seconds: number;
+        moonAnomalisticMonthDays: number;
+        moonNodalMonthDays: number;
+        moonSynodicMonthDays: number;
+        moonTropicalMonthDays: number;
+        moonFullMoonCycleEarthDays: number;
+        moonFullMoonCycleICRFDays: number;
+        moonApsidalMeetsNodalDays: number;
+        moonLunarLevelingCycleDays: number;
+        moonDraconicYearICRFDays: number;
+        moonDraconicYearEarthDays: number;
+    };
+    certifiedAxialPrecessionJ2000Years: () => number;
+    precessionRateJ2000ArcsecPerYr: () => number;
+    EPOCH_PARAMS: import("../layer0/index.js").EpochParams;
+    ageAnchorYear: any;
+    alphaGia: {
+        alphaAt: (year: number) => number;
+        j2RatioAt: (year: number) => number;
+        laggedL1PermilAt: (year: number) => number;
+        kPerPermille: number;
+        relaxationKyr: number;
+        laggedL1Terms: Array<{
+            period_kyr: number;
+            a: number;
+            b: number;
+        }>;
+    };
+    earthMoiFactorAtAge: (tMa: number) => number;
+    withLatticeAlpha: (build: () => void) => void;
+    j2RatioAtYear: (year: number) => number;
+    iEarthAtAge: (tMa: number) => number;
+    recession: {
+        distanceMetresAtAge: (tMa: number) => number | null;
+        slopeMetresPerMaAtAge: (tMa: number) => number | null;
+        ratioAtAge: (tMa: number) => number | null;
+        spanMa: number;
+    };
+    meanObliquity: {
+        meanObliquityDegAtAge: (tMa: number) => number | null;
+        cosRatioAtAge: (tMa: number) => number | null;
+        meanObliquityJ2000Deg: number;
+        spanMa: number;
+    };
+    dtCycles: {
+        taperAt: (year: number) => number;
+        taperDerivativeAt: (year: number) => number;
+        cycleDeltaTSecondsAt: (cycleKey: string, year: number) => number;
+        cycleLodSecondsAt: (cycleKey: string, year: number) => number;
+        swingDeltaTSecondsAt: (year: number) => number;
+        swingLodSecondsAt: (year: number) => number;
+        swingLodRateAt: (year: number) => number;
+    };
+    s3ArcsecPerYr: number;
+    deepLod: {
+        lodSecondsAtAge: (t_Ma: number) => number | null;
+        lodSecondsAtAgeWithAlpha: (t_Ma: number, alpha: number) => number | null;
+        lodHoursAtAge: (t_Ma: number) => number | null;
+        lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+        lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
+        lunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+        lunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+        precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
+        lunarTorqueFactorAtAge: (tMa: number) => (number | null);
+        precessionTorqueTermAtAge: (tMa: number) => (number | null);
+        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;
+        dLodDtDecompositionAtAge: (t_Ma: number) => {
+            tidal: number | null;
+            gia: number | null;
+            stack: number | null;
+            resonator: number | null;
+            net_L2: number | null;
+            net_L3: number | null;
+            net_L4: number | null;
+        };
+    };
+    l0: import("../layer0/index.js").EpochPrimitives;
+    deepEcc: {
+        eccAt: (tYr: number) => number;
+        eccRateAt: (tYr: number) => number;
+        modulation: (tMa: number, s: number) => number;
+        channelIntegral: (T: number, s: number) => number;
+        eFactorAt: (tYr: number) => number;
+        e0: number;
+        g0: number;
+    };
+    eCompModulation: (tMa: number, s: number) => number;
+    moonChain: {
+        distanceMetresAtAge: (tMa: number) => number | null;
+        distanceKmAtAge: (tMa: number) => number | null;
+        solarDeltaAKmAtAge: (tMa: number, aApparentKm: number) => number | null;
+        distanceCorrectedKmAtAge: (tMa: number) => number | null;
+        siderealMonthSecondsAtAge: (tMa: number) => number | null;
+        synodicMonthSecondsAtAge: (tMa: number) => number | null;
+        tropicalMonthSecondsAtAge: (tMa: number) => number | null;
+        perigeePrecessionSecondsAtAge: (tMa: number) => number | null;
+        nodePrecessionSecondsAtAge: (tMa: number) => number | null;
+        anomalisticMonthSecondsAtAge: (tMa: number) => number | null;
+        nodalMonthSecondsAtAge: (tMa: number) => number | null;
+        apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
+        lunarLevelingSecondsAtAge: (tMa: number) => number | null;
+    };
+    meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
+    meanMoonDistanceAtAge: (tMa: number) => number | null;
+    meanLodSecondsAtAge: (tMa: number) => number | null;
+    meanLodHoursAtAge: (t_Ma: number) => number | null;
+    meanLunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    meanLunisolarPrecessionPeriodYearsAtAge: (tMa: number) => number | null;
+    meanLunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    meanLunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+    precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
+    meanObliquityCitedDegAtAge: (tMa: number) => number | null;
+    meanAuAtAge: (tMa: number) => number;
+    meanSiderealYearSecondsAtAge: (tMa: number) => number;
+    meanTropicalYearSecondsAtAge: (t_Ma: number) => number;
+    meanTropicalYearDaysAtAge: (t_Ma: number) => number | null;
+    meanYearInDaysAtAge: (t_Ma: number) => number | null;
+    meanDeltaTSecondsAtAge: (tMa: number) => number;
+    bondCycleDeltaTCorrection: (year: number) => number;
+    hallstattCycleDeltaTCorrection: (year: number) => number;
+    jose5CycleDeltaTCorrection: (year: number) => number;
+    jose4CycleDeltaTCorrection: (year: number) => number;
+    bondCycleLodCorrection: (year: number) => number;
+    hallstattCycleLodCorrection: (year: number) => number;
+    jose5CycleLodCorrection: (year: number) => number;
+    jose4CycleLodCorrection: (year: number) => number;
+    dtCycleLodCorrectionSum: (year: number) => number;
+    resonatorSwingDeltaTCorrection: (year: number) => number;
+    resonatorSwingLodCorrection: (year: number) => number;
+    resonatorSwingLodRate: (year: number) => number;
+    meanLodSecondsWithCorrectionsAtAge: (t_Ma: number) => number | null;
+    dLodDtDecompositionAtAge: (t_Ma: number) => {
+        tidal: number | null;
+        gia: number | null;
+        stack: number | null;
+        resonator: number | null;
+        net_L2: number | null;
+        net_L3: number | null;
+        net_L4: number | null;
+    };
+    h5Correction: (year: number) => number;
+    computeLodKinematicSecondsAtEpoch: (year: number) => number;
+    computeLodRealSecondsAtEpoch: (year: number) => number;
+    meanSolarDeltaAAtAge: (tMa: number, aApparentKm: number) => number | null;
+    meanMoonDistanceCorrectedAtAge: (tMa: number) => number | null;
+    meanMoonSiderealMonthAtAge: (tMa: number) => number | null;
+    meanSynodicMonthAtAge: (tMa: number) => number | null;
+    meanTropicalMonthAtAge: (tMa: number) => number | null;
+    meanLunarPerigeePrecessionAtAge: (tMa: number) => number | null;
+    meanLunarNodePrecessionAtAge: (tMa: number) => number | null;
+    meanAnomalisticMonthAtAge: (tMa: number) => number | null;
+    meanNodalMonthAtAge: (tMa: number) => number | null;
+    meanApsidalMeetsNodalAtAge: (tMa: number) => number | null;
+    meanLunarLevelingCycleAtAge: (tMa: number) => number | null;
+    meanStellarDayAtAge: (tMa: number) => number | null;
+    meanStellarDayProjectedAtAge: (tMa: number) => number | null;
+    meanSiderealDayAtAge: (tMa: number) => number | null;
+    meanPlanetSemiMajorAxisAtAge: (planetName: string, tMa: number) => number | null;
+    /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
+    meanPlanetOrbitalPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
+    /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
+    meanEarthPlanetSynodicPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
+}>;
+/**
+ * Moon-channel eccentricity sensitivities (perigee/node), the [g/g₀]^s
+ * exponents of the factored deep-time law — Meeus-effective (S_N moved
+ * 1.0 → 1.018 with the v4 frame-attribution batch). ONE home since B2 (they
+ * were literals in model.js, tools/lib/deep-time.js and script.js `_FW_MOON`).
+ */
+export const MOON_ECC_SENSITIVITY_PERIGEE: 2.407;
+export const MOON_ECC_SENSITIVITY_NODE: 1.018;
+/**
+ * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
+ * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin
+ * gates it on BOTH env vars, DT_CORRECTIONS and DT_RESONATOR), else by
+ * `resonator`. The four cycles' δLOD terms follow their flags only when
+ * `cycleLodFollowsFlags` is true (the browser's research toggles switch a
+ * cycle's ΔT and δLOD together); the Node twin sums them unconditionally.
+ * Their ΔT terms always follow their flags. Identical with every flag ON.
+ */
+export type DeltaTFlags = {
+    bond: boolean;
+    hallstatt: boolean;
+    jose5: boolean;
+    jose4: boolean;
+    resonator: boolean;
+    resonatorDeltaT?: boolean;
+    cycleLodFollowsFlags?: boolean;
+};

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

@@ -94,6 +94,13 @@ export type DeepTimeLodDeps = {
      * difference on moonDistanceMetresAtAge (the LOD-rate decomposition's tidal term).
      */
     moonDistanceSlopeMetresPerMaAtAge?: ((tMa: number) => number | null) | undefined;
+    /**
+     * -
+     * OPTIONAL cos ε̄(t)/cos ε̄(0) of the cited mean obliquity (earth/mean-
+     * obliquity-history) — the factor the precession RATE of date carries on the
+     * constant (plan 12 §8 option (a)); absent → 1, the constant alone, bit-identical.
+     */
+    meanObliquityCosRatioAtAge?: ((tMa: number) => number | null) | undefined;
 };
 /**
  * @typedef {Object} DeepTimeLodConstants
@@ -142,6 +149,10 @@ export type DeepTimeLodDeps = {
  *   OPTIONAL analytic slope of the distance history (metres per Ma, positive
  *   into the past means a growing a... sign as da/dt_Ma); absent → a central
  *   difference on moonDistanceMetresAtAge (the LOD-rate decomposition's tidal term).
+ * @property {(tMa: number) => number|null} [meanObliquityCosRatioAtAge] -
+ *   OPTIONAL cos ε̄(t)/cos ε̄(0) of the cited mean obliquity (earth/mean-
+ *   obliquity-history) — the factor the precession RATE of date carries on the
+ *   constant (plan 12 §8 option (a)); absent → 1, the constant alone, bit-identical.
  */
 /** @param {DeepTimeLodDeps} deps */
 export function createDeepTimeLod(deps: DeepTimeLodDeps): {
@@ -150,6 +161,9 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
     lodHoursAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
+    lunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    lunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+    precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
     lunarTorqueFactorAtAge: (tMa: number) => (number | null);
     precessionTorqueTermAtAge: (tMa: number) => (number | null);
     siderealYearSecondsAtAge: (t_Ma: number) => number;

package/types/deltat/moon-distance-artifact.d.cts

@@ -15,5 +15,8 @@ export const MOON_DISTANCE_ARTIFACT: Readonly<{
     aEarthRadii: number[];
     aMinEarthRadii: number[];
     aMaxEarthRadii: number[];
+    obliquityDeg: number[];
+    obliquityMinDeg: number[];
+    obliquityMaxDeg: number[];
 }>;
 export const MOON_DISTANCE_ARTIFACT_HASH: "850cd0b266d46a71";

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

@@ -25,9 +25,11 @@
  *   measured at the 500-yr artifact cadence: 0.16″/0.18″ rms vs IAU-2006
  *   over 1900–2100/1600–2400 vs the era tier's 0.31″/0.63″, and 0.0396°
  *   vs La2004 over −200 kyr vs the deep tier's 0.069° — the C-1 verdict,
- *   plan 02) and zetaModes serves only as the TAIL beyond the span (the
- *   seam at the span edge is the extraction residual, far outside every
- *   certified window). solarTorqueShareJ2000 and
+ *   plan 02) and zetaModes serves only as the TAIL beyond the span, joined
+ *   to the series at each edge so n̂ is continuous there (plan 12 §8.3 — the
+ *   J2000-point anchor left a 1.6° step at −10 Myr that kicked the hybrid's
+ *   mean obliquity +0.9° for the whole tail; far outside every certified window).
+ *   solarTorqueShareJ2000 and
  *   dynamicalEllipticityRatioAtYearFn: THE PRECESSION OF DATE — the two
  *   factors on the precession constant (see alphaAtGeneral below); absent,
  *   the constant is the secular law alone, bit-identical.

package/types/earth/mass-chain.d.cts

@@ -0,0 +1,78 @@
+/**
+ * The mass chain in §9 order: Δa (the Sun's tidal coupling to the Earth–Moon
+ * barycentric wobble, leading order), the corrected Moon distance, Kepler on
+ * the Moon's orbit → G(M_E + M_M), the DE440 split → GM_E and GM_M, Kepler on
+ * Earth's orbit → G(M_☉ + M_E), the subtraction → GM_☉ and M_☉.
+ * @param {{
+ *   moonDistanceKm: number, moonSiderealMonthDays: number, meanSiderealYearDays: number,
+ *   meanLengthOfDaySeconds: number, massRatioEarthMoon: number, currentAUDistanceKm: number,
+ *   gravitationalConstantKm3KgS2: number,
+ * }} c
+ * @returns {{
+ *   meanSiderealYearSeconds: number, moonOrbitalShiftKm: number, moonDistanceCorrectedKm: number,
+ *   gmEarthMoonSystemKm3S2: number, gmEarthAloneKm3S2: number, gmMoonAloneKm3S2: number,
+ *   gmSunPlusEarthKm3S2: number, gmSunKm3S2: number, sunMassKg: number,
+ * }}
+ */
+export function computeMassChainJ2000(c: {
+    moonDistanceKm: number;
+    moonSiderealMonthDays: number;
+    meanSiderealYearDays: number;
+    meanLengthOfDaySeconds: number;
+    massRatioEarthMoon: number;
+    currentAUDistanceKm: number;
+    gravitationalConstantKm3KgS2: number;
+}): {
+    meanSiderealYearSeconds: number;
+    moonOrbitalShiftKm: number;
+    moonDistanceCorrectedKm: number;
+    gmEarthMoonSystemKm3S2: number;
+    gmEarthAloneKm3S2: number;
+    gmMoonAloneKm3S2: number;
+    gmSunPlusEarthKm3S2: number;
+    gmSunKm3S2: number;
+    sunMassKg: number;
+};
+/**
+ * The Moon's J2000 periods on the IAU inputs. The two precession inputs are
+ * equinox-OF-DATE observables (3231.50 d / 6798.38 d; the legacy "ICRF" name
+ * is kept); the star-referenced periods remove the equinox's own precession,
+ * 1/T_E = 1/T_I ∓ 1/T_p (the perigee advances, the node regresses), T_p the
+ * J2000 precession anchor in days. The months follow kinematically.
+ * @param {{
+ *   moonSiderealMonthDays: number, moonApsidalPrecessionDaysInput: number, moonNodalPrecessionDaysInput: number,
+ *   axialPrecessionJ2000Days: number, meanSiderealYearDays: number, meanSolarYearDays: number,
+ * }} c
+ * @returns {{
+ *   moonApsidalPrecessionDaysEarth: number, moonNodalPrecessionDaysEarth: number,
+ *   moonApsidalJ2000Seconds: number, moonNodalJ2000Seconds: number, moonSiderealMonthJ2000Seconds: number,
+ *   moonAnomalisticMonthDays: number, moonNodalMonthDays: number, moonSynodicMonthDays: number, moonTropicalMonthDays: number,
+ *   moonFullMoonCycleEarthDays: number, moonFullMoonCycleICRFDays: number,
+ *   moonApsidalMeetsNodalDays: number, moonLunarLevelingCycleDays: number,
+ *   moonDraconicYearICRFDays: number, moonDraconicYearEarthDays: number,
+ * }}
+ */
+export function computeMoonPeriodsJ2000(c: {
+    moonSiderealMonthDays: number;
+    moonApsidalPrecessionDaysInput: number;
+    moonNodalPrecessionDaysInput: number;
+    axialPrecessionJ2000Days: number;
+    meanSiderealYearDays: number;
+    meanSolarYearDays: number;
+}): {
+    moonApsidalPrecessionDaysEarth: number;
+    moonNodalPrecessionDaysEarth: number;
+    moonApsidalJ2000Seconds: number;
+    moonNodalJ2000Seconds: number;
+    moonSiderealMonthJ2000Seconds: number;
+    moonAnomalisticMonthDays: number;
+    moonNodalMonthDays: number;
+    moonSynodicMonthDays: number;
+    moonTropicalMonthDays: number;
+    moonFullMoonCycleEarthDays: number;
+    moonFullMoonCycleICRFDays: number;
+    moonApsidalMeetsNodalDays: number;
+    moonLunarLevelingCycleDays: number;
+    moonDraconicYearICRFDays: number;
+    moonDraconicYearEarthDays: number;
+};

package/types/earth/mean-obliquity-history.d.cts

@@ -0,0 +1,32 @@
+export type MeanObliquityTable = {
+    /**
+     * - ascending, 0 first, uniform 1-Myr step
+     */
+    ageMa: ReadonlyArray<number>;
+    /**
+     * - the cited mean obliquity column, degrees
+     */
+    obliquityDeg: ReadonlyArray<number>;
+};
+/**
+ * @typedef {Object} MeanObliquityTable
+ * @property {ReadonlyArray<number>} ageMa - ascending, 0 first, uniform 1-Myr step
+ * @property {ReadonlyArray<number>} obliquityDeg - the cited mean obliquity column, degrees
+ */
+/**
+ * @param {{ table: MeanObliquityTable }} deps
+ * @returns {{
+ *   meanObliquityDegAtAge: (tMa: number) => number|null,
+ *   cosRatioAtAge: (tMa: number) => number|null,
+ *   meanObliquityJ2000Deg: number,
+ *   spanMa: number,
+ * }}
+ */
+export function createMeanObliquityHistory(deps: {
+    table: MeanObliquityTable;
+}): {
+    meanObliquityDegAtAge: (tMa: number) => number | null;
+    cosRatioAtAge: (tMa: number) => number | null;
+    meanObliquityJ2000Deg: number;
+    spanMa: number;
+};

package/types/earth/precession-composed.d.cts

@@ -24,14 +24,20 @@ export function computeSolarTorqueShare(c: {
  *   moonDistanceMetresAtAge: (tMa: number) => (number | null),
  *   moonDistanceJ2000Metres: number,
  *   yearToTMa: (year: number) => number,
+ *   meanObliquityCosRatioAtAge?: (tMa: number) => (number | null),
  * }} deps - p0 = the model's derived J2000 rate, 1,296,000/axial0 (S5; a
  *   number or a LAZY getter — the callers' certified year laws read the
  *   deep-time factory this is built inside, so the anchor is resolved on
  *   the first RATE use, never at construction; the torque term needs no
- *   anchor); lodJ2000Seconds the SAME day basis lodSecondsAtAge(0) returns.
+ *   anchor); lodJ2000Seconds the SAME day basis lodSecondsAtAge(0) returns;
+ *   meanObliquityCosRatioAtAge the cited mean obliquity's cos ε̄(t)/cos ε̄(0)
+ *   (OPTIONAL — absent, the rate of date is the constant alone, bit-identical).
  * @returns {{
  *   composedRateArcsecPerYrAtAge: (tMa: number) => (number | null),
  *   composedPeriodYearsAtAge: (tMa: number) => (number | null),
+ *   torqueRateArcsecPerYrAtAge: (tMa: number) => (number | null),
+ *   torquePeriodYearsAtAge: (tMa: number) => (number | null),
+ *   obliquityFactorAtAge: (tMa: number) => (number | null),
  *   composedRateArcsecPerYrAtYear: (year: number) => (number | null),
  *   composedPeriodYearsAtYear: (year: number) => (number | null),
  *   composedRateRatioAtAge: (tMa: number) => (number | null),
@@ -49,9 +55,13 @@ export function createComposedPrecession(deps: {
     moonDistanceMetresAtAge: (tMa: number) => (number | null);
     moonDistanceJ2000Metres: number;
     yearToTMa: (year: number) => number;
+    meanObliquityCosRatioAtAge?: (tMa: number) => (number | null);
 }): {
     composedRateArcsecPerYrAtAge: (tMa: number) => (number | null);
     composedPeriodYearsAtAge: (tMa: number) => (number | null);
+    torqueRateArcsecPerYrAtAge: (tMa: number) => (number | null);
+    torquePeriodYearsAtAge: (tMa: number) => (number | null);
+    obliquityFactorAtAge: (tMa: number) => (number | null);
     composedRateArcsecPerYrAtYear: (year: number) => (number | null);
     composedPeriodYearsAtYear: (year: number) => (number | null);
     composedRateRatioAtAge: (tMa: number) => (number | null);

package/types/index.d.ts

@@ -14,6 +14,7 @@ export { createSecularSeriesOverride } from "./planets/secular-series.cjs";
 export { createDeepEccChannel } from "./moon/deep-ecc-channel.cjs";
 export { createDeepOrbitalHistory } from "./earth/deep-orbital-history.cjs";
 export { createSiderealYearChannel } from "./earth/sidereal-year-channel.cjs";
+export { createMeanObliquityHistory } from "./earth/mean-obliquity-history.cjs";
 export { createPlanetModel } from "./planets/model.cjs";
 export { PLANET_DISPLAY_NAMES } from "./planets/display-names.cjs";
 export { createDeltaTCycles } from "./deltat/cycles.cjs";
@@ -58,6 +59,11 @@ import { REFERENCE_DATA } from './constants/index.js';
 import { assembleModel } from './model.js';
 export { GENERATED as DEFAULT_CONSTANTS, CONSTANTS_HASH, MODEL_VERSION, PREPRINT_DOI, REFERENCE_DATA };
 export { FITTED_COEFFICIENTS, COEFFICIENTS_HASH } from "./constants/index.js";
+export { GREGORIAN_START, GREGORIAN_START_JD, calendarToJD, jdToCalendar, jdToDecimalYear, dayFractionToHMS, julianCalendarToJD, jdToJulianCalendar } from "./time/calendar.cjs";
+export { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wrap180Closed, wrap2pi, wrapPi } from "./math/angles.cjs";
+export { createDeepTimeChain, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE } from "./deltat/deep-time-chain.js";
+export { computeMassChainJ2000, computeMoonPeriodsJ2000 } from "./earth/mass-chain.cjs";
+export { SECONDS_PER_DAY, DAYS_PER_JULIAN_YEAR, JULIAN_CENTURY_DAYS, SECONDS_PER_JULIAN_YEAR } from "./constants/units.cjs";
 export { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG } from "./earth/frame-of-date.cjs";
 export { derivePlanetGeometry, scaffoldOrbitInputs } from "./planets/geometry.cjs";
 export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, computeApsidalSecularDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from "./planets/keplerian-chain.cjs";

package/types/layer0/index.d.ts

@@ -1,7 +1,8 @@
-export function createEpochPrimitives({ params: p, alphaAtAgeMa, moonDistanceMetresAtAge }: {
+export function createEpochPrimitives({ params: p, alphaAtAgeMa, moonDistanceMetresAtAge, meanObliquityCosRatioAtAge }: {
     params: EpochParams;
     alphaAtAgeMa: (tMa: number) => number;
     moonDistanceMetresAtAge: (tMa: number) => number | null;
+    meanObliquityCosRatioAtAge?: (tMa: number) => number | null;
 }): EpochPrimitives;
 export type EpochParams = {
     /**
@@ -53,7 +54,7 @@ export type EpochPrimitives = {
     lodSeconds: (year: number) => number | null;
     siderealYearSeconds: (year: number) => number;
     /**
-     * sidereal less one equinox turn per T_p(t) (the composed precession period)
+     * sidereal less one equinox turn per T_p(t) (the composed precession period on the RATE of date)
      */
     tropicalYearSeconds: (year: number) => number;
 };

package/types/math/angles.d.cts

@@ -0,0 +1,18 @@
+/** π/180 — degrees to radians. */
+export const DEG2RAD: number;
+/** 180/π — radians to degrees. */
+export const RAD2DEG: number;
+/** Arcseconds in a full turn. */
+export const ARCSEC_PER_TURN: 1296000;
+/** Arcseconds per radian, 206264.806…  */
+export const ARCSEC_PER_RAD: number;
+/** @param {number} x degrees @returns {number} the same angle in [0, 360) */
+export function wrap360(x: number): number;
+/** @param {number} x degrees @returns {number} the same angle in [−180, 180) */
+export function wrap180(x: number): number;
+/** @param {number} x degrees @returns {number} the same angle in (−180, 180] */
+export function wrap180Closed(x: number): number;
+/** @param {number} x radians @returns {number} the same angle in [0, 2π) */
+export function wrap2pi(x: number): number;
+/** @param {number} x radians @returns {number} the same angle in [−π, π) */
+export function wrapPi(x: number): number;

package/types/model.d.ts

@@ -1,3 +1,9 @@
+/**
+ * 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
@@ -41,7 +47,7 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         lunarSmallParameter: number;
         moonMassFraction: number;
         solarTidalShiftKm: number;
-        moonDistanceCorrectedKm: any;
+        moonDistanceCorrectedKm: number;
         gmEarthMoonSystemKm3S2: number;
         gmEarthKm3S2: number;
         gmMoonKm3S2: number;
@@ -69,10 +75,18 @@ 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 — the composed torque rate's period on the derived J2000 anchor (25,771.4 at J2000). @param {number} year @returns {number} */
+        /** 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 composed rate ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³], ″/yr (p₀ = 50.2883, derived). @param {number} year @returns {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. */
@@ -142,6 +156,10 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
     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;

package/types/moon/arguments.d.cts

@@ -121,12 +121,7 @@ export function createMoonArguments({ constants, fns }: {
         };
     };
 };
-/**
- * Astronomical JD → decimal year: Julian calendar before 1582-10-15,
- * Gregorian after (Meeus Ch. 7 inverse). The S3-aligned coordinate for the
- * Sun secular deviations. @param {number} jd @returns {number}
- */
-export function jdToDecimalYear(jd: number): number;
+import { jdToDecimalYear } from "../time/calendar.cjs";
 /**
  * @typedef {{Lp: number, D: number, M: number, Mp: number, F: number}} MoonArgsDeg
  */
@@ -137,3 +132,4 @@ export function jdToDecimalYear(jd: number): number;
  *  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). */
 export const SECULAR_COMPLETION_TAPER_YEARS: 27943.083333333332;
+export { jdToDecimalYear };