@essrt/physics 21.0.0 → 22.0.1

package/src/planets/model.cjs

@@ -7,8 +7,9 @@
  *   inputs     the chain's g-mode beat (wobble) and the spin channel's
  *              derived J2000 obliquity, carried on the record since plan 06
  *              Phase 7 commit 2
- *   geometry   (the type-branched ellipse family on the OBSERVED J2000
- *              eccentricity; minor bodies use their record's base)
+ *   geometry   (the ellipse family — geometry.cjs's own taxonomy since
+ *              plan 11 step 3 — on the OBSERVED J2000 eccentricity; minor
+ *              bodies use their record's base)
  *
  * The ψ and K amplitude laws that opened this chain (ψ constant →
  * inclination amplitude/mean; K constant → eccentricity amplitude/base/
@@ -35,7 +36,6 @@ const { derivePlanetGeometry, scaffoldOrbitInputs } = require('./geometry.cjs');
 
 /**
  * @typedef {Object} PlanetModelBody
- * @property {string} [type] - geometry type branch ('I' | 'II' | 'III')
  * @property {number} solarYearInput - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]
@@ -108,11 +108,12 @@ function createPlanetModel(env, bodies) {
     // eccentricityPhaseJ2000Deg stood here. With the law gone, the ellipse
     // geometry below reads the OBSERVED J2000 eccentricity where it used to
     // read the law's System-Reset base — the observation the chain is
-    // anchored on, not a construction. The stored base survives only for the
-    // additional bodies, which carry theirs in model-parameters.json.
+    // anchored on, not a construction. The additional bodies' base is their
+    // observed J2000 eccentricity too (astro-reference additionalBodiesReference;
+    // plan 11 step 2 deleted the stored copy).
 
     rec.geometry = derivePlanetGeometry({
-      key, type: b.type,
+      key,   // the ellipse family is the law's own taxonomy (geometry.cjs SCAFFOLD_ELLIPSE_FAMILY, plan 11 step 3)
       // plan 09 step 7b: the chain's J2000 elements where the artifact has the body, else the observed period input
       ...scaffoldOrbitInputs(key, env.chainArtifact, b.solarYearInput),
       orbitalEccentricityBase: b.orbitalEccentricityBase !== undefined

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/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,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/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: "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/types/earth/frame-of-date.d.cts

@@ -35,3 +35,12 @@ export function computeEarthFrameOfDate(sample: {
  * @returns {number} θ in radians, in (−π, π]
  */
 export function solveWheelAngleForLongitude(lambdaRad: number, e: number, periRad: number): number;
+/**
+ * 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
+ */
+export function computeSceneWorldFrameJ2000(): number[][];
+/** The J2000 equinox's azimuth in the scene's world XZ plane, degrees (atan2(x, z)) — the view convention: the equinox along world +X. */
+export const SCENE_EQUINOX_VIEW_AZIMUTH_DEG: 90;

package/types/index.d.ts

@@ -15,6 +15,7 @@ export { createDeepEccChannel } from "./moon/deep-ecc-channel.cjs";
 export { createDeepOrbitalHistory } from "./earth/deep-orbital-history.cjs";
 export { createSiderealYearChannel } from "./earth/sidereal-year-channel.cjs";
 export { createPlanetModel } from "./planets/model.cjs";
+export { PLANET_DISPLAY_NAMES } from "./planets/display-names.cjs";
 export { createDeltaTCycles } from "./deltat/cycles.cjs";
 export { createDeepTimeLod } from "./deltat/deep-time.cjs";
 export { deltaTEspenakMeeusCanonSeconds } from "./deltat/historical.cjs";
@@ -56,7 +57,7 @@ import { REFERENCE_DATA } from './constants/index.js';
 import { assembleModel } from './model.js';
 export { GENERATED as DEFAULT_CONSTANTS, CONSTANTS_HASH, MODEL_VERSION, PREPRINT_DOI, REFERENCE_DATA };
 export { FITTED_COEFFICIENTS, COEFFICIENTS_HASH } from "./constants/index.js";
-export { computeEarthFrameOfDate, solveWheelAngleForLongitude } from "./earth/frame-of-date.cjs";
+export { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG } from "./earth/frame-of-date.cjs";
 export { derivePlanetGeometry, scaffoldOrbitInputs } from "./planets/geometry.cjs";
 export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, computeApsidalSecularDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from "./planets/keplerian-chain.cjs";
 export { computeEquatorNodeOriginSFrameDeg, convertNodeSFrameToEquatorOriginDeg } from "./planets/inv-plane-frame.cjs";

package/types/model.d.ts

@@ -261,6 +261,8 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
                 meanDeg: number;
                 maxDeg: number;
             } | null;
+            spinAxisEclipticJ2000AtYear: (year: number) => number[] | null;
+            orbitNormalEclipticJ2000AtYear: (year: number) => number[];
             stepYr: number;
             maxSpanYr: 10000000;
         }>;

package/types/moon/series.d.cts

@@ -12,7 +12,6 @@ export type MeeusTermTable = Array<[number, number, number, number, number]>;
  *     moonDistanceJ2000Km: number,
  *     j2000JD: number,
  *     julianCenturyDays: number,
- *     moonMeeusLpCorrectionDeg: number,
  *     fwA2RateDegPerCy: number,
  *     fwA3RateDegPerCy: number,
  *   },
