@essrt/physics 1.16.0 → 2.0.0

package/src/constants/generated.d.ts

@@ -2,14 +2,14 @@
 // Gives the TypeScript website full type safety at the boundary (§2g) while
 // packages/physics stays JavaScript.
 
-export declare const CONSTANTS_HASH: "844b86db470b9937";
+export declare const CONSTANTS_HASH: "ccc3b87d48a05bcd";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "844b86db470b9937";
+  readonly hash: "ccc3b87d48a05bcd";
   readonly additionalBodies: {
     "pluto": {
       "name": string;

package/src/constants/generated.js

@@ -32,17 +32,17 @@
  * carry it so a counterfactual is reproducible (§2d).
  * @type {string}
  */
-export const CONSTANTS_HASH = "844b86db470b9937";
+export const CONSTANTS_HASH = "ccc3b87d48a05bcd";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
-export const MODEL_VERSION = "v11.0";
+export const MODEL_VERSION = "v12.0";
 
 /** Canonical preprint DOI — single source: public/input/model-version.json. */
 export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";
 
 /** @type {Readonly<Record<string, unknown>>} */
 export const DEFAULT_CONSTANTS = Object.freeze({
-  hash: "844b86db470b9937",
+  hash: "ccc3b87d48a05bcd",
   additionalBodies: {
     "pluto": {
       "name": "Pluto",
@@ -207,8 +207,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     }
   },
   earth: {
-    "earthtiltMean": 23.413527499101747,
-    "earthInvPlaneInclinationAmplitude": 0.6360475695866625,
+    "earthtiltMean": 23.413526608469965,
+    "earthInvPlaneInclinationAmplitude": 0.636069722432876,
     "eccentricityBase": 0.015386009686374918,
     "eccentricityAmplitude": 0.0013559440307290061
   },
@@ -228,7 +228,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "perihelionPassageJ2000_JD": 2451547.042,
     "earthInclinationJ2000_deg": 1.57869,
     "earthInclinationRate_arcsecPerCentury": -18,
-    "deltaTStart": 55.85417672145975,
+    "deltaTStart": 55.47104318258222,
     "sunTilt": 7.155
   },
   foundational: {
@@ -238,7 +238,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "startmodelJD": 2451716.5,
     "startmodelYear": 2000.5,
     "correctionDays": -0.828832119703292,
-    "correctionSun": 0.4968839022012559,
+    "correctionSun": 0.4971540149655503,
     "temperatureGraphMostLikely": 14.5,
     "startAngleModel": 89.91949879,
     "useVariableSpeed": true,
@@ -248,7 +248,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "moonStartposApsidal": 347.5476,
     "moonStartposNodal": 64.0435,
     "moonStartposMoon": 67.8443,
-    "moonMeeusLpCorrection": 0.009081766072295829
+    "moonMeeusLpCorrection": 0.00908174470925335
   },
   moonMeeus: {
     "elpW1T2Decomposition_arcsecPerCy2": {
@@ -432,8 +432,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "mercury": {
       "name": "Mercury",
       "eocFraction": -0.527,
-      "startpos": 83.65212772217342,
-      "angleCorrection": 0.971595956003533,
+      "startpos": 83.64651142941864,
+      "angleCorrection": 0.9727205912043573,
       "perihelionEclipticFraction": [
         8,
         11
@@ -462,7 +462,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "name": "Venus",
       "eocFraction": 0.436,
       "startpos": 249.28769121856186,
-      "angleCorrection": -2.7506208719906127,
+      "angleCorrection": -2.789778861448692,
       "perihelionEclipticFraction": [
         -8,
         6
@@ -488,7 +488,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "name": "Mars",
       "eocFraction": -0.066224,
       "startpos": 121.4634461571797,
-      "angleCorrection": -2.110262688821635,
+      "angleCorrection": -2.104652936098656,
       "perihelionEclipticFraction": [
         8,
         36
@@ -517,7 +517,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "name": "Jupiter",
       "eocFraction": 0.495,
       "startpos": 13.887946755371052,
-      "angleCorrection": 0.9306112039524912,
+      "angleCorrection": 0.9333205235053406,
       "perihelionEclipticFraction": [
         8,
         39
@@ -545,8 +545,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "saturn": {
       "name": "Saturn",
       "eocFraction": 0.54,
-      "startpos": 11.279046793476084,
-      "angleCorrection": -0.17887349563889998,
+      "startpos": 11.279308392461772,
+      "angleCorrection": -0.1784695052193328,
       "perihelionEclipticFraction": [
         -8,
         65
@@ -575,7 +575,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "name": "Uranus",
       "eocFraction": 0.53,
       "startpos": 44.900388945775106,
-      "angleCorrection": -0.7329070274012324,
+      "angleCorrection": -0.7343433510960402,
       "perihelionEclipticFraction": [
         1,
         3
@@ -605,7 +605,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "name": "Neptune",
       "eocFraction": 0.585,
       "startpos": 47.9552024382285,
-      "angleCorrection": 2.332347646694906,
+      "angleCorrection": 2.3327665835675617,
       "perihelionEclipticFraction": [
         2,
         1
@@ -764,7 +764,7 @@ export const REFERENCE_DATA = Object.freeze({
     "sunModelTrueErrorDeg": 0.003,
     "sunTropicalYearDiffSeconds": 0.1,
     "sunSiderealYearDiffSeconds": 0.02,
-    "frameworkSunVsJplRmsArcsec": 0.93,
+    "frameworkSunVsJplRmsArcsec": 0.8,
     "meeusCh25SunVsJplRmsArcsec": 1.22,
     "moonSeriesLonVsJplRmsArcsec": 2.84,
     "moonSeriesLatVsJplRmsArcsec": 0.35,

package/src/eclipse/sun-planetary-completion.cjs

@@ -110,79 +110,83 @@ const ARG_D0 = 297.8501921;
  * [[6 mean-longitude multipliers lMe,lV,lE,lM,lJ,lS],
  *  [6 mean-anomaly multipliers MMe,MV,ME,MMa,MJ,MS], cos″, sin″].
  * 70 terms ≥ 0.05″ from the N3 framework-carrier extraction (79-term
- * LSQ on the D2 derived signal, fidelity 0.61″).
+ * LSQ on the D2 derived signal, fidelity 0.616″) — re-derived under the
+ * ONE eccentricity law (unification; D2 combined residual 0.60″,
+ * composed difference vs the previous table ≤1.5″ at the ancient presets,
+ * detrended 0.033 min). Individual coefficients redistribute between the
+ * near-degenerate ±M sideband families; only the composed function ships.
  * @type {Array<[number[], number[], number, number]>}
  */
 const TERMS = [
-  [[0, 0, 1, 0, -2, 0], [0, 0, 1, 0, 0, 0], 23.4846, 2.3099],
-  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, 0, 0], -1.7235, -14.0062],
-  [[0, 0, 2, 0, -3, 0], [0, 0, 0, 0, 1, 0], 7.4894, -7.7252],
-  [[0, 3, -3, 0, 0, 0], [0, 0, -1, 0, 0, 0], 5.7813, 1.7078],
-  [[0, 3, -4, 0, 0, 0], [0, 0, 0, 0, 0, 0], 1.1877, 5.8649],
-  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, 0, 0], 1.6357, -4.9485],
-  [[0, 2, -3, 0, 0, 0], [0, 0, 1, 0, 0, 0], -4.9053, -0.3900],
-  [[0, 1, -1, 0, 0, 0], [0, 0, 0, 0, 0, 0], 0.1351, 4.8331],
-  [[0, 0, 2, -3, 0, 0], [0, 0, 0, 1, 0, 0], -0.2239, -3.1687],
-  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, 1, 0], -2.1749, -1.7705],
-  [[0, 0, 1, -1, 0, 0], [0, 0, 0, -1, 0, 0], 0.8107, -2.5127],
-  [[0, 0, 1, 0, -1, 0], [0, 0, -1, 0, 0, 0], -2.5403, -0.1821],
-  [[0, 0, -1, 2, 0, 0], [0, 0, 0, 0, 0, 0], -1.3152, -2.0575],
-  [[0, 2, -2, 0, 0, 0], [0, 0, 0, 0, 0, 0], -1.8073, -0.8142],
-  [[0, 3, -3, 0, 0, 0], [0, -1, 0, 0, 0, 0], -0.8611, 1.7708],
-  [[0, 0, -1, 2, 0, 0], [0, 0, -1, 0, 0, 0], 1.2659, 1.1195],
-  [[0, 0, 2, -2, 0, 0], [0, 0, -1, 0, 0, 0], 1.1646, 1.0865],
-  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, -1, 0], 0.4560, -1.2051],
-  [[0, 0, -2, 4, 0, 0], [0, 0, 0, -1, 0, 0], -0.2198, 1.2107],
-  [[0, 3, -3, 0, 0, 0], [0, 0, 0, 0, 0, 0], 0.5260, -1.0137],
-  [[0, 2, -2, 0, 0, 0], [0, 0, -1, 0, 0, 0], -1.1195, 0.1107],
-  [[0, 0, 2, -2, 0, 0], [0, 0, 0, 0, 0, 0], 0.8202, -0.7125],
-  [[0, 3, -4, 0, 0, 0], [0, 0, -1, 0, 0, 0], 0.9722, 0.3524],
-  [[0, 0, 2, -3, 0, 0], [0, 0, 0, 0, 0, 0], 0.9763, -0.2018],
-  [[0, 3, -4, 0, 0, 0], [0, 0, 1, 0, 0, 0], 0.1162, 0.9095],
-  [[0, 0, 2, -2, 0, 0], [0, 0, 0, -1, 0, 0], 0.2940, 0.8000],
-  [[0, 2, -3, 0, 0, 0], [0, 0, 0, 0, 0, 0], 0.2210, -0.6939],
-  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, -1, 0], -0.2188, 0.6572],
-  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, 0, 0], 0.6135, 0.3150],
-  [[0, 1, -1, 0, 0, 0], [0, -1, 0, 0, 0, 0], 0.6161, -0.1008],
-  [[0, 0, 1, -1, 0, 0], [0, 0, 0, 1, 0, 0], 0.4387, -0.4123],
-  [[0, 0, 2, 0, -2, 0], [0, 0, -1, 0, 0, 0], -0.1825, 0.5176],
-  [[0, 2, -3, 0, 0, 0], [0, 1, 0, 0, 0, 0], 0.1125, -0.5348],
-  [[0, 0, -2, 4, 0, 0], [0, 0, -1, 0, 0, 0], -0.4800, 0.2245],
-  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, 0], 0.0845, -0.4012],
-  [[0, 0, 1, 0, -1, 0], [0, 0, 1, 0, 0, 0], 0.1622, -0.3517],
-  [[0, 0, 2, 0, -3, 0], [0, 0, -1, 0, 0, 0], 0.3786, 0.0361],
-  [[0, 0, 2, -3, 0, 0], [0, 0, 0, -1, 0, 0], 0.0701, -0.3556],
-  [[0, 0, 1, -1, 0, 0], [0, 0, 0, 0, 0, 0], -0.2643, -0.2045],
-  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, 1, 0], 0.2650, -0.1462],
-  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, 1], 0.3016, 0.0060],
-  [[0, 0, -1, 2, 0, 0], [0, 0, 0, 1, 0, 0], -0.0582, -0.2912],
-  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, 1, 0], -0.2051, -0.2135],
-  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, -1, 0], 0.0858, 0.2762],
-  [[0, 0, 2, -3, 0, 0], [0, 0, -1, 0, 0, 0], 0.1911, 0.2163],
-  [[0, 0, -1, 2, 0, 0], [0, 0, 0, -1, 0, 0], 0.2815, -0.0450],
-  [[0, 2, -3, 0, 0, 0], [0, 0, -1, 0, 0, 0], 0.2490, 0.0551],
-  [[0, 0, 1, 0, 0, -1], [0, 0, -1, 0, 0, 0], -0.0456, 0.2371],
-  [[0, 0, 2, 0, -3, 0], [0, 0, 0, 0, 0, 0], 0.1560, -0.1469],
-  [[0, 3, -4, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.1412, 0.1568],
-  [[0, 3, -4, 0, 0, 0], [0, -1, 0, 0, 0, 0], -0.1616, -0.1314],
-  [[0, 0, -2, 4, 0, 0], [0, 0, 0, 0, 0, 0], 0.1728, -0.0849],
-  [[0, 1, -1, 0, 0, 0], [0, 0, -1, 0, 0, 0], -0.0449, 0.1649],
+  [[0, 3, -3, 0, 0, 0], [0, 0, -1, 0, 0, 0], -9.2558, -0.3232],
+  [[0, 3, -4, 0, 0, 0], [0, 0, 0, 0, 0, 0], -0.1950, -9.2459],
+  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, 1, 0], -1.4662, -6.4341],
+  [[0, 1, -1, 0, 0, 0], [0, 0, 0, 0, 0, 0], 0.1345, 4.8332],
+  [[0, 2, -2, 0, 0, 0], [0, 0, 0, 0, 0, 0], -2.1364, -2.9549],
+  [[0, 2, -3, 0, 0, 0], [0, 0, 1, 0, 0, 0], -2.8927, -1.1913],
+  [[0, 0, 1, 0, -2, 0], [0, 0, 1, 0, 0, 0], 2.8982, -1.0900],
+  [[0, 2, -2, 0, 0, 0], [0, 0, -1, 0, 0, 0], 0.7820, 2.8796],
+  [[0, 0, 1, 0, -1, 0], [0, 0, -1, 0, 0, 0], -2.5401, -0.1829],
+  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, 0, 0], -0.2456, -2.1163],
+  [[0, 0, 1, -1, 0, 0], [0, 0, 0, -1, 0, 0], 2.1028, 0.1572],
+  [[0, 0, 2, -2, 0, 0], [0, 0, -1, 0, 0, 0], 2.0062, 0.4493],
+  [[0, 0, 2, 0, -3, 0], [0, 0, 0, 0, 1, 0], 0.2164, 1.7839],
+  [[0, 2, -3, 0, 0, 0], [0, 0, 0, 0, 0, 0], -1.5284, 0.8126],
+  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, -1, 0], -0.9353, -1.1855],
+  [[0, 0, 2, -3, 0, 0], [0, 0, 0, 1, 0, 0], 0.4254, -1.4438],
+  [[0, 3, -3, 0, 0, 0], [0, -1, 0, 0, 0, 0], -1.2408, 0.7438],
+  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, 0, 0], 1.3145, -0.1985],
+  [[0, 0, 2, 0, -2, 0], [0, 0, -1, 0, 0, 0], -0.4728, 1.1130],
+  [[0, 3, -4, 0, 0, 0], [0, 0, -1, 0, 0, 0], 0.9746, 0.3524],
+  [[0, 0, -1, 2, 0, 0], [0, 0, -1, 0, 0, 0], -0.9141, -0.1697],
+  [[0, 3, -4, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.7605, 0.1501],
+  [[0, 0, -1, 2, 0, 0], [0, 0, 0, -1, 0, 0], 0.6794, -0.0282],
+  [[0, 0, -1, 2, 0, 0], [0, 0, 0, 0, 0, 0], -0.6021, 0.2300],
+  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, 0, 0], -0.2055, -0.6053],
+  [[0, 0, 1, -1, 0, 0], [0, 0, 0, 0, 0, 0], -0.6212, -0.0277],
+  [[0, 0, 2, 0, -3, 0], [0, 0, 0, 0, 0, 0], 0.0277, 0.5852],
+  [[0, 0, -2, 4, 0, 0], [0, 0, -1, 0, 0, 0], -0.4801, 0.2245],
+  [[0, 0, 2, -2, 0, 0], [0, 0, 0, -1, 0, 0], 0.0605, -0.5014],
+  [[0, 0, -2, 4, 0, 0], [0, 0, 0, 0, 0, 0], 0.4682, 0.1671],
+  [[0, 0, 2, -2, 0, 0], [0, 0, 0, 0, 0, 0], 0.2949, 0.3878],
+  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, 0], 0.0845, -0.4011],
+  [[0, 2, -3, 0, 0, 0], [0, 1, 0, 0, 0, 0], 0.0589, -0.3744],
+  [[0, 3, -3, 0, 0, 0], [0, 0, 0, 0, 0, 0], -0.3257, -0.1827],
+  [[0, 0, 2, 0, -3, 0], [0, 0, -1, 0, 0, 0], 0.0149, 0.3557],
+  [[0, 0, -2, 4, 0, 0], [0, 0, 0, -1, 0, 0], -0.3152, -0.0762],
+  [[0, 3, -4, 0, 0, 0], [0, -1, 0, 0, 0, 0], 0.0714, 0.3142],
+  [[0, 0, -1, 2, 0, 0], [0, 0, 0, 1, 0, 0], 0.1331, -0.2622],
+  [[0, 0, 2, -3, 0, 0], [0, 0, -1, 0, 0, 0], 0.2826, 0.0458],
+  [[0, 2, -3, 0, 0, 0], [0, 0, -1, 0, 0, 0], -0.1689, -0.2243],
+  [[0, 0, 1, 0, -1, 0], [0, 0, 1, 0, 0, 0], -0.0262, -0.2757],
+  [[0, 0, 2, -3, 0, 0], [0, 0, 0, 0, 0, 0], 0.2267, -0.1326],
+  [[0, 0, 1, 0, -2, 0], [0, 0, 0, 0, -1, 0], -0.2245, -0.0951],
+  [[0, 0, 1, 0, 0, -1], [0, 0, -1, 0, 0, 0], -0.0453, 0.2368],
+  [[0, 2, -2, 0, 0, 0], [0, -1, 0, 0, 0, 0], -0.2075, 0.0622],
+  [[0, 3, -4, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.1413, 0.1567],
+  [[0, 1, -1, 0, 0, 0], [0, 0, -1, 0, 0, 0], 0.0858, -0.1399],
   [[0, 0, 2, 0, -3, 0], [0, 0, 1, 0, 0, 0], 0.1632, -0.0164],
-  [[0, 2, -2, 0, 0, 0], [0, -1, 0, 0, 0, 0], 0.0534, -0.1430],
-  [[0, 0, 2, -3, 0, 0], [0, 0, 1, 0, 0, 0], 0.1316, -0.0275],
-  [[0, 0, -2, 4, 0, 0], [0, 0, 1, 0, 0, 0], 0.1216, 0.0451],
-  [[0, 0, 2, -2, 0, 0], [0, 0, 0, 1, 0, 0], -0.1046, 0.0620],
-  [[0, 0, 1, -1, 0, 0], [0, 0, 1, 0, 0, 0], -0.1075, -0.0488],
-  [[0, 0, 2, 0, 0, -2], [0, 0, 0, 0, 0, 0], -0.0385, 0.1030],
-  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, -1], -0.0308, 0.1049],
-  [[0, 1, -1, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.0875, -0.0250],
+  [[0, 0, 2, -3, 0, 0], [0, 0, 1, 0, 0, 0], 0.1316, -0.0276],
+  [[0, 0, -2, 4, 0, 0], [0, 0, 1, 0, 0, 0], 0.1219, 0.0452],
+  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, 1, 0], -0.1243, -0.0271],
+  [[0, 1, -1, 0, 0, 0], [0, -1, 0, 0, 0, 0], 0.1139, -0.0172],
+  [[0, 0, 2, 0, 0, -2], [0, 0, 0, 0, 0, 0], -0.0385, 0.1031],
+  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, -1], -0.0309, 0.1052],
+  [[0, 0, 1, 0, -1, 0], [0, 0, 0, 0, 1, 0], 0.0945, 0.0475],
+  [[0, 0, 2, -3, 0, 0], [0, 0, 0, -1, 0, 0], -0.0981, -0.0051],
+  [[0, 0, 1, 0, 0, -1], [0, 0, 0, 0, 0, 1], -0.0827, 0.0500],
+  [[0, 0, 2, 0, -2, 0], [0, 0, 0, 0, -1, 0], 0.0866, -0.0204],
+  [[0, 1, -1, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.0085, -0.0855],
   [[0, 0, 1, 0, -2, 0], [0, 0, -1, 0, 0, 0], -0.0738, 0.0295],
   [[0, 0, 2, 0, -3, 0], [0, 0, 0, 0, -1, 0], 0.0097, 0.0787],
-  [[0, 1, -1, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.0016, 0.0762],
-  [[0, 2, -2, 0, 0, 0], [0, 0, 1, 0, 0, 0], 0.0562, -0.0490],
-  [[0, 3, -3, 0, 0, 0], [0, 0, 1, 0, 0, 0], 0.0638, -0.0362],
-  [[0, 0, 2, 0, -2, 0], [0, 0, 1, 0, 0, 0], -0.0085, 0.0686],
-  [[0, 2, -2, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.0453, 0.0505],
+  [[0, 1, -1, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.0016, 0.0764],
+  [[0, 0, 1, -1, 0, 0], [0, 0, 0, 1, 0, 0], 0.0655, 0.0358],
+  [[0, 3, -3, 0, 0, 0], [0, 0, 1, 0, 0, 0], 0.0517, -0.0513],
+  [[0, 0, 2, 0, -2, 0], [0, 0, 1, 0, 0, 0], -0.0085, 0.0687],
+  [[0, 2, -2, 0, 0, 0], [0, 1, 0, 0, 0, 0], -0.0274, 0.0580],
+  [[0, 0, 2, -2, 0, 0], [0, 0, 0, 1, 0, 0], 0.0521, 0.0217],
+  [[0, 0, 1, -1, 0, 0], [0, 0, 1, 0, 0, 0], 0.0188, 0.0523],
+  [[0, 2, -2, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.0421, -0.0340],
   [[0, 0, 1, -1, 0, 0], [0, 0, -1, 0, 0, 0], -0.0353, -0.0379],
 ];
 
@@ -240,7 +244,7 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy })
  *  at derivation time — the matched-pair fingerprint asserted by test:model.
  *  Unchanged from v1/v2: SUN_HARMONICS did not move in the D2 or N3
  *  landings. */
-const PAIRED_SUN_HARMONICS_SHA256 = 'e2cf42e9770c9e0a';
+const PAIRED_SUN_HARMONICS_SHA256 = 'cbc189cea1c20292';   // eccentricity unification: Step-0 refit → N2/N3 re-derived (fidelity 0.616″, 70 terms)
 
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
  *  full-precision carrier rates (deg/cy TT) the TERMS table was extracted

package/src/model.js

@@ -86,7 +86,10 @@ export function assembleModel(C, F) {
 
   const earthtiltMean = C.earth.earthtiltMean;
   const earthInclAmplitude = C.earth.earthInvPlaneInclinationAmplitude;
-  const eccentricityBase = C.earth.eccentricityBase;
+  // eccentricityAmplitude: the v11 H/16-beat amplitude, RETAINED ONLY as the
+  // Law-4 K calibration input (plan IP-eccentricity-unification, decision D2:
+  // K keeps its value). Earth's eccentricity law no longer uses it — see
+  // eccentricityAt (base' derived) below.
   const eccentricityAmplitude = C.earth.eccentricityAmplitude;
   const earthRAAngle = 2 * earthInclAmplitude - (earthInclAmplitude * earthInclAmplitude) / earthtiltMean;
   const earthInclMean = C.earthOrbital.earthInclinationJ2000_deg
@@ -305,13 +308,46 @@ export function assembleModel(C, F) {
     }
     return obliq;
   };
+  // ECCENTRICITY UNIFICATION — ONE law for the whole model (plan
+  // IP-eccentricity-unification, decision D1 = form (e)). Earth's orbital
+  // eccentricity is the magnitude of a vector sum, e = |v₁ + v₂|: v₁ of
+  // length base' along the ICRF perihelion (rotating once per H/3), v₂ of
+  // length base'/2 FIXED in space along the inclination-cycle anchor
+  // direction. Their relative angle θ rides the System-Reset lattice phase
+  // (θ(t) = 3(t − balancedYear)/H·360° − 180°, θ(J2000) = 81.178°) — the
+  // SAME anchor the inclination law and the Moon channel ride (their
+  // extremes coincide: max −23,200, min −79,100). Then
+  //     e(t) = base' · (1 + cos θ(t) / 2),   base' = e(J2000) / (1 + cos θ(J2000) / 2)
+  // base' is DERIVED from the observed J2000 eccentricity and the anchor
+  // (0.015520; the Law-5 balance moves 99.8636 → 99.8645%, Earth's weight
+  // there being 0.05%). No new constants. Why H/3 and not the former H/16
+  // beat law: eccentricity is frame-invariant and may carry only fixed-frame
+  // lattice periods; H/16 is the OF-DATE perihelion period (13 + 3 = 16 —
+  // the H/3 rotation seen from the H/13 equinox) and belongs to ϖ_of-date;
+  // measured, the H/16 law's present ė (−0.84e-5/cy) is 5× below the
+  // observed −4.20e-5 while this law reads −4.31e-5 (Phase-0 record:
+  // tools/explore/fq7s-h3-law-candidate.mjs; JPL Sun 1.49″, registry 0.80″,
+  // syzygy 3.72″). Consumers: the eclipse Sun (equation of centre), the
+  // besselian Sun distance, the Moon channel (E-factor, perigee/node T²),
+  // and the cardinal-point braid — one eccentricity everywhere.
+  // ONE implementation for all three runtimes: moon/ecc-channel.cjs (the
+  // Node engine's deep-time.js and the browser's script.js instantiate the
+  // same channel with the same inputs). Its phase counter runs from J2000
+  // (θ₀ = ϖ_ICRF(J2000) − 21.77° = 81.178°, the System-Reset anchor in
+  // anchor form — doc 66 §1); base' is derived inside from the observed
+  // J2000 eccentricity.
+  const moonEcc = createMoonEccChannel({
+    cyclesBetween,
+    eccentricityBase: C.earth.eccentricityBase,
+    perihelionLongitudeJ2000Deg: C.earthOrbital.earthPerihelionLongitudeJ2000,
+    inclinationCycleAnchorDeg: C.earthOrbital.earthInclinationCycleAnchor,
+    eccentricityJ2000: C.earthOrbital.earthEccentricityJ2000,
+  });
   /** @param {number} year @returns {number} */
-  const eccentricityAt = (year) => {
-    const th = phaseRadians(balancedYear, year, 16);
-    return Math.sqrt(eccentricityBase * eccentricityBase
-      + eccentricityAmplitude * eccentricityAmplitude
-      - 2 * eccentricityBase * eccentricityAmplitude * Math.cos(th));
-  };
+  const eccentricityAt = (year) => moonEcc.eccAt(year - 2000);
+  /** de/dyear of the one law — the cardinal braid's equation-of-centre
+   *  derivative rides it. @param {number} year @returns {number} */
+  const eccentricityRateAt = (year) => moonEcc.eccRateAt(year - 2000);
   /** @param {number} year @returns {number} */
   const inclinationDeg = (year) => earthInclMean
     - earthInclAmplitude * Math.cos(phaseRadians(balancedYear, year, 3));
@@ -379,8 +415,6 @@ export function assembleModel(C, F) {
       balancedYear,
       meanSolarYearDays,
       hJ2000: H,
-      eccentricityBase,
-      eccentricityAmplitude,
       tiltMeanDeg: earthtiltMean,
       raAngleDeg: earthRAAngle,
       inclAmplitudeDeg: earthInclAmplitude,
@@ -396,6 +430,7 @@ export function assembleModel(C, F) {
       meanHAtAgeMa: (tMa) => deepLod.hAtAge(tMa),
       meanYearRealLodDays: (tMa) => deepLod.yearInDaysAtAge(tMa),
       eccentricityAt,
+      eccentricityRateAt,
     },
   });
   /** Tropical year: mean of the four cardinal intervals. @param {number} year @returns {number} */
@@ -612,14 +647,9 @@ export function assembleModel(C, F) {
   const moonAnomalisticMonthDays = totalDaysInH / (nSid - nApsidalEJ2000);
   const moonSynodicMonthDays = totalDaysInH / (nSid + 13 - H);
 
-  // The framework H/3 eccentricity line (the Moon channel's view of the
-  // wobble movement — NOT the H/16 orbital eccentricity above)
-  const moonEcc = createMoonEccChannel({
-    cyclesBetween,
-    eccentricityBase,
-    perihelionLongitudeJ2000Deg: C.earthOrbital.earthPerihelionLongitudeJ2000,
-    inclinationCycleAnchorDeg: C.earthOrbital.earthInclinationCycleAnchor,
-  });
+  // moonEcc — the model's ONE eccentricity law — is created above (with
+  // eccentricityAt); the Moon channel's E-factor e(t)/e(J2000) and its
+  // perigee/node T² channel ride it, e(J2000) being the observed value exactly.
 
   // Layer-2 month/precession chain (Brouwer-Clemence m² scaling × the
   // e_E-line modulation)
@@ -860,12 +890,11 @@ export function assembleModel(C, F) {
   // D2 completion table is unchanged (residual 2lE re-fit ≈ 0.1″ — noise).
   const sunL0Deg = C.earthOrbital.sunMeanLongitudeJ2000_deg;
   const sunTropicalRateDegPerCy = 360 * julianCenturyDays / meanSolarYearDays;
-  /** @param {number} year @returns {number} */
-  const sunEccentricityAt = (year) => {
-    const cos3 = -Math.cos(phaseRadians(balancedYear, year, 3));
-    const cos3J2000 = -Math.cos(phaseRadians(balancedYear, 2000, 3));
-    return eccentricityAt(year) + (eccentricityBase / 2) * (cos3 - cos3J2000);
-  };
+  // The eclipse Sun's eccentricity IS the model's one eccentricity law
+  // (eccentricityAt above) — the former sunEccentricityAt (FQ-7-Sun option
+  // C-small: the J2000-anchored H/3 line) was the additive-anchor form of the
+  // same movement and is retired by the unification (measured identical on
+  // every modern gate).
   const sunMeanLongitudeDegAt = (() => {
     // SW PHASE B — the CLOSED FORM on the integrated lattice phase
     // (supersedes the E5 ±3,000-yr trapezoid table; valid at EVERY epoch,
@@ -943,7 +972,7 @@ export function assembleModel(C, F) {
     frameworkSun: {
       sunMeanLongitudeJ2000Deg: sunL0Deg,
       tropicalRateDegPerCy: sunTropicalRateDegPerCy,
-      eccentricityAt: sunEccentricityAt,
+      eccentricityAt,
       perihelionLongitudeDegAt: earthPerihelionDeg,
       meanLongitudeDegAt: sunMeanLongitudeDegAt,
     },
@@ -1087,7 +1116,7 @@ export function assembleModel(C, F) {
       frameworkSunDeps: Object.freeze({
         sunMeanLongitudeJ2000Deg: sunL0Deg,
         tropicalRateDegPerCy: sunTropicalRateDegPerCy,
-        eccentricityAt: sunEccentricityAt,
+        eccentricityAt,
         perihelionLongitudeDegAt: earthPerihelionDeg,
         meanLongitudeDegAt: sunMeanLongitudeDegAt,
       }),

package/src/moon/ecc-channel.cjs

@@ -8,10 +8,14 @@
  * max. The observed J2000 e (−0.86%), ė (+1.7%) and the sign of ë are
  * PREDICTIONS of this line, not inputs.
  *
- * NOT the same quantity as the Sun machinery's eccentricity: Earth's own
- * orbital e for the cardinal-point/EoC path keeps the H/16 perihelion law
- * (`cardinal`'s injected eccentricityAt). This line is the MOON channel's
- * view of the H/3 movement — two different laws, deliberately.
+ * ONE law for the eclipse chain (FQ-7-Sun option C-small, doc 66 §5): the
+ * eclipse Sun's equation of centre rides this same H/3 line (anchored at
+ * its J2000 value, model.js sunEccentricityAt) — eccentricity is
+ * frame-invariant and may carry only fixed-frame lattice periods; the
+ * H/16 perihelion cycle is the H/3 rotation seen from the H/13 equinox
+ * (13 + 3 = 16) and belongs to ϖ, not e (doc 108). Only the
+ * cardinal-point/deep-time path still keeps the H/16 law (`cardinal`'s
+ * injected eccentricityAt) — a recorded, not yet resolved, split.
  *
  * Extracted VERBATIM from src/script.js (_FW_ECC, _fwEarthEcc,
  * _eCompModulation, _fwChannelIntegral, the framework branch of _fwEFactor),
@@ -39,6 +43,7 @@
  * @param {{
  *   cyclesBetween: (yearA: number, yearB: number, divisorN: number) => (number | null),
  *   eccentricityBase: number,
+ *   eccentricityJ2000?: number,
  *   perihelionLongitudeJ2000Deg: number,
  *   inclinationCycleAnchorDeg: number,
  * }} deps — cyclesBetween MUST be the engine's own H-lattice phase counter
@@ -51,10 +56,21 @@ function createMoonEccChannel({
   eccentricityBase,
   perihelionLongitudeJ2000Deg,
   inclinationCycleAnchorDeg,
+  eccentricityJ2000,
 }) {
   const th0 = (perihelionLongitudeJ2000Deg - inclinationCycleAnchorDeg) * Math.PI / 180;
-  const amplitude = eccentricityBase / 2;
-  const mean = eccentricityBase;
+  // ECCENTRICITY UNIFICATION (plan IP-eccentricity-unification, D1 = form
+  // (e)): this line is the model's ONE eccentricity law for every consumer
+  // (eclipse Sun, Sun distance, Moon E-factor and T² channel, cardinal
+  // braid, scene). When the observed J2000 eccentricity is supplied the
+  // mean is DERIVED from it and the shared anchor — base' = e0 / (1 +
+  // cos θ₀ / 2) = 0.015520 — so e(J2000) is exact by construction; the
+  // legacy form (mean = the Law-5 base, e(J2000) a −0.86% prediction) is
+  // kept for callers that do not supply it.
+  const mean = eccentricityJ2000 !== undefined
+    ? eccentricityJ2000 / (1 + Math.cos(th0) / 2)
+    : eccentricityBase;
+  const amplitude = mean / 2;
   const e0 = mean + amplitude * Math.cos(th0);        // J2000 anchor, exact by construction
   const g0 = Math.pow(1 - e0 * e0, -1.5);             // (1−e²)^(−3/2) at the anchor
 
@@ -102,7 +118,21 @@ function createMoonEccChannel({
     return eccAt(tYr) / e0;
   }
 
-  return { eccAt, modulation, channelIntegral, eFactorAt, th0, amplitude, mean, e0, g0 };
+  /** de/dyear of the line at tYr years from J2000 (the cardinal braid's
+   *  equation-of-centre derivative): −amplitude·sin(θ)·dθ/dt, with dθ/dt
+   *  from the phase counter's two-point difference over ±1 yr (exact for a
+   *  linear phase, ppm-class under the integrated one). Returns 0 past the
+   *  tidal-lock asymptote, where eccAt returns the mean.
+   *  @param {number} tYr @returns {number} */
+  function eccRateAt(tYr) {
+    const c = cyclesBetween(2000, 2000 + tYr, 3);
+    if (c === null) return 0;
+    const cp = cyclesBetween(2000, 2000 + tYr + 1, 3), cm = cyclesBetween(2000, 2000 + tYr - 1, 3);
+    const dth = 2 * Math.PI * ((cp === null || cm === null) ? 0 : (cp - cm) / 2);
+    return -amplitude * Math.sin(th0 + 2 * Math.PI * c) * dth;
+  }
+
+  return { eccAt, eccRateAt, modulation, channelIntegral, eFactorAt, th0, amplitude, mean, e0, g0 };
 }
 
 module.exports = { createMoonEccChannel };

package/src/planets/predict.cjs

@@ -87,7 +87,7 @@ function beatPair(t1, t2) {
  * @property {(year: number) => number} computeObliquityEarthDeg
  * @property {(year: number) => number} computeEccentricityEarth
  * @property {number} obliquityMeanDeg - the fitted solstice-obliquity mean
- * @property {number} eccentricityMean - sqrt(base² + amplitude²)
+ * @property {number} eccentricityMean - Earth's mean eccentricity (the one law's base′)
  */
 
 /**