@neaps/tide-predictor 0.3.0 → 0.4.0

package/dist/index.js

@@ -181,132 +181,6 @@ const astro = (time) => {
 };
 var astronomy_default = astro;
 
-//#endregion
-//#region src/harmonics/prediction.ts
-const modulus = (a, b) => {
-	return (a % b + b) % b;
-};
-const addExtremesOffsets = (extreme, offsets) => {
-	if (typeof offsets === "undefined" || !offsets) return extreme;
-	if (extreme.high && offsets.height?.high) if (offsets.height.type === "fixed") extreme.level += offsets.height.high;
-	else extreme.level *= offsets.height.high;
-	if (extreme.low && offsets.height?.low) if (offsets.height.type === "fixed") extreme.level += offsets.height.low;
-	else extreme.level *= offsets.height.low;
-	if (extreme.high && offsets.time?.high) extreme.time = new Date(extreme.time.getTime() + offsets.time.high * 60 * 1e3);
-	if (extreme.low && offsets.time?.low) extreme.time = new Date(extreme.time.getTime() + offsets.time.low * 60 * 1e3);
-	return extreme;
-};
-const getExtremeLabel = (label, highLowLabels) => {
-	if (typeof highLowLabels !== "undefined" && typeof highLowLabels[label] !== "undefined") return highLowLabels[label];
-	return {
-		high: "High",
-		low: "Low"
-	}[label];
-};
-const predictionFactory = ({ timeline, constituents: constituents$1, start }) => {
-	const getLevel = (hour, modelBaseSpeed, modelU, modelF, modelBaseValue) => {
-		const amplitudes = [];
-		let result = 0;
-		constituents$1.forEach((constituent) => {
-			const amplitude = constituent.amplitude;
-			const phase = constituent.phase;
-			const f = modelF[constituent.name];
-			const speed = modelBaseSpeed[constituent.name];
-			const u = modelU[constituent.name];
-			const V0 = modelBaseValue[constituent.name];
-			amplitudes.push(amplitude * f * Math.cos(speed * hour + (V0 + u) - phase));
-		});
-		amplitudes.forEach((item) => {
-			result += item;
-		});
-		return result;
-	};
-	const prediction = {};
-	prediction.getExtremesPrediction = (options) => {
-		const { labels, offsets } = typeof options !== "undefined" ? options : {};
-		const results = [];
-		const { baseSpeed, u, f, baseValue } = prepare();
-		let goingUp = false;
-		let goingDown = false;
-		let lastLevel = getLevel(0, baseSpeed, u[0], f[0], baseValue);
-		timeline.items.forEach((time, index) => {
-			const hour = timeline.hours[index];
-			const level = getLevel(hour, baseSpeed, u[index], f[index], baseValue);
-			if (level > lastLevel && goingDown) results.push(addExtremesOffsets({
-				time: timeline.items[index - 1],
-				level: lastLevel,
-				high: false,
-				low: true,
-				label: getExtremeLabel("low", labels)
-			}, offsets));
-			if (level < lastLevel && goingUp) results.push(addExtremesOffsets({
-				time: timeline.items[index - 1],
-				level: lastLevel,
-				high: true,
-				low: false,
-				label: getExtremeLabel("high", labels)
-			}, offsets));
-			if (level > lastLevel) {
-				goingUp = true;
-				goingDown = false;
-			}
-			if (level < lastLevel) {
-				goingUp = false;
-				goingDown = true;
-			}
-			lastLevel = level;
-		});
-		return results;
-	};
-	prediction.getTimelinePrediction = () => {
-		const results = [];
-		const { baseSpeed, u, f, baseValue } = prepare();
-		timeline.items.forEach((time, index) => {
-			const hour = timeline.hours[index];
-			const prediction$1 = {
-				time,
-				hour,
-				level: getLevel(hour, baseSpeed, u[index], f[index], baseValue)
-			};
-			results.push(prediction$1);
-		});
-		return results;
-	};
-	const prepare = () => {
-		const baseAstro = astronomy_default(start);
-		const baseValue = {};
-		const baseSpeed = {};
-		const u = [];
-		const f = [];
-		constituents$1.forEach((constituent) => {
-			const value = constituent._model.value(baseAstro);
-			const speed = constituent._model.speed(baseAstro);
-			baseValue[constituent.name] = d2r * value;
-			baseSpeed[constituent.name] = d2r * speed;
-		});
-		timeline.items.forEach((time) => {
-			const uItem = {};
-			const fItem = {};
-			const itemAstro = astronomy_default(time);
-			constituents$1.forEach((constituent) => {
-				const constituentU = modulus(constituent._model.u(itemAstro), 360);
-				uItem[constituent.name] = d2r * constituentU;
-				fItem[constituent.name] = modulus(constituent._model.f(itemAstro), 360);
-			});
-			u.push(uItem);
-			f.push(fItem);
-		});
-		return {
-			baseValue,
-			baseSpeed,
-			u,
-			f
-		};
-	};
-	return Object.freeze(prediction);
-};
-var prediction_default = predictionFactory;
-
 //#endregion
 //#region src/node-corrections/index.ts
 const corrections = {
@@ -429,45 +303,10 @@ const corrections = {
 var node_corrections_default = corrections;
 
 //#endregion
-//#region src/constituents/constituent.ts
-/**
-* Computes the dot notation of two arrays
-*/
-const dotArray = (a, b) => {
-	const results = [];
-	a.forEach((value, index) => {
-		results.push(value * b[index]);
-	});
-	return results.reduce((total, value) => total + value);
-};
-const astronimicDoodsonNumber = (astro$1) => {
-	return [
-		astro$1["T+h-s"],
-		astro$1.s,
-		astro$1.h,
-		astro$1.p,
-		astro$1.N,
-		astro$1.pp,
-		astro$1["90"]
-	];
-};
-const astronomicSpeed = (astro$1) => {
-	const results = [];
-	astronimicDoodsonNumber(astro$1).forEach((number) => {
-		results.push(number.speed);
-	});
-	return results;
-};
-const astronomicValues = (astro$1) => {
-	const results = [];
-	astronimicDoodsonNumber(astro$1).forEach((number) => {
-		results.push(number.value);
-	});
-	return results;
-};
-const constituentFactory = (name, coefficients$1, u, f) => {
+//#region src/constituents/definition.ts
+function defineConstituent(name, coefficients$1, u, f) {
 	if (!coefficients$1) throw new Error("Coefficient must be defined for a constituent");
-	const constituent = {
+	return Object.freeze({
 		name,
 		coefficients: coefficients$1,
 		value: (astro$1) => {
@@ -478,14 +317,9 @@ const constituentFactory = (name, coefficients$1, u, f) => {
 		},
 		u: typeof u !== "undefined" ? u : node_corrections_default.uZero,
 		f: typeof f !== "undefined" ? f : node_corrections_default.fUnity
-	};
-	return Object.freeze(constituent);
-};
-var constituent_default = constituentFactory;
-
-//#endregion
-//#region src/constituents/compound-constituent.ts
-const compoundConstituentFactory = (name, members) => {
+	});
+}
+function defineCompoundConstituent(name, members) {
 	const coefficients$1 = [];
 	members.forEach(({ constituent, factor }) => {
 		constituent.coefficients.forEach((coefficient, index) => {
@@ -493,7 +327,7 @@ const compoundConstituentFactory = (name, members) => {
 			coefficients$1[index] += coefficient * factor;
 		});
 	});
-	const compoundConstituent = {
+	return Object.freeze({
 		name,
 		coefficients: coefficients$1,
 		speed: (astro$1) => {
@@ -524,69 +358,71 @@ const compoundConstituentFactory = (name, members) => {
 			});
 			return f.reduce((previous, value) => previous * value);
 		}
-	};
-	return Object.freeze(compoundConstituent);
-};
-var compound_constituent_default = compoundConstituentFactory;
+	});
+}
+/**
+* Computes the dot notation of two arrays
+*/
+function dotArray(a, b) {
+	const results = [];
+	a.forEach((value, index) => {
+		results.push(value * b[index]);
+	});
+	return results.reduce((total, value) => total + value);
+}
+function astronimicDoodsonNumber(astro$1) {
+	return [
+		astro$1["T+h-s"],
+		astro$1.s,
+		astro$1.h,
+		astro$1.p,
+		astro$1.N,
+		astro$1.pp,
+		astro$1["90"]
+	];
+}
+function astronomicSpeed(astro$1) {
+	const results = [];
+	astronimicDoodsonNumber(astro$1).forEach((number) => {
+		results.push(number.speed);
+	});
+	return results;
+}
+function astronomicValues(astro$1) {
+	const results = [];
+	astronimicDoodsonNumber(astro$1).forEach((number) => {
+		results.push(number.value);
+	});
+	return results;
+}
+var definition_default = {};
 
 //#endregion
-//#region src/constituents/index.ts
-const constituents = {};
-constituents.Z0 = constituent_default("Z0", [
-	0,
-	0,
-	0,
-	0,
-	0,
-	0,
-	0
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.SA = constituent_default("Sa", [
-	0,
-	0,
-	1,
-	0,
-	0,
-	0,
-	0
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.SSA = constituent_default("Ssa", [
-	0,
-	0,
+//#region src/constituents/M2.ts
+/**
+* Lunar semi-diurnal (M2).
+* Primary principal lunar constituent; largest semi-diurnal tidal component globally.
+* Amplitude varies 10-20% over lunar node cycle; typically 0.2-0.5m in coastal areas.
+* Base constituent for many compound shallow-water constituents.
+*/
+var M2_default = defineConstituent("M2", [
 	2,
 	0,
 	0,
 	0,
-	0
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.MM = constituent_default("MM", [
-	0,
-	1,
-	0,
-	-1,
-	0,
-	0,
-	0
-], node_corrections_default.uZero, node_corrections_default.fMm);
-constituents.MF = constituent_default("MF", [
-	0,
-	2,
-	0,
-	0,
 	0,
 	0,
 	0
-], node_corrections_default.uMf, node_corrections_default.fMf);
-constituents.Q1 = constituent_default("Q1", [
-	1,
-	-2,
-	0,
-	1,
-	0,
-	0,
-	1
-], node_corrections_default.uO1, node_corrections_default.fO1);
-constituents.O1 = constituent_default("O1", [
+], node_corrections_default.uM2, node_corrections_default.fM2);
+
+//#endregion
+//#region src/constituents/O1.ts
+/**
+* Lunar diurnal (O1).
+* Primary lunar diurnal constituent; one-per-lunar-day oscillation.
+* Amplitude varies 10-20% due to lunar node effects.
+*/
+var O1_default = defineConstituent("O1", [
 	1,
 	-1,
 	0,
@@ -595,7 +431,30 @@ constituents.O1 = constituent_default("O1", [
 	0,
 	1
 ], node_corrections_default.uO1, node_corrections_default.fO1);
-constituents.K1 = constituent_default("K1", [
+
+//#endregion
+//#region src/constituents/2MK3.ts
+/**
+* Shallow-water terdiurnal (2MK3 = M2 + O1).
+* Lunar-lunar interaction from M2 semi-diurnal and O1 diurnal components.
+* Generated in shallow-water environments; typical amplitude <1 cm.
+*/
+var _2MK3_default = defineCompoundConstituent("2MK3", [{
+	constituent: M2_default,
+	factor: 1
+}, {
+	constituent: O1_default,
+	factor: 1
+}]);
+
+//#endregion
+//#region src/constituents/K1.ts
+/**
+* Lunisolar diurnal (K1).
+* Combined lunar and solar diurnal constituent; strongest diurnal tide in many regions.
+* Often comparable in amplitude to O1; amplitude ratio K1/O1 varies with latitude.
+*/
+var K1_default = defineConstituent("K1", [
 	1,
 	1,
 	0,
@@ -604,52 +463,88 @@ constituents.K1 = constituent_default("K1", [
 	0,
 	-1
 ], node_corrections_default.uK1, node_corrections_default.fK1);
-constituents.J1 = constituent_default("J1", [
-	1,
-	2,
-	0,
-	-1,
-	0,
-	0,
-	-1
-], node_corrections_default.uJ1, node_corrections_default.fJ1);
-constituents.M1 = constituent_default("M1", [
-	1,
-	0,
-	0,
-	0,
-	0,
-	0,
-	1
-], node_corrections_default.uM1, node_corrections_default.fM1);
-constituents.P1 = constituent_default("P1", [
-	1,
-	1,
+
+//#endregion
+//#region src/constituents/2MK5.ts
+/**
+* Shallow-water fifth-diurnal from M2-K1 interaction.
+*
+* Note: Found in coastal tide predictions, especially in enclosed bays and harbors.
+* Amplitude typically 0.1-0.5 cm depending on location and water depth.
+* Shallow-water constituent only; not present in deep ocean.
+*
+* @see NOAA CO-OPS
+*/
+var _2MK5_default = defineCompoundConstituent("2MK5", [{
+	constituent: M2_default,
+	factor: 2
+}, {
+	constituent: K1_default,
+	factor: 1
+}]);
+
+//#endregion
+//#region src/constituents/2MO5.ts
+/**
+* Shallow-water fifth-diurnal from M2-O1 interaction.
+*
+* Note: Primarily shallow-water coastal phenomenon, not present in deep ocean.
+* Amplitude typically <0.5 cm except in extreme shallow-water or enclosed basins.
+* Found in coastal tide predictions alongside other shallow-water constituents.
+*
+* @see NOAA CO-OPS
+*/
+var _2MO5_default = defineCompoundConstituent("2MO5", [{
+	constituent: M2_default,
+	factor: 2
+}, {
+	constituent: O1_default,
+	factor: 1
+}]);
+
+//#endregion
+//#region src/constituents/S2.ts
+/**
+* Solar semi-diurnal (S2).
+* Principal solar semi-diurnal constituent; largest solar tide component.
+* Amplitude typically 20-50% of M2; varies with geographic location and latitude.
+* Ratio M2/S2 determines tidal form factor (diurnal vs semi-diurnal regimes).
+* Base constituent for many compound shallow-water constituents.
+*/
+var S2_default = defineConstituent("S2", [
+	2,
+	2,
 	-2,
 	0,
 	0,
 	0,
-	1
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.S1 = constituent_default("S1", [
-	1,
-	1,
-	-1,
-	0,
-	0,
-	0,
 	0
 ], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.OO1 = constituent_default("OO1", [
-	1,
-	3,
-	0,
-	0,
-	0,
-	0,
-	-1
-], node_corrections_default.uOO1, node_corrections_default.fOO1);
-constituents["2N2"] = constituent_default("2N2", [
+
+//#endregion
+//#region src/constituents/2MS6.ts
+/**
+* Sixth-diurnal shallow-water interaction: 2×M2 + S2.
+*
+* Generated only in shallow water; typical amplitude 0.01-0.2 meters depending on depth.
+* Included in IHO shallow-water constituent tables and NOAA analysis (order #37).
+*/
+var _2MS6_default = defineCompoundConstituent("2MS6", [{
+	constituent: M2_default,
+	factor: 2
+}, {
+	constituent: S2_default,
+	factor: 1
+}]);
+
+//#endregion
+//#region src/constituents/2N2.ts
+/**
+* Lunar semi-diurnal (2N2).
+* Second-order lunar semi-diurnal from Moon's orbital ellipticity.
+* Amplitude typically 5-10% of M2; significant in semi-diurnal analysis.
+*/
+var _2N2_default = defineConstituent("2N2", [
 	2,
 	-2,
 	0,
@@ -658,7 +553,15 @@ constituents["2N2"] = constituent_default("2N2", [
 	0,
 	0
 ], node_corrections_default.uM2, node_corrections_default.fM2);
-constituents.N2 = constituent_default("N2", [
+
+//#endregion
+//#region src/constituents/N2.ts
+/**
+* Lunar elliptic semi-diurnal (N2).
+* Primary lunar elliptic semi-diurnal constituent from Moon's orbital variations.
+* Amplitude typically 5-15% of M2; third-largest semi-diurnal constituent.
+*/
+var N2_default = defineConstituent("N2", [
 	2,
 	-1,
 	0,
@@ -667,25 +570,149 @@ constituents.N2 = constituent_default("N2", [
 	0,
 	0
 ], node_corrections_default.uM2, node_corrections_default.fM2);
-constituents.NU2 = constituent_default("NU2", [
+
+//#endregion
+//#region src/constituents/J1.ts
+/**
+* Lunar diurnal (J1).
+* Shallow-water lunar diurnal constituent; generally much smaller than O1 and K1.
+* Typically important only in shallow-water systems and enclosed basins.
+* Amplitude usually <5% of O1; rarely significant in routine predictions.
+*/
+var J1_default = defineConstituent("J1", [
+	1,
 	2,
+	0,
 	-1,
+	0,
+	0,
+	-1
+], node_corrections_default.uJ1, node_corrections_default.fJ1);
+
+//#endregion
+//#region src/constituents/2Q1.ts
+/**
+* Shallow-water diurnal (2Q1 = N2-J1).
+* Derived from interaction of semi-diurnal N2 and diurnal J1 constituents.
+* Amplitude typically very small; rarely significant.
+* Found only in shallow-water regions with strong tidal distortion.
+*/
+var _2Q1_default = defineCompoundConstituent("2Q1", [{
+	constituent: N2_default,
+	factor: 1
+}, {
+	constituent: J1_default,
+	factor: -1
+}]);
+
+//#endregion
+//#region src/constituents/2SM2.ts
+/**
+* Shallow-water semi-diurnal (2SM2).
+* Compound constituent: 2×S2 - M2
+* Solar-lunar interaction constituent generated in shallow-water environments.
+* Amplitude typically <2% of M2; complementary to MU2.
+*/
+var _2SM2_default = defineCompoundConstituent("2SM2", [{
+	constituent: S2_default,
+	factor: 2
+}, {
+	constituent: M2_default,
+	factor: -1
+}]);
+
+//#endregion
+//#region src/constituents/L2.ts
+/**
+* Lunar elliptic semi-diurnal (L2).
+* Secondary lunar elliptic semi-diurnal from Moon's orbital variations and perigee effects.
+* Amplitude typically 1-3% of M2; often grouped with other lunar elliptic constituents.
+* Important in detailed harmonic analyses of semi-diurnal tides.
+*/
+var L2_default = defineConstituent("L2", [
 	2,
+	1,
+	0,
 	-1,
 	0,
 	0,
+	2
+], node_corrections_default.uL2, node_corrections_default.fL2);
+
+//#endregion
+//#region src/constituents/3L2.ts
+/**
+* Triple lunar elliptic; 3 times L2 interaction.
+*
+* Warning: Not in standard IHO constituent bank. Mainly historical/theoretical interest.
+* Very small amplitude (<0.1 cm); found only in extreme shallow water or enclosed basins.
+* Often ignored in routine tide predictions.
+*
+* @see Schureman, P. (1958). Manual of Harmonic Analysis and Prediction of Tides
+*/
+var _3L2_default = defineCompoundConstituent("3L2", [{
+	constituent: L2_default,
+	factor: 3
+}]);
+
+//#endregion
+//#region src/constituents/3N2.ts
+/**
+* Triple N2 shallow-water harmonic (3 × N2).
+*
+* Note: Shallow-water constituent with definition based on compound frequency estimates.
+* Typical amplitude <0.5 cm; often <0.1 cm except in extreme shallow-water or enclosed basins.
+* Rarely significant in routine tide predictions.
+*
+* @see NOAA CO-OPS shallow-water constituents
+*/
+var _3N2_default = defineCompoundConstituent("3N2", [{
+	constituent: N2_default,
+	factor: 3
+}]);
+
+//#endregion
+//#region src/constituents/EP2.ts
+/**
+* Lunar elliptic semi-diurnal constituent (ε2).
+*/
+var EP2_default = defineConstituent("EP2", [
+	2,
+	-3,
+	2,
+	1,
+	0,
+	0,
 	0
 ], node_corrections_default.uM2, node_corrections_default.fM2);
-constituents.M2 = constituent_default("M2", [
+
+//#endregion
+//#region src/constituents/K2.ts
+/**
+* Lunisolar semi-diurnal (K2).
+* Combined lunar and solar semi-diurnal constituent from declination effects.
+* Amplitude typically 10-30% of S2; second-largest solar-related semi-diurnal component.
+* Important in semi-diurnal tidal analysis, especially at mid-latitudes.
+*/
+var K2_default = defineConstituent("K2", [
+	2,
 	2,
-	0,
 	0,
 	0,
 	0,
 	0,
 	0
-], node_corrections_default.uM2, node_corrections_default.fM2);
-constituents.LAM2 = constituent_default("LAM2", [
+], node_corrections_default.uK2, node_corrections_default.fK2);
+
+//#endregion
+//#region src/constituents/LAMBDA2.ts
+/**
+* Lunar semi-diurnal (λ2, lambda2).
+* Lunar semi-diurnal constituent from Moon's perigee effects.
+* Amplitude typically <5% of M2; important in detailed constituent analysis.
+* IHO standard designation (previously abbreviated as LAM2).
+*/
+var LAMBDA2_default = defineConstituent("LAMBDA2", [
 	2,
 	1,
 	-2,
@@ -694,149 +721,835 @@ constituents.LAM2 = constituent_default("LAM2", [
 	0,
 	2
 ], node_corrections_default.uM2, node_corrections_default.fM2);
-constituents.L2 = constituent_default("L2", [
-	2,
+
+//#endregion
+//#region src/constituents/LAM2.ts
+/**
+* Alias for compatibility between NOAA and IHO
+*
+* @see LAMBDA2
+*/
+var LAM2_default = LAMBDA2_default;
+
+//#endregion
+//#region src/constituents/M1.ts
+/**
+* Lunar diurnal elliptic (M1).
+* Secondary lunar diurnal constituent from Moon's elliptical orbit.
+* Typically very small amplitude; usually <1% of O1.
+* Rarely significant except in detailed harmonic analyses.
+*/
+var M1_default = defineConstituent("M1", [
 	1,
 	0,
-	-1,
 	0,
 	0,
-	2
-], node_corrections_default.uL2, node_corrections_default.fL2);
-constituents.T2 = constituent_default("T2", [
-	2,
-	2,
-	-3,
 	0,
 	0,
-	1,
-	0
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.S2 = constituent_default("S2", [
-	2,
-	2,
-	-2,
+	1
+], node_corrections_default.uM1, node_corrections_default.fM1);
+
+//#endregion
+//#region src/constituents/M3.ts
+/**
+* Lunar terdiurnal (M3).
+* Third-diurnal lunar constituent from Moon's orbital motion.
+* Typically found in shallow water and resonant systems.
+* Amplitude usually <2% of M2; important in some shallow-water analyses.
+*/
+var M3_default = defineConstituent("M3", [
+	3,
+	0,
+	0,
 	0,
 	0,
 	0,
 	0
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.R2 = constituent_default("R2", [
-	2,
+], (a) => {
+	return node_corrections_default.uModd(a, 3);
+}, (a) => {
+	return node_corrections_default.fModd(a, 3);
+});
+
+//#endregion
+//#region src/constituents/M4.ts
+/**
+* Shallow-water quarter-diurnal (M4 = 2×M2).
+* First overtide of M2; generated in shallow-water environments.
+* Amplitude typically 2-10% of M2; largest of the quarter-diurnal constituents.
+* Important indicator of shallow-water tidal distortion.
+*/
+var M4_default = defineCompoundConstituent("M4", [{
+	constituent: M2_default,
+	factor: 2
+}]);
+
+//#endregion
+//#region src/constituents/M6.ts
+/**
+* Shallow-water sixth-diurnal (M6 = 3×M2).
+* Second overtide of M2; generated in shallow-water environments.
+* Amplitude typically 0.5-3% of M2; important in extreme shallow water.
+* Indicator of significant tidal distortion and non-linear effects.
+*/
+var M6_default = defineCompoundConstituent("M6", [{
+	constituent: M2_default,
+	factor: 3
+}]);
+
+//#endregion
+//#region src/constituents/M8.ts
+/**
+* Shallow-water eighth-diurnal (M8 = 4×M2).
+* Third overtide of M2; generated in extreme shallow-water environments.
+* Amplitude typically <0.5% of M2; very rarely significant.
+* Found only in highly distorted tidal regimes (extreme shallow water, enclosed basins).
+*/
+var M8_default = defineCompoundConstituent("M8", [{
+	constituent: M2_default,
+	factor: 4
+}]);
+
+//#endregion
+//#region src/constituents/MA2.ts
+/**
+* Lunar variational semi-diurnal constituent (μ2, mu2).
+* Derived from Moon's orbital parameter variations.
+*
+* Note: Often included with M2 family in modern analysis. Minor constituent with
+* location-dependent amplitude.
+*/
+var MA2_default = defineConstituent("MA2", [
 	2,
+	0,
 	-1,
 	0,
 	0,
-	-1,
-	2
-], node_corrections_default.uZero, node_corrections_default.fUnity);
-constituents.K2 = constituent_default("K2", [
-	2,
+	0,
+	0
+], node_corrections_default.uM2, node_corrections_default.fM2);
+
+//#endregion
+//#region src/constituents/MB2.ts
+/**
+* Lunar elliptic constituent from parameter variations.
+*
+* From https://iho.int/mtg_docs/com_wg/IHOTC/IHOTC_Misc/TWCWG_Constituent_list.pdf:
+*
+* > MB2 was originally called Ma2 but this became ambiguous when spoken, or typed on
+* > computers without lower case, and so it was initially changed to MA2*. However, this
+* > in turn was thought to be clumsy and hence MB2 was finally adopted. Although
+* > theoretically they should have the same values of u and f as M2, they are so small
+* > that they are commonly treated as having values of u = 0 and f = 1.
+*
+* @see Schureman, P. (1958). Manual of Harmonic Analysis and Prediction of Tides
+*/
+var MB2_default = defineConstituent("MB2", [
 	2,
 	0,
+	1,
 	0,
 	0,
 	0,
 	0
-], node_corrections_default.uK2, node_corrections_default.fK2);
-constituents.M3 = constituent_default("M3", [
-	3,
+], node_corrections_default.uM2, node_corrections_default.fM2);
+
+//#endregion
+//#region src/constituents/MF.ts
+/**
+* Lunar fortnightly (MF).
+* Long-period constituent from lunar inequality interactions.
+* Significant in long-term water level records and coastal resonances.
+*/
+var MF_default = defineConstituent("MF", [
 	0,
+	2,
 	0,
 	0,
 	0,
 	0,
 	0
-], (a) => {
-	return node_corrections_default.uModd(a, 3);
-}, (a) => {
-	return node_corrections_default.fModd(a, 3);
-});
-constituents.MSF = compound_constituent_default("MSF", [{
-	constituent: constituents.S2,
+], node_corrections_default.uMf, node_corrections_default.fMf);
+
+//#endregion
+//#region src/constituents/MK3.ts
+/**
+* Shallow-water terdiurnal (MK3).
+* Compound constituent: M2 + K1
+* Lunisolar interaction from M2 semi-diurnal and K1 diurnal components.
+* Generated in shallow-water environments; typical amplitude 0.5-2 cm.
+* Often paired with 2MK3 in terdiurnal tide analysis.
+*/
+var MK3_default = defineCompoundConstituent("MK3", [{
+	constituent: M2_default,
 	factor: 1
 }, {
-	constituent: constituents.M2,
-	factor: -1
+	constituent: K1_default,
+	factor: 1
 }]);
-constituents["2Q1"] = compound_constituent_default("2Q1", [{
-	constituent: constituents.N2,
+
+//#endregion
+//#region src/constituents/MKS2.ts
+/**
+* Three-way shallow-water interaction of M2, K1, and S2.
+*
+* Warning: Not in standard IHO constituent bank. Shallow-water specific constituent
+* with definition varying by application and water depth. Use with caution and document
+* the specific convention used in your analysis.
+*
+* @see NOAA CO-OPS shallow-water constituents
+*/
+var MKS2_default = defineCompoundConstituent("MKS2", [
+	{
+		constituent: M2_default,
+		factor: 1
+	},
+	{
+		constituent: K1_default,
+		factor: 1
+	},
+	{
+		constituent: S2_default,
+		factor: -1
+	}
+]);
+
+//#endregion
+//#region src/constituents/MM.ts
+/**
+* Lunar monthly (MM).
+* Long-period constituent from lunar declination variations.
+* Important for long-term water level studies.
+*/
+var MM_default = defineConstituent("MM", [
+	0,
+	1,
+	0,
+	-1,
+	0,
+	0,
+	0
+], node_corrections_default.uZero, node_corrections_default.fMm);
+
+//#endregion
+//#region src/constituents/MN4.ts
+/**
+* Shallow-water quarter-diurnal (MN4).
+* Compound constituent: M2 + N2
+* Lunar-lunar elliptic interaction from M2 and N2 semi-diurnal components.
+* Generated in shallow-water environments; typical amplitude 1-5 cm.
+* Often significant in shallow seas and estuaries.
+*/
+var MN4_default = defineCompoundConstituent("MN4", [{
+	constituent: M2_default,
 	factor: 1
 }, {
-	constituent: constituents.J1,
-	factor: -1
+	constituent: N2_default,
+	factor: 1
 }]);
-constituents.RHO = compound_constituent_default("RHO", [{
-	constituent: constituents.NU2,
+
+//#endregion
+//#region src/constituents/MS4.ts
+/**
+* Shallow-water quarter-diurnal (MS4).
+* Compound constituent: M2 + S2
+* Lunar-solar interaction from M2 and S2 semi-diurnal components.
+* Generated in shallow-water environments; typical amplitude 1-5 cm.
+* Often second-largest quarter-diurnal component after M4.
+*/
+var MS4_default = defineCompoundConstituent("MS4", [{
+	constituent: M2_default,
 	factor: 1
 }, {
-	constituent: constituents.K1,
-	factor: -1
-}]);
-constituents.MU2 = compound_constituent_default("MU2", [{
-	constituent: constituents.M2,
-	factor: 2
-}, {
-	constituent: constituents.S2,
-	factor: -1
-}]);
-constituents["2SM2"] = compound_constituent_default("2SM2", [{
-	constituent: constituents.S2,
-	factor: 2
-}, {
-	constituent: constituents.M2,
-	factor: -1
+	constituent: S2_default,
+	factor: 1
 }]);
-constituents["2MK3"] = compound_constituent_default("2MK3", [{
-	constituent: constituents.M2,
+
+//#endregion
+//#region src/constituents/MSF.ts
+/**
+* Lunisolar synodic fortnightly (MSF = S2-M2).
+* Long-period constituent representing beat frequency between solar S2 and lunar M2.
+* Manifests as fortnightly modulation of tidal range (spring-neap cycle).
+* Amplitude typically 5-15% of M2; represents the primary spring-neap variation.
+* Important for tidal range predictions and coastal flooding assessments.
+*/
+var MSF_default = defineCompoundConstituent("MSF", [{
+	constituent: S2_default,
 	factor: 1
 }, {
-	constituent: constituents.O1,
-	factor: 1
+	constituent: M2_default,
+	factor: -1
 }]);
-constituents.MK3 = compound_constituent_default("MK3", [{
-	constituent: constituents.M2,
+
+//#endregion
+//#region src/constituents/MSQM.ts
+/**
+* Lunar-solar interaction compound constituent.
+*
+* Warning: Non-standard constituent not in IHO or modern NOAA standard tables.
+* Definition varies significantly across sources. Rarely used in modern tide prediction.
+* Appears in Schureman's tables as variant shallow-water interaction.
+*
+* @see NOAA CO-OPS shallow-water constituents
+* @see Schureman shallow-water analysis tables
+*/
+var MSQM_default = defineCompoundConstituent("MSQM", [
+	{
+		constituent: M2_default,
+		factor: 1
+	},
+	{
+		constituent: S2_default,
+		factor: 1
+	},
+	{
+		constituent: K1_default,
+		factor: 1
+	}
+]);
+
+//#endregion
+//#region src/constituents/T2.ts
+/**
+* Solar elliptic semi-diurnal (T2).
+* Larger solar elliptic semi-diurnal constituent from solar declination effects.
+* Amplitude typically <5% of S2; increases at higher latitudes.
+* Forms part of solar semi-diurnal analysis alongside S2 and R2.
+*/
+var T2_default = defineConstituent("T2", [
+	2,
+	2,
+	-3,
+	0,
+	0,
+	1,
+	0
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/MTM.ts
+/**
+* Lunar-solar shallow-water interaction (M2 modulated by lunar orbit).
+*
+* Warning: Not in modern IHO standard constituents; mostly historical interest.
+* Often replaced by specific ν2, λ2, or other lunar elliptic terms in modern analysis.
+* Amplitude is location and depth-dependent.
+*
+* @see NOAA CO-OPS shallow-water variants
+* @see Schureman Manual
+*/
+var MTM_default = defineCompoundConstituent("MTM", [{
+	constituent: M2_default,
 	factor: 1
 }, {
-	constituent: constituents.K1,
+	constituent: T2_default,
 	factor: 1
 }]);
-constituents.MN4 = compound_constituent_default("MN4", [{
-	constituent: constituents.M2,
-	factor: 1
+
+//#endregion
+//#region src/constituents/MU2.ts
+/**
+* Shallow-water semi-diurnal (MU2 or μ2).
+* Compound constituent: 2×M2 - S2
+* Lunar variational constituent generated in shallow-water environments.
+* Amplitude typically <2% of M2; found in coastal shallow-water predictions.
+*/
+var MU2_default = defineCompoundConstituent("MU2", [{
+	constituent: M2_default,
+	factor: 2
 }, {
-	constituent: constituents.N2,
-	factor: 1
+	constituent: S2_default,
+	factor: -1
 }]);
-constituents.M4 = compound_constituent_default("M4", [{
-	constituent: constituents.M2,
+
+//#endregion
+//#region src/constituents/N4.ts
+/**
+* Second overtide of N2; shallow-water quarter-diurnal harmonic.
+* Amplitude ranges 0.25 to 2.25 times mean due to 18.613-year lunar node cycle.
+*
+* Note: Shallow-water constituent, typically found in tide predictions for coastal areas.
+*/
+var N4_default = defineCompoundConstituent("N4", [{
+	constituent: N2_default,
 	factor: 2
 }]);
-constituents.MS4 = compound_constituent_default("MS4", [{
-	constituent: constituents.M2,
+
+//#endregion
+//#region src/constituents/NU2.ts
+/**
+* Lunar elliptic semi-diurnal (NU2).
+* Secondary lunar elliptic semi-diurnal from Moon's orbital variations.
+* Amplitude typically 2-5% of M2; smaller than N2.
+* Important in detailed semi-diurnal constituent analysis.
+*/
+var NU2_default = defineConstituent("NU2", [
+	2,
+	-1,
+	2,
+	-1,
+	0,
+	0,
+	0
+], node_corrections_default.uM2, node_corrections_default.fM2);
+
+//#endregion
+//#region src/constituents/OO1.ts
+/**
+* Lunar diurnal (OO1).
+* Second-order lunar diurnal constituent from Moon's orbital eccentricity.
+* Typically very small amplitude; <2% of O1 in most locations.
+* Important for detailed harmonic analysis in some regions.
+*/
+var OO1_default = defineConstituent("OO1", [
+	1,
+	3,
+	0,
+	0,
+	0,
+	0,
+	-1
+], node_corrections_default.uOO1, node_corrections_default.fOO1);
+
+//#endregion
+//#region src/constituents/P1.ts
+/**
+* Solar diurnal (P1).
+* Principal solar diurnal constituent; one-per-solar-day oscillation.
+* Amplitude typically 1/3 of K1; varies with latitude.
+* Forms key component of diurnal tidal analysis alongside K1 and O1.
+*/
+var P1_default = defineConstituent("P1", [
+	1,
+	1,
+	-2,
+	0,
+	0,
+	0,
+	1
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/Q1.ts
+/**
+* Elliptic lunar diurnal (Q1).
+* One solar day and one lunar day modulation interaction.
+* Amplitude typically 2-5% of O1; important in diurnal constituent analysis.
+*/
+var Q1_default = defineConstituent("Q1", [
+	1,
+	-2,
+	0,
+	1,
+	0,
+	0,
+	1
+], node_corrections_default.uO1, node_corrections_default.fO1);
+
+//#endregion
+//#region src/constituents/R2.ts
+/**
+* Solar elliptic semi-diurnal (R2).
+* Smaller solar elliptic semi-diurnal constituent from solar declination effects.
+* Amplitude typically <2% of S2; smallest of the solar semi-diurnal group.
+* Rarely significant except in very detailed analyses.
+*/
+var R2_default = defineConstituent("R2", [
+	2,
+	2,
+	-1,
+	0,
+	0,
+	-1,
+	2
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/R3.ts
+/**
+* Solar elliptic terdiurnal; 1.5 times R2 (smaller solar elliptic semi-diurnal).
+*
+* Note: Generated only in shallow water; very small amplitude (<0.05 cm typical).
+* Usually negligible except in extreme shallow-water environments.
+*
+* @see NOAA CO-OPS
+*/
+var R3_default = defineCompoundConstituent("R3", [{
+	constituent: R2_default,
+	factor: 1.5
+}]);
+
+//#endregion
+//#region src/constituents/RHO1.ts
+/**
+* Shallow-water diurnal (RHO1, ρ1).
+* Compound constituent: NU2-K1 (ν2-K1)
+* Derived from interaction of semi-diurnal NU2 and diurnal K1 constituents.
+* Amplitude typically very small; rarely significant.
+* Found only in shallow-water regions with strong tidal distortion.
+*/
+var RHO1_default = defineCompoundConstituent("RHO1", [{
+	constituent: NU2_default,
 	factor: 1
 }, {
-	constituent: constituents.S2,
-	factor: 1
+	constituent: K1_default,
+	factor: -1
 }]);
-constituents.S4 = compound_constituent_default("S4", [{
-	constituent: constituents.S2,
-	factor: 2
+
+//#endregion
+//#region src/constituents/RHO.ts
+/**
+* Alias for compatibility between NOAA and IHO
+*
+* @see RHO1
+*/
+var RHO_default = RHO1_default;
+
+//#endregion
+//#region src/constituents/S1.ts
+/**
+* Solar diurnal (S1).
+* Secondary solar diurnal constituent related to solar declination inequality.
+* Usually very small amplitude; typically dominated by K1 and P1 in diurnal analysis.
+* Rarely included in routine harmonic analyses.
+*/
+var S1_default = defineConstituent("S1", [
+	1,
+	1,
+	-1,
+	0,
+	0,
+	0,
+	0
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/S3.ts
+/**
+* Solar terdiurnal overtide; shallow-water only.
+*
+* Note: Generated only in shallow water; not found in deep ocean.
+* Typical amplitude <1 cm except in extreme shallow-water environments.
+*
+* @see NOAA CO-OPS
+*/
+var S3_default = defineCompoundConstituent("S3", [{
+	constituent: S2_default,
+	factor: 1.5
 }]);
-constituents.M6 = compound_constituent_default("M6", [{
-	constituent: constituents.M2,
-	factor: 3
+
+//#endregion
+//#region src/constituents/S4.ts
+/**
+* Shallow-water quarter-diurnal (S4).
+* Compound constituent: 2×S2
+* First overtide of S2; generated in shallow-water environments.
+* Amplitude typically 1-3% of S2; smallest of the common quarter-diurnal constituents.
+* Purely solar, so has fixed amplitude and phase at any location.
+*/
+var S4_default = defineCompoundConstituent("S4", [{
+	constituent: S2_default,
+	factor: 2
 }]);
-constituents.S6 = compound_constituent_default("S6", [{
-	constituent: constituents.S2,
+
+//#endregion
+//#region src/constituents/S6.ts
+/**
+* Shallow-water sixth-diurnal (S6).
+* Compound constituent: 3×S2
+* Second overtide of S2; generated in shallow-water environments.
+* Amplitude typically 0.5-1% of S2; smallest of the common overtides.
+* Rarely significant except in extreme shallow water or resonant systems.
+*/
+var S6_default = defineCompoundConstituent("S6", [{
+	constituent: S2_default,
 	factor: 3
 }]);
-constituents.M8 = compound_constituent_default("M8", [{
-	constituent: constituents.M2,
-	factor: 4
+
+//#endregion
+//#region src/constituents/SA.ts
+/**
+* Solar annual (Sa).
+* Long-term constituent driven by solar declination variations over the year.
+*/
+var SA_default = defineConstituent("Sa", [
+	0,
+	0,
+	1,
+	0,
+	0,
+	0,
+	0
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/SGM.ts
+/**
+* Lunar diurnal variational constituent (σ1, sigma1).
+* Derived from Moon's declination variations.
+*
+* Note: Often has small amplitude; closely related to K1 and O1 variations.
+*/
+var SGM_default = defineConstituent("SGM", [
+	1,
+	-3,
+	2,
+	0,
+	0,
+	0,
+	-1
+], node_corrections_default.uO1, node_corrections_default.fO1);
+
+//#endregion
+//#region src/constituents/SSA.ts
+/**
+* Solar semi-annual (Ssa).
+* Semi-annual constituent from solar declination with twice-yearly periodicity.
+*/
+var SSA_default = defineConstituent("Ssa", [
+	0,
+	0,
+	2,
+	0,
+	0,
+	0,
+	0
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/T3.ts
+/**
+* Solar elliptic terdiurnal; 1.5 times T2 (larger solar elliptic semi-diurnal).
+*
+* Note: Generated only in shallow water; minimal amplitude (<0.1 cm typical).
+* Rarely significant except in extreme shallow-water or enclosed basins.
+*
+* @see NOAA CO-OPS
+*/
+var T3_default = defineCompoundConstituent("T3", [{
+	constituent: T2_default,
+	factor: 1.5
 }]);
+
+//#endregion
+//#region src/constituents/Z0.ts
+/**
+* Mean sea level (Z0).
+* Not a tidal constituent in the strict sense, but represents the mean sea level offset
+* or the "zero" reference level used in tidal predictions.
+* No astronomical forcing; typically determined from harmonic analysis of observed data.
+*/
+var Z0_default = defineConstituent("Z0", [
+	0,
+	0,
+	0,
+	0,
+	0,
+	0,
+	0
+], node_corrections_default.uZero, node_corrections_default.fUnity);
+
+//#endregion
+//#region src/constituents/index.ts
+const constituentModules = {
+	"./2MK3.ts": _2MK3_default,
+	"./2MK5.ts": _2MK5_default,
+	"./2MO5.ts": _2MO5_default,
+	"./2MS6.ts": _2MS6_default,
+	"./2N2.ts": _2N2_default,
+	"./2Q1.ts": _2Q1_default,
+	"./2SM2.ts": _2SM2_default,
+	"./3L2.ts": _3L2_default,
+	"./3N2.ts": _3N2_default,
+	"./EP2.ts": EP2_default,
+	"./J1.ts": J1_default,
+	"./K1.ts": K1_default,
+	"./K2.ts": K2_default,
+	"./L2.ts": L2_default,
+	"./LAM2.ts": LAM2_default,
+	"./LAMBDA2.ts": LAMBDA2_default,
+	"./M1.ts": M1_default,
+	"./M2.ts": M2_default,
+	"./M3.ts": M3_default,
+	"./M4.ts": M4_default,
+	"./M6.ts": M6_default,
+	"./M8.ts": M8_default,
+	"./MA2.ts": MA2_default,
+	"./MB2.ts": MB2_default,
+	"./MF.ts": MF_default,
+	"./MK3.ts": MK3_default,
+	"./MKS2.ts": MKS2_default,
+	"./MM.ts": MM_default,
+	"./MN4.ts": MN4_default,
+	"./MS4.ts": MS4_default,
+	"./MSF.ts": MSF_default,
+	"./MSQM.ts": MSQM_default,
+	"./MTM.ts": MTM_default,
+	"./MU2.ts": MU2_default,
+	"./N2.ts": N2_default,
+	"./N4.ts": N4_default,
+	"./NU2.ts": NU2_default,
+	"./O1.ts": O1_default,
+	"./OO1.ts": OO1_default,
+	"./P1.ts": P1_default,
+	"./Q1.ts": Q1_default,
+	"./R2.ts": R2_default,
+	"./R3.ts": R3_default,
+	"./RHO.ts": RHO_default,
+	"./RHO1.ts": RHO1_default,
+	"./S1.ts": S1_default,
+	"./S2.ts": S2_default,
+	"./S3.ts": S3_default,
+	"./S4.ts": S4_default,
+	"./S6.ts": S6_default,
+	"./SA.ts": SA_default,
+	"./SGM.ts": SGM_default,
+	"./SSA.ts": SSA_default,
+	"./T2.ts": T2_default,
+	"./T3.ts": T3_default,
+	"./Z0.ts": Z0_default,
+	"./definition.ts": definition_default
+};
+const constituents = {};
+for (const [path, module] of Object.entries(constituentModules)) {
+	if (path.includes("index.ts")) continue;
+	const name = path.split("/").pop()?.replace(/\..*$/, "") ?? "";
+	if (name === "definition") continue;
+	constituents[name] = module;
+}
 var constituents_default = constituents;
 
+//#endregion
+//#region src/harmonics/prediction.ts
+const modulus = (a, b) => {
+	return (a % b + b) % b;
+};
+const addExtremesOffsets = (extreme, offsets) => {
+	if (typeof offsets === "undefined" || !offsets) return extreme;
+	if (extreme.high && offsets.height?.high) if (offsets.height.type === "fixed") extreme.level += offsets.height.high;
+	else extreme.level *= offsets.height.high;
+	if (extreme.low && offsets.height?.low) if (offsets.height.type === "fixed") extreme.level += offsets.height.low;
+	else extreme.level *= offsets.height.low;
+	if (extreme.high && offsets.time?.high) extreme.time = new Date(extreme.time.getTime() + offsets.time.high * 60 * 1e3);
+	if (extreme.low && offsets.time?.low) extreme.time = new Date(extreme.time.getTime() + offsets.time.low * 60 * 1e3);
+	return extreme;
+};
+const getExtremeLabel = (label, highLowLabels) => {
+	if (typeof highLowLabels !== "undefined" && typeof highLowLabels[label] !== "undefined") return highLowLabels[label];
+	return {
+		high: "High",
+		low: "Low"
+	}[label];
+};
+const predictionFactory = ({ timeline, constituents: constituents$1, start }) => {
+	const getLevel = (hour, modelBaseSpeed, modelU, modelF, modelBaseValue) => {
+		const amplitudes = [];
+		let result = 0;
+		constituents$1.forEach((constituent) => {
+			const amplitude = constituent.amplitude;
+			const phase = constituent.phase;
+			const f = modelF[constituent.name];
+			const speed = modelBaseSpeed[constituent.name];
+			const u = modelU[constituent.name];
+			const V0 = modelBaseValue[constituent.name];
+			amplitudes.push(amplitude * f * Math.cos(speed * hour + (V0 + u) - phase));
+		});
+		amplitudes.forEach((item) => {
+			result += item;
+		});
+		return result;
+	};
+	const prediction = {};
+	prediction.getExtremesPrediction = (options) => {
+		const { labels, offsets } = typeof options !== "undefined" ? options : {};
+		const results = [];
+		const { baseSpeed, u, f, baseValue } = prepare();
+		let goingUp = false;
+		let goingDown = false;
+		let lastLevel = getLevel(0, baseSpeed, u[0], f[0], baseValue);
+		timeline.items.forEach((time, index) => {
+			const hour = timeline.hours[index];
+			const level = getLevel(hour, baseSpeed, u[index], f[index], baseValue);
+			if (level > lastLevel && goingDown) results.push(addExtremesOffsets({
+				time: timeline.items[index - 1],
+				level: lastLevel,
+				high: false,
+				low: true,
+				label: getExtremeLabel("low", labels)
+			}, offsets));
+			if (level < lastLevel && goingUp) results.push(addExtremesOffsets({
+				time: timeline.items[index - 1],
+				level: lastLevel,
+				high: true,
+				low: false,
+				label: getExtremeLabel("high", labels)
+			}, offsets));
+			if (level > lastLevel) {
+				goingUp = true;
+				goingDown = false;
+			}
+			if (level < lastLevel) {
+				goingUp = false;
+				goingDown = true;
+			}
+			lastLevel = level;
+		});
+		return results;
+	};
+	prediction.getTimelinePrediction = () => {
+		const results = [];
+		const { baseSpeed, u, f, baseValue } = prepare();
+		timeline.items.forEach((time, index) => {
+			const hour = timeline.hours[index];
+			const prediction$1 = {
+				time,
+				hour,
+				level: getLevel(hour, baseSpeed, u[index], f[index], baseValue)
+			};
+			results.push(prediction$1);
+		});
+		return results;
+	};
+	const prepare = () => {
+		const baseAstro = astronomy_default(start);
+		const baseValue = {};
+		const baseSpeed = {};
+		const u = [];
+		const f = [];
+		constituents$1.forEach((constituent) => {
+			const model = constituents_default[constituent.name];
+			if (!model) return;
+			const value = model.value(baseAstro);
+			const speed = model.speed(baseAstro);
+			baseValue[constituent.name] = d2r * value;
+			baseSpeed[constituent.name] = d2r * speed;
+		});
+		timeline.items.forEach((time) => {
+			const uItem = {};
+			const fItem = {};
+			const itemAstro = astronomy_default(time);
+			constituents$1.forEach((constituent) => {
+				const model = constituents_default[constituent.name];
+				if (!model) return;
+				const constituentU = modulus(model.u(itemAstro), 360);
+				uItem[constituent.name] = d2r * constituentU;
+				fItem[constituent.name] = modulus(model.f(itemAstro), 360);
+			});
+			u.push(uItem);
+			f.push(fItem);
+		});
+		return {
+			baseValue,
+			baseSpeed,
+			u,
+			f
+		};
+	};
+	return Object.freeze(prediction);
+};
+var prediction_default = predictionFactory;
+
 //#endregion
 //#region src/harmonics/index.ts
 const getDate = (time) => {
@@ -867,13 +1580,11 @@ const harmonicsFactory = ({ harmonicConstituents, offset }) => {
 		if (typeof constituent.name === "undefined") throw new Error("Harmonic constituents must have a name property");
 		if (constituents_default[constituent.name] !== void 0) constituents$1.push({
 			...constituent,
-			_model: constituents_default[constituent.name],
 			phase: d2r * constituent.phase
 		});
 	});
 	if (offset !== false) constituents$1.push({
 		name: "Z0",
-		_model: constituents_default.Z0,
 		phase: 0,
 		amplitude: offset
 	});