@@ -46,7 +45,6 @@ export function createMoonSeries({ constants, fns }: {
         moonDistanceJ2000Km: number;
         j2000JD: number;
         julianCenturyDays: number;
-        moonMeeusLpCorrectionDeg: number;
         fwA2RateDegPerCy: number;
         fwA3RateDegPerCy: number;
     };

package/types/planets/display-names.d.cts

@@ -0,0 +1,2 @@
+/** @type {Readonly<Record<string, string>>} */
+export const PLANET_DISPLAY_NAMES: Readonly<Record<string, string>>;

package/types/planets/geometry.d.cts

@@ -1,6 +1,5 @@
 export type PlanetGeometryBody = {
     key: string;
-    type?: string;
     periodDays: number;
     semiMajorAxisAu?: number;
     orbitalEccentricityBase?: number;
@@ -12,7 +11,6 @@ export type PlanetGeometryBody = {
 /**
  * @typedef {{
  *   key: string,
- *   type?: string,
  *   periodDays: number,
  *   semiMajorAxisAu?: number,
  *   orbitalEccentricityBase?: number,
@@ -40,6 +38,7 @@ export type PlanetGeometryBody = {
  *   elipticOrbit: (number | undefined), perihelionDistance: (number | undefined),
  *   speedKmh: number, rotationPeriodHours: (number | undefined),
  *   eccentricityPerihelion: (number | undefined), lowestPoint: (number | undefined),
+ *   ellipseFamily: (string | undefined),
  * }}
  */
 export function derivePlanetGeometry(body: PlanetGeometryBody, env: {
@@ -58,6 +57,7 @@ export function derivePlanetGeometry(body: PlanetGeometryBody, env: {
     rotationPeriodHours: (number | undefined);
     eccentricityPerihelion: (number | undefined);
     lowestPoint: (number | undefined);
+    ellipseFamily: (string | undefined);
 };
 /**
  * The scaffold's orbit inputs for one body — ONE home for the rule (plan 09
@@ -74,3 +74,10 @@ export function scaffoldOrbitInputs(key: string, chainArtifact: any, fallbackPer
     periodDays: number;
     semiMajorAxisAu: (number | undefined);
 };
+/**
+ * Each body's scaffold ellipse family (the survey's measured taxonomy above).
+ * Pluto, Halley and Ceres take their body-name branches; a key absent here
+ * has no ellipse family (geometry-only).
+ * @type {Readonly<Record<string, string>>}
+ */
+export const SCAFFOLD_ELLIPSE_FAMILY: Readonly<Record<string, string>>;

package/types/planets/model.d.cts

@@ -1,8 +1,4 @@
 export type PlanetModelBody = {
-    /**
-     * - geometry type branch ('I' | 'II' | 'III')
-     */
-    type?: string | undefined;
     /**
      * - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
      */
@@ -47,7 +43,6 @@ export type PlanetModelRecord = {
 };
 /**
  * @typedef {Object} PlanetModelBody
- * @property {string} [type] - geometry type branch ('I' | 'II' | 'III')
  * @property {number} solarYearInput - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]

package/types/planets/spin-channel.d.cts

@@ -99,6 +99,8 @@ export function createPlanetSpinChannel(deps: {
         meanDeg: number;
         maxDeg: number;
     } | null;
+    spinAxisEclipticJ2000AtYear: (year: number) => number[] | null;
+    orbitNormalEclipticJ2000AtYear: (year: number) => number[];
     stepYr: number;
     maxSpanYr: 10000000;
 }>;
@@ -152,6 +154,8 @@ export function createPlanetSpinChannelFromArtifacts(deps: {
         meanDeg: number;
         maxDeg: number;
     } | null;
+    spinAxisEclipticJ2000AtYear: (year: number) => number[] | null;
+    orbitNormalEclipticJ2000AtYear: (year: number) => number[];
     stepYr: number;
     maxSpanYr: 10000000;
 }>;