@essrt/physics 20.0.1 → 22.0.0

package/src/planets/spin-channel.cjs

@@ -46,6 +46,16 @@
  * `computeObliquityJ2000Deg`), and — through the chain's own g-mode beat
  * (keplerian-chain `computeSecularShape`) — the "wobble period" the K
  * device beat from its integer axial and obliquity fractions.
+ *
+ * Since plan 11 step 3 the channel also serves the POLE of date itself —
+ * `spinAxisEclipticJ2000AtYear` (the integrated ŝ(t), angular-momentum
+ * sense) beside `orbitNormalEclipticJ2000AtYear` (n̂(t)) — so both scene
+ * twins place each planet's rotation axis in world through the frame bridge
+ * exactly as Earth's axis is placed (R4): the pole's azimuth used to ride the
+ * retired scaffold's start angle (one turn of the pole per orbit), only its
+ * tilt magnitude was read from here. Mercury's Cassini-locked pole is the
+ * J2000 pole carried with its orbit normal (the minimal rotation n̂₀ → n̂(t)),
+ * the lock's constant obliquity.
  */
 
 'use strict';
@@ -124,6 +134,19 @@ const dot = (a, b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
 const unit = (a) => { const r = Math.hypot(a[0], a[1], a[2]); return [a[0] / r, a[1] / r, a[2] / r]; };
 /** @param {number[]} a @param {number[]} b */
 const angleDeg = (a, b) => Math.acos(Math.max(-1, Math.min(1, dot(a, b)))) * R2D;
+/**
+ * v carried by the minimal rotation that takes unit a onto unit b (Rodrigues
+ * about a × b); v itself when a and b coincide.
+ * @param {number[]} a @param {number[]} b @param {number[]} v
+ */
+const rotateAligned = (a, b, v) => {
+  const k = cross(a, b);
+  const s = Math.hypot(k[0], k[1], k[2]), c = dot(a, b);
+  if (s < 1e-15) return v.slice();
+  const kh = [k[0] / s, k[1] / s, k[2] / s];
+  const kv = cross(kh, v), kd = dot(kh, v) * (1 - c);
+  return [v[0] * c + kv[0] * s + kh[0] * kd, v[1] * c + kv[1] * s + kh[1] * kd, v[2] * c + kv[2] * s + kh[2] * kd];
+};
 
 /**
  * The IAU pole (ICRF equatorial RA/Dec) as an ecliptic-J2000 unit vector.
@@ -272,9 +295,9 @@ function createPlanetSpinChannel(deps) {
   const stepYr = deps.stepYr ?? (cassiniLocked ? STEP_MAX_YR
     : Math.min(STEP_MAX_YR, Math.max(STEP_MIN_YR, axialPrecessionPeriodYearsJ2000 / STEPS_PER_PRECESSION_PERIOD)));
   // append-only trajectories from J2000, one per direction: sample i sits at t = sign · i · step
-  /** @typedef {{ eps: number[], last: number[] }} Track */
+  /** @typedef {{ eps: number[], sv: number[][], last: number[] }} Track — the obliquity and the unit spin vector per sample */
   /** @returns {Track} */
-  const mkTrack = () => ({ eps: [obliquityJ2000], last: s0.slice() });
+  const mkTrack = () => ({ eps: [obliquityJ2000], sv: [s0.slice()], last: s0.slice() });
   const fwd = mkTrack(), bwd = mkTrack();
   /** @param {Track} tr @param {number} sign @param {number} count */
   const grow = (tr, sign, count) => {
@@ -284,6 +307,7 @@ function createPlanetSpinChannel(deps) {
       const h = sign * stepYr;
       const s = rk4(tr.last, t, h);
       tr.eps.push(angleDeg(s, orbitNormalAt(t + h)));
+      tr.sv.push(s);
       tr.last = s;
     }
   };
@@ -333,6 +357,24 @@ function createPlanetSpinChannel(deps) {
     }
     return { minDeg: mn, meanDeg: sum / cnt, maxDeg: mx };
   };
