@essrt/physics 21.0.0 → 22.0.0

package/package.json

@@ -1,9 +1,9 @@
 {
   "name": "@essrt/physics",
-  "version": "21.0.0",
+  "version": "22.0.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core — 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": "v24.0",
+    "modelVersion": "v25.0",
     "_note": "The model identity this package version ships (IP-unified-architecture §10: package semver tracks the API axis; the model version + CONSTANTS_HASH/COEFFICIENTS_HASH in src/constants/ carry the scientific identity). The version-pinning gate enforces this pairing against public/input/model-version.json."
   },
   "type": "module",
@@ -54,6 +54,7 @@
     "./planets/inv-plane-frame": { "types": "./types/planets/inv-plane-frame.d.cts", "default": "./src/planets/inv-plane-frame.cjs" },
     "./planets/secular-series": { "types": "./types/planets/secular-series.d.cts", "default": "./src/planets/secular-series.cjs" },
     "./planets/model": { "types": "./types/planets/model.d.cts", "default": "./src/planets/model.cjs" },
+    "./planets/display-names": { "types": "./types/planets/display-names.d.cts", "default": "./src/planets/display-names.cjs" },
     "./deltat/cycles": { "types": "./types/deltat/cycles.d.cts", "default": "./src/deltat/cycles.cjs" },
     "./deltat/deep-time": { "types": "./types/deltat/deep-time.d.cts", "default": "./src/deltat/deep-time.cjs" },
     "./deltat/historical": { "types": "./types/deltat/historical.d.cts", "default": "./src/deltat/historical.cjs" },

package/src/climate/l1-orbital.cjs

@@ -88,7 +88,9 @@ function laggedL1Terms(l1Terms, relaxationKyr) {
  *   j2RateJ2000PerYr?: number,
  *   j2J2000?: number,
  * }} cfg - relaxationKyr, alphaGiaRateJ2000PerYr and j2RateJ2000PerYr are
