@essrt/physics 23.6.0 → 24.1.0

package/types/earth/deep-orbital-history.d.cts

@@ -98,3 +98,26 @@ export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSerie
         ellipticity: number;
     };
 };
+/**
+ * The CO-ROTATING anchored mode sum — the z-side evaluator of the history
+ * (the ϖ̇ fix, measured): the plain DC anchor residual is a spurious
+ * zero-frequency mode that preserves the J2000 VALUE but dilutes the local
+ * argument rate by |sum(0)|/|anchor| (z: 11.47 → 7.93″/yr, a 92 s-class
+ * anomalistic-year bias in the mode tier). Here the residual ROTATES at the
+ * raw sum's own local arg-rate g* = Im(ż·z̄)/|z|² (analytic), so at t = 0
+ * both the value AND the rate are exact, and at depth the residual is one
+ * |R|-amplitude term at a physical frequency instead of a DC offset (deep
+ * envelope measured unchanged: [0.0001, 0.0672] → [0.0003, 0.0678] over
+ * 10–50 Myr). Exported: e(t) = hypot(sum(t)) needs no spin integration,
+ * so a fixture can read a deep eccentricity node here (bit-identical at the
+ * grid nodes; the published surface is the grid's linear interpolation and
+ * stays on the history).
+ * @param {ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>} modes
+ * @param {[number, number]} anchor the J2000 value (q, p) the sum is pinned to
+ * @returns {(t: number) => [number, number]} t in years from J2000
+ */
+export function createCoRotatingAnchoredModeSum(modes: ReadonlyArray<{
+    omegaRadPerYr: number;
+    re: number;
+    im: number;
+}>, anchor: [number, number]): (t: number) => [number, number];

package/types/eclipse/besselian.d.cts

@@ -54,24 +54,6 @@ export type BesselianDeps = {
         speedOfLightKmS: number;
     };
 };
-/**
- * @typedef {Object} BesselianDeps
- * @property {(dDaysTT: number) => {lonDeg: number, latDeg: number, distKm: number}} moonFullAtDaysTT
- *   full-series Moon: ecliptic-of-date longitude/latitude (deg) + distance (km), TT axis (days since J2000)
- * @property {(jdUT: number) => number} sunLonDegAt - geometric mean sun longitude (deg), JD(UT) axis with the finder's internal ΔT
- * @property {(T: number) => number} sunCompletionDeg - planetary completion (deg) at T centuries TT, SUBTRACTED from the finder sun (see eclipse/sun-planetary-completion.cjs; the finders deliberately stay without it — their fitted anchors and certified canon statistics were produced on the bare form, and elongation-class timing absorbs the omission into the fitted phases)
- * @property {(year: number) => number} sunAberrationDegAt - the Sun's derived planetary aberration κ (deg) at calendar year, SUBTRACTED from the geometric Sun longitude (the apparent Sun; see the header)
- * @property {(jd: number) => number} deltaTSecondsAt - framework ΔT (J2000-zeroed convention)
- * @property {(year: number) => number} obliquityDegAt - framework obliquity (deg) at calendar year
- * @property {(year: number) => number} eccentricityAt - framework Earth-orbit eccentricity at calendar year
- * @property {(year: number) => number} perihelionLongitudeDegAt - framework perihelion longitude (deg) at calendar year
- * @property {(jd: number) => number} yearFromJD - the model's own JD → calendar-year conversion
- * @property {{ j2000JD: number, julianCenturyDays: number, earthDiameterKm: number,
- *   moonDiameterKm: number, sunDiameterKm: number, sunDistanceKm: number,
- *   earthFlatteningInverse: number, ttBridgeSeconds: number,
- *   gmstMeanSiderealT0Deg: number, gmstMeanSiderealRateDegPerDay: number,
- *   gmstMeanSiderealT2Deg: number, speedOfLightKmS: number }} constants
- */
 /** @param {BesselianDeps} deps */
 export function createBesselian(deps: BesselianDeps): {
     umbraGroundAt: (jdUT: number) => {
@@ -96,4 +78,41 @@ export function createBesselian(deps: BesselianDeps): {
         };
         centralDurationSeconds: number | null;
     };
+    gmstDeg: (jdUT: number) => number;
 };
