@essrt/physics 4.3.0 → 4.4.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "4.3.0",
+  "version": "4.4.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": "v13.0",

package/src/cardinal/one-source-structure.cjs

@@ -101,6 +101,31 @@ function createCardinalStructure({ sampleAt, tropicalYearSecondsAtYearFn }) {
       + eocOffsetSeconds(year + 1, type) - eocOffsetSeconds(year, type);
   }
 
+  /**
+   * The ANOMALISTIC year — time between successive perihelion passages,
+   * SI seconds. At perihelion the equation of center is ZERO by definition
+   * (M = 0), so perihelion passages are pure mean-anomaly events and the
+   * closed form needs no EoC term at all:
+   *
+   *   T_anom(year) = T_mean(year+½) · 360 / (360 − Δϖ_yr)
+   *
+   * with Δϖ_yr = ϖ(year+1) − ϖ(year) (wrapped; the year-over-year apsidal
+   * advance RELATIVE TO THE EQUINOX — periOfDateDeg is equinox-referenced,
+   * so this is the climatic-precession rate, ~0.017°/yr ≈ +25 min at
+   * J2000). Exact difference form; driven by the ENGINE's ϖ(t) — this is
+   * the model's physics, replacing the retired K-family anomalistic
+   * harmonic fit (the family that measured worst against the one-source
+   * movement: 0.41 s broadband).
+   * @param {number} year the interval starts at this year's perihelion
+   * @returns {number} seconds
+   */
+  function anomalisticYearSeconds(year) {
+    const p0 = sampleAt(year).periOfDateDeg;
+    const p1 = sampleAt(year + 1).periOfDateDeg;
+    const dPeriDeg = ((p1 - p0 + 540) % 360) - 180;   // wrapped, prograde +
+    return tropicalYearSecondsAtYearFn(year + 0.5) * 360 / (360 - dPeriDeg);
+  }
+
   /**
    * The four year lengths minus the mean — the e(t)-proportional spread.
    * @param {number} year
@@ -116,7 +141,7 @@ function createCardinalStructure({ sampleAt, tropicalYearSecondsAtYearFn }) {
     return out;
   }
 
-  return { eocOffsetSeconds, yearLengthSeconds, spreadSeconds };
+  return { eocOffsetSeconds, yearLengthSeconds, spreadSeconds, anomalisticYearSeconds };
 }
 
 module.exports = { createCardinalStructure, equationOfCenterDeg, CARDINAL_LONGITUDE_DEG };

package/src/model.js

@@ -1175,6 +1175,7 @@ export function assembleModel(C, F, laws = {}) {
       yearLengthSeconds: /** @param {number} year @param {'VE'|'SS'|'AE'|'WS'} type @returns {number} */ (year, type) => cardinalStructureM.yearLengthSeconds(year, type),
       eocOffsetSeconds: /** @param {number} year @param {'VE'|'SS'|'AE'|'WS'} type @returns {number} */ (year, type) => cardinalStructureM.eocOffsetSeconds(year, type),
       spreadSeconds: /** @param {number} year */ (year) => cardinalStructureM.spreadSeconds(year),
+      anomalisticYearSeconds: /** @param {number} year @returns {number} */ (year) => cardinalStructureM.anomalisticYearSeconds(year),
     }),
     moon: Object.freeze({
       distanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceMetresAtAge(yearToTMa(year)) / 1000,