@essrt/physics 1.5.0 → 1.6.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.5.0",
+  "version": "1.6.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core \u2014 pure, injectable-constants, no I/O, no globals, no DOM. The complete model incl. fitted coefficients; every published version immutably ships one recorded model identity.",
   "essrt": {
     "modelVersion": "v11.0",

package/src/model.js

@@ -783,7 +783,20 @@ export function assembleModel(C, F) {
       argsAt: /** @param {number} jdTT */ (jdTT) => moonArgs.argsAt(jdTT),
       eFactorForD: fwEFactor,
       eFactorAtJdTT: /** @param {number} jdTT */ (jdTT) => fwEFactor(jdTT - j2000JD),
-      getMoonDistanceKm: () => moonDistanceKm,
+      // 20.3d(i): the Driver-1 ratio at the EVALUATED epoch (per-jd pure
+      // evaluator). MATCHED TRIPLE with the tools-lib and browser getters:
+      // identical age arithmetic — the SI-linear year (yearFromJD, never
+      // the calendar year: the recorded linear-vs-calendar mirror trap)
+      // against the MODEL-START anchor (startModelYearWithCorrection, the
+      // browser's J2000_CALENDAR_YEAR — NOT this model's 2000.0 yearToTMa
+      // convention: the 0.4977-yr difference is 4.95e-11 of distance and
+      // failed the bit-exact parity gate). Falls back to the J2000
+      // constant without a jd (legacy call shape) or past the domain.
+      getMoonDistanceKm: /** @param {number} [jdTT] */ (jdTT) => {
+        if (jdTT === undefined) return moonDistanceKm;
+        const d = moonDistanceMetresAtAge((startModelYearWithCorrection - yearFromJD(jdTT)) / 1e6);
+        return d === null ? moonDistanceKm : d / 1000;
+      },
       getEccentricityBase: () => C.moonReference.moonOrbitalEccentricityBase,
       deltaTSeconds: /** @param {number} jd */ (jd) => (jdTTFromUT(jd) - jd) * 86400,
       jdToSIyear: yearFromJD,

package/src/moon/series.cjs

@@ -48,7 +48,7 @@
  *     argsAt: (jdTT: number) => {Lp: number, D: number, M: number, Mp: number, F: number},
  *     eFactorForD: (dDaysTT: number, T: number, T2: number) => number,
  *     eFactorAtJdTT: (jdTT: number, T: number, T2: number) => number,
- *     getMoonDistanceKm: () => number,
+ *     getMoonDistanceKm: (jdTT?: number) => number,
  *     getEccentricityBase: () => number,
  *     deltaTSeconds: (jdUT: number) => number,
  *     jdToSIyear: (jd: number) => number,
@@ -147,9 +147,17 @@ function createMoonSeries({ constants, fns }) {
    *  @param {number} Dr @param {number} Mr @param {number} Mpr
    *  @param {number} Fr @param {number} E @param {number} E2
    *  @returns {number} km */
-  function fullDistanceKm(Dr, Mr, Mpr, Fr, E, E2) {
+  // 20.3d(i): the Driver-1 ratio is evaluated AT THE EVALUATED EPOCH — the
+  // getter receives the evaluation jdTT so each engine can answer with its
+  // pure a_M(t) evaluator (per-jd deterministic; previously the browser
+  // answered with whatever epoch the render loop last set — a scene-state
+  // timing accident — while the Node engines pinned the J2000 constant).
+  // Engines whose getter ignores the argument keep their old behavior.
+  /** @param {number} Dr @param {number} Mr @param {number} Mpr @param {number} Fr
+   *  @param {number} E @param {number} E2 @param {number} [jdTT] @returns {number} */
+  function fullDistanceKm(Dr, Mr, Mpr, Fr, E, E2, jdTT) {
     const Sr = sumTableCos(moonR, Dr, Mr, Mpr, Fr, E, E2);
-    return (moonRMeanKm + Sr * 1e-3) * (getMoonDistanceKm() / moonDistanceJ2000Km);
+    return (moonRMeanKm + Sr * 1e-3) * (getMoonDistanceKm(jdTT) / moonDistanceJ2000Km);
   }
 
   /** The PRODUCTION scene evaluation at d days TT from J2000 (the caller
@@ -192,7 +200,7 @@ function createMoonSeries({ constants, fns }) {
     const fullSl = Sl + (2 * eocHalf / D2R * 1e6) * Math.sin(Mpr)
                      + (1.25 * eocHalf * eocHalf / D2R * 1e6) * Math.sin(2 * Mpr);
     const lonDeg = Lp / D2R + fullSl * 1e-6 + moonMeeusLpCorrectionDeg;
-    const distKm = fullDistanceKm(Dr, Mr, Mpr, Fr, E, E2);
+    const distKm = fullDistanceKm(Dr, Mr, Mpr, Fr, E, E2, j2000JD + d);
     return { thetaAddRad, lonDeg, latRad, latDeg, distKm, T };
   }
 
@@ -242,7 +250,7 @@ function createMoonSeries({ constants, fns }) {
     const args = argsAt(jdTT);
     const Dr = args.D * D2R, Mr = args.M * D2R, Mpr = args.Mp * D2R, Fr = args.F * D2R;
     const E = eFactorAtJdTT(jdTT, T, T * T);
-    return fullDistanceKm(Dr, Mr, Mpr, Fr, E, E * E);
+    return fullDistanceKm(Dr, Mr, Mpr, Fr, E, E * E, jdTT);
   }
 
   return { sceneEvalAt, truncatedLonDeg, truncatedBetaDeg, truncatedDistanceKm, additionalArgs };