+/**
+ * @typedef {Object} BesselianDeps
+ * @property {(dDaysTT: number) => {lonDeg: number, latDeg: number, distKm: number}} moonFullAtDaysTT
+ *   full-series Moon: ecliptic-of-date longitude/latitude (deg) + distance (km), TT axis (days since J2000)
+ * @property {(jdUT: number) => number} sunLonDegAt - geometric mean sun longitude (deg), JD(UT) axis with the finder's internal ΔT
+ * @property {(T: number) => number} sunCompletionDeg - planetary completion (deg) at T centuries TT, SUBTRACTED from the finder sun (see eclipse/sun-planetary-completion.cjs; the finders deliberately stay without it — their fitted anchors and certified canon statistics were produced on the bare form, and elongation-class timing absorbs the omission into the fitted phases)
+ * @property {(year: number) => number} sunAberrationDegAt - the Sun's derived planetary aberration κ (deg) at calendar year, SUBTRACTED from the geometric Sun longitude (the apparent Sun; see the header)
+ * @property {(jd: number) => number} deltaTSecondsAt - framework ΔT (J2000-zeroed convention)
+ * @property {(year: number) => number} obliquityDegAt - framework obliquity (deg) at calendar year
+ * @property {(year: number) => number} eccentricityAt - framework Earth-orbit eccentricity at calendar year
+ * @property {(year: number) => number} perihelionLongitudeDegAt - framework perihelion longitude (deg) at calendar year
+ * @property {(jd: number) => number} yearFromJD - the model's own JD → calendar-year conversion
+ * @property {{ j2000JD: number, julianCenturyDays: number, earthDiameterKm: number,
+ *   moonDiameterKm: number, sunDiameterKm: number, sunDistanceKm: number,
+ *   earthFlatteningInverse: number, ttBridgeSeconds: number,
+ *   gmstMeanSiderealT0Deg: number, gmstMeanSiderealRateDegPerDay: number,
+ *   gmstMeanSiderealT2Deg: number, speedOfLightKmS: number }} constants
+ */
+/** Mean sidereal time at Greenwich, degrees (mean equinox — pairs with the
+ *  mean-of-date body longitudes; coefficients single-sourced from
+ *  astro-reference.json physicalConstants). THE ONE GMST of the model: the
+ *  eclipse location tier, the API's observer tier, the eclipse and lunar
+ *  audits and the browser's diagnostics all read it. Meeus 12.4's
+ *  cubic term (−T³/38,710,000) is deliberately absent: it is 0.9″ in the
+ *  eclipse canon and at most 93″ (6 s of time) at the edge of the observer
+ *  era, and every recorded eclipse and lunar value was measured without it.
+ *  @param {number} jdUT
+ *  @param {{ j2000JD: number, julianCenturyDays: number, gmstMeanSiderealT0Deg: number, gmstMeanSiderealRateDegPerDay: number, gmstMeanSiderealT2Deg: number }} K
+ *  @returns {number} */
+export function computeGmstDeg(jdUT: number, K: {
+    j2000JD: number;
+    julianCenturyDays: number;
+    gmstMeanSiderealT0Deg: number;
+    gmstMeanSiderealRateDegPerDay: number;
+    gmstMeanSiderealT2Deg: number;
+}): number;

package/types/index.d.ts

@@ -66,7 +66,7 @@ export { computeMassChainJ2000, computeMoonPeriodsJ2000 } from "./earth/mass-cha
 export { SECONDS_PER_DAY, DAYS_PER_JULIAN_YEAR, JULIAN_CENTURY_DAYS, SECONDS_PER_JULIAN_YEAR } from "./constants/units.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, computeApsidalSecularDegPerYr as computeApsidalTangentDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from "./planets/keplerian-chain.cjs";
