@essrt/physics 7.0.0 → 9.1.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: "7940c564f7b5f0a6";
+export declare const CONSTANTS_HASH: "b9c4d78647577e67";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "7940c564f7b5f0a6";
+  readonly hash: "b9c4d78647577e67";
   readonly additionalBodies: {
     "pluto": {
       "name": string;
@@ -372,6 +372,113 @@ export declare const DEFAULT_CONSTANTS: {
       "rotationPeriodDays": number;
     };
   };
+  readonly planetSpinPhysical: {
+    "mercury": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": unknown[];
+    };
+    "venus": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": unknown[];
+    };
+    "mars": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "jupiter": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "momentOfInertiaFactorAlternative": number;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "saturn": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "momentOfInertiaFactorAlternative": [number, number];
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "uranus": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+      }>;
+    };
+    "neptune": {
+      "j2": number;
+      "j2ReferenceRadiusKm": number;
+      "momentOfInertiaFactor": number;
+      "momentOfInertiaFactorClass": string;
+      "rotationRateDegPerDay": number;
+      "poleRaJ2000Deg": number;
+      "poleDecJ2000Deg": number;
+      "cassiniLocked": boolean;
+      "satellites": Array<{
+        "name": string;
+        "gmKm3S2": number;
+        "semiMajorAxisKm": number;
+        "orbitalPeriodDays": number;
+        "retrograde": boolean;
+      }>;
+    };
+  };
   readonly planets: {
     "mercury": {
       "name": string;
@@ -386,9 +493,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "venus": {
       "name": string;
@@ -403,9 +507,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": null;
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "mars": {
       "name": string;
@@ -420,9 +521,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "jupiter": {
       "name": string;
@@ -437,9 +535,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "saturn": {
       "name": string;
@@ -454,9 +549,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
     };
     "uranus": {
       "name": string;
@@ -471,10 +563,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": [number, number];
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
-      "eccentricityCycleApprox": boolean;
     };
     "neptune": {
       "name": string;
@@ -489,10 +577,6 @@ export declare const DEFAULT_CONSTANTS: {
       "inclinationCycleAnchor": number;
       "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
-      "obliquityCycleFraction": null;
-      "axialPrecessionFraction": [number, number];
-      "eccentricityCycleFraction": [number, number];
-      "eccentricityCycleApprox": boolean;
     };
   };
   readonly timeReference: {

package/src/constants/generated.js

@@ -33,7 +33,7 @@
  * carry it so a counterfactual is reproducible (§2d).
  * @type {string}
  */
-export const CONSTANTS_HASH = "7940c564f7b5f0a6";
+export const CONSTANTS_HASH = "b9c4d78647577e67";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
 export const MODEL_VERSION = "v14.0";
@@ -43,7 +43,7 @@ export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";
 
 /** @type {Readonly<Record<string, unknown>>} */
 export const DEFAULT_CONSTANTS = Object.freeze({
-  hash: "7940c564f7b5f0a6",
+  hash: "b9c4d78647577e67",
   additionalBodies: {
     "pluto": {
       "name": "Pluto",
@@ -171,10 +171,10 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "eros": 16.84
   },
   cardinalPointAnchors: {
-    "SS": 2451716.5748,
-    "WS": 2451900.06782,
-    "VE": 2451623.738137,
-    "AE": 2451810.304796
+    "SS": 2451716.575,
+    "WS": 2451900.0674,
+    "VE": 2451623.816,
+    "AE": 2451810.2271
   },
   deepTime: {
     "alpha1PerMa": -0.000099375895103,
@@ -250,7 +250,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
     "moonStartposApsidal": 347.5476,
     "moonStartposNodal": 64.0435,
     "moonStartposMoon": 67.8443,
-    "moonMeeusLpCorrection": 0.009097354859355077
+    "moonMeeusLpCorrection": 0
   },
   moonMeeus: {
     "elpW1T2Decomposition_arcsecPerCy2": {
@@ -430,6 +430,216 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "rotationPeriodDays": 0.671300001591743
     }
   },
+  planetSpinPhysical: {
+    "mercury": {
+      "j2": 0.0000503,
+      "j2ReferenceRadiusKm": 2440,
+      "momentOfInertiaFactor": 0.346,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": 6.1385108,
+      "poleRaJ2000Deg": 281.0103,
+      "poleDecJ2000Deg": 61.4155,
+      "cassiniLocked": true,
+      "satellites": []
+    },
+    "venus": {
+      "j2": 0.0000044044,
+      "j2ReferenceRadiusKm": 6051,
+      "momentOfInertiaFactor": 0.337,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": -1.4813688,
+      "poleRaJ2000Deg": 272.76,
+      "poleDecJ2000Deg": 67.16,
+      "cassiniLocked": false,
+      "satellites": []
+    },
+    "mars": {
+      "j2": 0.001956608,
+      "j2ReferenceRadiusKm": 3396,
+      "momentOfInertiaFactor": 0.3644,
+      "momentOfInertiaFactorClass": "spin-inferred",
+      "rotationRateDegPerDay": 350.891982443297,
+      "poleRaJ2000Deg": 317.68143,
+      "poleDecJ2000Deg": 52.8865,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Phobos",
+          "gmKm3S2": 0.0007087,
+          "semiMajorAxisKm": 9375,
+          "orbitalPeriodDays": 0.3187
+        },
+        {
+          "name": "Deimos",
+          "gmKm3S2": 0.0000962,
+          "semiMajorAxisKm": 23457,
+          "orbitalPeriodDays": 1.2625
+        }
+      ]
+    },
+    "jupiter": {
+      "j2": 0.014696572,
+      "j2ReferenceRadiusKm": 71492,
+      "momentOfInertiaFactor": 0.26393,
+      "momentOfInertiaFactorClass": "gravity-constrained",
+      "momentOfInertiaFactorAlternative": 0.2756,
+      "rotationRateDegPerDay": 870.536,
+      "poleRaJ2000Deg": 268.056595,
+      "poleDecJ2000Deg": 64.495303,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Io",
+          "gmKm3S2": 5959.91547,
+          "semiMajorAxisKm": 421800,
+          "orbitalPeriodDays": 1.762732
+        },
+        {
+          "name": "Europa",
+          "gmKm3S2": 3202.7121,
+          "semiMajorAxisKm": 671100,
+          "orbitalPeriodDays": 3.525463
+        },
+        {
+          "name": "Ganymede",
+          "gmKm3S2": 9887.83275,
+          "semiMajorAxisKm": 1070400,
+          "orbitalPeriodDays": 7.155588
+        },
+        {
+          "name": "Callisto",
+          "gmKm3S2": 7179.2834,
+          "semiMajorAxisKm": 1882700,
+          "orbitalPeriodDays": 16.69044
+        }
+      ]
+    },
+    "saturn": {
+      "j2": 0.016290573,
+      "j2ReferenceRadiusKm": 60330,
+      "momentOfInertiaFactor": 0.2181,
+      "momentOfInertiaFactorClass": "gravity-constrained",
+      "momentOfInertiaFactorAlternative": [
+        0.2257,
+        0.2438
+      ],
+      "rotationRateDegPerDay": 810.7939024,
+      "poleRaJ2000Deg": 40.589,
+      "poleDecJ2000Deg": 83.537,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Mimas",
+          "gmKm3S2": 2.50349,
+          "semiMajorAxisKm": 186000,
+          "orbitalPeriodDays": 0.942422
+        },
+        {
+          "name": "Enceladus",
+          "gmKm3S2": 7.21037,
+          "semiMajorAxisKm": 238400,
+          "orbitalPeriodDays": 1.370218
+        },
+        {
+          "name": "Tethys",
+          "gmKm3S2": 41.21353,
+          "semiMajorAxisKm": 295000,
+          "orbitalPeriodDays": 1.887802
+        },
+        {
+          "name": "Dione",
+          "gmKm3S2": 73.11607,
+          "semiMajorAxisKm": 377700,
+          "orbitalPeriodDays": 2.736916
+        },
+        {
+          "name": "Rhea",
+          "gmKm3S2": 153.94175,
+          "semiMajorAxisKm": 527200,
+          "orbitalPeriodDays": 4.517503
+        },
+        {
+          "name": "Titan",
+          "gmKm3S2": 8978.1371,
+          "semiMajorAxisKm": 1221900,
+          "orbitalPeriodDays": 15.945448
+        },
+        {
+          "name": "Iapetus",
+          "gmKm3S2": 120.51511,
+          "semiMajorAxisKm": 3561700,
+          "orbitalPeriodDays": 79.331002
+        }
+      ]
+    },
+    "uranus": {
+      "j2": 0.00351068,
+      "j2ReferenceRadiusKm": 25559,
+      "momentOfInertiaFactor": 0.225,
+      "momentOfInertiaFactorClass": "interior-model",
+      "rotationRateDegPerDay": -501.1600928,
+      "poleRaJ2000Deg": 257.311,
+      "poleDecJ2000Deg": -15.175,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Miranda",
+          "gmKm3S2": 4.3,
+          "semiMajorAxisKm": 129846,
+          "orbitalPeriodDays": 1.413479
+        },
+        {
+          "name": "Ariel",
+          "gmKm3S2": 83.5,
+          "semiMajorAxisKm": 190929,
+          "orbitalPeriodDays": 2.520379
+        },
+        {
+          "name": "Umbriel",
+          "gmKm3S2": 85.1,
+          "semiMajorAxisKm": 265986,
+          "orbitalPeriodDays": 4.144177
+        },
+        {
+          "name": "Titania",
+          "gmKm3S2": 226.9,
+          "semiMajorAxisKm": 436298,
+          "orbitalPeriodDays": 8.705869
+        },
+        {
+          "name": "Oberon",
+          "gmKm3S2": 205.3,
+          "semiMajorAxisKm": 583511,
+          "orbitalPeriodDays": 13.463237
+        }
+      ]
+    },
+    "neptune": {
+      "j2": 0.00340843,
+      "j2ReferenceRadiusKm": 25225,
+      "momentOfInertiaFactor": 0.23,
+      "momentOfInertiaFactorClass": "interior-model",
+      "rotationRateDegPerDay": 541.1397757,
+      "poleRaJ2000Deg": 299.36,
+      "poleDecJ2000Deg": 43.46,
+      "cassiniLocked": false,
+      "satellites": [
+        {
+          "name": "Triton",
+          "gmKm3S2": 1428.49546,
+          "semiMajorAxisKm": 354800,
+          "orbitalPeriodDays": 5.876994,
+          "retrograde": true
+        },
+        {
+          "name": "Proteus",
+          "gmKm3S2": 2.58342,
+          "semiMajorAxisKm": 117600,
+          "orbitalPeriodDays": 1.122315
+        }
+      ]
+    }
+  },
   planets: {
     "mercury": {
       "name": "Mercury",
@@ -446,19 +656,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 32.83,
       "inclinationCycleAnchor": 234.52,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 9,
-      "obliquityCycleFraction": [
-        8,
-        3
-      ],
-      "axialPrecessionFraction": [
-        -8,
-        9
-      ],
-      "eccentricityCycleFraction": [
-        2,
-        21
-      ]
+      "ascendingNodeCyclesIn8H": 9
     },
     "venus": {
       "name": "Venus",
@@ -475,16 +673,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 54.7,
       "inclinationCycleAnchor": 218.64,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 1,
-      "obliquityCycleFraction": null,
-      "axialPrecessionFraction": [
-        8,
-        91
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        19
-      ]
+      "ascendingNodeCyclesIn8H": 1
     },
     "mars": {
       "name": "Mars",
@@ -501,19 +690,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 354.87,
       "inclinationCycleAnchor": 236.07,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 64,
-      "obliquityCycleFraction": [
-        8,
-        21
-      ],
-      "axialPrecessionFraction": [
-        -1,
-        2
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        52
-      ]
+      "ascendingNodeCyclesIn8H": 64
     },
     "jupiter": {
       "name": "Jupiter",
@@ -530,19 +707,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 312.89,
       "inclinationCycleAnchor": 287.06,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 36,
-      "obliquityCycleFraction": [
-        1,
-        2
-      ],
-      "axialPrecessionFraction": [
-        -8,
-        21
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        44
-      ]
+      "ascendingNodeCyclesIn8H": 36
     },
     "saturn": {
       "name": "Saturn",
@@ -559,19 +724,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 118.81,
       "inclinationCycleAnchor": 116.26,
       "antiPhase": true,
-      "ascendingNodeCyclesIn8H": 36,
-      "obliquityCycleFraction": [
-        1,
-        3
-      ],
-      "axialPrecessionFraction": [
-        -4,
-        3
-      ],
-      "eccentricityCycleFraction": [
-        8,
-        163
-      ]
+      "ascendingNodeCyclesIn8H": 36
     },
     "uranus": {
       "name": "Uranus",
@@ -588,20 +741,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 307.8,
       "inclinationCycleAnchor": 21.33,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 11,
-      "obliquityCycleFraction": [
-        1,
-        2
-      ],
-      "axialPrecessionFraction": [
-        610,
-        1
-      ],
-      "eccentricityCycleFraction": [
-        1,
-        10
-      ],
-      "eccentricityCycleApprox": true
+      "ascendingNodeCyclesIn8H": 11
     },
     "neptune": {
       "name": "Neptune",
@@ -618,17 +758,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "ascendingNodeInvPlane": 192.04,
       "inclinationCycleAnchor": 174.04,
       "antiPhase": false,
-      "ascendingNodeCyclesIn8H": 3,
-      "obliquityCycleFraction": null,
-      "axialPrecessionFraction": [
-        -68,
-        1
-      ],
-      "eccentricityCycleFraction": [
-        2,
-        25
-      ],
-      "eccentricityCycleApprox": true
+      "ascendingNodeCyclesIn8H": 3
     }
   },
   timeReference: {
@@ -766,8 +896,8 @@ export const REFERENCE_DATA = Object.freeze({
     "sunModelTrueErrorDeg": 0.003,
     "sunTropicalYearDiffSeconds": 0.1,
     "sunSiderealYearDiffSeconds": 0.02,
-    "frameworkSunVsJplRmsArcsec": 0.8,
-    "meeusCh25SunVsJplRmsArcsec": 1.22,
+    "frameworkSunVsJplRmsArcsec": 1.03,
+    "meeusCh25SunVsJplRmsArcsec": 1.02,
     "moonSeriesLonVsJplRmsArcsec": 2.84,
     "moonSeriesLatVsJplRmsArcsec": 0.35,
     "mercuryNewtonianArcsecCy": 532,

package/src/earth/deep-orbital-history.cjs

@@ -320,6 +320,14 @@ function createDeepOrbitalHistory({
       eSinPeri: e * Math.sin(periOfDateDeg * D2R),
       eCosPeri: e * Math.cos(periOfDateDeg * D2R),
       inclEclDeg: i * R2D,
+      // The orbit normal's in-plane components in the J2000 ecliptic frame
+      // (n̂ = (sin i sin Ω, −sin i cos Ω, cos i); the third component follows).
+      // With epsDeg and equinoxLonJ2000Deg these reconstruct the FULL Earth
+      // frame of date — n̂, the equinox ĝ and the spin axis ŝ — for the scene
+      // (earth/frame-of-date.cjs). Smooth and small (47″/cy): linear
+      // interpolation on the grid is exact to the ζ series' own resolution.
+      orbitNormalX: n[0],
+      orbitNormalY: n[1],
       // The equinox node (ŝ×n̂) longitude in the J2000 ecliptic frame —
       // its year-over-year retrograde advance IS the general precession of
       // date, wobble included (the n̂(t) geometry generates the equinox
@@ -420,6 +428,8 @@ function createDeepOrbitalHistory({
           eSinPeri: lerp(a.eSinPeri, b.eSinPeri),
           eCosPeri: lerp(a.eCosPeri, b.eCosPeri),
           inclEclDeg: lerp(a.inclEclDeg, b.inclEclDeg),
+          orbitNormalX: lerp(a.orbitNormalX, b.orbitNormalX),
+          orbitNormalY: lerp(a.orbitNormalY, b.orbitNormalY),
           // C1 (cubic Hermite) — linear interpolation here made the
           // realized equinox RATE piecewise-constant per grid cell (the
           // measured per-year precession staircase). Node rates come from