+  /** the planet's orbit normal of date n̂(t), a unit vector in the J2000 ecliptic (the ζ table at the chain's anchor) @param {number} year */
+  const orbitNormalEclipticJ2000AtYear = (year) => orbitNormalAt(year - 2000);
+  /**
+   * the spin axis of date ŝ(t), a unit vector in the J2000 ecliptic in the
+   * angular-momentum sense (Venus and Uranus point south of their orbits):
+   * the integrated pole, interpolated linearly between the track's samples
+   * and renormalized (≤ 1.2° of pole motion per sample — the chord error is
+   * arcsecond-class, below what any renderer resolves); Mercury (Cassini
+   * lock) carries the J2000 pole with its orbit normal; null beyond ±10 Myr
+   * @param {number} year
+   */
+  const spinAxisEclipticJ2000AtYear = (year) => {
+    if (cassiniLocked) return rotateAligned(n0, orbitNormalAt(year - 2000), s0);
+    const b = bracket(year - 2000);
+    if (!b) return null;
+    const a = b.tr.sv[b.i], c = b.tr.sv[b.i + 1];
+    return unit([a[0] + b.f * (c[0] - a[0]), a[1] + b.f * (c[1] - a[1]), a[2] + b.f * (c[2] - a[2])]);
+  };
 
   return Object.freeze({
     key,
@@ -350,6 +392,8 @@ function createPlanetSpinChannel(deps) {
     obliquityDegAtYear,
     spinPrecessionRateArcsecPerYrAtYear,
     obliquityEnvelopeDeg,
+    spinAxisEclipticJ2000AtYear,
+    orbitNormalEclipticJ2000AtYear,
     stepYr,
     maxSpanYr: MAX_SPAN_YR,
   });

package/types/chain-cycles/index.d.cts

@@ -5,7 +5,6 @@ export type PeriodSecondsAtAge = (tMa: number) => (number | null);
 /**
  * @param {{
  *   ageAnchorYear: number,
- *   tropicalYearSecondsAtAge: (tMa: number) => (number | null),
  *   tropicalYearJ2000Seconds: number,
  *   isDeepTime: () => boolean,
  *   gridMinYear?: number,
@@ -15,12 +14,18 @@ export type PeriodSecondsAtAge = (tMa: number) => (number | null);
  *   maxCacheEntries?: number,
  * }} deps — ageAnchorYear is the year whose age is 0 in the period
  *   functions' t_Ma coordinate (startmodelYear = 2000.5, the certified
- *   convention); tropicalYearJ2000Seconds feeds the snapshot branch;
+ *   convention); tropicalYearJ2000Seconds is the length of ONE tick of the
+ *   year axis the integrator is fed — the linear SI-year label, each tick
+ *   exactly one J2000 tropical year of seconds — so it is the integrand's
+ *   weight in BOTH branches (plan 09 step 7c′: the deep branch had weighted
+ *   each tick by the tropical year OF DATE, adding that year's deep-time
+ *   drift to every lunar count — ~5° of Moon at ±100 kyr, ~1.3 orbits at
+ *   1 Ma — although the tick itself never changes length; the owner's call:
+ *   a plain correctness fix, the tick already is an SI tick);
  *   isDeepTime is read PER CALL so the engine's runtime toggle works.
  */
-export function createChainCycleIntegrator({ ageAnchorYear, tropicalYearSecondsAtAge, tropicalYearJ2000Seconds, isDeepTime, gridMinYear, gridMaxYear, gridStepYears, gridAnchorYear, maxCacheEntries, }: {
+export function createChainCycleIntegrator({ ageAnchorYear, tropicalYearJ2000Seconds, isDeepTime, gridMinYear, gridMaxYear, gridStepYears, gridAnchorYear, maxCacheEntries, }: {
     ageAnchorYear: number;
-    tropicalYearSecondsAtAge: (tMa: number) => (number | null);
     tropicalYearJ2000Seconds: number;
     isDeepTime: () => boolean;
     gridMinYear?: number;

package/types/climate/l1-orbital.d.cts

@@ -73,7 +73,9 @@ export function laggedL1Terms(l1Terms: Array<{
  *   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/types/constants/coefficients.d.ts

@@ -1,74 +1,11 @@
 // 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: "82c86b078b32fff8";
+export declare const COEFFICIENTS_HASH: "91e2d041c7b6e038";
 
 export declare const FITTED_COEFFICIENTS: {
-  readonly ANOMALISTIC_YEAR_HARMONICS: Array<[number, number, number]>;
   readonly BALANCE_PRESETS: Array<Array<number | string>>;
-  readonly CARDINAL_POINT_ANCHORS: {
-    "SS": number;
-    "WS": number;
-    "VE": number;
-    "AE": number;
-  };
-  readonly CARDINAL_POINT_ANCHORS_ADJUSTED: {
-    "SS": number;
-    "WS": number;
-    "VE": number;
-    "AE": number;
-  };
-  readonly CARDINAL_POINT_DERIVED: {
-    "lincoef": number;
-    "h0": number;
-    "h1": number;
-    "eccOrderDivisors": [number, number];
-    "note": string;
-  };
-  readonly CARDINAL_POINT_ECC_TERMS: {
-    "SS": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "WS": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "VE": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "AE": Array<{
-      "order": number;
-      "sin": number;
-      "cos": number;
-    }>;
-  };
-  readonly CARDINAL_POINT_HARMONICS: {
-    "SS": Array<[number, number, number]>;
-    "WS": Array<[number, number, number]>;
-    "VE": Array<[number, number, number]>;
-    "AE": Array<[number, number, number]>;
-  };
-  readonly CARDINAL_POINT_JOINT_TERMS: {
-    "quadratureDeg": {
-      "SS": number;
-      "AE": number;
-      "WS": number;
-      "VE": number;
-    };
-    "terms": Array<{
-      "order": number;
-      "div": number;
-      "sin": number;
-      "cos": number;
-    }>;
-    "note": string;
-  };
   readonly CLIMATE_FORMULA_COEFFS: {
     "config": {
       "L1_lines": Array<{
@@ -303,105 +240,105 @@ export declare const FITTED_COEFFICIENTS: {
           "delta_l3": number;
         };
         "carbon_amplification_ratios": {
-          "19.0014": {
+          "19.0015": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "19.1554": {
+          "19.1555": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "22.4452": {
+          "22.4453": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "23.1931": {
+          "23.1933": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "23.7599": {
+          "23.7601": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "27.2689": {
+          "27.2691": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "28.8424": {
+          "28.8426": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.0047": {
+          "29.0050": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.6172": {
+          "29.6174": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.8091": {
+          "29.8094": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "39.5560": {
+          "39.5564": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.2292": {
+          "40.2297": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.9239": {
+          "40.9244": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "41.6445": {
+          "41.6450": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "53.6174": {
+          "53.6182": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
@@ -634,17 +571,4 @@ export declare const FITTED_COEFFICIENTS: {
   };
   readonly MEEUS_LATITUDE_TERMS: Array<number[]>;
   readonly MEEUS_LONGITUDE_TERMS: Array<number[]>;
-  readonly PERI_HARMONICS_RAW: Array<[number, number, number]>;
-  readonly PERI_OFFSET: number;
-  readonly RA_DAY_OFFSET_MS: {
-    "mean": number;
-    "terms": Array<[number, number, number]>;
-    "note": string;
-  };
-  readonly SIDEREAL_YEAR_HARMONICS: Array<[number, number, number]>;
-  readonly SOLSTICE_OBLIQUITY_HARMONICS: Array<[number, number, number]>;
-  readonly SOLSTICE_OBLIQUITY_MEAN_FITTED: number;
-  readonly SUN_LONGITUDE_HARMONICS: Array<[number, number, number]>;
-  readonly SUN_LONGITUDE_MEAN: number;
-  readonly TROPICAL_YEAR_HARMONICS: Array<[number, number, number]>;
 };

package/types/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: "6c2beb70bb3d5000";
+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: "6c2beb70bb3d5000";
-  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,48 +102,28 @@ export declare const DEFAULT_CONSTANTS: {
       "thermalPumpFactor": number;
     };
   };
-  readonly earth: {
-    "earthtiltMean": number;
-    "earthInvPlaneInclinationAmplitude": number;
-    "eccentricityBase": number;
-    "eccentricityAmplitude": number;
-    "kDeviceWheelPhaseAtJ2000Cycles": number;
-  };
   readonly earthOrbital: {
     "j2000EpochYear": number;
     "perihelionalignmentYear": number;
     "juneSolstice2000_JD": number;
     "obliquityJ2000_deg": number;
     "obliquityRate_arcsecPerCentury": number;
+    "generalPrecessionJ2000_arcsecPerYr": number;
     "earthEccentricityJ2000": number;
     "earthEccentricityDotJ2000": number;
     "earthEccentricityDotDotJ2000": number;
     "earthPerihelionLongitudeJ2000": number;
     "sunMeanLongitudeJ2000_deg": number;
     "earthAscendingNodeInvPlane": number;
-    "earthInclinationCycleAnchor": number;
     "perihelionPassageJ2000_JD": number;
     "earthInclinationJ2000_deg": number;
     "earthInclinationRate_arcsecPerCentury": number;
     "deltaTStart": number;
     "sunTilt": number;
   };
-  readonly foundational: {
-    "holisticyearLength": number;
-    "stepYears": number;
-    "inputmeanlengthsolaryearindays": number;
-    "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": {
@@ -211,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;
@@ -290,6 +221,7 @@ export declare const DEFAULT_CONSTANTS: {
       "invPlaneInclinationJ2000": number;
       "meanAnomaly": number;
       "trueAnomaly": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "venus": {
       "solarYearInput": number;
@@ -302,6 +234,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "mars": {
       "solarYearInput": number;
@@ -314,6 +247,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "jupiter": {
       "solarYearInput": number;
@@ -326,6 +260,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "saturn": {
       "solarYearInput": number;
@@ -338,6 +273,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "uranus": {
       "solarYearInput": number;
@@ -350,6 +286,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
     "neptune": {
       "solarYearInput": number;
@@ -362,6 +299,7 @@ export declare const DEFAULT_CONSTANTS: {
       "meanAnomaly": number;
       "trueAnomaly": number;
       "rotationPeriodDays": number;
+      "ascendingNodeInvPlaneClosedForm": number;
     };
   };
   readonly planetSpinPhysical: {
@@ -471,79 +409,9 @@ export declare const DEFAULT_CONSTANTS: {
       }>;
     };
   };
-  readonly planets: {
-    "mercury": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "venus": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "mars": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "jupiter": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "saturn": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "uranus": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-    "neptune": {
-      "name": string;
-      "eocFraction": number;
-      "startpos": number;
-      "angleCorrection": number;
-      "type": string;
-      "mirrorPair": string;
-      "ascendingNodeInvPlane": number;
-      "ascendingNodeCyclesIn8H": number;
-    };
-  };
   readonly timeReference: {
+    "sceneEpochJD": number;
+    "sceneEpochYear": number;
     "j2000JD": number;
     "julianCenturyDays": number;
     "gregorianStartJD": number;
@@ -556,10 +424,7 @@ export declare const DEFAULT_CONSTANTS: {
     "tropicalYearSS": number;
     "tropicalYearAE": number;
     "tropicalYearWS": number;
-    "tropicalYearMean": number;
-    "anomalisticYear": number;
     "siderealYear": number;
-    "iauPrecessionJ2000": number;
     "tropicalYearRateSecPerCentury": number;
     "solarDay": number;
     "siderealDay": number;
@@ -603,10 +468,6 @@ export declare const REFERENCE_DATA: {
     "greatattractorDistance": number;
     "milkywaySpeed": number;
   };
-  readonly giaCoxChaoPeltier: {
-    "dJ2DtPerYr": number;
-    "j2ToAlphaFactor": number;
-  };
   readonly jplEclipticInclinationTrends: {
     "mercury": number;
     "venus": number;
@@ -716,6 +577,7 @@ export declare const REFERENCE_DATA: {
     "meeusPeriLong2000AD": number;
     "meeusPeriLong2500AD": number;
     "meeusPeriLong3000AD": number;
+    "simon1994AnomalisticYearDays": number;
     "mainstreamAxialPrecKyr": number;
     "mainstreamPeriPrecKyr": number;
     "mainstreamObliqCycleKyr": number;