+export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, computeApsidalTangentDegPerYr, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from "./planets/keplerian-chain.cjs";
 export { computeEquatorNodeOriginSFrameDeg, convertNodeSFrameToEquatorOriginDeg } from "./planets/inv-plane-frame.cjs";
 export { computePlanetPrecessionConstant, createPlanetSpinChannel, createPlanetSpinChannelFromArtifacts, computeObliquityJ2000Deg } from "./planets/spin-channel.cjs";
 export { CHAIN_ARTIFACT, CHAIN_ARTIFACT_HASH } from "./planets/chain-artifact.js";

package/types/model/sun-and-eclipses.d.ts

@@ -23,6 +23,7 @@ export function assembleSunAndEclipses(ctx: import("./lunar-theory.js").Ctx): {
             };
             centralDurationSeconds: number | null;
         };
+        gmstDeg: (jdUT: number) => number;
     };
     eclipseFinders: {
         sunLonDegAt: (jd: number) => number;

package/types/model/surface.d.ts

@@ -227,6 +227,8 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
             latDeg: number;
             lonDeg: number;
         } | null;
+        /** THE ONE GMST — mean sidereal time at Greenwich, degrees, at a UT JD (eclipse/besselian.cjs computeGmstDeg; the observer tier's hour angle). @param {number} jdUT @returns {number} */
+        gmstDegAtJD: (jdUT: number) => number;
         solarLocalCircumstances: (jdGreatest: number, latDeg: number, lonDeg: number) => {
             kind: "none" | "partial" | "annular" | "total";
             magnitude: number;

package/types/model.d.ts

@@ -183,6 +183,7 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
             latDeg: number;
             lonDeg: number;
         } | null;
+        gmstDegAtJD: (jdUT: number) => number;
         solarLocalCircumstances: (jdGreatest: number, latDeg: number, lonDeg: number) => {
             kind: "none" | "partial" | "annular" | "total";
             magnitude: number;

package/types/planets/keplerian-chain.d.cts

@@ -155,7 +155,7 @@ export function computePlanetElementsAtYear(year: number, planetChain: KcPlanetC
  * @param {Object<string,KcPlanetChain>=} allChains the full chain map
  * @returns {number} deg/yr, prograde positive
  */
-export function computeApsidalSecularDegPerYr(year: number, planetChain: KcPlanetChain, allChains?: {
+export function computeApsidalTangentDegPerYr(year: number, planetChain: KcPlanetChain, allChains?: {
     [x: string]: KcPlanetChain;
 } | undefined): number;
 /** The chain's canonical constructor from a PARSED governed artifact (K2.1 →

package/types/planets/model.d.cts

@@ -1,8 +1,8 @@
 export type PlanetModelBody = {
     /**
-     * - the observed sidereal period input, SI days (the scaffold period of the no-chain bodies; the chain planets ride the artifact's elements)
+     * - the body's orbital-period input, SI days — for the planets the period OF DATE (measured in sun-and-eclipses.js: their rates sit 1.4–1.6°/cy above the VSOP sidereal mean motions, the precession plus input rounding), not the sidereal period the field was once documented as; the scaffold period of the no-chain bodies — the chain planets ride the artifact's elements
      */
-    solarYearInput: number;
+    orbitalPeriodInputDays: number;
     invPlaneInclinationJ2000?: number | undefined;
     longitudePerihelion?: number | undefined;
     /**
@@ -43,7 +43,7 @@ export type PlanetModelRecord = {
 };
 /**
  * @typedef {Object} PlanetModelBody
- * @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} orbitalPeriodInputDays - the body's orbital-period input, SI days — for the planets the period OF DATE (measured in sun-and-eclipses.js: their rates sit 1.4–1.6°/cy above the VSOP sidereal mean motions, the precession plus input rounding), not the sidereal period the field was once documented as; the scaffold period of the no-chain bodies — the chain planets ride the artifact's elements
  * @property {number} [invPlaneInclinationJ2000]
  * @property {number} [longitudePerihelion]
  * @property {number} [wobblePeriodYears] - carriers: the chain's g-mode beat (input since Phase 7 commit 2)