- *   model-parameters deepTime constants; alphaJ2000 the IERS
+ *   assembled deepTime block (DEFAULT_CONSTANTS.deepTime: τ the banked record
+ *   data/deep-time-fitted.json, the rate derived from astro-reference
+ *   giaCoxChaoPeltier — plan 11 step 1); alphaJ2000 the IERS
  *   moment-of-inertia factor; j2J2000 Earth's J₂ (astro-reference)
  * @returns {{
  *   alphaAt: (year: number) => number,

package/src/constants/coefficients.d.ts

@@ -1,6 +1,6 @@
 // GENERATED — do not edit. Regenerate: node tools/constants/generate.mjs --write
 // The fitted coefficients' shapes for TypeScript consumers (§2g); values live
-// in coefficients.js, emitted VERBATIM from fitted-coefficients.json.
+// in coefficients.js, emitted VERBATIM from the fit artifacts under data/.
 
 export declare const COEFFICIENTS_HASH: "91e2d041c7b6e038";
 

package/src/constants/coefficients.js

@@ -2,7 +2,8 @@
  * GENERATED — do not edit. Regenerate:
  *   node tools/constants/generate.mjs --write
  *
- * Source: public/input/fitted-coefficients.json, emitted VERBATIM.
+ * Source: the shipped fit artifacts under data/ (EXTRA_COEFFICIENT_SOURCES in
+ * the generator), each emitted VERBATIM.
  *
  * Separate from DEFAULT_CONSTANTS because the lifecycle differs (§2j): these are
  * fitting-pipeline output, regenerated wholesale when a step runs, never

package/src/constants/generated.d.ts

@@ -2,56 +2,14 @@
 // Gives the TypeScript website full type safety at the boundary (§2g) while
 // packages/physics stays JavaScript.
 
-export declare const CONSTANTS_HASH: "c161ec19e5984b9e";
+export declare const CONSTANTS_HASH: "f4f43386df2781b8";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "c161ec19e5984b9e";
-  readonly additionalBodies: {
-    "pluto": {
-      "name": string;
-      "orbitalEccentricityBase": number;
-      "type": string;
-      "angleCorrection": number;
-      "startpos": number;
-      "ascendingNodeInvPlane": number;
-      "invPlaneInclinationMean": number;
-      "invPlaneInclinationAmplitude": number;
-    };
-    "halleys": {
-      "name": string;
-      "orbitalEccentricityBase": number;
-      "type": string;
-      "angleCorrection": number;
-      "startpos": number;
-      "ascendingNodeInvPlane": number;
-      "invPlaneInclinationMean": number;
-      "invPlaneInclinationAmplitude": number;
-    };
-    "eros": {
-      "name": string;
-      "orbitalEccentricityBase": number;
-      "type": string;
-      "angleCorrection": number;
-      "startpos": number;
-      "ascendingNodeInvPlane": number;
-      "invPlaneInclinationMean": number;
-      "invPlaneInclinationAmplitude": number;
-    };
-    "ceres": {
-      "name": string;
-      "orbitalEccentricityBase": number;
-      "orbitDistanceOverride": number;
-      "angleCorrection": number;
-      "startpos": number;
-      "ascendingNodeInvPlane": number;
-      "invPlaneInclinationMean": number;
-      "invPlaneInclinationAmplitude": number;
-    };
-  };
+  readonly hash: "f4f43386df2781b8";
   readonly additionalBodiesReference: {
     "pluto": {
       "solarYearInput": number;
@@ -64,6 +22,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "halleys": {
       "solarYearInput": number;
@@ -143,21 +102,6 @@ export declare const DEFAULT_CONSTANTS: {
       "thermalPumpFactor": number;
     };
   };
-  readonly earth: {
-    "earthtiltMean": number;
-    "earthInvPlaneInclinationAmplitude": number;
-    "eccentricityBase": number;
-    "eccentricityAmplitude": number;
-    "kDeviceInclinationPhaseJ2000Deg": number;
-    "kDeviceWheelPeriodYears": {
-      "3": number;
-      "5": number;
-      "8": number;
-      "13": number;
-      "16": number;
-    };
-    "kDeviceWheelPhaseAtJ2000Cycles": number;
-  };
   readonly earthOrbital: {
     "j2000EpochYear": number;
     "perihelionalignmentYear": number;
@@ -177,19 +121,9 @@ export declare const DEFAULT_CONSTANTS: {
     "deltaTStart": number;
     "sunTilt": number;
   };
-  readonly foundational: {
-    "startmodelJD": number;
-    "startmodelYear": number;
-    "correctionDays": number;
-    "correctionSun": number;
-    "startAngleModel": number;
-    "useVariableSpeed": boolean;
-  };
-  readonly moon: {
-    "moonStartposApsidal": number;
-    "moonStartposNodal": number;
-    "moonStartposMoon": number;
-    "moonMeeusLpCorrection": number;
+  readonly giaCoxChaoPeltier: {
+    "dJ2DtPerYr": number;
+    "j2ToAlphaFactor": number;
   };
   readonly moonMeeus: {
     "elpW1T2Decomposition_arcsecPerCy2": {
@@ -216,21 +150,13 @@ export declare const DEFAULT_CONSTANTS: {
     "moonApsidalPrecessionDaysInputICRF": number;
     "moonNodalPrecessionDaysInputICRF": number;
     "moonDistance": number;
+    "moonRecessionAlpha1PerMa": number;
     "moonEclipticInclinationJ2000": number;
     "moonInclinationConstantBrownELP": number;
     "moonOrbitalEccentricityBase": number;
     "moonObliquityEclipticJ2000": number;
     "moonTilt": number;
   };
-  readonly perihelionPassageRef: {
-    "mercury": number;
-    "venus": number;
-    "mars": number;
-    "jupiter": number;
-    "saturn": number;
-    "uranus": number;
-    "neptune": number;
-  };
   readonly physicalConstants: {
     "currentAUDistance": number;
     "G_CONSTANT": number;
@@ -295,6 +221,7 @@ export declare const DEFAULT_CONSTANTS: {
       "invPlaneInclinationJ2000": number;
       "meanAnomaly": number;
       "trueAnomaly": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "venus": {
       "solarYearInput": number;
@@ -307,6 +234,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "mars": {
       "solarYearInput": number;
@@ -319,6 +247,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "jupiter": {
       "solarYearInput": number;
@@ -331,6 +260,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "saturn": {
       "solarYearInput": number;
@@ -343,6 +273,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "uranus": {
       "solarYearInput": number;
@@ -355,6 +286,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "neptune": {
       "solarYearInput": number;
@@ -367,6 +299,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
   };
   readonly planetSpinPhysical: {
@@ -476,72 +409,9 @@ export declare const DEFAULT_CONSTANTS: {
       }>;
     };
   };
-  readonly planets: {
-    "mercury": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "venus": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "mars": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "jupiter": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "saturn": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "uranus": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-    "neptune": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-    };
-  };
   readonly timeReference: {
+    "sceneEpochJD": number;
+    "sceneEpochYear": number;
     "j2000JD": number;
     "julianCenturyDays": number;
     "gregorianStartJD": number;
@@ -598,10 +468,6 @@ export declare const REFERENCE_DATA: {
     "greatattractorDistance": number;
     "milkywaySpeed": number;
   };
-  readonly giaCoxChaoPeltier: {
-    "dJ2DtPerYr": number;
-    "j2ToAlphaFactor": number;
-  };
   readonly jplEclipticInclinationTrends: {
     "mercury": number;
     "venus": number;

package/src/constants/generated.js

@@ -2,7 +2,11 @@
  * GENERATED — do not edit. Regenerate:
  *   node tools/constants/generate.mjs --write
  *
- * Source of truth: public/input/{model-parameters,astro-reference}.json
+ * Source of truth: public/input/astro-reference.json (model-parameters.json is
+ * DELETED, plan 11 step 6 — its leaves have their homes); the deepTime block
+ * is ASSEMBLED from astro-reference.json (moonReference.moonRecessionAlpha1PerMa,
+ * giaCoxChaoPeltier), data/deep-time-fitted.json and data/deltaT-4flag-fit.json
+ * shipped_taper (plan 11 step 1).
  * Classification and rationale: tools/constants/generate.mjs
  *
  * Two exports, deliberately separate (§2d):
@@ -12,7 +16,6 @@
  *   REFERENCE_DATA     validation targets + presentation. Single-sourced so
  *                      nothing duplicates them, but NOT injectable:
  *   moonGrailWilliams2014            target
- *   giaCoxChaoPeltier                target
  *   obliquityChapront2002            target
  *   planetSpinObserved               target
  *   ascendingNodesSouamiSouchay      target
@@ -33,59 +36,17 @@
  * carry it so a counterfactual is reproducible (§2d).
  * @type {string}
  */
-export const CONSTANTS_HASH = "c161ec19e5984b9e";
+export const CONSTANTS_HASH = "f4f43386df2781b8";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
-export const MODEL_VERSION = "v24.0";
+export const MODEL_VERSION = "v25.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: "c161ec19e5984b9e",
-  additionalBodies: {
-    "pluto": {
-      "name": "Pluto",
-      "orbitalEccentricityBase": 0.2488273,
-      "type": "I",
-      "angleCorrection": 2.469281,
-      "startpos": 71.555,
-      "ascendingNodeInvPlane": 101.06,
-      "invPlaneInclinationMean": 15.7162,
-      "invPlaneInclinationAmplitude": 0.717024
-    },
-    "halleys": {
-      "name": "Halley's",
-      "orbitalEccentricityBase": 0.96714291,
-      "type": "III",
-      "angleCorrection": -1.619816,
-      "startpos": 80,
-      "ascendingNodeInvPlane": 59.56,
-      "invPlaneInclinationMean": 150,
-      "invPlaneInclinationAmplitude": 0.1
-    },
-    "eros": {
-      "name": "Eros",
-      "orbitalEccentricityBase": 0.2229512,
-      "type": "II",
-      "angleCorrection": 0.047888,
-      "startpos": 57.402,
-      "ascendingNodeInvPlane": 10.36,
-      "invPlaneInclinationMean": 9.25,
-      "invPlaneInclinationAmplitude": 0.5
-    },
-    "ceres": {
-      "name": "Ceres",
-      "orbitalEccentricityBase": 0.0755347,
-      "orbitDistanceOverride": 2.76596,
-      "angleCorrection": 0,
-      "startpos": 0,
-      "ascendingNodeInvPlane": 80.89,
-      "invPlaneInclinationMean": 0.43,
-      "invPlaneInclinationAmplitude": 0.05
-    }
-  },
+  hash: "f4f43386df2781b8",
   additionalBodiesReference: {
     "pluto": {
       "solarYearInput": 90465,
@@ -97,7 +58,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNode": 110.30393,
       "meanAnomaly": 15.55009,
       "trueAnomaly": 26.31965048,
-      "rotationPeriodDays": 6.38720012152536
+      "rotationPeriodDays": 6.38720012152536,
+      "ascendingNodeInvPlaneClosedForm": 101.06
     },
     "halleys": {
       "solarYearInput": 27503,
@@ -189,21 +151,6 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "thermalPumpFactor": 0.454
     }
   },
-  earth: {
-    "earthtiltMean": 23.413526608469965,
-    "earthInvPlaneInclinationAmplitude": 0.636069722432876,
-    "eccentricityBase": 0.015386009686374918,
-    "eccentricityAmplitude": 0.0013559440307290061,
-    "kDeviceInclinationPhaseJ2000Deg": 81.17719000000001,
-    "kDeviceWheelPeriodYears": {
-      "3": 111772.33333333333,
-      "5": 67063.4,
-      "8": 41914.625,
-      "13": 25793.615384615383,
-      "16": 20957.3125
-    },
-    "kDeviceWheelPhaseAtJ2000Cycles": 0.908498525276082
-  },
   earthOrbital: {
     "j2000EpochYear": 2000,
     "perihelionalignmentYear": 1246.03125,
@@ -223,19 +170,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "deltaTStart": 56.07959443353602,
     "sunTilt": 7.155
   },
-  foundational: {
-    "startmodelJD": 2451716.5,
-    "startmodelYear": 2000.5,
-    "correctionDays": -0.828832119703292,
-    "correctionSun": 0.4971540149655503,
-    "startAngleModel": 89.91949879,
-    "useVariableSpeed": true
-  },
-  moon: {
-    "moonStartposApsidal": 347.5476,
-    "moonStartposNodal": 64.0435,
-    "moonStartposMoon": 67.8443,
-    "moonMeeusLpCorrection": 0
+  giaCoxChaoPeltier: {
+    "dJ2DtPerYr": -2.7e-11,
+    "j2ToAlphaFactor": 2
   },
   moonMeeus: {
     "elpW1T2Decomposition_arcsecPerCy2": {
@@ -262,21 +199,13 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "moonApsidalPrecessionDaysInputICRF": 3231.493,
     "moonNodalPrecessionDaysInputICRF": 6798.38,
     "moonDistance": 384399.07,
+    "moonRecessionAlpha1PerMa": -0.000099375895103,
     "moonEclipticInclinationJ2000": 5.1573,
     "moonInclinationConstantBrownELP": 5.1453964,
     "moonOrbitalEccentricityBase": 0.054900489,
     "moonObliquityEclipticJ2000": 1.5424,
     "moonTilt": 6.687
   },
-  perihelionPassageRef: {
-    "mercury": 2460335.9,
-    "venus": 2455464.42,
-    "mars": 2456499.441,
-    "jupiter": 2464224.5,
-    "saturn": 2452875.9,
-    "uranus": 2439699.8,
-    "neptune": 2409432.4
-  },
   physicalConstants: {
     "currentAUDistance": 149597870.698828,
     "G_CONSTANT": 6.6743e-20,
@@ -340,7 +269,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNode": 48.33033155,
       "invPlaneInclinationJ2000": 6.3472858,
       "meanAnomaly": 156.6364301,
-      "trueAnomaly": 164.1669319
+      "trueAnomaly": 164.1669319,
+      "ascendingNodeInvPlaneClosedForm": 32.83
     },
     "venus": {
       "solarYearInput": 224.695,
@@ -352,7 +282,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 2.1545441,
       "meanAnomaly": 324.9668371,
       "trueAnomaly": 324.5198504,
-      "rotationPeriodDays": 243.022699230302
+      "rotationPeriodDays": 243.022699230302,
+      "ascendingNodeInvPlaneClosedForm": 54.7
     },
     "mars": {
       "solarYearInput": 686.93,
@@ -364,7 +295,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 1.6311858,
       "meanAnomaly": 109.2630844,
       "trueAnomaly": 118.9501056,
-      "rotationPeriodDays": 1.02595659586635
+      "rotationPeriodDays": 1.02595659586635,
+      "ascendingNodeInvPlaneClosedForm": 354.87
     },
     "jupiter": {
       "solarYearInput": 4330.53,
@@ -376,7 +308,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 0.3219652,
       "meanAnomaly": 32.47179744,
       "trueAnomaly": 35.69428061,
-      "rotationPeriodDays": 0.413541666975253
+      "rotationPeriodDays": 0.413541666975253,
+      "ascendingNodeInvPlaneClosedForm": 312.89
     },
     "saturn": {
       "solarYearInput": 10747,
@@ -388,7 +321,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 0.9254704,
       "meanAnomaly": 325.663876,
       "trueAnomaly": 321.7910116,
-      "rotationPeriodDays": 0.440023148755863
+      "rotationPeriodDays": 0.440023148755863,
+      "ascendingNodeInvPlaneClosedForm": 118.81
     },
     "uranus": {
       "solarYearInput": 30586,
@@ -400,7 +334,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 0.9946692,
       "meanAnomaly": 145.7292678,
       "trueAnomaly": 148.5142459,
-      "rotationPeriodDays": 0.718329998141018
+      "rotationPeriodDays": 0.718329998141018,
+      "ascendingNodeInvPlaneClosedForm": 307.8
     },
     "neptune": {
       "solarYearInput": 59800,
@@ -412,7 +347,8 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "invPlaneInclinationJ2000": 0.7354155,
       "meanAnomaly": 262.5003424,
       "trueAnomaly": 261.2242728,
-      "rotationPeriodDays": 0.671300001591743
+      "rotationPeriodDays": 0.671300001591743,
+      "ascendingNodeInvPlaneClosedForm": 192.04
     }
   },
   planetSpinPhysical: {
@@ -625,72 +561,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       ]
     }
   },
-  planets: {
-    "mercury": {
-      "name": "Mercury",
-      "eocFraction": -0.527,
-      "startpos": 83.64651142941864,
-      "angleCorrection": 0.9727205912043573,
-      "type": "I",
-      "mirrorPair": "uranus",
-      "ascendingNodeInvPlane": 32.83
-    },
-    "venus": {
-      "name": "Venus",
-      "eocFraction": 0.436,
-      "startpos": 249.28769121856186,
-      "angleCorrection": -2.789778861448692,
-      "type": "I",
-      "mirrorPair": "neptune",
-      "ascendingNodeInvPlane": 54.7
-    },
-    "mars": {
-      "name": "Mars",
-      "eocFraction": -0.066224,
-      "startpos": 121.4634461571797,
-      "angleCorrection": -2.104652936098656,
-      "type": "II",
-      "mirrorPair": "jupiter",
-      "ascendingNodeInvPlane": 354.87
-    },
-    "jupiter": {
-      "name": "Jupiter",
-      "eocFraction": 0.495,
-      "startpos": 13.887946755371052,
-      "angleCorrection": 0.9333205235053406,
-      "type": "III",
-      "mirrorPair": "mars",
-      "ascendingNodeInvPlane": 312.89
-    },
-    "saturn": {
-      "name": "Saturn",
-      "eocFraction": 0.54,
-      "startpos": 11.279308392461772,
-      "angleCorrection": -0.1784695052193328,
-      "type": "III",
-      "mirrorPair": "earth",
-      "ascendingNodeInvPlane": 118.81
-    },
-    "uranus": {
-      "name": "Uranus",
-      "eocFraction": 0.53,
-      "startpos": 44.900388945775106,
-      "angleCorrection": -0.7343433510960402,
-      "type": "III",
-      "mirrorPair": "mercury",
-      "ascendingNodeInvPlane": 307.8
-    },
-    "neptune": {
-      "name": "Neptune",
-      "eocFraction": 0.585,
-      "startpos": 47.9552024382285,
-      "angleCorrection": 2.3327665835675617,
-      "type": "III",
-      "mirrorPair": "venus",
-      "ascendingNodeInvPlane": 192.04
-    }
-  },
   timeReference: {
+    "sceneEpochJD": 2451716.5,
+    "sceneEpochYear": 2000.5,
     "j2000JD": 2451545,
     "julianCenturyDays": 36525,
     "gregorianStartJD": 2299160.5,
@@ -760,10 +633,6 @@ export const REFERENCE_DATA = Object.freeze({
     "greatattractorDistance": 200000000,
     "milkywaySpeed": 2160000
   },
-  giaCoxChaoPeltier: {
-    "dJ2DtPerYr": -2.7e-11,
-    "j2ToAlphaFactor": 2
-  },
   jplEclipticInclinationTrends: {
     "mercury": -0.00595,
     "venus": -0.00079,

package/src/earth/frame-of-date.cjs

@@ -26,11 +26,35 @@
  * the K sun plane and the RA frame's ecliptic parted by 20.5″ at J2000 and 10′
  * at −3000, the rendered Sun −19″ in declination against Horizons at the
  * March 2000 equinox (+19″ in September; 0.1–1″ at the solstice).
+ *
+ * THE SCENE'S WORLD FRAME (plan 11 step 4) is the J2000 ecliptic itself:
+ * world Y is the J2000 ecliptic pole and the J2000 equinox lies in the world
+ * XZ plane at the view azimuth `SCENE_EQUINOX_VIEW_AZIMUTH_DEG` (atan2(x, z)).
+ * `computeSceneWorldFrameJ2000` is the bridge R (ecliptic J2000 → world) both
+ * twins place every frame of date through. The azimuth is the ONE scene
+ * literal the retired K-device scaffold leaves behind — the pose the scaffold
+ * encoded at J2000, measured on the Node twin (its 0.026° pole tilt dropped),
+ * kept for visual continuity; no observable rides it (every physical reading
+ * is taken in the placed Earth frame, which is R-invariant by construction).
  */
 'use strict';
 
 const D2R = Math.PI / 180;
 
+/** The J2000 equinox's azimuth in the scene's world XZ plane, degrees (atan2(x, z)) — the view convention. */
+const SCENE_EQUINOX_VIEW_AZIMUTH_DEG = 158.2265109779131;
+
+/**
+ * The scene's world frame: R maps a J2000-ecliptic vector to world, columns
+ * [ĝ, n̂ × ĝ, n̂] = [the equinox, longitude 90°, the ecliptic pole] with
+ * n̂ = world Y exactly and ĝ at the view azimuth.
+ * @returns {number[][]} a 3×3 rotation matrix, rows × columns
+ */
+function computeSceneWorldFrameJ2000() {
+  const a = SCENE_EQUINOX_VIEW_AZIMUTH_DEG * D2R, s = Math.sin(a), c = Math.cos(a);
+  return [[s, c, 0], [0, 0, 1], [c, -s, 0]];
+}
+
 /**
  * @typedef {{x: number, y: number, z: number}} Vec3
  */
@@ -102,4 +126,4 @@ function solveWheelAngleForLongitude(lambdaRad, e, periRad) {
   return Math.atan2(Math.sin(th), Math.cos(th));
 }
 
-module.exports = { computeEarthFrameOfDate, solveWheelAngleForLongitude };
+module.exports = { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG };